The Periodic Response of Tidal Flat Sediments to Runoff Variation of Upstream Main River: A Case Study in the Liaohe Estuary Wetland, China

(1) Background: To reveal the intrinsic relationship between the tidal flat sediments in an estuary wetland and the runoff from the upstream river. This research was conducted in the tidal flats of the Liaohe estuary wetland. (2) Methods: The 137Cs and 210Pb dating technique was used to reconstruct the time correspondence between tidal flat sediments and runoff, and the periodic response was explored between the changes in the tidal flat sediments and runoff based on the spectrum analysis method. (3) Results: The average sedimentation rate in the tidal flat was 2.24 cm·year−1 during the past 50 years. The amount of fine sediment particles deposited on the estuary tidal flat was directly related to the amount of sediments transported by the river and inversely proportional to the ability of rivers to transport fine matter. The high frequency reproduction cycle of 14–15 years in the flood season flow and 5–6 years in the annual sediment discharge of the Liaohe River correspond to the high and low frequency reproduction cycles of the median size of sediments in Liaohe estuarine wetland tidal flats. (4) Conclusions: The research clarified the hydrological constraints of the action law between Liaohe River runoff and the estuarine sediments. The periodic response between Liaohe River runoff and the sediment was established.


Introduction
The ocean is the final destination of fine-grained matter from land, and estuarine wetlands are an important land-sea transition zone. Frontal tidal flats of estuarine wetlands are the interface of the land-sea exchange of materials, where sediment deposition is the most sensitive [1]. The sediment characters of estuary tidal flats are affected by the changes in runoff and sediment transportation directly from upstream rivers. Regarding the relationship between sediment deposition in estuaries and the hydrological processes of rivers, it is very important to establish the time correspondence between the tidal flat sediment deposition parameters and the time series data of rivers, such as runoff and sediment discharge. Therefore, establishing this correspondence requires the identification of the environmental information carried by the fine sediments deposited in the tidal flats of the estuary [2].

Sample Collection
To obtain information on both the sedimentary history and sediment deposition rate of the tidal flats in the Liaohe estuarine wetland, three vertical sediment columns (C1, C2, C3) were sampled randomly in the central area of the tidal flats using a gravity sediment sampler on 16 May 2017 ( Figure  1). The sampling points were located in the core area of the Liaohe Estuary National Nature Reserve, China. Human activity is not allowed across the entire region, therefore the tidal flat was formed naturally and there has been no human disturbance in the past 50 years or more.
According to the deposition rate of the northern part of Liaodong Bay and Gaizhou Beach, over nearly 30 years [33,34] the sampling depth of the sediment column was determined to be 140 cm. Due to the high moisture content of sediments, it was not appropriate to split the samples on site. Therefore, the sample columns were sent to the laboratory for freezing immediately after they were collected. The samples were cut at regular intervals after 24 h of storage. The surface layer was

Sample Collection
To obtain information on both the sedimentary history and sediment deposition rate of the tidal flats in the Liaohe estuarine wetland, three vertical sediment columns (C1, C2, C3) were sampled randomly in the central area of the tidal flats using a gravity sediment sampler on 16 May 2017 ( Figure 1). The sampling points were located in the core area of the Liaohe Estuary National Nature Reserve, China. Human activity is not allowed across the entire region, therefore the tidal flat was formed naturally and there has been no human disturbance in the past 50 years or more.
According to the deposition rate of the northern part of Liaodong Bay and Gaizhou Beach, over nearly 30 years [33,34] the sampling depth of the sediment column was determined to be Water 2020, 12, 61 4 of 15 140 cm. Due to the high moisture content of sediments, it was not appropriate to split the samples on site. Therefore, the sample columns were sent to the laboratory for freezing immediately after they were collected. The samples were cut at regular intervals after 24 h of storage. The surface layer was subdivided at a resolution of 1.5 cm from 0-60 cm and a resolution of 2 cm from 61-140 cm, which referenced the sample segmentation methods from related studies [20,25]. The divided samples were air-dried according to the corresponding number of layers and the air-drying period was 90 days. In the last 30 days, the sample weights were measured every 10 days, and the constant weight of each sample was regarded as the weight that did not change, and the weight was accurate to 0.01 g. After air-drying, the samples were sieved using a 1.50 mm sieve and then sealed for 137 Cs and 210 Pb activity testing.
The testing of the sediment particle size used the same column that was used to test 137 Cs and 210 Pb. The preparation process was the same as the preparation of the 137 Cs and 210 Pb sample testing.
The samples were sieved by using a 1.50 mm sieve after air-drying and then dispersed for particle size testing. The particle sizes were measured with a Matersizer 3000 laser particle sizer (Malvern, UK) with a particle size range of 0.20-3000.00 µm. The sediment components and particle size parameters were statistically analyzed and categorized using the Folk-Ward method [35].

Sedimentary Dating
The sedimentary history was determined by the 137 Cs and 210 Pb dating techniques. The 137 Cs and 210 Pb activities were performed using a high purity germanium detector gamma spectrometer (BE5030, HPGe, San Ramon, CA, USA) with a system energy resolution of 2.25 KeV. The test time of each sample was 249,456 s to 56,342 s to ensure that the measurement error was controlled below 6%.
(1) 137 Cs dating technique The year was calibrated by analyzing the 137 Cs-specific activity peak. The sediment's deposition time information was determined based on the 137 Cs nuclear accumulation precipitation peak in the Northern Hemisphere, which appeared in 1963, and the Chernobyl nuclear leakage secondary accumulation peak, which appeared in 1986 [36,37].
The deposition rate was calculated using the following formula: In the formula, D r is the deposition rate (cm·year −1 ), T 1 and T 2 are the years specified by the half-life dating method (year); and H is the time slice depth or depth difference between two time horizons (cm).
(2) 210 Pb dating technique The specific activity of 210 Pb ex at the bottom of the sediment sample did not reach the base value. In addition, the study area was an estuary tidal flat, where the sediment mainly came from upstream erosion sediment transportation. Therefore, the CIC model was used to conduct a dating analysis based on 210 Pb ex . The constant initial concentration model of 210 Pb (CIC) assumed that the sediment 210Pb ex has a constant initial concentration in the model. The initial specific activity of 210Pb ex defined at the sediment-water interface was A 0 , and the specific activity of 210 Pb ex at a certain sediment mass depth m was A m . The A m has an exponential decay relationship with the increase of mass depth [38]. The formulas were as follows: Water 2020, 12, 61

of 15
In the formulas, A 0 is the initial specific activity of 210 Pb ex at the sediment-water interface; A m is the specific activity of 210 Pb ex at mass depth m (Bq·g −1 ); λ is the decay coefficient of 210 Pb, 0.03114·year −1 ; m is the mass depth which refers to the cumulative value of sediment above a certain depth (g·cm −2 ); t is deposition time (y); s is the mass deposition rate in time t (g·(cm 2 ·year) −1 ). For the verification of the deposition rate determined by the 137 Cs, the mass deposition rate (g (cm 2 ·year) −1 ) was converted to the deposition rate (cm·y −1 ) according to the sediment bulk density and porosity correction.

Relationship between the Sediment Deposition Characteristics and Runoff
Based on the dating results, the age-related information of the tidal flat sediments in the Liaohe estuary wetland was determined through sedimentation rate inversion, and the time-correspondence relationship between the estuary tidal flat sediments and the Liaohe River runoff was reconstructed. The relationship between the tidal flat sediments and runoff was explored using the time-correspondence relationship. The mean size, median size, clay content and 137 Cs specific activity were selected to characterize the sedimentary features of the estuarine tidal flat sediments. The specific surface area of the sediments was selected as an auxiliary index. The main index of the flow in the flood season, the mean annual flow, and the annual sediment discharge, which are close to the sediment transportation of the Liaohe River, were selected as the key indicators. The above runoff indicator data were obtained from the "Liaohe River Hydrological Yearbook" during the period from 1987-2015. Data from the Liujianfang Hydrological Station were cited in this study, which is the last hydrological station before the Liaohe River inflow to the estuary [39], located in Taian County, Anshan City, China, 63 km from the Liaohe estuary wetland.
The key hydrological indicators were measured and calculated as follows: (1) The flow in the flood season: This refers to the river flow in the relatively concentrated rainfall season of each year, which was calculated according to the daily flow data of the river hydrological observation section. The flood season of Liaohe River Basin ranged from June to September in each year. The formula was as follows: where F f is the flow in the flood season (m 3 ·s −1 ), F i is daily average flow rate (m 3 ·s −1 ), and n is the days of the flood season (D). (2) The mean annual flow rate: This refers to the average value of river flow during each year, which was calculated based on the daily flow data of river hydrological observation sections.
The formula was as follow: where F m is the mean annual flow rate (m 3 ·s −1 ), F i is daily average flow rate (m 3 ·s −1 ), and N is the days of each year (D). (3) The annual sediment discharge: This refers to the total sediment transported through the river hydrological observation section in one year, which was calculated from the daily sediment concentration data of river hydrological observation sections. The formula was as follows: where S is the annual sediment discharge (t·year −1 ), A i is daily average sediment concentration (kg·m −3 ), Q i is daily runoff (m 3 ), and N is the number of days in a year (D).
All the daily data were measured by automatic observation device, and then the mean daily data were automatically generated.
The correlation matrix of the Pearson correlation coefficient was used to analyze the correlations between hydrological indicators of Liaohe River and sediment particle size index. The analysis was completed in Matlab2010b software.

Periodic Relationship between Sediment and Runoff
The cycle reproducibility on the flow in the flood season, the mean annual flow, the annual sediment discharge of the Liaohe River runoff index and the mean and median sizes of the deposited sediments in the Liaohe estuary tidal flats were determined by spectrum analysis, which used a fast Fourier transform to explore the periodicity of the data. According to the spectral analysis results, the frequency corresponding to the peak values in the maximum power, the maximum secondary power and the maximum third power were selected to obtain the periodicity of each index.
The main calculation process is as follows: Step 1. Centralize the original data and perform a fast Fourier transform (FFT) to turn the original data x t into X T ( f ).
Step 2. Calculate the power spectral density.
According to X T ( f ) = A + JB (A is the real part, B is the imaginary part), the power spectral density calculation formula can be simplified as: where f i is the line frequency and T is the length of the data interval.
The line frequency f i = i/T = i/N, i = 0, 1, 2, · · · , N − 1; N is the amount of the data.
Step 3. Drawing the spectrogram The spectrogram will be plotted based on the calculated power spectral density P( f i ) and line frequency f i . The period will be calculated according to the frequency corresponding to the peak value of the spectrum, and the period is the reciprocal of the frequency.

Sedimentary History of the Tidal Flats in the Liaohe Estuary Wetland
(1) 137 Cs dating technique According to the linear quasi-determination of sediment bulk density and sedimentation depth of sediment column C1, C2 and C3, the sediment bulk density shows a linear increasing trend with the increase of sediment depth (R 2 C1 = 0.8391, R 2 C2 = 0.8543, R 2 C3 = 0.8812,). So the 137 Cs profile peak value was suitable for the dating analysis. The average 137 Cs-specific activity of the sediment sample column was 1.10 Bq·kg −1 (Figure 2(C1)). There were two complete accumulation peaks with increasing depth increases: the maximum peak was at 120.00 cm with a 137 Cs-specific activity of 2.47 Bq·kg −1 , and the sub-maximum peak was at 78.00 cm with a 137 Cs-specific activity of 1.52 Bq·kg −1 . Consequently, it was determined that the 137 Cs accumulation peak at 120 cm was 1963, and the 137 Cs accumulation peak at 78 cm was 1986. The sediment deposition rate was 1.83 cm·y −1 from 1963 to 1986, the deposition rate was 2.52 cm·year −1 from 1986 to 2017, and the average deposition rate was 2.22 cm·year −1 from 1963 to 2017.
124.00 cm and the sub-maximum peak was at 84.00 cm. The sediment deposition rate was 1.74 cm·year −1 from 1963 to 1986, the deposition rate was 2.71 cm·year −1 from 1986 to 2017, and the average deposition rate was 2.29 cm·year −1 from 1963 to 2017. Similarly known, the maximum peak of 137 Csspecific activity was at 124.00 cm and the sub-maximum peak was at 84.00 cm based on the data analysis of Figure 2 C3. The sediment deposition rate was 1.65 cm·year −1 from 1963 to 1986, the deposition rate was 2.71 cm·year −1 from 1986 to 2017, and the average deposition rate was 2.26 cm·y −1 from 1963 to 2017.  Likewise, we can see from (Figure 2(C2)) that the maximum peak of 137 Cs-specific activity was at 124.00 cm and the sub-maximum peak was at 84.00 cm. The sediment deposition rate was 1.74 cm·year −1 from 1963 to 1986, the deposition rate was 2.71 cm·year −1 from 1986 to 2017, and the average deposition rate was 2.29 cm·year −1 from 1963 to 2017. Similarly known, the maximum peak of 137 Cs-specific activity was at 124.00 cm and the sub-maximum peak was at 84.00 cm based on the data analysis of Figure 2(C3). The sediment deposition rate was 1.65 cm·year −1 from 1963 to 1986, the deposition rate was 2.71 cm·year −1 from 1986 to 2017, and the average deposition rate was 2.26 cm·y −1 from 1963 to 2017.
(2) 210 Pb dating technique According to the linear quasi-determination of the ln( 210 Pb ex ) value and sedimentation depth of sediment column C1, C2 and C3, the ln( 210 Pb ex ) shows a linear decreasing trend with the increase of sediment depth (R 2 C1 = 0.8985, R 2 C2 = 0.7871, R 2 C3 = 0.8048,). So, the CIC model was suitable for the 210 Pb ex dating analysis. The deposition time of the entire sediment sample column C1 from the surface layer to the bottom layer was 2017 to 1953, which has a 64-year sediment record (Figure 3(C1)). The sediment deposition rate was 2.21cm·y −1 from 1963 to 1986, the deposition rate was 2.50 cm·y −1 from 1986 to 2017, and the average deposition rate was 2.22 cm·year −1 from 1963 to 2017 based on the 210 Pb dating result.
Water 2020, 12, x FOR PEER REVIEW 8 of 15 According to the linear quasi-determination of the ln( 210 Pbex) value and sedimentation depth of sediment column C1, C2 and C3, the ln( 210 Pbex) shows a linear decreasing trend with the increase of sediment depth (R 2 C1 = 0.8985, R 2 C2 = 0.7871, R 2 C3 = 0.8048,). So, the CIC model was suitable for the 210 Pbex dating analysis. The deposition time of the entire sediment sample column C1 from the surface layer to the bottom layer was 2017 to 1953, which has a 64-year sediment record (Figure 3 C1). The sediment deposition rate was 2.21cm·y −1 from 1963 to 1986, the deposition rate was 2.50 cm·y −1 from 1986 to 2017, and the average deposition rate was 2.22 cm·year −1 from 1963 to 2017 based on the 210 Pb dating result.

Relationship between Sediment Deposition in the Tidal Flats of the Liaohe Estuary Wetland and the Runoff from the Liaohe River
There is a consistent correlation between the mean size, median size, clay content, specific surface area and 137 Cs-specific activity of the sediment granularity parameters and the flow in the flood season, the mean annual flow, and the annual sediment discharge from the Liaohe runoff index (Figure 4).

Relationship between Sediment Deposition in the Tidal Flats of the Liaohe Estuary Wetland and the Runoff from the Liaohe River
There is a consistent correlation between the mean size, median size, clay content, specific surface area and 137 Cs-specific activity of the sediment granularity parameters and the flow in the flood season, the mean annual flow, and the annual sediment discharge from the Liaohe runoff index (Figure 4).  According to the measured hydrological data from Liujianfang Hydrological Station on the Liaohe River, 2010, 1998, 1995, 1993, 1991 and 1987 had higher amounts of runoff and sediment transportation. The minimum values of the flow in the flood season, the mean annual flow and the annual sediment discharge were more than 500.00 m 3 ·s −1 , 180.00 m 3 ·s −1 and 5.00 million t·year −1 respectively among these years. The corresponding depth of the deposited sediment in these years had mean sizes greater than 43.00 μm, with an average value of 46.07 μm and a maximum value of 51.13 μm. The corresponding median sizes were greater than 41.00 μm, with an average value of 43.63 μm and a maximum value of 49.03 μm. The corresponding clay contents were less than 3.00%, except in 1991, and the average clay content was 2.97%. The corresponding 137 Cs-specific activities were less than 1.00 Bq·kg −1 , except in 2010, and the average value was 0.85 Bq·kg −1 .
The years 2013, 2008 and 2005 had a medium amount of runoff and sediment discharge. The values of the flow in flood season, the mean annual flow and the annual sediment discharge were in the range of 400.00 to 600.00 m 3 ·s −1 , 100.00 to 180.00 m 3 ·s −1 and 2.00 to 4.30 million t·y −1 respectively among these years. The mean sizes of the corresponding sediment depths in these years were between 42 μm and 46.00 μm, and the average value was 44.03 μm. The corresponding median sizes were between 40.00 μm and 44.00 μm, and the average value was 41.53 μm. The average mean size and median size decreased by 2.04 μm and 2.10 μm, respectively, compared to the average values in the higher runoff and sediment discharge years. The corresponding clay contents were between 3.00% and 4.00%, and the average value was 3.61%, which increased 0.63% compared to the average value in the higher runoff and sediment discharge years. The corresponding 137 Cs-specific activities showed increasing trends in 2013 and 2005, except in 2008, with an average value of 1.29 Bq·kg −1 , which increased 0.32 Bq·kg −1 compared to the average value in the higher runoff and sediment discharge years.
The other years had low flows and sediment discharge. The maximum values of the flow in the flood season, the mean annual flow and the annual sediment discharge were less than 400.00 m 3 ·s −1 , 100.00 m 3 ·s −1 and 2.00 million t·y −1 respectively among those years. The average mean sizes of the corresponding sediment depths in those years was 46.69 μm, and the average median size was 43.83 μm, which increased 0.43 μm and 0.19 μm compared to the average values in the higher runoff and sediment discharge years and increased 2.46 μm and 2.29 μm compared to the average values in the medium runoff and sediment discharge years. The corresponding average clay content was 2.80%, which decreased 0.18% compared to the average value in the higher runoff and sediment discharge years and decreased 0.81% compared to the average value in the medium runoff and sediment discharge years. The corresponding average 137 Cs-specific activity was 1.02 Bq·kg −1 , which increased 0.05 Bq·kg −1 , which is close to the average value in the higher runoff and sediment discharge years, and decreased 0.11 Bq·kg −1 compared to the average value in the medium runoff and sediment discharge years. According to the measured hydrological data from Liujianfang Hydrological Station on the Liaohe River, 2010, 1998, 1995, 1993, 1991 and 1987 had higher amounts of runoff and sediment transportation. The minimum values of the flow in the flood season, the mean annual flow and the annual sediment discharge were more than 500.00 m 3 ·s −1 , 180.00 m 3 ·s −1 and 5.00 million t·year −1 respectively among these years. The corresponding depth of the deposited sediment in these years had mean sizes greater than 43.00 µm, with an average value of 46.07 µm and a maximum value of 51.13 µm. The corresponding median sizes were greater than 41.00 µm, with an average value of 43.63 µm and a maximum value of 49.03 µm. The corresponding clay contents were less than 3.00%, except in 1991, and the average clay content was 2.97%. The corresponding 137 Cs-specific activities were less than 1.00 Bq·kg −1 , except in 2010, and the average value was 0.85 Bq·kg −1 .
The years 2013, 2008 and 2005 had a medium amount of runoff and sediment discharge. The values of the flow in flood season, the mean annual flow and the annual sediment discharge were in the range of 400.00 to 600.00 m 3 ·s −1 , 100.00 to 180.00 m 3 ·s −1 and 2.00 to 4.30 million t·y −1 respectively among these years. The mean sizes of the corresponding sediment depths in these years were between 42 µm and 46.00 µm, and the average value was 44.03 µm. The corresponding median sizes were between 40.00 µm and 44.00 µm, and the average value was 41.53 µm. The average mean size and median size decreased by 2.04 µm and 2.10 µm, respectively, compared to the average values in the higher runoff and sediment discharge years. The corresponding clay contents were between 3.00% and 4.00%, and the average value was 3.61%, which increased 0.63% compared to the average value in the higher runoff and sediment discharge years. The corresponding 137 Cs-specific activities showed increasing trends in 2013 and 2005, except in 2008, with an average value of 1.29 Bq·kg −1 , which increased 0.32 Bq·kg −1 compared to the average value in the higher runoff and sediment discharge years.
The other years had low flows and sediment discharge. The maximum values of the flow in the flood season, the mean annual flow and the annual sediment discharge were less than 400.00 m 3 ·s −1 , 100.00 m 3 ·s −1 and 2.00 million t·y −1 respectively among those years. The average mean sizes of the corresponding sediment depths in those years was 46.69 µm, and the average median size was 43.83 µm, which increased 0.43 µm and 0.19 µm compared to the average values in the higher runoff and sediment discharge years and increased 2.46 µm and 2.29 µm compared to the average values in the medium runoff and sediment discharge years. The corresponding average clay content was 2.80%, which decreased 0.18% compared to the average value in the higher runoff and sediment discharge years and decreased 0.81% compared to the average value in the medium runoff and sediment discharge years. The corresponding average 137 Cs-specific activity was 1.02 Bq·kg −1 , which increased 0.05 Bq·kg −1 , which is close to the average value in the higher runoff and sediment discharge years, and decreased 0.11 Bq·kg −1 compared to the average value in the medium runoff and sediment discharge years.

Sediment Deposition and Runoff Cycles
Based on the previous analysis of the corresponding relationship between the tidal flat sediments and the upstream river runoff and sediment discharge, the cyclical relationship between them will be deeply explored. There is clear cycle reproducibility in all the indexes of flow in the flood season, the mean annual flow and annual sediment discharge of the Liaohe River and the mean particle size and median size of the sediments in the Liaohe estuary wetland tidal flats ( Figure 5 and Table 1). On a 30 year time scale, there was a 14 to 15 year recurrence cycle, a 3 to 4 year recurrence cycle and a 2 to 3 year recurrence cycle of the flow in the flood season. The cycle reproducibility of the mean annual flow and the annual sediment discharge were consistent, with a recurrence cycle of 5 to 6 years, a recurrence cycle of 3 to 4 years and a recurrence cycle of 2 to 3 years. Within a sedimentary depth range of 30 years, there was a recurrence cycle of 14 to 15 years, a recurrence cycle of 7 to 8 years and a recurrence cycle of 4 to 5 years for the mean particle size. There was a recurrence cycle of 14 to 15 years, a recurrence cycle of 5 to 6 years, and a recurrence cycle of 3 to 4 years for the median size. Overall, there was a correspondence cycle of 14 to 15 years between the mean particle size and the median size of the tidal flat sediments in the Liaohe estuary wetland and the flow in the flood season in the Liaohe River. There were 5 to 6 year and 3 to 4 year correspondence cycles between the median size of the tidal flat sediments, the mean annual flow and the annual sediment discharge of the Liaohe River. The cycle of 4 to 5 years in the mean particle size of the tidal flat sediments was close to the cycle of 5 to 6 years in the mean annual flow and the annual sediment discharge of the Liaohe River.

Discussion
In the present study, the deposition rate of the three sampling points was different, and this may be due to the fact that the dating results of 137 Cs and 210 Pb were not absolute dates. Besides, the corresponding sedimentary time at the same mass depth was not the same, and resulted in different  Table 1. Cycle analysis of Liaohe runoff index and sediment particle size parameters of tidal flat sediments in Liaohe estuary wetland.

Indexes of Runoff Discharge
Frequency Meanwhile, according to the analysis of the spectrum power ( Figure 5), the cycle corresponding to the high-power spectrum peak value (5.8 × 10 5 a 2 ·cpd −1 ) of the flow in the flood season was consistent with the cycle corresponding to the high-power spectral peak values (132.75 a 2 ·cpd −1 and 130.41 a 2 ·cpd −1 ) of the mean particle size and median size, which was a 14 to 15 year recurrence period. The cycles corresponding to the high-power spectrum peak values (2.59 × 10 4 a 2 ·cpd −1 and 7.62 × 10 5 a 2 ·cpd −1 ) of the mean annual flow and the annual sediment discharge were consistent with the cycles corresponding to the high-power spectral peak values (30.69 a 2 ·cpd −1 ) of median size, which was a 5 to 6 year recurrence period.

Discussion
In the present study, the deposition rate of the three sampling points was different, and this may be due to the fact that the dating results of 137 Cs and 210 Pb were not absolute dates. Besides, the corresponding sedimentary time at the same mass depth was not the same, and resulted in different calculation results [40,41]. Furthermore, there was a one year difference between the dating results of 137 Cs and 210 Pb, which implied that the dating results of the two methods were basically consistent and credible. Additionally, the average deposition rate of the three sampling points ranged from 2.22 to 2.27 cm·y −1 , with the maximum difference of 0.05 cm·y −1 and an average output value of 0.03 cm·y −1 , which indicated that there was little difference in the deposition rates, and the values were almost the same in the area. Thus, it was reasonable and accurate to use 137 Cs and 210 Pb for dating analysis and sediment deposition rate calculation. Therefore, the average deposition rates of the Liaohe estuary tidal flat determined by two methods were as follows: 2.24 cm·y −1 during the period of 1963 to 2017 and 1.80 cm·y −1 during the period of 1986 to 1963. Furthermore, the sedimentation rate was 2.62 cm·y −1 in the Liaohe estuary wetland tidal flats from 1986 to 2017, which was close to results from the Gaizhou Beach area of the Liaohe estuary and the Liaohe estuary region of Northern Liaodong Bay [33,34]. The results indicated that the deposition rate in the tidal flats in the Liaohe estuary wetland identified in this study was consistent with the deposition rate in the larger regional environment of the Liaohe Delta. In addition, the deposition rate C2 was greater than C1 and C3, which may be caused by extreme rainfall storm events in the upstream basin [42][43][44]. For example, there were extreme flood events of the Liaohe River in 1998 and 1993. The greater runoff capacity could transport a large amount of sediment to the relatively far distance of sampling point C2, and resulted in little deposition in C1, while the larger distance of C3 lead to relatively little sediment. Therefore, there existed differences in the spatial distribution of the deposition rate within the area.
The sediment discharge from the Liaohe River was negatively correlated with the clay content of the tidal flat sediments, and positively correlated with the mean particle size and median size of the sediment, and significantly negatively correlated with the clay content and 137 Cs-specific activity of the sediments on a nearly 30-year time scale (Table 2). According to the results of related research, it was believed that the amount of fine sediment particles deposited on the estuary tidal flat was directly related to the amount of sediments transported by the river and inversely proportional to the ability of rivers to transport fine matter [26]. In the Evrotas River research, in Greece, the sediment transportation to the downstream estuary was significantly reduced due to the decrease in runoff from the upstream river during the dry season [45]. It was shown that the sediments deposited in the Yangtze River estuary, China are significantly reduced in the context of the decline of the upper reaches of the river runoff and sediment discharge [46]. The studies have shown that the changes of river runoff directly affected the amount of sediment transported to the estuary, which was consistent with the findings of this study. Other research has shown that the proportion of coarse-grained sand in sediment carried by river runoff increased, and the ratio of clay content decreased as the flow of the upstream river increased, resulting in the increase of sediment particle size in an estuary in Fujian, Southeastern China [24]. In the Yellow River estuary research, in China, the estuary sediment particle size has an obvious coarsening trend due to the increase in upstream runoff and sediment transport by the implementation of the water and sediment adjustment project in the Yellow River [47,48]. In general, the fine-grained clay content tended to decrease, and the sediment particle size showed a coarsening trend of the estuary wetland tidal flats with the increase in sediment discharge of rivers, which were consistent with the findings of this study. There are few studies discussing the cyclical relationship between tidal flat sediments in estuary wetlands, river runoff and sediment transport. Some studies have shown that the sedimentary volume in the Yangtze River estuary, in China, and the runoff in the basin have a 9 to 11 year cycle correspondence [49]. The mean particle size of the sediments in the Pearl River estuary had a 2 to 4 year cycle correspondence with the Pearl River sediment transport in China [50]. However, due to notable differences between the hydrological characteristics of the Yangtze River, the Pearl River and the Liaohe River, the relevant results have little significance for the study of sediment deposition in the tidal flats of the Liaohe estuary wetland. Other studies have shown that river runoff has a significant cyclical effect on the tidal constant of the estuary. There are 2.50, 5.50 and 6.50 years of common movement and interdecadal cycle relationships in the La Plata estuary in Spain [51]. The study in Spain has not yet determined whether the river runoff and sediment transportation have a cyclical effect on the fine-sediment deposition in the estuary. The results in this research indicated that the long recurrence times of the flow in the flood season in the Liaohe River affected the high frequency recurrence of the grain size parameters of the estuary tidal flats. The long recurrence times of the mean annual flow and the annual sediment discharge affected the low frequency reproduction of the sediment particles. There is a 14 to 15 year cycle correspondence between the mean particle size and the median size of the tidal flat sediments in the Liaohe estuary wetland and the flow in the flood season in the Liaohe River. There is a 5 to 6 year cyclical and a 3 to 4 year cyclical correspondence between the median size of the tidal flat sediments and the mean annual flow and the annual sediment discharge of the Liaohe River.

Conclusions
In the present study, the dating results of 210 Pb ex were consistent with that of 137 Cs, which indicated that the determined sediment dating results were credible. The average sedimentation rate in the center of the Liaohe estuary wetland tidal flats was 2.24 cm·y −1 over the past 50 years. The sedimentation rates from 1963 to 2017 and 1963 to 1986 were 2.62 cm·y −1 and 1.80 cm·y −1 , respectively. Besides, this study clarified that the amount of fine-sediment particles deposited on the estuary tidal flat was directly related to the amount of sediments transported by the river and inversely proportional to the ability of rivers to transport fine matter. After the runoff sediment transport capacity was reduced to a certain extent, the flow in flood season, the mean annual flow and the annual sediment discharge ranged from 400.00 to 600.00 m 3 ·s −1 , 100.00 to 180.00 m 3 ·s −1 and 2.00 to 4.30 million t·y −1 , respectively. Furthermore, soil particle size and median size decreased with the decreasing of annual sediment discharge, while the clay content of sediment fine particles increased with decreasing of annual sediment discharge. The sediment carrying capacity was limited when the flow in the flood season was less than 400.00 m 3 ·s −1 and the mean annual flow was less than 100.00 m 3 ·s −1 . The total amount of fine sediments transported to the estuary tidal flats was reduced due to the lack of sediment transportation power of the runoff. In addition, there was a 14 to 15 year corresponding periodic relationship between the mean particle size and the median size of the tidal flat sediments and the flow in the flood season in the Liaohe River. There was a 5 to 6 year and a 3 to 4 year corresponding periodic relationship between the median size of the tidal flat sediments and the mean annual flow and the annual sediment discharge of the Liaohe River. The results will provide theoretical support for revealing the sedimentation and erosion balance of estuarine sediments and the stable maintenance of tidal flats.