Several aspects of global climate change, e.g., the rise of sea level and water temperature anomalies, suggest the advantages of studying wave distributions. In this study, WAVEWATCH-III (WW3) (version 6.07), which is a well-known numerical wave model, was employed for simulating waves over global seas from 1993–2020. The European Centre for Medium-Range Weather Forecasts (ECMWF), Copernicus Marine Environment Monitoring Service (CMEMS), current and sea level were used as the forcing fields in the WW3 model. The validation of modelling simulations against the measurements from the National Data Buoy Center (NDBC) buoys and Haiyang-2B (HY-2B) altimeter yielded a root mean square error (RMSE) of 0.49 m and 0.63 m, with a correlation (COR) of 0.89 and 0.90, respectively. The terms calculated by WW3-simulated waves, i.e., breaking waves, nonbreaking waves, radiation stress, and Stokes drift, were included in the water temperature simulation by a numerical circulation model named the Stony Brook Parallel Ocean Model (sbPOM). The water temperature was simulated in 2005–2015 using the high-quality Simple Ocean Data Assimilation (SODA) data. The validation of sbPOM-simulated results against the measurements obtained from the Array for Real-time Geostrophic Oceanography (Argo) buoys yielded a RMSE of 1.12 °C and a COR of 0.99. By the seasonal variation, the interrelation of the currents, sea level anomaly, and significant wave heights (SWHs) were strong in the Indian Ocean. In the strong current areas, the distribution of the sea level was consistent with the SWHs. The monthly variation of SWHs, currents, sea surface elevation, and sea level anomalies revealed that the upward trends of SWHs and sea level anomalies were consistent from 1993–2015 over the global ocean. In the Indian Ocean, the SWHs were obviously influenced by the SST and sea surface wind stress. The rise of wind stress intensity and sea level enlarges the growth of waves, and the wave-induced terms strengthen the heat exchange at the air–sea layer. It was assumed that the SST oscillation had a negative response to the SWHs in the global ocean from 2005–2015. This feedback indicates that the growth of waves could slow down the amplitude of water warming.

The purpose of our work was to assess the feasibility of hindcasting waves using WAVEWATCH-III (WW3) in a typhoon by assembling winds from multiple remote-sensed products. During the typhoon season in 2021–2022, the swath wind products in the Western Pacific Ocean were collected from scatterometers and radiometers. Cyclonic winds with a spatial resolution of 0.125° at intervals of 6 h were obtained by assembling the remote-sensed winds from those satellites. The maximum wind speeds, V_max, were verified using the reanalysis data from the National Hurricane Center (NHC), yielding a root-mean-squared error (RMSE) of 4.79 m/s and a scatter index (SI) value of 0.2. The simulated wave spectrum was compared with the measurements from Surface Waves Investigation and Monitoring (SWIM) carried out on the Chinese–French Oceanography Satellite (CFOSAT), yielding a correlation coefficient (Cor) of 0.80, squared error (Err) of 0.49, RMSE of significant wave height (SWH) of 0.48 m with an SI of 0.25, and an RMSE of the peak wave period (PWP) of 0.95 s with an SI of 0.10. The bias of wave (WW3 minus European Centre for Medium-Range Weather Forecasts (ECMWFs) reanalysis (ERA-5)) concerning the bias of wind (assembling minus ERA-5) showed that the WW3-simulated SWH with the assembling wind forcing was significantly higher than that with the ERA-5 wind forcing. Moreover, the bias of SWH gradually increased with an increasing bias of wind speed; i.e., the bias of SWH increased up to 4 m as the bias of wind speed reached 30 m/s. It was concluded that the assembling wind from multiple scatterometers and radiometers is a promising source for wave simulations via WW3 in typhoons.

The strong wind-induced current and sea level have influences on the wave distribution in a tropical cyclone (TC). In particular, the wave–current interaction is significant in the period in which the TC passed the mesoscale eddy. In this study, the wave fields of Typhoon Chan-hom (2015) are hindcastly simulated using a coupled oceanic model that utilizes a nested triangle grid, i.e., the finite-volume community ocean model-simulating waves nearshore (FVCOM-SWAVE) model. The forcing wind field is composited from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis data and the simulation using a parametric Holland model, denoted as H-E. The open boundary fields include tide data from TPOX.5 and the hybrid coordinate ocean model (HYCOM) global datasets, including sea surface temperature (SST), sea surface salinity, sea surface current, and sea level data. The simulated oceanic parameters (e.g., the significant wave height, SWH) are validated against the measurements from the Jason-2 altimeter, yielding a root mean square error (RMSE) of 0.58 m for the SWH, a correlation (COR) coefficient of 0.94, and a scatter index (SI) of 0.23. Similarly, the simulated SSTs are compared with the remote sensing products of the remote sensing system (REMSS) and the measurements from Argos, yielding an RMSE of 0.95, and an SI of 80% near the typhoon center, and the ratio is reduced to approximately <20% away from the typhoon center, indicating that Stokes transport is an essential aspect in the water mixing during a TC. The mesoscale eddies are detected by the sea level anomalies (SLA) fusion data from AVISO. It is found that the significant wave heights, Stokes drift, and Stokes transport inside the eddy area were higher than those outside the eddy area. These parameters inside the cold mesoscale eddies were higher than t inside the warm mesoscale eddies. Otherwise, the SST mainly increased within the cold mesoscale eddies area, while decreased within the warm mesoscale eddies area. The influence of mesoscale eddies on the SST was in proportion to the eddy radius and eddy EKE.