Earth’s response to climate change can be observed from the Sea level (SL) variation, which is one of the important parameters . SL rising is receiving more attention for global researchers as one of the catastrophic consequences of climate change. SL rise is a matter of urgent concern, especially because many low-lying parts of the world undergoes consequent threats    . Since 1980，an increasing trend in SL observed Globally. However, SL rise over China’s coastal sea level is higher than that of the world  . Since 2000, China’s coastal SL continued to be increasing and 2012 reached its highest level since 1980, however 2005 year shows decline . From 2010 to 2012 China’s coastal SL have high periodical oscillation, several high superimposed periodical oscillations had obvious effect on the sea level. In 2012, coastal temperature and SST were higher 0.4˚C. The anomaly of the wind leads to water accumulation is one of the causes of rising sea levels. In 2012, it was the tropical cyclones landing time and had broad influence, especially in August 2012 has six tropical cyclones affected China’s coastal area and had made obvious effect on sea level rise in that month. Subtropical high in 2012, northerly wind has influence in rising in SL over the east China sea and south China . In the context of global warming, sea level rise has become a major global environmental problem. Impact of rising SL on the economic development of coastal areas, marine ecological environment and people’s production and life is increasingly significant.
Global sea level rise is caused by climate warming, swelling of the sea, continental glaciers and melting of polar ice caps. Under the background of global climate change, China’s coastal have higher temperatures and SST and lower pressure, then leading to increase in sea level. SL rise will lead to threat in the coastal cities flood control and water supply system. The change of offshore SL is related to the human living environment and its development prospects, so SL gets more attentions. In recent decades, many experts have made lot of investigations to get consistent understanding. In recent 100 years, the global sea level is on the rise and the average rise rate is about 1.0 - 2.0 mm/a according to previous studies  - . Average annual variability of the sea level along the coast of China is in the same range however, there are significant differences in SL in different regions.
Mean global SL and mean global sea surface temperature (SST) has strong correlation , which suggests that the steric effect impact sea level signature. In addition to global warming factors, in the long-term scale, SL changes are affected by structural changes in the earth, the activities of the great glacial period and the effects of the atmosphere, the ocean itself, etc. Numerous studies have shown that the change in SL is closely related to the water thermal expansion .     . Further this is supported by the strong similarity of the spatial distribution of satellite-observed in sea level variations and surface temperature   . By 2100 rise in temperatures is predicted by 2˚C - 4.4˚C leading to global average sea level rise of 2 - 6 mm per year up to 2100  which impacts the protected coastal habitats. Local SST variations of about 0.58˚C/yr observed and would be adequate for sea level rise. Temperature to sea level relationship may differ from the surface temperature to sea level relationship on this long timescale, which is caused by the ocean circulation    . In this paper we find out the relation between the SL and SST over the Zhejiang province coastal areas.
2. Study Area and Data
Zhejiang is a province with the most islands in China. Zhejiang coastline extending 6486 kilometers and a total ocean area of 220,000 square kilometers. In addition, the province has a large number of bays with over 60 natural ports. Zhejiang is located in the southern part of the Yangtze River Delta on the southeast coast of China, which accounts for sediment transport to sea. The Study area we chosen 26˚N - 32˚N: 119˚E - 124˚E, which in the East China. Further the coastal areas have been chosen to find out the variations of SL over Zhejiang Province given in the Figure 1 and they are as follows: (a) Hangzhou Bay (HZB) 30˚N - 31˚N; 120˚E - 122˚E, (b) Ningbo (NB) 28˚N - 31˚N; 120˚E - 123˚E, and (c) Wenzhou (WZ) 27˚N - 29˚N; 119˚E - 122˚E.
Figure 1. Study area which includes the coastal areas of Zhejiang province.
Sea Surface Height Relative to Geoid over the area has been obtained from the operational ORAS4 data used for this analysis. The European Centre for Medium-Range Weather Forecasting (ECMWF) estimates the state of the global ocean via the operational system Ocean-S4. Ocean-S4 gives an estimate of the history of the ocean from September 1957 to present (with a few days delay) via the Ocean Reanalysis System 4 (ORAS4), as well as the latest ocean conditions, provided by the real-time extension Ocean Real Time Analysis System 4 (ORTA4). In this study, the data over the period of January 1958 to December 2015 monthly SL data used to find the variations in the study areas.
To validate the relation between the SL and SST relation we have used the dataset HadISST (Hadley Centre Global Sea Ice and Sea Surface Temperature) is a combination of monthly globally complete fields of SST and sea ice concentration for 1871-present. HadISST uses reduced space optimal interpolation applied to SSTs from the Marine Data Bank (mainly ship tracks) and ICOADS through 1981 and a blend of in-situ and adjusted satellite-derived SSTs for 1982-onwards. The “bucket correction” was applied to SSTs for 1871-1941. HadISST is primarily intended to be used as boundary conditions for atmospheric models. With above mentioned data sets SL variations computed and presented in following sections and given the relation between SL and SST.
3.1. Monthly Variations of SL
Monthly SL variations over Zhejiang province areas are presented in the Figures 2(a)-(c). One can clearly observe that SL is showing an increasing trend over Zhejiang province areas (HZB, NB, WZ). SL is varying between −0.26 m and 0.29 m over HZB, −0.24 m and 0.28 m over NB, −0.24 m and 0.29 m over WZ respectively. Lower SL monthly variations can be observed in 1961-62, 1967-68. Higher SL anomalies observed in 2013, however over HZB and NB 2014 also showing higher SL. Higher values of SL can be explained as the more rainfall leading to more river discharge and higher SST values over the study area. Monthly plot is revealing there is a seasonal variation, which are discussing in the next section.
3.2. Seasonal Variations of SL at Different Seasons
Figure 3 depicts the seasonal variation SL over the Zhejiang province areas. SL variations are indicating increasing trend in every season. Autumn season SL in higher and spring season is lower than the other seasons. Higher SL in autumn is due to higher rainfall, more river runoff, however it is found that higher typhoon activity observed in the summer especially in the August month . From 1975 to 1987 seasonal SL is showing decreasing trend, then have an increasing trend.
Over the HZB (Figure 3(a)), the seasonal ranges in different seasons are given as follows: in the spring SL varies between −3.8 and −19.1 cm, in summer −5.09 and 20.8 cm, in autumn 7.2 - 20.8 cm and in winter −16.7 and 3.1 cm respectively.
Figure 2. Monthly variation of Sea level anomalies over different areas of Zhejiang Province ((a) Hangzhou Bay (HZB), (b) Ningbo (NB) and (c) Wenzhou (WZ) for the period 1958-2015.
In NB (Figure 3(b)), spring SL varies between −2.7 and −18.5 cm, in summer −5.9 and 20.5 cm, in autumn 6.96 and 20.54 cm and in winter −15.7 and 4.4 cm respectively. In WZ (Figure 3(c)), the seasonal ranges in in the spring SL varies between −0.7 and −19.1 cm, in summer −7.4 and 20.2 cm, in autumn 6.3 and 20.2 cm and in winter −15.5 and 6.5 cm respectively. These ranges clearly indicating that SL of WZ higher that other two areas and is due to the impact of lower latitude as well the intrusion of Kuroshio current. The steric effect of the SST is higher in the lower latitudes.
3.3. Annual Variations of SL
The interannual SL variations of the Zhejiang provincial coastal areas (HZB, NB and WZ) are depicted in Figure 4 and showing an increasing trend. The SL in the WZ is higher than other two areas (NB and HZB). This variation is SL can be explained as the latitudinal effect. Previous studies concluded that large deflection of the Kuroshio can lead to the rise in SL . Because during the great bend of the Kuroshio, it carries a large amount of high temperature and high salt water into the East China Sea, It has a great influence on the marine environment in the East China Sea  , Sea level change is the appearance of this effect. At the same time, we can also see some years of low sea level, for instance
Figure 3. Seasonal variation of SL over different areas of Zhejiang Province ((a) Hangzhou Bay (HZB), (b) Ningbo (NB) and (c) Wenzhou (WZ)) for the period 1958-2015.
Figure 4. Annual variation of SL over different areas of Zhejiang Province ((a) Hangzhou Bay (HZB), (b) Ningbo (NB) and (c) Wenzhou (WZ)) for the period 1958-2015.
1963, 1968, 1981, 1985, 1988, 1995, 2002, 2007 and 2015. SL decreased in the study period in the El Nino, Lower SL clearly indicating that El Nino events have a great influence on Zhejiang provincial areas to modify SL. Western Pacific sea water cooling, so that the sea level dropped and the East Pacific Ocean warming, sea level rise during the formation and development of El Nino .
Figure 5 clearly indicates the increasing trend of SL over HZB, NB and WZ areas of Zhejiang province. In the study period, lower anomalies of SL observed till 1973, there is an increase in SL for in 1976, then decrease again and fluctuations in SL can be seen till 1995 and an increase trend till 2015. However, 2005 is special case which has lower SL. Year 2005 is a special ENSO events in recent decades, and have the distinctive evolution features compared to the typical ENSO events. The typical El Nino has a fast decaying after its peak . Weak warm phase in the first half of 2005, mean SST anomalies were observed in the central equatorial Pacific Ocean are higher, which peaks in May. In the western equatorial Pacific remained slightly above or near-normal in the last few months of 2005, however recent monthly mean anomaly in December 2005 decreasing below the long-term average.
3.4. Decadal Variations
Decadal anomalies clearly illustrate an increasing trend of SL (Figure 6). The decade 1968-1977, SL increased however in 1978-1987-decade SL declined, after
Figure 5. Annual SL anomaly variations over different areas of Zhejiang Province ((a) Hangzhou Bay (HZB), (b) Ningbo (NB) and (c) Wenzhou (WZ)) for the period 1958-2015.
Figure 6. Decadal variation of SL over different areas of Zhejiang Province (((a) Hangzhou Bay (HZB), (b) Ningbo (NB) and (c) Wenzhou (WZ)).
that there is an increase trend in SL. Increase in SL after 1987 may be connect with the climate change, global warming and increase in temperature over the tropical Pacific Ocean.
3.5. SL Variations Related to SST
SST variations will inevitably lead to the change of sea level. This article uses the sea surface temperature (HadISST) data to find out the variations of sea surface temperature in Zhejiang province areas. We use HadISST (1958-2015) SST data and sea surface height data to find out the relation over the Zhejiang provincial areas and presented in Figure 7.
In order to further study the correlation between the SL and SST, we analyzed the correlation between SL and SST for the study period (1958-2015). The correlation coefficients are as follows: over the HZB the correlation coefficient (CC) between SL and SST is 0.72, at NB the CC is 0.69 and at WZ the CC is 0.61. The CC clearly indicating there is a significant relation between SL and SST. However, WZ is lower than HZB, this shows over the HZB region the influence of SST on SL is higher. The southern city of WZ having an influence of Kuroshio warm currents during the summer and in the winter, it was influenced by the cooler current coming from the northern side, also due to the rainfall and river runoff.
Global warming leads to expansion of the ocean, polar ice caps and terrigenous Melting glaciers are the main reason of global sea level rise . From the global sea level rise distribution Sea level changes in the world oceans has obvious regional characteristics . China’s coastal areas are located in a relatively high level of sea level rise rate is higher . Regional SL changes in addition to the affected by the global sea level change, also by the local SST, current, wind, temperature, pressure, precipitation and so on hydrologic meteorological elements , higher SST and lower pressure, sea level rise  . Since 1990, with extreme weather and climate events occur frequently, resulting in the number of abnormal seasonal sea level increase obvious, anomaly values increased . The affection of the East Asian monsoon changes also is the important cause of the change of China’s coast and the neighboring seas Marine environmental elements, such as Seawater temperature and salinity, sea level,
Figure 7. Relation between Sea level (SL) and Sea surface temperature (SST) during the study period (1958-2015) over Zhejiang province areas ((a) Hangzhou Bay (HZB), (b) Ningbo (NB) and (c) Wenzhou (WZ)).
etc  . The relation between the sea temperature and SL is given in the Figures 7(a)-(c). There is a positive relation can be seen between SST and SL with the correlation of 0.7284. And Figure 7 illustrating relation between SST and SL and indicating an increasing trend.
Over the Zhejiang province areas SL is indicating and increasing trend from monthly, seasonal, annual and decadal variations. SL is increased 0.17 cm per year from 1958 to 2015. As previous studies indicating that SST is having an impact on SL, we found there is a positive relation between Zhejiang areas SST and SL. HZB area illustrating higher relation between SL and SST, however the seasonal annual plots showing latitudinal effect on the SL. WZ indicating higher SL than HZB. Local SST, river discharge, sediment in the water are having sufficient impact on the SL. However, Zhejiang provincial areas experience several typhoons, which may produce higher rainfall and changes in the sea surface heights, which may also have sufficient effect on SL.
Authors express their gratitude to ECMWF for providing sea surface height data freely and Hadley Centre Global Sea Ice and Sea Surface Temperature for providing SST data.
 Church, J., Gregory, J.M., Huybrechts, P., Kuhn, M., Lambeck, K., Nhuan, M.T., Qin, D. and Woodworth, P.L. (2001) Changes in Sea Level. In: Houghton, J.T., et al., Eds., Climate Change 2001: The Scientific Basis, Cambridge University Press, Cambridge, 639-693.
 Woodworth, P., Gregory, J.M. and Nicholls, R.J. (2004) Long-Term Sea Level Changes and Their Impacts. In: Robinson, A.R. and Brink, K., Eds., The Sea, Biological-Physical Interactions in the Sea, Vol. 12, Harvard University Press, Cambridge, 717-752.
 Li, C.Y., Sun, S.Q. and Mu, M.Q. (2001) Origin of the TBO-Interaction between Anomalous East-Asian Winter Monsoon and ENSO Cycle. Advances in Atmospheric Sciences, 18, 555-566. https://doi.org/10.1007/s00376-001-0044-y
 Chen, Z.Y., Huang, Y.H., Zhou, T.H., Tang, E.X., Yu, Y.F. and Tian, H. (1991) A Preliminary Study on Mean Sea Level of the Changjiang River Estuary. Oceanologia et Limnologia Sinica, 22, 315-320.
 Cazenave, A., Dominh, K., Gennero, M.C. and Ferret, B. (1998) Global Mean Sea Level Changes Observed by Topex-Poseidon Altimetry and ERS-1. Physics and Chemistry of the Earth, 23, 1069-1075. https://doi.org/10.1016/S0079-1946(98)00146-3
 Repert, J.P., Donguy, J.R., Elden, G., et al. (1985) Relations between Sea Level, Thermocline Depth, Heat Content, and Dynamic Height in the Tropical Pacific Ocean. Journal of Geophysical Research, 90, 11719-11725. https://doi.org/10.1029/JC090iC06p11719
 Chambers, D.P., Tapley, B.D. and Stewart, R.H. (1997) Long-Period Ocean Heat Storage Rates and Basin-Scale Heat Fluxes from TOPEX. Journal of Geophysical Research, 102, 10525-10533. https://doi.org/10.1029/96JC03644
 Chen, J.L., Shum, C.K., Wilson, C.R., et al. (2000) Seasonal Sea Level Change from TOPEX/Poseidon Observation and Thermal Contribution. Journal of Geodesy, 73, 638-647. https://doi.org/10.1007/s001900050002
 Cabanes, C., Cazenave, A. and Le Provost, C. (2001) Sea Level Rise during Past 40 Years Determined from Satellite and in Situ Observations. Science, 294, 840-842. https://doi.org/10.1126/science.1063556
 Leuliette, E.W. and Wahr, J.M. (1999) Coupled Pattern Analysis of Sea Surface Temperature and TOPEX/POSEIDON Sea Surface Height. Journal of Physical Oceanography, 29, 599-611. https://doi.org/10.1175/1520-0485(1999)029<0599:CPAOSS>2.0.CO;2
 Nerem, R.S. and Mitchum, G.T. (2001) Sea Level Change. In: Fu, L.-L. and Cazenave, A., Eds., Satellite Altimetry and Earth Sciences: A Handbook of Techniques and Applications, Academic Press, Cambridge, 329-347. https://doi.org/10.1016/S0074-6142(01)80153-4
 Nerem, R.S. and Mitchum, G.T. (2001) Observations of Sea Level Change from Satellite Altimetry. In: Douglas, B.C., Kearney, M.S. and Leatherman, S.P., Eds., Sea Level Rise: History And Consequences, Academic Press, Cambridge, 121-159. https://doi.org/10.1016/S0074-6142(01)80009-7
 IPCC (2007) Climate Change 2007: The Physical Science Basis, Summary for Policy makers. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, Intergovernmental Panel on Climate Change.
 Nong, G.T., Najjar, R.G., Seidov, D. and Peterson, W.H. (2000) Simulation of Ocean Temperature Change Due to the Opening of Drake Passage. Geophysical Research Letters, 27, 2689-2692. https://doi.org/10.1029/1999GL011072
 Cramer, B.S., Toggweiler, J.R., Wright, J.D., Katz, M.E. and Miller, K.G. (2009) Ocean Over-Turning since the Late Cretaceous: Inferences from a New Benthic Foraminiferal Isotope Compilation. Paleoceanography, 24, 1-14. https://doi.org/10.1029/2008PA001683
 Katz, M.E., Cramer, B.S., Toggweiler, J.R., Esmay, G., Liu, C., Miller, K.G., Rosenthal, Y., Wade, B.S. and Wright, J.D. (2011) Impact of Antarctic Circumpolar Current Development on Late Paleogene Ocean Structure. Science, 332, 1076-1079. https://doi.org/10.1126/science.1202122
 Gao, S., Zhu, L., Zhang, W., et al. (2020) Western North Pacific Tropical Cyclone Activity in 2018: A Season of Extremes. Scientific Reports, 10, Article No. 5610. https://doi.org/10.1038/s41598-020-62632-5
 Zhang, J.L., Chen, M.C., Chen, M.X., et al. (2009) Characteristics of SST Variation/Change in the Coastal Region of the East China Sea. Proceeding of the Nineteen (2009) International offshore and Polar Engineering Conference, Osaka, 21-26 June 2009, 483-489.