Long-Term Trends of Sea Surface Wind in the Northern South China Sea under the Background of Climate Change

: The long-term trends of sea surface wind are of great importance to our understanding of the effects of climate change on the marine environment. In the northern South China Sea (SCS), the long-term changes in coastal sea surface wind are not well-understood. Based on the latest reanalysis (ERA5) data from 1979 to 2019, our analysis showed a decreasing trend in the annual mean wind speed in the coastal area and an increasing trend in the open sea. There was a signiﬁcant weakening trend in the easterly wind component in the coastal and continental shelf areas, whereas there was an increasing trend in the northerly wind component in the open sea. The Mann–Kendall mutation analysis suggested that there were signiﬁcant changes in the wind speed and frequency of strong wind. Signiﬁcant correlations were found between the variation of the wind ﬁeld and El Niño–Southern Oscillation by wave coherence analysis. The strengthening of the wind stress curl was an important factor for the enhancement of coastal upwelling along the coast of the northern SCS. The wind ﬁeld plays an important role in modulating the climatic change of signiﬁcant wave height.


Introduction
Climate change has become a prominent focus for both marine and atmospheric research worldwide. Global warming has resulted in different degrees of temperature increase in the oceans and on land. Such a temperature change alters the monsoon system and thus causes changes in the sea surface wind, ocean currents, sea surface temperature (SST) and precipitation, among other factors, over different timescales [1][2][3][4][5][6][7]. The South China Sea (SCS) is a marginal sea of the western Pacific Ocean. It resides at the junction of the Eurasian, Pacific and Indo-Australian plates ( Figure 1). The East Asian monsoon prevails over the SCS. The location of the SCS makes it highly sensitive to climatic changes in atmospheric-marine processes [7,8]. Investigating the long-term trends of sea surface wind in the SCS is important for understanding the effect of climate change on the local marine environment.
Many studies have focused on the change in the East Asian monsoon in response to climate change. Using the monsoon intensity index, Shi [9] investigated the trend of winter monsoon intensity in East Asia from 1950 to 1989 and found that the winter monsoon showed a decreasing trend while the temperature in mainland China showed a clear increasing trend. Analyses of the daily ground data from 194 stations and the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR)'s reanalysis data indicated that the southwest monsoon intensity in China decreased rapidly from 1957 to 1978 and then increased from 1979 to 2000, whereas the southeast monsoon intensity decreased consistently from 1957 to 2000 [10]. After analyzing the observed data at meteorological stations, Jiang et al. [11] found that the temperature and pressure difference between the Asian continent and the Pacific Ocean decreased significantly owing to global warming. Therefore, the East Asian monsoon weakens in both winter and summer, which leads to the weakening of the average wind speed. The results also indicated decreasing trends, albeit different in different regions, in the annual mean wind speed, days of wind gusts and maximum wind speed in mainland China. Jiang et al. [12] used the global climate model to predict the future trends of wind in China and indicated that, compared with that in 1981-2000, the maximum wind speed will decrease in 2046-2065 and 2080-2099. After exploring the meteorological observation data from 1961 to 2012, Shi et al. [13] also demonstrated a weakening trend in wind speed in Guangdong Province (adjacent to the southern coast of China; Figure 2a). In addition, the analysis of 51 years of data  from the Guangzhou meteorological observation station showed a weakening trend in the local annual maximum wind speed [14]. Using NCEP/NCAR and ERA40 data from 1948 to 2007, Wang et al. [1] analyzed the winter monsoon wind speed at 850 hPa and found no significant changes in the long-term trend of the meridional wind component over the central SCS. However, both reanalysis datasets showed a significant increasing trend in the easterly wind component. The wind indices defined in the lower troposphere during the summer monsoon period showed a coherent decreasing trend from 1960 to 1998, thereby suggesting a weakening SCS summer monsoon [15][16][17].  From the perspective of global change, the sea surface wind showed significant spatial and temporal variations. Ward et al. [18] found no clear trend of sea surface wind speed in the second half of the 20th century but found a decreasing trend in the equatorial waters, South Atlantic tropical waters and North Pacific tropical waters. Gulev et al. [19,20] found a strong increasing trend in wind speed in the middle latitudes of the North Atlantic from 1964 to 1993 but found no significant trend in the Northeast Atlantic. After analyzing the buoy data provided by the Meteorological Service of Canada and the National Data Buoy Center of the United States, Gower [21] revealed an increasing trend in the sea surface wind speed in the high latitude area of the North Pacific from 1977 to 1999. Using satellite assimilation data from 1991 to 2008, Young et al. [22] studied the trend of global sea surface wind speed and found an increasing trend in many areas worldwide, with the increasing trend in extreme wind speed being more significant. Islek et al. [6] used 40 years (1979-2019) of ERA-Interim and CFSR data to analyze the temporal and spatial variation characteristics of the wind field in the Black Sea. Their results showed that the wind fields in the western and northern Black Sea were stable and continuous, whereas those in the eastern Black Sea were complex and volatile owing to climate change.
Most studies have focused on the temporal and spatial variations of sea surface wind and the dynamic mechanisms of the SCS monsoon outbreak. However, the long-term trends of wind speed and direction in the northern SCS, especially in the Pearl River Estuary (PRE) and the coastal upwelling area, are unclear. Some physical processes, such as coastal upwelling, river plume transport and estuarine circulation, are controlled by local wind forcing. The long-term trends of sea surface wind in the local hotspot areas should be thoroughly investigated. Changes in the wind field due to climate change could significantly affect the estuarine ecosystem. Wilson et al. [23] indicated that climatic processes, namely the long-term changes in both the direction and directional constancy of summertime winds, control the ventilation of bottom waters and thus the seasonal development of hypoxia in the western Long Island Sound. Scully [24] also indicated that the decadal-scale variability of climate forcing caused the change in wind direction over the Chesapeake Bay, which was responsible for the increased hypoxic volume from the early 1980s to the present. Assessing the trends of sea surface wind is important for understanding the response of the marine environment to climate change.
In this study, we used the latest European Centre for Medium-Range Weather Forecasts (ECMWF) fifth-generation reanalysis (ERA5) data to analyze the long-term trends of the wind field in the northern SCS over the past 41 years (1979-2019). The results of this study will help to reveal the long-term trends of sea surface wind and their effect on the local marine environment in the northern SCS against the background of climate change.
The remainder of this paper is organized as follows: Section 2 introduces the data and methods. Section 3 presents the results and discussion. Section 4 summarizes the conclusions.

Study Area
The northern SCS has a coastline in the northeast-southwest direction ( Figure 1). It is connected to the Pacific Ocean and the East China Sea through the Luzon Strait and the Taiwan Strait, respectively. The wind field in the northern SCS has strong regional and seasonal characteristics. Wind forcing is an important factor affecting the marine environment. As the most important estuary in the northern SCS, the PRE discharges large amounts of fresh water, nutrients, sediment and terrestrial pollutants into the northern SCS [25][26][27][28]. The numerical model results obtained by Pan et al. [29] indicated that the wind direction affects the generation of plumes in the PRE. Therefore, the trend of the wind field in the PRE was the main focus of this study. In addition, there is significant upwelling along the east coast of Hainan Island and eastern Guangdong [30][31][32][33][34]. Coastal upwelling can transport nutrients to the sea surface and promote the growth of phytoplankton, zooplankton and other organisms, thereby forming a better fishing ground. Studies have indicated that the summer southwesterly monsoon is the main factor affecting the upwelling intensity along the northern coast of the SCS [30,31]. Analyzing the wind trends in the coastal upwelling area was another focus of this study.

Data
The trends of sea surface wind in the northern SCS were investigated using the latest ERA5 hourly and monthly wind data at 10 m above the sea surface. The SST and significant wave height data from ERA5 were also used to evaluate the responses of the marine environment to the changes in the wind regime. ERA5 is the global reanalysis product produced by ECMWF. The ECMWF reanalysis began in 1979, followed by the production of ERA15 in the mid-1990s, ERA40 from 2001 to 2003, and ERA-Interim from 2006 to 2019 [35]. ERA5 is based on a hybrid incremental 4D-Var system [36], which uses observations from over 200 satellite instruments or types of conventional data. The spatial and temporal resolutions of ERA5 are 0.25 • × 0.25 • and 1 h, respectively. ERA5 has higher spatial and temporal resolutions than the previous four generations of datasets.

Mann-Kendall Test
The Mann-Kendall (M-K) test was used to detect the change point of the wind field. The M-K test is a climatic change-point analysis method that is widely used in meteorological and hydrological time series [37][38][39][40]. It is a nonparametric test that requires the data to be independent and to tolerate outliers. UF and UB are two M-K statistics used to indicate the trend and the change point in the data. Given a time series x 1 , x 2 , . . . , x n , where: UF k is calculated as follows: when k = 1 and UF 1 = 0. UB k is calculated with the reverse sequence. The change point of the time series is defined as the point where the curves of UF and UB intersect within the confidence interval. UF and UB obey the normal distribution. With the significance level α = 0.05, the critical value of UF and UB is ±1.96 [41]. A negative UF value indicates a decreasing trend, whereas a positive UF value indicates an increasing trend. The trend changes significantly if UF and UB exceed the critical value.

Cross-Wavelet Transform
This study investigated the correlation between the wind field and the El Niño-Southern Oscillation (ENSO). The Oceanic Niño Index (ONI) provided by the National Oceanic and Atmospheric Administration is used to indicate the intensity of ENSO. Two major steps of our analysis were (1) empirical orthogonal function (EOF) analysis and (2) cross-wavelet transform (XWT) analysis. First, the major spatiotemporal patterns of the wind field were identified using the EOF method, which were used to extract significant space-time signals from the wind data. XWT was then used on the EOF time coefficients and ONI to determine the correlation between the wind field and ENSO.
XWT analysis is a signal analysis method that combines cross-spectrum analysis with wavelet transform [42]. It can effectively reflect the time-frequency correlations between two time series and reveal their resonance period and phase relationship at multiple timescales [43]. Given two time series, X i and Y i , XWT is calculated as follows: where W y * i (s) is the complex conjugate of W y i (s) and s is the wavelet scale. The crosswavelet power is defined as W XY i , which represents the power resonance information of X n and Y n in the time-frequency domain. A higher power value represents a stronger correlation between the high energy region of X n and Y n .
Wavelet coherence (WTC) is based on the continuous wavelet transform and XWT of two time series to reflect the time-frequency correlation degree. It is calculated as follows: where S is the smooth operator. In this study, the Morlet wavelet was selected as the wavelet function. Morlet wavelets have non-orthogonality, provide a good balance between time and frequency localization and contain more vibration information [44].

Spatial Variation of the Wind Field in the Northern SCS
The annual mean wind fields averaged from 1979 to 2019 are shown in Figure 2a. Owing to the morphologic effect of the Taiwan Strait and the Luzon Strait, the wind speed in the northeastern SCS was significantly stronger than that in other areas [45]. The wind speed decreased gradually from the northeast to the southwest. The northeasterly wind occupied the eastern portion of the northern SCS, whereas the easterly wind appeared in the western portion, especially in the offshore area of western Guangdong Province and eastern Hainan Island.
The trend of the annual mean wind speed over the past 41 years in the northern SCS was calculated using least-squares regression. The trend of the annual mean wind speed ( Figure 2b) showed significant spatial variations. A clear decreasing trend was observed in the coastal area. In the PRE and its adjacent sea, the maximum rate of decrease was −0.015 m·s −1 ·y −1 . A slightly weaker rate of decrease of approximately 0.01 m·s −1 ·y −1 appeared over the coastal waters of the Leizhou Peninsula and eastern Guangdong. The rate of decrease southwest of Hainan Island was also −0.015 m·s −1 ·y −1 . The decreasing trend in wind speed gradually weakened from the coastal area to the open sea. An increasing trend gradually appeared around 18 • N-19 • N, and the increasing trend offshore of southeast Hainan Island was stronger than that in other areas.
To assess the long-term changes in the annual mean wind direction, the trends of wind vectors were obtained by calculating the trends of the zonal wind component (u) and the meridional wind component (v). The resulting trend vectors represent the change in wind direction and also help to illustrate the contribution of each wind component to the long-term trends of wind speed. As shown in Figure 2a, the annual mean wind was generally northeasterly or easterly (u > 0). Thus, the trends of wind vectors shown in Figure 2b indicate that the zonal wind component decreased, which meant that the easterly wind component weakened over the past 41 years. The weakening of the easterly wind was more notable in the coastal and continental shelf areas than in the open sea. However, the meridional wind increased, which meant that the northerly wind component strengthened. The strengthening of northerly wind was more significant in the open sea than in the coastal area. Thus, the decreasing trend in wind speed in the coastal sea was due to the weakening of easterly wind, whereas the increasing trend in wind in the open sea was due to the strengthening of the northerly wind.
The monthly mean wind fields are shown in Figure 3a. Dominated by the East Asian monsoon, the sea surface wind over the northern SCS showed clear seasonal variations.
A strong northeasterly wind appeared in the winter months (December, January and February) and a strong southerly/southwesterly wind appeared in the summer months (June, July and August). The spring and autumn months are the monsoon transition period. The trends of the monthly mean wind fields showed remarkable seasonal variations (Figure 3b). The wind speed over the continental shelf area showed an increasing trend in December and January and a decreasing trend in February. The trends in the spring (March to May) and summer (June to August) months showed complex spatiotemporal variations. In March, the wind speed showed an increasing trend in almost the entire northern SCS. The highest growth rate of 0.035 m·s −1 ·y −1 appeared in the southeast corner. The wind vectors showed a remarkable increasing trend in the northerly wind component, indicating the strengthening of the northerly wind. In April and May, the wind speed along the east coast of Guangdong showed a significant decreasing trend, the maximum rate of which reached −0.035 m·s −1 ·y −1 . In June, the wind speed showed an increasing trend in the entire northern SCS. According to the wind vectors, both the meridional and zonal wind components showed increasing wind, indicating the strengthening of the southwesterly wind in June. In July and August, the wind speed in the continental shelf showed a decreasing trend. The decreasing trend in July was more significant (with a rate of −0.025 m·s −1 ·y −1 ) in the coastal upwelling area of eastern Guangdong and eastern Hainan Island. Such long-term trends could cause changes in the coastal upwelling process, which in turn could affect fishing in the northern SCS. The trends of the wind vectors indicated a decrease in the meridional wind, suggesting the weakening of the southerly wind in July and August. From September to November, the wind speed showed a decreasing trend over a large domain. The weakening of the easterly wind contributed to the decrease in wind speed in September and October, and the weakening of the northerly wind contributed to the decrease in wind speed in November.

Wind Field Trends in Typical Areas
According to the spatial distribution of the wind speed trend shown in Figure 2b, the decreasing trends in the PRE and the adjacent sea and coastal waters of the Leizhou Peninsula and eastern Guangdong were more significant than those in other areas. These hotspots were selected as representative areas to conduct further analyses. To verify the reliability of our analysis, data recorded at Macau Airport (22.22 • N, 113.50 • E) during 1997-2015 were compared with ERA5 data at station A (located in the lower PRE, close to the Macao Airport). The monthly mean results were plotted in Figure 4. It can be seen that both the wind direction and speed are in good agreement between ERA5 data and airport observations. The annual mean wind speed and its trend from these two data sources also showed good consistency in both magnitude and trend (as shown in Figure 5a, subplot for station A). Therefore, the ERA5 dataset used in this study is robust in the coastal area. The annual mean wind speed at each station is shown in Figure 5a. The annual mean wind speed at each station showed a decreasing trend accompanied by significant interannual oscillation. The rates of decrease at stations A, B, C and D were −0.0106, −0.0130, −0.0102 and −0.0068 m·s −1 ·y −1 , respectively. The PRE and its adjacent water had the greatest rate of decrease in the annual mean wind field. Wilson et al. (2008) indicated that the weakening of the wind field will weaken the mixing in the estuary, thereby enhancing water stratification. Such a condition prevents the transport of dissolved oxygen to the bottom of the estuary, thereby making the estuary more prone to hypoxic events. Hypoxia has frequently been reported during the summer in the PRE [26,46]. Our analyses also indicated a notable decreasing trend in wind speed in July in the PRE and adjacent waters (Figure 3b). The weakening wind in July is assumed to contribute to hypoxia in the PRE. Further research is needed to reveal the related dynamic mechanisms.  The quartile values (75%) of the daily mean wind speeds at stations A, B, C and D over 41 years were used as the strong wind standard at each station ( Table 1). The strong wind days at each station were calculated as the days in each year with a daily mean wind speed greater than the standard value. The results are shown in Figure 5b. The strong wind days also showed decreasing trends at stations A, B and C. There was no significant trend at station D. The rates of decrease in strong wind days are listed in Table 1. The highest rate of decrease was at station A, which is located in the PRE. Strong winds play an important role in modulating the synoptic events in coastal hydrodynamic processes. The decrease in strong wind days in the estuary may significantly affect vertical mixing, lateral circulation, nutrient retention and sediment transport. Changes in the wind field must be taken into account when examining the long-term changes in coastal circulation and related environmental processes.

M-K Test of the Wind Fields in Representative Areas
M-K mutation analyses were conducted using the monthly mean wind speed at stations A, B, C and D. The algorithms for UF and UB were provided in Section 3, and the results were shown in Figure 6. At station A, the UF curve fluctuated from 1979 to 2002 and, since 2003, UF has been less than 0, indicating that the wind speed continuously decreased. UF passed the significance test after 2010, indicating that the decreasing trend was more significant after 2010. The UF and UB curves intersect around 2004, indicating that the wind speed change occurred around this time. The annual mean wind speed was 2.07 m·s −1 before the change point. After the change point, the wind speed decreased by 13% to 1.81 m·s −1 . At station B, UF < 0 after 2004, indicating a decrease in wind speed. Similar to that at station A, the change point at station B also occurred around 2004. At station C, the wind speed decreased continuously after 1986 (UF < 0), and the change point occurred around 1994. There was no clear change point of wind speed at station D. Further analyses of the wind field change before and after the change point were conducted using a wind rose diagram. The results for stations A, B and C are shown in Figure 7. The dominant wind directions at each station showed minor changes before and after the change point. Decreased easterly wind and increased northeasterly wind were discerned. However, the comparison of the frequency of strong winds (wind speed > 8 m·s −1 ) before and after the change point indicated that the frequency of strong winds after the change point was lower than that before the change point. The strong wind frequency at stations A, B and C decreased by 1.25%, 4.10% and 2.86%, respectively. These results indicate that climate change has a significant effect on the wind speed and frequency of strong wind in the northern coastal area of the SCS but not on the wind direction.

Correlation between the Interannual Variation of the Wind Field and ENSO
The general 3-6 year periodicity of the SCS monsoon tends to match that of ENSO [1]. ENSO has a significant effect on the Pacific subtropical high and the Southeast Asian monsoon circulation. The EOF analysis decomposes the spatiotemporal variations of the geophysical field into a combination of orthogonal spatial patterns with corresponding principal components (PCs) in a linear fashion [47]. It serves as a powerful tool for identifying major modes of climate variability such as ENSO [48]. Figure 8 shows the EOF analysis of the normalized monthly mean sea surface wind speed in the northern SCS from 1979 to 2019. The leading mode accounted for 65.8% of the total variance, representing the spatial pattern corresponding to the distribution of the mean wind speed, which had a clear land-sea difference. The second and third modes accounted for 12.31% and 8.77% of the total variance, respectively. The spatial pattern of the second eigenvector reflects the phase difference between the coastal sea and open sea, whereas the spatial pattern of the third eigenvector reflects the phase difference between the northeastern and southwestern portions of the northern SCS. The wind speed in the eastern Guangdong coastal sea and the eastern Hainan coastal sea showed the opposite phase pattern.  The wind speed in the northern SCS exhibits large year-to-year variations, which can be clearly seen from the PCs (Figure 8). Previous studies indicated that ENSO was the dominant factor that controls the interannual variation of the SCS monsoon [1,49]. In order to identify the dominant frequency of northern SCS wind speed variations that are related to the ENSO, the corresponding PCs of each mode were used to conduct the correlation analysis with ONI using XWT and WTC in the time and frequency domains. Generally, the results of XWT revealed the frequency of wind speed change in the northern SCS that could be accounted for by ENSO, whereas the results of WTC identified the period during which the wind speed variations were significantly correlated with ENSO. For the first mode (PC1), the result of XWT (Figure 9a) indicated that the wind speed and ENSO had significant common power in the frequency of the 2-6 year band from 1984 to 2002 and in the 1-4 year band after 2002. The PC1 and ENSO showed significant negative correlation in the period of 3-6 years during the 41 years (Figure 9b). This is consistent with the results of Wang et al. [1]. The XWT of the second mode (PC2) and ONI showed significant common power in the frequency of the 2-6 year band from 1981 to 1993 with the wind field change ahead of ENSO for approximately 3/8 periods. After 1993, the wind field change was ahead of ENSO for approximately 1/8 periods (Figure 9c). The corresponding WTC analysis (Figure 9d

Long-Term Trends of Wind Stress Curl and Effects on Coastal Upwelling
Jing et al. [30] indicated that both the alongshore wind stress and the wind stress curl play significant roles in coastal upwelling in the eastern Hainan and eastern Guangdong coastal seas. Coastal upwelling is a ubiquitous, monsoon-driven physical process that occurs along the continental shelf of the northwestern SCS [31]. It is important to examine the long-term trends of wind stress curl in the northern SCS to assess its effect on coastal upwelling. The spatial distribution of the summer wind stress curl, averaged from 1979 to 2019, is shown in Figure 10a. The wind stress curl was positive in both the eastern Hainan and eastern Guangdong coastal seas, in which typical coastal upwelling exists. The maximum wind stress curl was 3 × 10 −7 N·× m −3 in eastern Hainan. The high wind stress curl is expected to induce strong Ekman pumping in this area, which is conducive to upwelling. The trend of summer wind stress curl, shown in Figure 10b, indicated a strong increasing trend (with a rate of approximately 1 × 10 −9 N·× m −3 ·y −1 ) in the coastal area, especially in the eastern Hainan and eastern Guangdong coastal seas. Water that rises to the surface as a consequence of upwelling is usually colder than the surface water. The summer mean SST averaged from 1979 to 2019 is shown in Figure 10c. Two typical coastal upwelling regions can be discerned. One is the eastern Hainan coastal upwelling area, where the SST was approximately 28.0 • C-approximately 1.5 • C colder than the surrounding water. The other one is the eastern Guangdong coastal upwelling area, where the SST was approximately 27 • C-approximately 1 • C colder than the surrounding water. Figure 10d presents the long-term trends of SST. The summer SST in the northern SCS showed a warming trend against the background of global warming. Spatial variations can be found in the SST trend. The warming trend was stronger in the Guangdong coastal sea. The maximum warming rate was 0.35 • C·y −1 . The rates of increase in SST in the two typical upwelling regions were lower than those in other areas, which indicated that the SST difference between the upwelling area and the ambient area increased. The coastal upwelling strengthened from 1979 to 2019.
Coastal upwelling is an important driving force of the spatial distribution of zooplankton abundance [50,51]. It occurs when alongshore wind stress produces offshore Ekman transport [= τ alongshore /(ρ f ), where τ alongshore is the alongshore component of surface wind stress, ρ is the seawater density and f is the Coriolis parameter] or when the wind stress curl produces upward motion at the bottom of the Ekman layer. Figure 10 indicated that the strengthening of coastal upwelling in the two typical coastal upwelling regions of the northern SCS was in good agreement with the increasing trend of the positive wind stress curl at the same sites. Although the increasing trend of wind stress curl can also be found around the PRE and its adjacent coastal sea, the corresponding wind stress curl is not positive and there are no obvious upwelling processes. Wind stress, which is proportional to the Ekman transport [52], showed a decreasing trend in July and August along the coast of the northern SCS (Figure 3b). Our results indicated that the strengthening of the positive wind stress curl was responsible for the strengthening of coastal upwelling at eastern Hainan Island. Further numerical model research is needed to evaluate the relative importance of long-term changes in alongshore wind stress and wind stress curl to the long-term changes in coastal upwelling in eastern Guangdong.

Effect of the Trend of Wind Stress on Significant Wave Height
The spatial distribution of the annual mean significant wave height (averaged from 1979 to 2019) in the northern SCS is shown in Figure 11a. The spatial distribution of the significant wave height was closely related to the spatial variation of the seabed topography. The wave height in the coastal shallow waters was lower than that in the offshore deep waters owing to energy dissipation in the shallow area. The spatial distribution of the long-term trends of annual mean wave height (Figure 11b) showed an increasing trend in the open sea, which gradually weakened towards the coast of the northern SCS. A decreasing trend was seen in the coastal area, consistent with the decreasing trend in wind speed shown in Figure 2b. Waves in the SCS are influenced by the East Asian monsoon and thus show a significant seasonal difference [53,54]. The spatial distributions of seasonal mean wind speed ( Figure 12a) and significant wave height (Figure 12b) showed similar seasonal variations. The high wave height in autumn and winter was consistent with the high wind speed during these seasons. The long-term trends of the seasonal mean wind speed (Figure 12c) showed a similar pattern to that of the seasonal mean significant wave height (Figure 12d). In spring and summer, the wave height showed a decreasing trend in the coastal area and an increasing trend in the offshore area. The decreasing trend was present in almost the entire northern SCS in autumn. In winter, a significant increasing trend was observed throughout the domain. Similar spatial patterns were found for the long-term trend of seasonal mean wind speed, which plays an important role in the change in wave height, which in turn may influence ocean engineering in the long term. The long-term reduction in wave height could also be related to the oceanic conditions [55] (Shimura et al., 2016) and extreme weather [4]. Studies in other oceans/seas have also indicated that the long-term trend of wave height showed clear spatial [56] (Aydogan and Ayat, 2018) and seasonal [57] (Zheng and Li, 2015) variations. Wind and wave energy are important for renewable energy, ocean engineering and coastal protection. Long-term changes in wind and waves also affect the coastal topography and the shape of beaches. Their variability and responses to climate change require further exploration.

Conclusions
In this study, the latest ERA5 data were used to investigate the long-term trends of the wind field in the northern SCS over the past 41 years (1979-2019). The main conclusions are as follows.
The long-term trends of the wind field in the northern SCS showed significant spatial and seasonal variations. The wind speed in the coastal area has mainly decreased, whereas that in the open sea has mainly increased. The PRE and its adjacent sea, the coastal waters of the Leizhou Peninsula and the coastal waters of eastern Guangdong are three representative areas with significant decreasing trends in annual mean wind speed and strong wind days. The easterly wind component showed a significant weakening trend in the coastal and continental shelf area, whereas the northerly wind component showed an increasing trend in the open sea. The M-K mutation analysis of the wind field in the representative area indicated significant changes in the wind speed and frequency of strong winds in the coastal area of the northern SCS. The changes in the wind direction are minor.
The wind field in the northern SCS was significantly correlated with ENSO. The PCs of the dominant EOF modes and ONI were analyzed statistically using XWT and WTC. We found a significant negative correlation between the first mode PC and ENSO in the 3-6 year band. The second mode PC had a weaker correlation with ENSO, whereas both the phase and period were different at different timescales. A significant positive correlation was observed between the third mode PC and ENSO in the 4-6 year band.
The effect of climate change on the wind field has a significant effect on marine environmental factors. During 1979-2019, the upwelling along the east coast of Hainan Island and Guangdong showed an increasing trend. The strengthening of the wind stress curl was a major cause of the increase in coastal upwelling. The long-term trends of the significant wave height and wind showed similar spatial and seasonal variations. The long-term trends of the wind field play an important role in modulating the effect of climate change on wave height in the northern SCS.  Data Availability Statement: Data available on request due to restrictions. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.