In this study, Version 5.16 of the WAVEWATCH-III (WW3) model is used to simulate parameters of typhoon-generated wave fields in the Western North Pacific Ocean during the period 1998–2017. From a database of more than 300 typhoons, typhoon tracks are partitioned into six groups by their direction of motion and longitude of recurvature track. For typhoons that recurve east of 140° E, or track toward mainland Asia, regions of high significant wave height (SWH) values are separated by a minimum in SWH near 30° N. Partitioning SWH into wind sea and swell components demonstrates that variations in typhoon tracks produce a much stronger signal in the wind sea component of the wave system. Empirical orthogonal function (EOF) analysis is used to compute the four leading modes of variation in average SWH simulated by the WW3 model. The first EOF mode contributes to 17.3% of the total variance; all other modes contribute less than 10%. The first EOF mode also oscillates on an approximately 1-year cycle during the period 1998–2017. Overall, typhoon-induced wave energy dominates north of 30° N. Temporal analysis of the leading principal component of SWH indicates that (a) the intensity of the wave pattern produced by westward-tracking typhoons decreased during the last 20 years, and (b) typhoons that recurve east of 140° E and those that track westward toward southeast Asia are largely responsible for the decadal variability of typhoon-induced wave distribution.

In this study, an empirical algorithm is proposed to retrieve significant wave height (SWH) from dual-polarization Sentinel-1 (S-1) synthetic aperture radar (SAR) imagery collected under cyclonic conditions. The retrieval scheme is based on the well-known CWAVE empirical function that is here updated to deal with multi-polarization S-1 SAR measurements collected using the interferometric wide (IW) and the Extra Wide-Swath (EW) imaging modes, under cyclonic conditions. First, a training dataset that consists of six S-1 SAR images collected under cyclonic conditions is exploited to both tune the retrieval function and to check the soundness of the retrievals against the co-located WAVEWATCH-III (WW3) numerical simulations. The comparison of simulation from the WW3 model and measurements from altimeter Jason-2 shows a 0.29m root mean square error (RMSE) of significant wave height (SWH). Then, a testing data-set that consists of two S-1 SAR images is exploited to provide a preliminary validation. The results, verified against both WW3 and European Centre for Medium-Range Weather Forecasts (ECMWF) data, show the soundness of the herein approach.

In this work, three types of machine learning algorithms are applied for synthetic aperture radar (SAR) wind retrieval in tropical cyclones (TCs), and the optimal method is confirmed. In total, 30 Sentinel-1 (S-1) images in dual-polarization (vertical–vertical [VV] and vertical–horizontal [VH] were collected during the period from 2016 to 2021, which were acquired in interferometric-wide and extra-wide modes with pixels of 10 m and 40 m, respectively. More than 100,000 sub-scenes with a spatial coverage of 3 km are extracted from these images. The dependences of variables estimated from sub-scenes, i.e., VV-polarized and VH-polarized normalized radar cross-section (NRCS), as well as the azimuthal wave cutoff wavelength, on wind speeds from the stepped-frequency microwave radiometer (SFMR) and the soil moisture active passive (SMAP) radiometer are studied, showing the linear relations between wind speed and these three parameters; however, the saturation of VV-polarized NRCS and the azimuthal wave cutoff wavelength is observed. This is the foundation of selecting input variables in machine learning algorithms. Two-thirds of the collocated dataset (20 images) are used for training the process using three machine learning algorithms, i.e., eXtreme Gradient Boosting (XGBoost), Multi-layer Perceptron, and K-Nearest Neighbor, and the coefficients are fitted after training completion through 20 images collocated with SFMR and SMAP data. Another 10 images are taken for validation up to 70 m/s, yielding a 2.53 m/s root mean square error (RMSE) with a 0.96 correlation and 0.12 scatter index (SI) using XGBoost. The result is better than the >5 m/s error achieved using the existing cross-polarized geophysical model function and the other two machine learning algorithms; moreover, the comparison between wind retrievals using XGBoost and Level-2 CyclObs products shows about 4 m/s RMSE and 0.18 SI. This suggests that the machine learning algorithm XGBoost is an effective method for inverting the TC wind field utilizing SAR measurements in dual-polarization.

The Chinese-French Oceanography SATellite (CFOSAT) jointly developed by the Chinese National Space Agency (CNSA) and the Centre National d’Etudes Spatiales (CNES) of France carries a wave spectrometer (Surface Waves Investigation and Monitoring, SWIM). SWIM has one nadir and five off-nadir beams to measure ocean surface waves. These near-nadir beams range from 0° to 10° at an interval of 2°. In this work, we investigated the performance of wave parameters derived from wave spectra measured by SWIM at off-nadir beams during the period 2020 to December 2022, e.g., incidence angles of 6°, 8° and 10°, which were collocated with the wave simulated by Simulating Waves Nearshore (SWAN). The validation of SWAN-simulated significant wave heights (SWHs) against National Data Buoy Center (NDBC) buoys of National Oceanic and Atmospheric Administration (NOAA) exhibited a 0.42 m root mean square error (RMSE) in the SWH. Our results revealed a RMSE of 1.02 m for the SWIM-measured SWH in the East Pacific Ocean compared with the SWH simulated by SWAN, as well as a 0.79 correlation coefficient (Cor) and a 1.17 squared error (Err) for the wave spectrum at an incidence angle of 10°, which are better than those (i.e., the RMSEs were > 1.1 m with Cors < 0.76 and Errs > 1.2) achieved at other incidence angles of SWH up to 14 m. This analysis indicates that the SWIM product is a relevant resource for wave monitoring over global seas. The collocated wave retrievals for more than 300 cases from Gaofen-3 (GF-3) synthetic aperture radar (SAR) images in China Seas were also used to verify the accuracy of SWIM-measured wave spectra. The energy of the SWIM-measured wave spectra represented by SWH was found to decrease with an increasing incidence angle in a case study. Moreover, the SWIM-measured wave spectra were most consistent with the SAR-derived wave spectra at an incidence angle of 10°, yielding a 0.77 Cor and 1.98 Err between SAR-derived and SWIM wave spectra under regular sea state conditions (SWH < 2 m). The error analysis indicates that the difference in SWH between SWIM at an incidence angle of 10° and SWAN has an increasing tendency with the growth in sea surface wind and sea state and it stabilizes to be 0.6 m at SWH > 4 m; however, the current and sea level have less influence on the uncertainties of the SWIM product.

Marine oil spills occur suddenly and pose a serious threat to ecosystems in coastal waters. Oil spills continuously affect the ocean environment for years. In this study, the oil spill caused by the accident of the Sanchi ship (2018) in the East China Sea was hindcast simulated using the oil particle-tracing method. Sea-surface winds from the European Centre for Medium-Range Weather Forecasts (ECMWF), currents simulated from the Finite-Volume Community Ocean Model (FVCOM), and waves simulated from the Simulating WAves Nearshore (SWAN) were employed as background marine dynamics fields. In particular, the oil spill simulation was compared with the detection from Chinese Gaofen-3 (GF-3) synthetic aperture radar (SAR) images. The validation of the SWAN-simulated significant wave height (SWH) against measurements from the Jason-2 altimeter showed a 0.58 m root mean square error (RMSE) with a 0.93 correlation (COR). Further, the sea-surface current was compared with that from the National Centers for Environmental Prediction (NCEP) Climate Forecast System Version 2 (CFSv2), yielding a 0.08 m/s RMSE and a 0.71 COR. Under these circumstances, we think the model-simulated sea-surface currents and waves are reliable for this work. A hindcast simulation of the tracks of oil slicks spilled from the Sanchi shipwreck was conducted during the period of 14–17 January 2018. It was found that the general track of the simulated oil slicks was consistent with the observations from the collected GF-3 SAR images. However, the details from the GF-3 SAR images were more obvious. The spatial coverage of oil slicks between the SAR-detected and simulated results was about 1 km². In summary, we conclude that combining numerical simulation and SAR remote sensing is a promising technique for real-time oil spill monitoring and the prediction of oil spreading.

Extreme surface winds and wave heights of tropical cyclones (TCs)—pose serious threats to coastal community, infrastructure and environments. In recent decades, progress in numerical wave modeling has significantly enhanced the ability to reconstruct and predict wave behavior. This review offers an in-depth overview of TC-related wave modeling utilizing different computational schemes, with a special attention to WAVEWATCH III (WW3) and Simulating Waves Nearshore (SWAN). Due to the complex air–sea interactions during TCs, it is challenging to obtain accurate wind input data and optimize the parameterizations. Substantial spatial and temporal variations in water levels and current patterns occurs when coastal circulation is modulated by varying underwater topography. To explore their influence on waves, this study employs a coupled SWAN and Finite-Volume Community Ocean Model (FVCOM) modeling approach. Additionally, the interplay between wave and sea surface temperature (SST) is investigated by incorporating four key wave-induced forcing through breaking and non-breaking waves, radiation stress, and Stokes drift from WW3 into the Stony Brook Parallel Ocean Model (sbPOM). 20 TC events were analyzed to evaluate the performance of the selected parameterizations of external forcings in WW3 and SWAN. Among different nonlinear wave interaction schemes, Generalized Multiple Discrete Interaction Approximation (GMD) Discrete Interaction Approximation (DIA) and the computationally expensive Wave-Ray Tracing (WRT) A refined drag coefficient (C_d) equation, applied within an upgraded ST6 configuration, reduce significant wave height (SWH) prediction errors and the root mean square error (RMSE) for both SWAN and WW3 wave models. Surface currents and sea level variations notably altered the wave energy and wave height distributions, especially in the area with strong TC-induced oceanic current. Finally, coupling four wave-induced forcings into sbPOM enhanced SST simulation by refining heat flux estimates and promoting vertical mixing. Validation against Argo data showed that the updated sbPOM model achieved an RMSE as low as 1.39 m, with correlation coefficients nearing 0.9881.