Analyzing the Effects of Short-Term Persistence and Shift in Sea Level Records along the US Coast

: Is sea level change affected by the presence of autocorrelation and abrupt shift? This question reﬂects the importance of trend and shift detection analysis in sea level. The primary factor driving the global sea level rise is often related to climate change. The current study investigates the changes in sea level along the US coast. The sea level records of 59 tide gauge data were used to evaluate the trend, shift, and persistence using non-parametric statistical tests. Mann-Kendall and Pettitt’s tests were utilized to estimate gradual trends and abrupt shifts, respectively. The study also assessed the presence of autocorrelation in sea level records and its effect on both trend and shift was examined along the US coast. The presence of short-term persistence was found in 57 stations and the trend signiﬁcance of most stations was not changed at a 95% conﬁdence level. Total of 25 stations showed increasing shift between 1990–2000 that was evaluated from annual sea level records. Results from the current study may contribute to understanding sea level variability across the contiguous US. This study extends an elaborative understanding of sea level trends and shifts which might be useful for water managers.


Introduction
According to the Intergovernmental Panel on Climate Change's (IPCC) [1] Fifth Assessment Report, the global sea level is estimated to rise 0.28-0.98 m by 2100  in different representative concentration pathways scenario. Sea level in the Atlantic and Gulfs along US coasts have risen in the 20th century [2]. The local sea-level change may vary considerably from global sea level change due to factors such as land subsidence and emergence, global sea-level fluctuations, ocean circulation, and other localized factors. Sea level change in Sitka, Alaska is decreasing due to land emergence whereas Louisiana and Atlantic City are experiencing rising levels due to subsidence, which is unrelated to global sea level change according to the National Ocean and Atmospheric Administration website (NOAA, https://oceantoday.noaa.gov/globalvslocalsealevel/). Therefore, depending on the locations, and the rate over time it can cause spatial and temporal variability in sea level change [3].
Most of the studies found an acceleration in sea level while some focused on sea-level deceleration. Houston and Dean [4] found that during the 20th century the 57 US tide gage records did not show any acceleration. However, Rahmstorf and Vermeer [5] disagreed with Houston and Dean [4] and illustrated that using different start years for trend shows acceleration in sea level. The trend literatures corroborate that the hydroclimatological time series exhibit trends and the presence of autocorrelation [6] might alter the significance of the trend. Furthermore, Meyssignac et al. [7] state that internal climate variability has a dominant signal in the tropical Pacific than the anthropogenic forcing on regional sea level trends. This had led to an intensive debate in sea level studies and is also documented by Visser et al. [8]. To address these issues, detection of long-term sea level variability within the earliest record is essential to shaping our perspective of climate change as well as for policymakers to implement proper management plans to minimize the threat to the coastal ecosystems [9,10].
The increase in the rate of mean sea level is one of the main facets of probable climate variability and change [11]. The sea level natural variability is a non-stationary process that also drives sea-level fluctuations [3,12,13]. Trend detection has been used widely for quantifying changes in hydroclimatological data. The trend detection of hydroclimatic data, which are non-stationary in nature, requires a robust statistical method that accounts for non-stationarity in time series [14]. Trends are estimated using both parametric and non-parametric tests. For most of the parametric tests, the time series are required to be normally distributed, which does not always apply to extreme events [15].
Non-parametric tests are the most widely used and recommended technique due to their advantages that it does not assume any initial probability distribution and accounts for non-stationarity in the time series [16,17]. Therefore, the non-parametric Mann-Kendall (MK) test that requires independent data is highly beneficial for the trend estimation of hydroclimatic data. MK test is a rank-based robust approach, which is not sensitive to the distribution of the data, or the outliers, and is also efficient with missing data. However, it has some limitations on autocorrelated datasets. Yue et al. [18] and Hamed and Rao [19] corroborated that the trend significance was overestimated with the presence of positive autocorrelation whereas underestimated by negative autocorrelation in the time series. The hydroclimatic variables have a tendency to exist in clusters in a course of time, i.e., during extreme events; this is termed 'scaling' or 'persistence' [20]. These results are due to a lack of independence and are termed as short-term persistence (STP). The presence of one significant autocorrelation (lag-1) results in STP. Teegavarapu and Schmidt [21] found strong autocorrelation in monthly sea level observations in US tide gages. In addition, Boon et al. [22] also found the presence of autocorrelation in Washington DC tide gage and address the importance of removing significant autocorrelation for sea level trend estimation in Chesapeake Bay station in the United States. In addition, Becker et al. [23] also used power law and corroborated that sea level change exhibits long-term correlations. Therefore, to evaluate the sea level trend within a time-series, the effects of STP should be carefully considered.
In geosciences, climate change and variability have created a need to detect persisting trends in gradual changes over a long historic period [24]. Similarly, studies have also detected sudden changes, which represent the abrupt changes in the hydroclimatological data [25]. Trends (gradual) are often observed throughout the time series and abrupt changes (shift) are the sudden changes that are due to non-homogeneity in the time series. Previously, Kerr [26] and Trenberth and Hurrell [27] demonstrated a sudden change during 1976-1977 in the Northern Pacific Ocean-atmosphere system which was later termed the climate regime shift of 1977 [28]. According to Villarini et al. [29], categorizing trend (linear) and shift (non-linear) is essential as trend estimation can be affected by the shift in the data. Shift estimation is essential to understand whether the observed variation in the data is an effect of long-term variability, or a consequence of an abrupt shift [25]. Change-point analysis is done by utilizing a non-parametric Pettitt's test that can detect the location of abrupt changes (shift) in the mean and variance in the time series. McCabe and Wolock [30] also argued the presence of a trend might continue in the future. A trend typically follows a shift and time series may remain the same until another shift occur. Thus, a simultaneous study of both trend and shift is essential for the improved evaluation of the sea level variability. Moreover, evaluating the persistence while estimating the shift in the sea level would be an addition to the research gap.
The goal of the current study is to evaluate long-term sea level variability across the US coast by evaluating both gradual and abrupt changes in decadal scale. The study was motivated by the need to analyze whether the sea level is rising, falling, or exhibits no trend utilizing purely statistical approaches. The study addressed the following key questions (1) What is the overall sea level variability along the coastal US? (2) What is the significance of autocorrelation in both trend and shift? and (3) Does STP affect the local sea level trends around the USA? These research questions will assess the sea level change and the effects of autocorrelation on trend and shift. The major objectives of the current study are to (1) detect gradual changes by utilizing the MK test by considering autocorrelation for more realistic error statistics for the estimated trend parameters, and (2) detect abrupt shifts and confirm trend significance effects due to the presence of autocorrelation. Mann-Kendall test was performed in three ways to address the objective: (i) Mann-Kendall without autocorrelation (MK1), (ii) Mann-Kendall with lag-1 autocorrelation and trend-free prewhitening (MK2), and (iii) Pettitt's test to detect and evaluate the significance of the shift in the time series.

Study Area and Data
The latitude, longitude, and time-series data used in the current study are the insitu water level measurements along the US coast acquired from the NOAA (https:// tidesandcurrents.noaa.gov/stations.html?type=Water+Levels). The dataset used in this study represents the local relative sea level change. Out of 95 stations from NOAA, 59 stations were selected for this study on basis of the existing continuous data as shown in Figure 1. The monthly and annual raw data were obtained and the long-term data with the station having more than 30 continuous years of time-series data with no omission were used for the current study. For the analysis, the stations with maximum continuous time-series datasets were obtained. In total, 59 time-series were obtained with the station that had 70 years of data having a maximum record, while each gauge had at least 30 years of continuous time series. In this study, variable record length was used to cover maximum number of stations with more than 30 years of continuous records. The data from March to May is referred to as spring season, June to August as summer season, September to November as fall season, and December to February as winter season.

Trend Estimation Techniques
Variations in sea level can be estimated using a robust non-parametric MK test to determine the trend. The key advantages of the MK test are due to its accuracy working with a skewed distribution [18,[31][32][33], and that it accounts for non-stationarity. In this test, each time series is ranked with reference to time in each station and test statistics S is computed. However, the MK test has limitations against autocorrelation resulting in spurious trends [19,34]. To account for this limitation, significant lag-1 autocorrelation significant at 90% was removed from the dataset to evaluate trends by utilizing pre-whitening to account for STP [18]. Persistence in a time-series results from an overestimation of a trend as an effect of an underestimated variance [19,35,36]. Thus, considering the STP from the literature review is important, and was underlined in this study. A brief description of MK tests extracted from Yue et al. [18] and Hamed [37] is presented below: The null hypothesis for the MK test states:

Hypothesis 2. Upward or Downward trend exists.
The test statistics S was obtained by using Equation (1) and function sgn(x) by Equation (2) below, x k and x i are observations at a total length of time n. The variance of S is calculated using Equation (3).
The z test statistics is calculated using Equation (4).
The MK test is based on the null hypothesis that the time series has no trend. The direction (upward or downward) of the trend was determined by the MK test based on the obtained test statistics S sign (positive or negative). The significance level is given by test statistics z and determines the chances of rejecting the null hypothesis whereas confidence interval was calculated by subtracting significance level from 1.
A robust non-parametric method, Thiel-Sen approach (TSA) [38,39], is used in the current study to estimate the magnitude of trend. According to Sen [39], it is effective with missing data and is also sensitive to outliers. It was computed using the β slope given in Equation (5). The estimate of the beta slope gives the rate of change to time. MK test without accounting for autocorrelation is referred to as MK1 hereafter.
The MK1 test has limitation to the presence of autocorrelation. To account this, the lag-1 autocorrelation coefficient (r k ) was calculated using Equation (6).
The autocorrelation was removed by utilizing pre-whitening technique if the data lies outside the confidence level i.e., if data are found to be serially dependent. Else, trend free pre-whitening (TFPW) was applied as defined by Yue et al. [18]. The β slope was removed from the series to get detrended series by using Equation (7). Then, the lag-1 autocorrelation was removed from the detrended series by using Equation (8), Then the trend is added back to the detrended series to get new series by using Equation (9). The MK test statistics S and z are again computed for the new series. MK test after accounting for autocorrelation is referred to as MK2 hereafter.
Separate MK1 test on each station was computed for annual and monthly time series and then MK2 were utilized independently to analyze the effects of short term persistence. Further, β slope on each station was also estimated separately for the different study period on the longest continuous available data.

Shift Detection
In addition to assessing the trend, abrupt changes were determined using the nonparametric Pettitt's test [40]. It is a widely used test to analyze the non-homogeneity in the data. It is a rank-based approach and is less sensitive to outliers and skewness of the data. The Pettitt test evaluates the anomaly in mean within a time series at n length of time considering a shift at j years. It uses Mann-Whitney statistics to estimate the probability that the trends on either side of the shift at j years are significantly different by trials from the same population.
The test statistics are calculated using Equation (10).
where sgn(x) is a function defined in Equation (2) above. The significant change point is obtained if test statistics |Uj,n| is maximum and is calculated using Equation (11), K T is the test statistics associated with the confidence level P.
At maximum |Uj,n|, the probability of shift is given by Equation (12), If the significance level is more than the p-value, the null hypothesis is rejected, and a shift is detected. The direction of change also depends on the considered significance level; negative change indicates the minimum value, and the positive change indicates the maximum value.
Also, serial correlations were accounted for computing with and without TFPW in Pettitt's test as proposed by Serinaldi and Kilsby [41] to evaluate the existence of significant shifts in the dataset. Like trend, the time series were detrended by deducting the trend slope [18] and removing the autocorrelation from the detrended dataset. This is the novel approach of the current study that makes the shift detection in sea level more reliable. Hence, with TFPW, the underlying trends can be detected, making Pettitt's test more dependable. A 95% significance level (or p ≤ 0.05) was utilized to assess the significance of Pettitt's test.

Sea Level Trend and Persistence
The annual sea level trend analysis was performed to evaluate trends and persistence in sea level records. A total of 52 stations showed a positive trend, one showed negative and 6 showed no trend resulting from the MK1 test as shown in Table 1. The annual trend and slope were also obtained from the NOAA website (https://tidesandcurrents. noaa.gov/sltrends/mslUSTrendsTable.html) for comparison with the current study and shown in Figure 2a. The total number of stations obtained from NOAA website is 95 and in the current study, 59 stations with continuous time series were used. The time period of two stations used in this study is same as NOAA. NOAA records showed a decreasing trend in three stations in the Northwest. Of those, two stations showed no trend resulting from the MK1 test as shown in Figure 2b. Similarly, one station on the east coast and one station on the west coast showed no trend resulting from MK1 trend analysis whereas most of the stations showed a positive trend (Table 1). However, the magnitude of trends obtained from NOAA had similar range as that obtained from the MK1 test and is shown in Figure 2a,b. The STP was also accounted for annual sea level records to estimate the trend and analyze the effects of the presence of serial correlation on the trend estimation. The sea level trend evaluated after accounting for STP showed results like MK1 where most of the stations showed an increasing trend. Figure 2c shows the spatial plot of sea level trend after accounting for STP. From this, it can be observed that the significance of the trend of no station was changed to a non-significant trend. Instead, 52 stations showed a positive annual sea level trend, even after accounting for STP.    Table 1 shows the total number of stations with significant MK1 trends at a 95% confidence level for all months. Through the MK1 test, the sea level at 38 stations showed a significant positive trend. Furthermore, as observed, most of the stations showed an increasing trend for all months (Jan-Dec). The results indicate that the sea level trends of 40 stations exhibited a positive trend in March to November and 1 to 4 stations exhibited a negative trend in February, and May to September. July had the greatest number of positive trends (52, or 88%) followed by August having 51 stations, 50 stations in September, and 49 stations in October, whereas February had the fewest number of positive trends (33). The greatest number of significant negative monthly trends were in September (4), and the absence of negative trends was in January, March, April, October, November, and December. A significant positive trend with more than 40 stations was found in 9 months mostly located along the East Coast; however, more than 32 stations showed increasing trends every month. Table 1 shows the temporal sea level variation as increasing/positive, decreasing/negative, and no trend/no change along the US coast, at the 95% confidence level.
The magnitude of the trend was computed by utilizing the TSA slope for significant stations showing a trend at the level of significance, p ≤ 0.05. The maximum and minimum trend magnitude was 0.85 cm/year for October and 0.62 cm/year for January month, respectively as shown in Figure 2d using the entire dataset. Each whisker shows the outliers, each red line shows the median, and each box represents the values within the 25th and 75th percentile of the sea level trend magnitude for each month. The results of sea level indicate that the magnitude of sea-level starts increasing from May with slight fluctuations through September. It was noteworthy that the annual and monthly median was above zero. The maximum annual trend slope along the US coast was 0.64 cm/year whereas the minimum trend magnitude was −0.20 cm/year using the entire data. Figure 2e shows the distribution of TSA slope by using the latest 30 years of available data from each station at different temporal scale (annual and monthly). The latest 30 years of data were used to observe the recent 30-year trend pattern in sea level because a minimum 30-year record is recommended to capture climate trends. It showed a similar range of trend magnitude that was obtained by using the entire dataset as shown in Figure 2d,e. The maximum annual trend slope along the US coast using the latest 30 years of available data was found to be 0.82 cm/year whereas the minimum trend magnitude was −0.11 cm/year. The 30-year trend magnitude was found to be 0.96 cm/year for September and 0.61 cm/year for January month.
The STP was accounted to estimate the monthly sea level trend and to assess the influence of presence of serial correlation on the significance of the trend estimation. The sea level trend evaluated with MK2 test showed similar results to MK1 as shown in Table 1. More than 28 stations showed significant trends after accounting for STP for all months as shown in Figure 3. Table 1 summarizes the total number of stations showing trends after accounting for STP. The highest lag-1 autocorrelation was observed in September and October having maximum stations with a positive trend (49 stations), followed by 39 stations in November. Out of 59 stations along the US coasts, the East Coast of the US had significant increasing trends throughout the year as observed from the spatiotemporal sea level variation from the MK1 test. Comparatively, the results accounting for STP indicated an increasing trend from summer through fall and decreasing from winter through spring with 48 stations showing an increasing trend in May. A maximum station showing a positive trend in sea level over the study period was observed for the spring season compared to the other seasons. The number of stations with an increasing trend increased from early spring-March (1 station)-through mid-fall-October (21 stations). The negative trend slowly reduced in winter with the least number of stations observed in February (1 station). Few stations showed a negative trend. From this, it can be observed that the trend significance of only a few stations was changed to a non-significant trend, whereas most of the stations showed a positive monthly sea level trend, even after accounting for STP, as shown in Table 1 and Figure 3.

Shift Detection in Sea Level
Thirty-eight out of 59 stations showed a significant positive abrupt change in Janu and 21 stations showed no shift in the data. The maximum number of stations showin positive shift was 52 (88%) in May, before accounting for autocorrelation, and Aug both before and after accounting for autocorrelation (Table 1). However, for May on stations showed a negative shift and 4 stations in September both before and after counting autocorrelation. The Pettitt test results identified that the shift years cluste from the 1960s to 2012. The maximum significant negative shift was found in four stati in September. Most of the increasing shifts were found between 1991 and 2000 as sho in Table 2. A maximum number of stations showed a positive shift from May thro October. A few stations in some months experienced a shift earlier in the year. Along w the trend, it showed an increase in the presence of abrupt changes. Furthermore, Pett test was also performed in annual sea level record and it was observed that 25 stati showed a shift between 1991-2000 as shown in Table 2. Figure 4a The step with pre-whitening was also utilized to consider the effects of serial co lation on the Pettitt test on the monthly and annual sea level records. The shift was inv tigated after accounting for serial correlation and the overall similar pattern of the s was observed when compared with shift before accounting for autocorrelation. Howe few stations had different shifts years after accounting for serial correlation in the Pe test in a few months as shown in Figure 4b. Overall, the abrupt shift analysis confirms the effects of STP on the stations and the significance of only a few stations changed a

Shift Detection in Sea Level
Thirty-eight out of 59 stations showed a significant positive abrupt change in January and 21 stations showed no shift in the data. The maximum number of stations showing a positive shift was 52 (88%) in May, before accounting for autocorrelation, and August, both before and after accounting for autocorrelation (Table 1). However, for May only 3 stations showed a negative shift and 4 stations in September both before and after accounting autocorrelation. The Pettitt test results identified that the shift years clustered from the 1960s to 2012. The maximum significant negative shift was found in four stations in September. Most of the increasing shifts were found between 1991 and 2000 as shown in Table 2. A maximum number of stations showed a positive shift from May through October. A few stations in some months experienced a shift earlier in the year. Along with the trend, it showed an increase in the presence of abrupt changes. Furthermore, Pettitt's test was also performed in annual sea level record and it was observed that 25 stations showed a shift between 1991-2000 as shown in Table 2. Figure 4a,b represents the spatial maps showing shift year before and after accounting autocorrelation in the shift. It can be noted that the nearby stations have a similar shift year. From Figure 4c it is evident that 57 stations out of 59 stations show the presence of lag-1 autocorrelation.
The step with pre-whitening was also utilized to consider the effects of serial correlation on the Pettitt test on the monthly and annual sea level records. The shift was investigated after accounting for serial correlation and the overall similar pattern of the shift was observed when compared with shift before accounting for autocorrelation. However, few stations had different shifts years after accounting for serial correlation in the Pettitt test in a few months as shown in Figure 4b. Overall, the abrupt shift analysis confirms that the effects of STP on the stations and the significance of only a few stations changed after applying the TFPW in the shift test. In addition, 53 out of 59 stations showed a positive shift before accounting for autocorrelation whereas 52 stations after accounting for autocorrelation resulted from shift analysis in annual sea level record as shown in Table 1. The spatial plot of shift significance before and after accounting autocorrelation at 95% confidence level is shown in Figure 4d,e respectively. applying the TFPW in the shift test. In addition, 53 out of 59 stations showed a positive shift before accounting for autocorrelation whereas 52 stations after accounting for autocorrelation resulted from shift analysis in annual sea level record as shown in Table 1. The spatial plot of shift significance before and after accounting autocorrelation at 95% confidence level is shown in Figure 4d,e respectively.  (Table 2). Only 11 out of 59 showed the presence of shift in spring, respectively. From this, we can infer that during spring the minimum number of stations exhibiting a shift was observed (only 1) with 11 stations being the maximum during that interval.   (Table 2). Only 11 out of 59 showed the presence of shift in spring, respectively. From this, we can infer that during spring the minimum number of stations exhibiting a shift was observed (only 1) with 11 stations being the maximum during that interval.

Discussion
Both MK and Pettitt's test works well on distribution-free time series thus, autocorrelation was tested in the current study, and if present, autocorrelation was removed with the TFPW technique as suggested by Yue et al. [18] and Hamed [37]. The presence of autocorrelation was identified before evaluating trend and shift for precise estimation of sea level change. The presence of autocorrelation alters the variance of the distribution. Serinaldi and Kilsby [41] emphasized the importance of incorporating the serial correlation for shift detection. The characteristics and statistical significance of these serial correlations and their effects on the sea level variability have not been identified in previous studies. In the current study, the long-term trend was analyzed to examine the persistence as it can pose significant deviations and affect the local trends. The inclusion of persistence was also useful to address evidence to different real and spurious trends [42].
This study also provides strong scientific evidence that the assumption of stationarity while assessing sea level trends may not be valid. The MK and Pettitt's tests are widely used as it accounts for non-stationary behavior in time series. The sea level was observed to have STP at most of the stations. In this study, both STP and the abrupt shift were detected which signifies the presence of non-stationarities in the dataset. The lag-1 autocorrelation was considered to remove significant autocorrelation in the sea level to correctly evaluate the variation in the sea level. The presence of STP in the data does not alter the magnitude of trends, but it significantly reduced the statistical significance of the trends in Kumar et al.'s research [35]. Previous studies have shown that the z-value of MK tests was significantly reduced after incorporating STP for annual flow statistics [35] temperature, precipitation, and streamflow [43]. However, in this study, it was observed that statistically significant sea level trends may not be statistically significant after accounting for the STP in only a few stations as shown in Table 1. Most stations on the East Coast exhibited an increasing trend, even after considering the STP.
The significance of autocorrelation was tested in the trend and shift detection and it was observed that the significance of trend did not change in most of the stations, even after accounting for the autocorrelation in annual and monthly sea level record as shown in Table 1, and Figure 2. A possible mode of non-stationarity was presented by using the shift analysis. Teegavarapu [44] also states the importance of considering the abrupt shift in sea level trend analysis. The time-series detecting shift and autocorrelation was the main focus of the current study to estimate the temporal trend as suggested by Teegavarapu [44]. The detection of shift showed that the variation in the sea level is also affected mostly by the shift in the data in addition to the long-term gradual changes as trends. Further, the analysis shows that the presence of autocorrelation in a significant shift does not change the trend significance in most of the stations as shown in Figure 3 and Table 1. In the future, another trend in the sea level may be observed with another abrupt shift in the data.
The 1970s shift years detected in this study are well documented in previous studies [45]. Further, Woodsworth et al. [46] also found positive acceleration related to an 'inflexion' point in the US tide gages in 1930. The determination of time of occurrence of sudden changes in time series aids researchers to find the evidence of natural or anthropogenic changes in climatic, hydrologic, or land use. In addition, the abrupt shift could also be due to non-homogeneity in the dataset which is often affected by the moving of the stations. The significant shift years also represent non-stationary behavior in sea level records. Therefore, it is very important to identify the period from where significant change has occurred in sea level record. In the current study, summer season showed a maximum of the station showing a shift in all categorized years with a maximum of 40 stations in 1991-2000. However, few nearby stations are observed in the different shift year range which could be due to the variable record length used in this study.
The World Meteorological Organization (WMO) [47] noted the unprecedented decadal rise of temperature in 1991-2000 and 2001-2010. The sea surface temperature is high during the summer and 1991-2000 was warm due to increased temperatures [47][48][49]. This may have increased the positive sea level trend as the current study observed a higher number of stations showing a shift during 1991-2000. Reager et al. [50] observed that sea-level rise has been slowed due to climate-driven processes and human activities. It is also discussed by Chao et al. [51] that water management in the reservoir could be one of the factors of slow sea level rise and found that the trend since 1930 would be almost linear than deceleration. The results of abrupt changes could be due to variation in gauging practices, land-use change, and regulations of rivers by reservoir operations, water management implementations as well as being an indication of the severity of a climate event. Also, the shift could be due to anthropogenic processes like water pumping, urbanization, terrestrial water storage changes, and storm surges which have not been considered in the current study and could be a potential area for future research.

Conclusions
Previous studies on sea level changes suggest both an increase and a decrease in sea level, whereas in this study the trend was found to be significant at 95% confidence level, even after accounting for autocorrelation. The current study intends to analyze both trends and shift at both spatial and temporal scales along the contiguous US coasts. The abrupt and gradual changes in sea level utilizing the available tide gauge (30-70-year record) dataset were analyzed in this study. The current study focused on long-term changes as trend and abrupt shift in monthly and annual sea level along the US coast.
Persistence at short-term temporal scales was considered to estimate its effects in both trend and shift. The study accounted for lag-1 autocorrelation to confirm STP. The key conclusions from the present analysis can be summarized as:

1.
The gradual trend analysis revealed the presence of STP in most of the stations and it is evident that most stations show a positive sea level trend, even after accounting autocorrelation. The maximum and minimum trend magnitude that was evaluated by using the entire dataset was found to be 0.85 cm/year for October and 0.62 cm/year for January, respectively.

2.
Previous studies focused on gradual trends while the current study is unprecedented as it considers the presence of both gradual and abrupt changes in sea level.

3.
The autocorrelation in sea level can significantly change the detection of a sea level trend and shift. However, the most of stations showed a positive trend before and after accounting for the autocorrelation. Accounting for autocorrelation yielded more reliable results. 4.
The number of stations showing significant shifts was similar after accounting for autocorrelation in sea level records. In the annual time series, 53 stations showed an increasing trend before autocorrelation, and 52 stations showed an increasing shift, even after the accounting of the autocorrelation in the shift.

5.
This study concludes that the presence of autocorrelation did not change the significance of trend and shift in most stations. Most stations showed increasing trend and shift after accounting for autocorrelation.
This study identified the shift years and influential seasons within the spatial distribution of stations undergoing significant shifts. These results could be useful to understand the behavioral change in each station. Additionally, this study is intended for future research in how we understand sea level changes, and is also beneficial for coastal water management.