Changes in Means and Extreme Events of Sea Surface Temperature in the East China Seas Based on Satellite Data from 1982 to 2017

Marginal seas are fundamental to humans for their importance in mariculture resources and commerce. Based on the NOAA 0.25 degree daily Optimum Interpolation (OI) sea surface temperature (SST) data set, spatiotemporal changes in mean and extreme SST in the East China Seas (ECSs) were examined for from 1982 to 2017. As a regional average, the annual mean SST has notably increased at a rate of 0.21 ± 0.08 ◦C per decade. The warming SST during 1982–2017 is probably related to the influence from a recent strengthening and westward extension of the WPSH. There are also notable warming trends in annual minimum and maximum SST. Spatially, the rapid warming of annual mean SSTs are located in the vicinity of the Yangtze Estuary, exceeding 0.2 ◦C per decade and part of the ECS-Kuroshio. This pattern may be largely affected by the spatial changes of minimum SST. Rapid warming of maximum SST can be found across the region, from the northern East China Sea (ECS) to the Bohai Sea. Since 1982, extreme hot days (EHDs) have undergone an obvious increasing trend, at a rate of 15.2 days per decade. Conversely, extreme cold days (ECDs) have been decreasing. Notably, the largest increase of EHDs appears in the western ECS and the Bohai Sea, which both have rich marine ecosystems. The trend of EHDs has a significant relationship to mean SST, suggesting that there will be a further increase in EHDs under continued warming in the ECSs. These findings emphasize the importance and urgency of strategies which should be planned for the adaptation and mitigation of specific types of extreme hot events in this region.


Introduction
Sea surface temperature (SST) is considered to be one of the most important indicators in quantifying climate change [1].A rapid warming is expected to have a large impact on regional marine biodiversity, both at the ecosystem level and at the population level [2].The Fifth Assessment Report of the Intergovernmental Panel on Climate Change affirmed that our warming of the climate system is unequivocal and its extremes are also changing [3].To better adapt to climate warming and decision-making, an improved and detailed understanding of regional climate changes is required.
Part of the Northwest Pacific Ocean, the Bohai Sea, Yellow Sea, and East China Sea (ECSs) is located at one of the largest continental shelves in the world, with a large number of estuaries, gulfs, freshwater inputs, and maricultures.There have been much literature focused on the SST variability over the China Seas [4][5][6][7][8].They found that the mean warming in the China Seas has been considerably faster than the global average.In addition, the ECSs are vulnerable to climate change, especially changes in climate extremes [4,6,9].However, the previous studies covering the China Seas typically used reconstructed or reanalyzed monthly mean SST data sets.The detailed extreme hot and cold SST changes and their relationship with changes in SST annual ranges (ARs) remain unknown, owing to the lower spatial and temporal resolution of data sets.At present, extreme water temperature events can cause shifts in species ranges, local extinctions, and economic loss on aquaculture and seafood industries [9][10][11].For instance, the rapid warming of the ECSs has been linked to frequent occurrences of harmful algae blooms and northward shifts of marine fishes [12].A Chinese Marine Disaster Bulletin from 2012 [13] indicated that severely harmful algae occurred at the coast of the southern ECS, causing a direct economic loss of two billion RMB.Changes in the SST at shorter than seasonal scales and variations in the frequency of SST extreme events in sub-region scales should be taken into account when planning for strategies to mitigate these changes.Thus, the present study adds some new insights into understanding the inter-annual variations of annual mean, annual minimum, annual maximum SSTs, and AR averaged of the ECSs and spatial distribution of trends in the ECSs from 1982 to 2017.Meanwhile, extreme hot and cold events are also detected in detail during the study period.Finally, we will discuss the possible preliminary impacts on marine fisheries.The information through our work will improve the understanding of the hydrography and climate trends in this region and benefit fisheries and marine biological research.
The remainder of the paper is organized as follows.Following the introduction, Section 2 briefly introduces the study region, data sources, and methods.Section 3 describes the changes of mean SST and extreme SST.Changes in extreme events are conducted in the following section.A summary and discussion are given in Section 5.

Study Region
The study region including the Bohai Sea, the Yellow Sea, and the East China Sea (ECS) is bounded by China, the Ryukyu (Nansei) Islands, Kyushu, and the Korean Peninsula (Figure 1).The study region also covered a vicinity area in order to show a better perspective for the SST variability, such as the SST variability in the ECS-Kuroshio.Thus, we approximated the study region in a range of 116-132

Data Source
A daily SST data set was acquired from the NOAA National Climatic Data Center (NCDC),

Data Source
A daily SST data set was acquired from the NOAA National Climatic Data Center (NCDC), available at https://www.ncdc.noaa.gov/oisstwith a high spatial resolution of 0.25 • × 0.25 • from 1982-2017 [14].This data set is derived from the Optimum Interpolation (OI) SST Analysis Product, referred herein as OISST v2, which uses SST data from an advanced very high-resolution radiometer infrared satellite from the Pathfinder satellite combined with buoy data, ship data, and sea ice data SST data sets.In order to apply the correction for bias in OISST, the satellite data have been classified into daytime and nighttime bins and corrected separately using the patterns of 15 days averaged by in situ SSTs using NOAA's OI algorithm.The bias-corrected daytime and nighttime satellite SST, ship, and buoy SSTs are merged based on noise-to-signal ratio maps for each data type, which have averaged weights of 15.1, 15.1, 1.0, and 15.1, respectively.Therefore, it can be interpreted as the bulk SST at about 0.5 m depth [15].
In this paper, we use the Taylor diagram to provide a statistical summary of how well the OISST v2 represents the 19 stations-based observations along the coastal ECSs.In order to do the comparison and evaluation, OISST v2 has been interpolated to the stations' locations.The correlation and standard deviation on the annual basis are computed and shown in a Taylor diagram (Figure 2).As seen in Figure 2, OISST v2 are clustered around the reference values (REF), except for 4 points (points 3, 13, 14, and 18).Most of them have positive correlations with observations in the range of 0.5~0.9.A total of 79% of these points has a lower standard deviation and smaller root mean square error, which suggests that OISST and in situ water temperatures in the coast area of China agree reasonably well.Additional analysis (not shown here) indicated that there is no statistically significant trend in the differences between the long-term water temperature records.OISST is currently the longest satellite data record that can be used to study long-term SST variability.This data set has been widely used in scientific researches [16][17][18][19].The suitability of the data base to identify extreme events was previously shown in work of Lima and Wethdy and Liao et al. [16,20].

Annual Frequency of Extremely Cold and Hot Days
Extreme days were carried out for each pixel as follows: daily SST anomalies were calculated to remove the seasonal cycle.This was achieved by subtracting from the SST of a certain day (e.g., January 1, 1982) the mean temperature of that day (January 1) of the climatic baseline period of 1982-2011 [17,21].Therefore, an extreme day can occur at any time of the year.An EHD (or ECD) is defined as a day with its SST anomaly above (or below for ECDs) the 90th (or 10th for ECDs) percentile of standardized anomaly of the climatic baseline period.Then, the number of EHDs (or

Annual Frequency of Extremely Cold and Hot Days
Extreme days were carried out for each pixel as follows: daily SST anomalies were calculated to remove the seasonal cycle.This was achieved by subtracting from the SST of a certain day (e.g., January 1, 1982) the mean temperature of that day (January 1) of the climatic baseline period of 1982-2011 [17,21].Therefore, an extreme day can occur at any time of the year.An EHD (or ECD) is defined as a day with its SST anomaly above (or below for ECDs) the 90th (or 10th for ECDs) percentile of standardized anomaly of the climatic baseline period.Then, the number of EHDs (or ECDs) was calculated as the number of days per year with extremely high (low for ECDs) SST anomalies.

Testing for Significance of Correlations and Trends
When test the significance of correlation between two time series and the significance of the presence of trends in a time series, an important question requires consideration: how large sample correlations and sample trends need to be, even if the stochastic processes, which generate the series, are not correlated at all and exhibit no trends.Firstly, we have to make an assumption, namely the processes X and Y share no correlation, or segments of length L of the process have no trend.Standard procedures are available in the literature, namely p values for correlations and Mann-Kendall for trends [22,23] when there are "no correlations" between the underlying processes and trends can hardly appear in limited segments of an infinite stationary time series.
In the case of correlations, the assumption is that the underlying processes are stationary (free of systematic trends) and serially independent, that is, X t and X t+1 for any t are independent.In the case of trends, the assumption is the independence of X t .However, in geophysical cases, these assumptions are not satisfied-the result is that the null hypotheses are more often falsely rejected (i.e., in cases where there are no correlations or no trends [22] than stipulated by the significance level (normally 5%).
A practical remedy for avoiding such errors is to deal with normalized series (mean = 0, standard deviation = 1) X t (and Y t ) as follows: (1) "detrend" the time series before testing for correlations between two time series X t and Y t .Firstly, determining the linear fit f X t and f Y t , and then do the hypothesis testing with (2) "prewhiten" the time series, by first determining the sample autocorrelation α = 1/L∑ t X t X t+1 of the time series X t of length L, and forming a series X t = X t − αX t , and then testing for the null hypothesis of no trend.
For both cases, the standard routines are applied.If the null hypothesis is rejected at the stipulated significance level of 5%, then the sample trend f X t , or the sample correlation 1/L∑ t X t X t+1 , is considered "significant".

Annual SST Changes
Figure 3 shows the time series of regional averaged annual mean SST, maximum SST and minimum SST anomalies in the ECS from 1982 to 2017, as well as AR SST, which was defined as the difference between the annual maximum SST and minimum SST.The annual mean SST increased at a rate of 0.21 ± 0.08 • C per decade, significant at the 95% level (p < 0.05) (Figure 3a).The warming trend is higher than the global average SST warming rate of 0.15 ± 0.03 • C per decade from 1982 to 2017.Moreover, we found that the largest annual mean SST anomaly appeared in 1998, followed by the second largest SST anomaly in 2016.Both of the 1997/1998 and 2015/2016 El Niño events are by far the strongest El Niño events since the beginning of the 20th century [24,25].The SST peaks were probably caused by the low-frequency El Niño/Southern Oscillation (ENSO) variabilities [26,27].When ENSO is closely related to the East Asian Monsoon [28,29], it can lead to weaker than normal East Asian Winter Monsoon (EAWM) and a stronger west Pacific subtropical high (WPSH) [8,[30][31][32], and also favors an abnormally high SST in the marginal seas.
The general warming may have been mainly caused by anthropogenic increase in atmospheric CO 2 concentration [33], and partly by the urbanization effect around the observational stations.In this subsection, however, we further analyzed the possible relationship between the warming and the large-scale atmosphere circulation modes.As shown in Figure 4a-c [34,35].During this study period, the increasing rate of the annual minimum SST anomalies is 0.21 ± 0.18 • C per decade, significant at the 95% level (p < 0.05) (Figure 3b).The annual maximum SST increased at a rate of 0.20 ± 0.13 • C per decade, significant at the 95% level (p < 0.05) (Figure 3c).However, during the studied period, regional averaged AR SST showed an insignificant decreasing range trend of −0.02 ± 0.02 • C per decade (Figure 3d), which may due to the similar warming amplitudes of annual maximum and minimum SST.
Atmosphere 2018, 9, x FOR PEER REVIEW 6 of 13 range trend of -0.02 ± 0.02 °C per decade (Figure 3d), which may due to the similar warming amplitudes of annual maximum and minimum SST.Spatially, the significant high warming rates of the annual mean SST appear in the western ECS and part of the ECS-Kuroshio, exceeding 0.2 °C per decade (Figure 5a).In them, warming rates of the Yangtze River estuary are especially large, up to 0.4 °C per decade.The rapid warming in the western ECS coincides with the previous findings [7,8].Studies pointed that the stream temperature in the vicinity of the Yangtze Estuary had increased by ~2 °C since 1986, thereby contributing to the extremely rapid warming observed in the western ECS [4,36].However, the other significant warming at the ECS-Kuroshio has not been captured using the coarse SST grid dataset in the previous studies.There are weak cooling trends occurs in the west of the Korea Island in Figure 5a which were not clearly shown in the previous studies and needs a further investigation.
We also estimate the trends distribution of annual minimum and maximum SST in the study region (Figure 5b and 5c).Comparing to the rapid warming of minimum SST in the range of the ECS-Kuroshio and the Yangtze River estuary, the maximum SST shows clear increasing in the most part from the northern ECS to the Bohai Sea.The inter-annual SST variability in the ECSs is not only Spatially, the significant high warming rates of the annual mean SST appear in the western ECS and part of the ECS-Kuroshio, exceeding 0.2 • C per decade (Figure 5a).In them, warming rates of the Yangtze River estuary are especially large, up to 0.4 • C per decade.The rapid warming in the western ECS coincides with the previous findings [7,8].Studies pointed that the stream temperature in the vicinity of the Yangtze Estuary had increased by ~2 • C since 1986, thereby contributing to the extremely rapid warming observed in the western ECS [4,36].However, the other significant warming at the ECS-Kuroshio has not been captured using the coarse SST grid dataset in the previous studies.There are weak cooling trends occurs in the west of the Korea Island in Figure 5a which were not clearly shown in the previous studies and needs a further investigation.

Changes in the Extreme Events and Their Relationship with Mean SST
In this section, we examine the changes of extreme SSTs based on the daily OISST v2. Figure 6 exhibits the time series of regional averaged annual EHDs and ECDs in the study region.The number of annual EHDs significantly increased, at a rate of 15.2 days per decade (significant at the 95% level, p < 0.05).In contrast, the number of annual ECDs remarkably decreased, at a rate of 10.5 days per decade (significant at the 95% level, p < 0.05).Annual EHDs and ECDs are positively (0.80) and negatively (-0.92) correlated with annual mean SST, respectively.In general, the increase of We also estimate the trends distribution of annual minimum and maximum SST in the study region (Figure 5b,c).Comparing to the rapid warming of minimum SST in the range of the ECS-Kuroshio and the Yangtze River estuary, the maximum SST shows clear increasing in the most part from the northern ECS to the Bohai Sea.The inter-annual SST variability in the ECSs is not only associated with large-scale climate oscillations in the Pacific (for instance, Pacific Decadal Oscillation, North Pacific Oscillation, East Asia Monsoon, ENSO, etc.) but also affect by many local processes, including the Yangtze River discharge, Kuroshio, and East Asian Monsoon, etc.The complex local processes may lead to regional differences in trends.Wu et al. [6] related the phenomenon of high SST increasing trend along the western Pacific marginal seas to the poleward shift and the intensification of the Kuroshio.In boreal winters, weakening northeasterly winds after 1982 would strengthen the northward advection of the Kuroshio and its cross-shelf currents into the ECS, due to a strong southward upper Ekman drift [31].It is actually coincident with the finding that a rapid warming of annual minimum SST in the range of the ECS-Kuroshio.As shown in Figure 5d, there is an out-of-phase change on both sides of 30 • N (Figure 5d).The AR SST is notably weakened in the band of the ECS-Kuroshio, with the decreasing center located at (127-128 • E, 28-30 • N).The decrease of the AR SST is associated with a rapid warming of the minimum SST and the slow warming of the maximum SST.AR SSTs in the Bohai Sea and the Yellow Sea are increasing, except some much small areas.

Changes in the Extreme Events and Their Relationship with Mean SST
In this section, we examine the changes of extreme SSTs based on the daily OISST v2. Figure 6 exhibits the time series of regional averaged annual EHDs and ECDs in the study region.The number of annual EHDs significantly increased, at a rate of 15.2 days per decade (significant at the 95% level, p < 0.05).In contrast, the number of annual ECDs remarkably decreased, at a rate of 10.5 days per decade (significant at the 95% level, p < 0.05).Annual EHDs and ECDs are positively (0.80) and negatively (-0.92) correlated with annual mean SST, respectively.In general, the increase of EHDs is accompanied by a decrease of ECDs (r = −0.76,p < 0.01).Therefore, the region with higher warming rates tends to have more EHDs and less ECDs.Spatially, significant increases in the frequency of EHDs were observed over the western ECS and the Bohai Sea (near shore areas even exceeding 20 days per decade), in association with a decrease in the frequency of ECDs simultaneously in Figure 7.This high rate is consistent with a rapid warming of annual mean SST in the western ECS.Lima and Wethey [17] studied trends of coastal warm and cold days in the world's coastlines from 1982 to 2010 using OISST v2.They detected an increase in EHDs on the range of 10~20 days per decade and the decrease in ECDs on the range of 5~15 days per decade in the coastline of the China Seas, consistent with our finding of the rapid increase of EHDs in the coastal areas.
Atmosphere 2018, 9, x FOR PEER REVIEW 10 of 13 EHDs is accompanied by a decrease of ECDs (r = -0.76,p < 0.01).Therefore, the region with higher warming rates tends to have more EHDs and less ECDs.Spatially, significant increases in the frequency of EHDs were observed over the western ECS and the Bohai Sea (near shore areas even exceeding 20 days per decade), in association with a decrease in the frequency of ECDs simultaneously in Figure 7.This high rate is consistent with a rapid warming of annual mean SST in the western ECS.Lima and Wethey [17] studied trends of coastal warm and cold days in the world's coastlines from 1982 to 2010 using OISST v2.They detected an increase in EHDs on the range of 10~20 days per decade and the decrease in ECDs on the range of 5~15 days per decade in the coastline of the China Seas, consistent with our finding of the rapid increase of EHDs in the coastal areas.EHDs is accompanied by a decrease of ECDs (r = -0.76,p < 0.01).Therefore, the region with higher warming rates tends to have more EHDs and less ECDs.Spatially, significant increases in the frequency of EHDs were observed over the western ECS and the Bohai Sea (near shore areas even exceeding 20 days per decade), in association with a decrease in the frequency of ECDs simultaneously in Figure 7.This high rate is consistent with a rapid warming of annual mean SST in the western ECS.Lima and Wethey [17] studied trends of coastal warm and cold days in the world's coastlines from 1982 to 2010 using OISST v2.They detected an increase in EHDs on the range of 10~20 days per decade and the decrease in ECDs on the range of 5~15 days per decade in the coastline of the China Seas, consistent with our finding of the rapid increase of EHDs in the coastal areas.

Summary and Discussions
Our main goal in this paper was to characterize the changes in mean and extreme SST in the ECSs.For this purpose, we used the NOAA OISST v2, which has high spatial and temporal resolutions.The analysis of SST variability in the ECSs successfully lead to a better understanding of this marine environment, including climate related hydrographic and ecosystem trends.The main findings of the present study can be summarized as follows: (1) As a regional average, annual mean SSTA over the ECSs is subject to more notable warming trend, with a rate of 0.21 ± 0.08 Significantly higher SSTs above average in some area are implicated in dramatic changes to the physical, chemical and biological state of the marine environment.It is worth noting that our analysis indicates that the warming trend in the ECSs is greater than the trend of the global averaged SST.This finding is consistent with previous results based on different SST data sets [6,8].The increase in SST makes the abundance of warm water species and the decrease of warm-temperate species in the Yangtze Estuary area [37][38][39].Also, it can cause latitudinal shifts in species distributions [39].Moreover, exposure and vulnerability to extreme events can destroy marine fishery assets and infrastructure.For instances, extreme higher water temperature can lead to abnormal metabolism of scallop, a cold-water species [40].In the boreal summer of 2017, Zhangzi Island (122.7 • E and 39 • N) of the Yellow Sea had experienced more than 40 EHDs, with the maximum daily SST anomaly of 4.12 • C.Under such conditions, scallop could not absorb nutrients normally, resulting in malnutrition, bacterial infection and eventually inducing widespread mortality of scallops in this area.
Currently, the frequency and intensity of extreme events are increasing globally as a consequence of anthropogenic and natural climate change influence [3,33].The latest research points out that the human-induced global warming has never stopped and the human influence is dominant in long-term warming [41].In the present study, we found the likelihood of extreme hot events increases from 1982 to 2017 in the ECSs, especially in the sub-regions of the Bohai Sea and the western ECS.Given the continued global warming, there will be more frequency and intense extreme hot events in the study region.Given the information mentioned above, contemporary warming events are possibly superimposed onto the warming trends, increasing the risk of extreme SST events on marine ecosystem.Thus, better targeted strategies should be designed for an adaptation and mitigation against specific types of hot extremes as early as possible.

13 Figure 1 .
Figure 1.Map showing the geographical location of the study region, including the Bohai Sea, the Yellow Sea, the East China Sea, and adjacent seas.

Figure 1 .
Figure 1.Map showing the geographical location of the study region, including the Bohai Sea, the Yellow Sea, the East China Sea, and adjacent seas.

Figure 2 .
Figure 2. Taylor diagram of the annual mean SSTs from OISST v2 compared with 19 stations-based observations, REF stands for reference values of SST observations.The radial distance from the origin is proportional to the standard deviation ratio.The triangles and circle denote the bias between OISST and stations-based observations.

Figure 2 .
Figure 2. Taylor diagram of the annual mean SSTs from OISST v2 compared with 19 stations-based observations, REF stands for reference values of SST observations.The radial distance from the origin is proportional to the standard deviation ratio.The triangles and circle denote the bias between OISST and stations-based observations.
Figure3shows the time series of regional averaged annual mean SST, maximum SST and minimum SST anomalies in the ECS from 1982 to 2017, as well as AR SST, which was defined as the difference between the annual maximum SST and minimum SST.The annual mean SST increased at a rate of 0.21 ± 0.08 • C per decade, significant at the 95% level (p < 0.05) (Figure3a).The warming trend is higher than the global average SST warming rate of 0.15 ± 0.03 • C per decade from 1982 to 2017.Moreover, we found that the largest annual mean SST anomaly appeared in 1998, followed by the second largest SST anomaly in 2016.Both of the 1997/1998 and 2015/2016 El Niño events are by far the strongest El Niño events since the beginning of the 20th century[24,25].The SST peaks were probably caused by the low-frequency El Niño/Southern Oscillation (ENSO) variabilities[26,27].When ENSO is closely related to the East Asian Monsoon[28,29], it can lead to weaker than normal East Asian Winter Monsoon (EAWM) and a stronger west Pacific subtropical high (WPSH)[8,[30][31][32], and also favors an abnormally high SST in the marginal seas.The general warming may have been mainly caused by anthropogenic increase in atmospheric CO 2 concentration[33], and partly by the urbanization effect around the observational stations.In this subsection, however, we further analyzed the possible relationship between the warming and the large-scale atmosphere circulation modes.As shown in Figure4a-c, the WPSH indices of

Figure 3 .
Figure 3.Time series of the regional averaged annual mean SST (a), annual minimum SST (b), annual maximum SST (c) anomalies and AR SST (d) in the ECSs from 1982-2017 (Units: °C, 1982-2011 base period).The black dashed lines are the linear trend lines based on the least squares estimator.

Figure 3 .
Figure 3.Time series of the regional averaged annual mean SST (a), annual minimum SST (b), annual maximum SST (c) anomalies and AR SST (d) in the ECSs from 1982-2017 (Units: • C, 1982-2011 base period).The black dashed lines are the linear trend lines based on the least squares estimator.

Figure 4 .
Figure 4. Time series of summer mean (black lines) and 9 year running mean (shaded lines) indices from 1982-2017 of the west Pacific subtropical high (WPSH): area (a), intensity (b), westward extension (c), and time series of the East Asia winter monsoon Index (d).Data are all from the Chinese National Climate Center (CNCC).

Figure 4 .
Figure 4. Time series of summer mean (black lines) and 9 year running mean (shaded lines) indices from 1982-2017 of the west Pacific subtropical high (WPSH): area (a), intensity (b), westward extension (c), and time series of the East Asia winter monsoon Index (d).Data are all from the Chinese National Climate Center (CNCC).

Figure 5 .
Figure 5. Linear trends spatial distribution of (a) annual mean SST, (b) annual maximum SST, (c) annual minimum SST, (d) AR SST in the ECSs (Units: °C per decade).

Figure 5 .
Figure 5. Linear trends spatial distribution of (a) annual mean SST, (b) annual maximum SST, (c) annual minimum SST, (d) AR SST in the ECSs (Units: • C per decade).

Figure 6 .
Figure 6.Regional averaged annual EHDs (red line) and ECDs (blue line) from 1982 to 2017 in the study region (units: days).

Figure 6 .
Figure 6.Regional averaged annual EHDs (red line) and ECDs (blue line) from 1982 to 2017 in the study region (units: days).

Figure 6 .
Figure 6.Regional averaged annual EHDs (red line) and ECDs (blue line) from 1982 to 2017 in the study region (units: days).
• C per decade which is higher than the global mean trend based on the same data set during 1982-2017.The most significant warming was observed at the western ECS and the Bohai Sea.The warming SST during 1982-2017 is probably related to the influence from the recent strengthening and westward extension of the WPSH.(2)Our results indicate that there are significant warming tendencies in annual minimum SST and maximum SST.Compared to the rapid warming of maximum SST from the northern ECS to the Bohai Sea, the rapid warming of minimum SST mainly covers the western ECS and part of the ECS-Kuroshio.AR SST shows an out-of-phase change on both sides of 30 • N during 1982-2017.(3) Rapid decreases in EHDs and increases in ECDs were detected in the ECSs from 1982 to 2017, at a rate of 15.2 days per decade and 10.5 days per decade, respectively.Spatially, the notable increases in the frequency of EHDs were observed in almost all of the ECSs, especially in the Bohai Sea and the western ECS, exceeding 20 days per decade.