Improvement for the Multi-scale Periodic Characteristics Revealing of Precipitation Signals and Its Impact Assessment on Soil Hydrological Process by Combining Hht and Cwt Approaches

This study conducts a detailed analysis of multi-scale periodic features involved in the annual and seasonal precipitation signals at the typical coastal reclamation region in China by selecting the suitable continuous wavelet transform (CWT) and innovatively combining the improved Hilbert–Huang transform (HHT), and further deduces the precipitation trend and its impacts on the future soil hydrological process. The Morlet wavelet transform is proved suitable in revealing the precipitation signals broad-scale period-icities, however, the critical mode mixing problem in CWT causes the poor significance in the fine-scale periodicities, which cannot well match the climate background. By combining the HHT approach, the fine-scale mode mixing drawback in CWT is effectively eliminated, and the the studied precipitation signals multi-scale periodicities are accurately revealed. Consequently, an overall decreasing trend of annual and seasonal precipitation in future years is demonstrated. Furthermore, by novelly using the cross wavelet transform (XWT) and wavelet transform coherence (WTC) approaches, the prominent correlations between the precipitation dynamics and soil and groundwater salinities dynamics, it is demonstrated that the precipitation increase can effectively leach the surface soil salt downwards into the deeper soil layers and groundwater with 5–7-day lag in the new cultivated tidal land. The revealed future decreasing trend of precipitation , especially in spring and summer, may aggravate the soil salinization at the coastal reclamation region, thus some reasonable salt leaching and evaporation suppression measures need to be taken to prevent the possible soil secondary salin-ization process.


Introduction
The precipitation dynamics always plays a prominent role in the soil hydrological process, which impacts significantly on the health of a water-soil-crop system.In particular, at the studied coastal region in northern Jiangsu, China, a huge amount of tidal lands are being enclosed for agricultural utilization, which has aroused international attention, e.g., the planned average annual reclamation amount of about 7000 hectares tidal lands during the year 2010-2020 at Dongtai County takes up about 40 % of the total amounts in China (Wang and Ke, 1997;Meng and Tang, 2010).These new areas commonly need 5-10 years of soil amelioration and desalting before they can be cultivated for wheat and corn in a rain-fed farming system.Thus, the precipitation dynamics at the coastal region play a key role in both the new-land soil desalting process and the cultivated-land soil secondary salinization process (Zhao et al., 2010), that makes the precise study of the precipitation dynamics very important for the agriculture development and environment health in these new areas.However, the uneven precipitation distribution during seasons and years, and the special and stochastic climatic and environmental factors, such as typhoon and storm tide, at the coastal region make the study hard to be precise enough and usually bring some environmental risks, such as the waterlog Published by Copernicus Publications on behalf of the European Geosciences Union.and salinization.In view of this, more effective approaches need to be introduced into the precise analysis of precipitation multi-scale characteristics and the accurate forecast of precipitation trend and its impacts on the soil hydrological process.
Since the hydrometeorological signals are highly nonstationary and their physical processes often operate under a large range of scales varying from 1 day to several decades, the continuous wavelet transform (CWT) has been introduced and developed to reveal the multiple time scales characteristics of these nonstationary signals (Mandelbrot and Wallis, 1968;Grossmann and Morlet, 1984;Mallat, 1989;Daubechies, 1994;Torrence and Compo, 1998;Labat et al., 2000).Because the CWT can provide information about both time and frequency simultaneously, and enables a separation to be made between features associated with different characteristic length scales, it has some advantages over the traditional Fourier transform and has been widely used in revealing the periodic characteristics of hydrometeorological signals at multiple time scales (Walker, 1997;Nener et al., 1999;Coulibaly and Bum, 2004;Labat et al., 2005;Shark and Yu, 2006;Liu et al., 2009;Yu et al., 2013).For hydrometeorological signals, the Morlet and Mexican hat wavelets are two types of the most popularly used continuous wavelets (Foufoula-Georgiou and Kumar, 1995).To well and truly assess the statistical significance of hydrometeorological signals multi-scale periodicities, the wavelet power spectrum and its confidence test are usually requested to be combined with the wavelet transform analysis (Ge, 2007;Zhang et al., 2010).In view of the above background, it is feasible to utilize the CWT approaches to effectively reveal the precipitation's multi-scale characteristics at the studied coastal reclamation region, although little similar research has been carried out at this area.
One major problem in CWT is the frequent appearance of mode mixing (Gong et al., 2005;Huang and Wu, 2008), which is defined as the CWT decomposed information of one scale (or frequency) either consisting of information of other widely disparate scales (or frequencies) or the decomposed information of a similar scale residing in different scales information.Mode mixing is often a consequence of signal intermittency, which may cause serious aliasing in the timefrequency distribution and make the signal's multi-scale periodicities unclear.To overcome the mode mixing problem in CWT, the presented study innovatively introduces the improved Hilbert-Huang transform (HHT) approach, based on the Ensemble Empirical Mode Decomposition (EEMD) algorithm (Wu and Huang, 2009), into the revealing capability improvement of hydrometeorological signals multi-scale periodic characteristics.As a new noise-assisted data analysis method, the EEMD algorithm is proposed based on the original HHT, which is based on the original Empirical Mode Decomposition (EMD) method (Huang et al., 1998(Huang et al., , 1999)), to overcome the similar mode mixing drawback in EMD as that in CWT.The improved EEMD method defines the true intrinsic mode function (IMF) components as the mean of an ensemble of trials, each consisting of the signal plus a white noise of finite amplitude.Although adding noise may result in smaller signal-to-noise ratio, the added white noise will provide a relatively uniform reference scale distribution to facilitate EMD.Therefore, the low signal-to-noise ratio does not affect the decomposition method but actually enhances it to avoid the mode mixing problems.The HHT is an adaptive time-frequency data analysis method and has been proved quite versatile in a broad range of applications for extracting signals from data generated in noisy nonlinear and nonstationary processes.
One crucial issue in assessing the impact of precipitation dynamics on the future soil hydrological process is to precisely reveal the correlation between the precipitation dynamics and soil and groundwater salinities dynamics.However, the classic regression analysis cannot well reveal the correlation, because the large numbers of zero value data in the daily hydrometeorological signals highly increase the signals noisy and nonlinear dynamics (Yu et al., 2013).Thus, the cross wavelet transform (XWT) and wavelet transform coherence (WTC) approaches are novelly introduced into the study to discuss the resonance feature between the precipitation and salinity signals at multiple time scales.By combining the continuous wavelet analysis and cross spectrum analysis (Roth, 1971;Liu, 1994), the XWT and WTC approaches have been developed and utilized to expose the common power, relative phase and significant coherence of two signals in time-frequency space (Torrence and Webster, 1999;Grinsted et al., 2004).Because of the embedded wavelet theory, the XWT and WTC approaches have advantages over the single cross spectrum analysis in revealing the cross information of nonstationary signals about both time and frequency simultaneously and have been popularly used for meteorological, geophysical, mechanical and electronic signals (Gurley et al., 2003;Jevrejeva et al., 2003;Yang et al., 2009;Yu et al., 2014).
In the remainder of this paper, the studied sites and data details and the methodologies of CWT, XWT, WTC and HHT are presented.In the next section, the results of multi-scale periodic characteristics of different precipitation signals are revealed by selecting the suitable CWT type, and the future precipitation trends at the corresponding dominant scales are deduced and forecasted with the comparison with the traditional linear fitting results.In the discussion section, the mode mixing problems in CWT are discussed by analyzing the climate background of precipitation's multi-scale periodicities and improved by innovatively combining the EEMD based HHT approach.Furthermore, the impacts of the forecasted future precipitation dynamics on the future soil hydrological process are assessed by quantitatively analyzing the correlation between the precipitation and salinization dynamics, where the XWT and WTC approaches are novelly used to reveal the oscillation coherence and time lag property between the precipitation and groundwater salinity signals at daily scales.Finally, the conclusions are drawn.

Studied sites and data
Three meteorological stations in the coastal region of northern Jiangsu, China are selected as the study sites in this paper (Fig. 1).The Ganyu station (34.84The monitoring site (32.65 • N, 120.90 • E) is located at a newly cultivated farmland at the typical reclamation area Dongtai County, which is about 5 km from the Yellow Sea.This tidal land was enclosed from the sea in 2004 and first cultivated for wheat and corn in June 2011.The soil electrical conductivity (EC) data in the different soil layers (10, 50 and 100 cm) were monitored 6 times each month from June 2011 to May 2013.The daily groundwater EC data and the corresponding daily precipitation data covering 1 whole year (365 days) from 11 May 2012 to 10 May 2013 were monitored and collected.

Continuous wavelet transform (CWT) and wavelet power spectrum test
The CWT is a mathematical tool which allows the decomposition of the signal f (t) in terms of elementary contributions called wavelets, which can be thought of as a packet of sine waves of varying amplitude and wavelength (Sadowskey, 1996;Torrence and Compo, 1998;Labat et al., 2000).For time series f (t) ∈ L 2 (R), R is the domain of real numbers, or finite energy signal, the CWT of signal f (t) with the analyzing wavelet φ is the convolution of f (t) with a set of dilated and translated wavelets: where a is the dilation (scale or frequency) parameter, b is the translation (position or time) parameter, φ(t) is the complex conjugate functions of φ(t), and δt is the time interval of data series.In this study, two popularly used wavelets, the Mexican hat and Morlet wavelets (Eq.2), were selected as the φ (Mallat, 1989;Daubechies, 1994;Torrence and Compo, 1998).The Mexican hat wavelet, as a real-value wavelet function, is the second derivative of the Gaussian function.
The global wavelet power spectrum is defined as where N is the length of data.The statistical significance of wavelet power can be assessed relative to the null hypotheses that the signal is generated by a stationary process with a given background power spectrum (P ) (Torrence and Compo, 1998).For the global wavelet power spectrum test, the theoretical spectrum P is defined as P = σ 2 P k x 2 v /v, where σ 2 is the variance of data series, x 2 v is the chi square when degrees of freedom is v at the requested confidence level, and P k is the Fourier white noise or red noise power spectrum.P k is defined as Eq. ( 3) where r(1) is the autocorrelation coefficient under a first order lag of data series, and T is the period of data series.Theoretically, for Morlet wavelet, T ≈ 1.033a, while for Mexican hat wavelet, T ≈ 3.974a (Torrence and Compo 1998).When r(1) > 0.1, P k is red noise spectrum, and when r(1) ≤ 0.1, P k is white noise spectrum and equals the default value 1.0 (Allen and Smith, 1996;Ge, 2007).The v values in Mexican hat and Morlet wavelet power spectrums are calculated by

Cross wavelet transform (XWT) and wavelet transform coherence (WTC)
The XWT of two signals x(t) and y(t) with finite energy is defined as where Z v (p) is the confidence level associated with the probability p for a probability density function (PDF) defined by the square root of the product of two chi-square distributions.In this study 0.05 significance level was selected and Z 2 (95 %)= 3.999 (Grinsted et al., 2004).Because the wavelet is not completely localized in time, the XWT has edge artifacts and a cone of influence (COI) in which edge effects cannot be ignored is introduced.Here the COI is taken as the area in which the wavelet power caused by a discontinuity at the edge has dropped to e −2 of the value at the edge (Grinsted et al., 2004).The phase difference between the components of two signals, revealed from the complex angles of W xy (a, b), can be estimated by calculating the mean and confidence interval of the phase difference.The circular mean of the phase over regions with higher than 5 % statistical significance that are outside the COI are used to quantify the phase relationship.The circular mean of a set of angles (α i , i = 1...n) is defined as (Zar, 1999) The circular standard deviation, reflecting the scatter of angles around the mean, is defined as The circular standard deviation is analogous to the linear standard deviation, in that it varies from zero to infinity.In some cases there might be reasons for calculating the mean phase angle for each scale, and then the phase angle can be quantified as a number of days.
Analysis of the wavelet transform coherence (WTC) of two signals x(t) and y(t) is another useful measure to reveal the high common wavelet power and the localized correlation coefficient of the two signals in time frequency space.The wavelet coherence is defined by Torrence and Webster (1999) as where S is a smoothing operator.For the Morlet wavelet a suitable smoothing operator is given by Torrence and Webster (1999): where S scale denotes smoothing along the wavelet scale axis, S time denotes smoothing in time, c 1 and c 2 are normalization constants and is the rectangle function.The statistical significance level of WTC is estimated using Monte Carlo methods (Grinsted et al., 2004).

Hilbert-Huang transform (HHT)
In the original HHT, based on the Empirical Mode Decomposition (EMD) approach, the data series x(t) is decomposed in terms of truth, intrinsic mode functions (IMFs), and trend (R) components (Huang et al., 1998(Huang et al., , 1999)): where c j is the IMF component, and R n is the residue of data x(t) after n number of IMFs are extracted.
As proposed by Wu and Huang (2009), the first step of Ensemble Empirical Mode Decomposition (EEMD) process is to add a white noise series w(t) to the targeted data.The added white noise is treated as the possible random noise that would be encountered in the measurement process.Under such conditions, the ith "artificial" observation will be (10) In the case of only one observation, each multipleobservation ensemble is mimicked by adding not arbitrary but different realizations of white noise (w i (t)) to that single observation as given in Eq. ( 11).The second step is to decompose the data with added white noise into IMFs.Repeat step 1 and step 2 again and again, but with different white noise series each time, and at last obtain the ensemble means of corresponding IMFs of the decompositions as the final result.As the ensemble number approaches infinity, the truth, c j (t), as defined by EEMD, is where R j k (t) is the contribution to the j th IMF from the added white noise of the kth trial of the j th IMF in the noiseadded signal.The amplitude of noise w i (t) is not necessarily small, usually setting an amplitude of 0.2 of that of the standard deviation of the corresponding data.But the ensemble number of the trials N has to be large, usually not lower than 100.
The HHT marginal spectrum represents the accumulated amplitude (energy) over the entire data span in a probabilistic sense and offers a measure of the total amplitude (energy) contribution from each frequency (period) value, serving as an alternative spectrum expression of the traditional Fourier spectrum.The marginal spectrum is defined as (Huang and Wu, 2008) where [0, T ] is the temporal domain within which the data is defined, H expresses the amplitude in terms of a function of time and frequency, and ω is the instantaneous frequency.All the CWT, XWT, WTC and HHT in this paper were carried out in MATLAB.

Precipitation trend forecast based on the conventional linear fitting method
The average coefficients of variation for the monthly precipitation during each year reached 122.1, 95.4 and 75.6 % at the coastal Ganyu, Dongtai and Shanghai meteorological stations, respectively, and for the annual precipitation among years reached 24.5, 24.0 and 18.4 %, respectively.The prominent uneven precipitation distribution during seasons and years in the studied area and a slightly decreasing trend of this kind of unevenness with the latitude dropping were revealed.The summer precipitation made up 61.0, 51.0 and 41.4 % of the annual precipitation at Ganyu, Dongtai and Shanghai, respectively, with the significant correlation coefficients of 0.88, 0.84 and 0.77, respectively, between the summer and annual precipitation at the 0.05 confidence levels.Further more, the 10-year running means of summer precipitation signals showed the most similar dynamic curves with the annual precipitation signals (Fig. 2).The above analysis indicated that the summer precipitation dynamics played a dominant role in the annual precipitation changes in the studied area.As one conventional method for assessing the precipitation dynamics, the original precipitation data series are usually treated by linear fitting, in which the 10 times of the slope of fitting line, i.e., the climate tendency rate, denotes the precipitation changes of every 10 years (Xu et al., 2004;Wang et al., 2006).A negative climate tendency rate denotes a decreasing precipitation trend, while a positive one denotes an increasing trend.The fitting equations for the original annual and seasonal precipitation data series, together with the fitting equations for the 10-year running means data, at the three studied stations were computed and shown in Table 1.Results showed that the climate tendency rates of annual precipitation at Ganyu, Dongtai and Shanghai were −8.83, −13.21 and 12.88 mm/10 years, respectively, which indicated the decreasing trends of annual precipitation at Ganyu and Dongtai and the increasing trend at Shanghai.The climate tendency rates of seasonal precipitation at different meteorological stations were calculated in the same way.However, it should be noted that most of the fitting equations based on the original precipitation data did not pass the statistical test at 0.05 confidence level.Although the statistic significance of fitting equations based on the 10-year running means data increased slightly, the overall performance of precipitation trend forecast based on the conventional analysis was too poor to be used for the accurate precipitation dynamics assessment.In view of this, the popular wavelet analysis method continuous wavelet transform (CWT) was introduced into the precipitation trend forecast by revealing the periodic characteristics of precipitation signals at multiple temporal-frequency scales.

Multi-scale periodic characteristics of annual and seasonal precipitation signals based on CWT
The Morlet and Mexican hat wavelet transform coefficients of annual and seasonal precipitation signals at different meteorological stations were calculated by MATLAB.Taking the annual precipitation at Dongtai as an example (Fig. 3), the Morlet and Mexican hat wavelet transform coefficients clearly indicated the distribution of the annual precipitation signal at multiple time scales.The solid isograms indicated Results showed that the Morlet wavelet transform coefficients of Dongtai annual precipitation signal generated five kinds of quasi periodic oscillations (QPOs), namely at 20-, 14-, 7-, 5-and 2-year time scales (Figs.3a and 4a).The corresponding average changing periodicities (T ) of Dongtai annual precipitation signal (i.e., the average cycle years between each two time domains with positive wavelet coefficients) at the five time scales were 24, 16, 9, 5 and 2 years, respectively, which were obtained by calculating and averaging the year numbers of each two neighboring high and low precipitation periods.Among the five QPOs, the 24-and 16year periods passed the 0.005 confidence test, and the 9-year period passed the 0.2 confidence test, but the 5-and 2-year periods showed very poor significances.With respect to the Mexican hat wavelet transform, only one QPO at an 8-year time scale was revealed and its corresponding period T was 31 years, which passed the 0.2 confidence test.The multiscale periodic characteristics for the Dongtai seasonal precipitation signals and for the annual and seasonal precipitation signals at the Ganyu and Shanghai stations based on the Morlet and Mexican hat wavelet transforms were obtained in the same way, and were shown in Tables 2-4.It should be noted that the global wavelet power spectra of the Shanghai autumn and winter precipitation series were tested using the red noise model because their r(1) were 0.13 and 0.16, respectively.
Comparing the two types of CWTs performances, the Morlet wavelet transform showed more significant advantage than the Mexican hat in capturing the oscillatory behaviors and revealing the local features of the studied precipitation signals periodicities, especially prominent in increasing the confidence level of the wavelet power spectrum test.In view of this, the dominant precipitation periods with the highest confidence level revealed from the Morlet wavelet transform were selected to further deduce and forecast the annual and seasonal precipitation trends in the studied area.

Precipitation trend forecast based on CWT
According to the dominant periods of annual and seasonal precipitation signals revealed by the Morlet wavelet transforms, the precipitation trends in future years were forecasted by deducing the future fluctuating of the wavelet transform coefficient curves at the corresponding dominant time scales.Taking the Dongtai precipitation forecast as an example, two dominant periods of 24 and 16 years, both of which passed 0.005 significance tests, were revealed in both the annual and summer precipitation signals, and the dominant periods of 16, 23 and 17 years were revealed in the spring, autumn and winter precipitation signals, respectively.The corresponding wavelet coefficient curves for the annual and seasonal precipitation at different dominant time scales were shown in Fig. 5.Each single cycle in the wavelet coefficient curve was constituted by a high-precipitation half-period, in which the wavelet coefficients were positive, and a low-precipitation half-period, in which the coefficients were negative.Results in Fig. 5 showed that all the wavelet coefficients of annual and seasonal precipitation in 2011 located at the end of each high-precipitation half-period, which indicated that the annual and summer precipitation in the next 8-12 years, the spring precipitation in the next 8 years, the autumn precipitation in the next 12 years, and the winter precipitation in the next 8 years would probably be in the relatively low periods.Similarly, the future precipitation trends at the Ganyu and Shanghai stations were deduced and forecasted.Results showed that the annual and summer precipitation in the next 8-10 years, the autumn precipitation in the next 6 years and the winter precipitation in the next 7 years at the Ganyu station would probably be in the relatively low periods, because all their wavelet coefficients at the corresponding dominant periods in 2011 located at the end of each high-precipitation half-period.With respect to the Ganyu spring precipitation, it would still be in a relatively low period in the next 2 years and then revert to a relatively high period in the next 3 years, because the wavelet coefficient in 2011 located at the first year of one low-precipitation half-period.At the Shanghai station, all the annual and summer precipitation in the next 5-12 years, the spring precipitation in the next 12 years and the winter precipitation in the next 8 years would probably be in the relatively low periods.With respect to the Shanghai autumn precipitation, it would still be in a relatively high period in the next 2 years and then revert to a relatively low period in the next 12 years.The above forecast showed an overall decreasing trend of precipitation in the next future years in the studied coastal region.

Improvement of multi-scale periodic characteristics revealing capabilities by combining continuous wavelet transform (CWT) and Hilbert-Huang transform (HHT) methods
The previous studies (Zhu and Wang, 2001;Wei and Zhang, 2009;Zhu et al., 2011) have reported the significant decadal quasi periodic oscillation (QPO) of summer precipitation at about 30 years and interannual QPOs of summer precipitation at about 2 and 8 years in both the Huaihe River valley and the Yangtze River lower reach in China.The climate background is the corresponding decadal and interannual QPOs characteristics in the East Asian summer monsoon circulation (Chang et al., 2000).Among the three studied stations in the coastal region, the Ganyu station is located at the north part of Huaihe River Valley, the Shanghai station is located at the south part of the Yangtze River Lower Reach, and the Dongtai station is located at the crossing part of the two valleys.Based on the popular CWT method, the summer precipitation signals in the studied stations show significant broad-scale periodicities at approximate 9 years and 25 years, which match the approximate 8-and 30-year scale QPO characteristics in the two valleys precipitation.However, the summer precipitation periodicity at fine scale in the studied stations, revealed by the CWT, shows very poor significances; this does not match the significant interannual QPO characteristic at about a 2-year scale in the two valleys precipitation.The main reason causing the poor performance of CWT method on revealing the fine-scale periodicity in the coastal precipitation is the mode mixing problem in the CWT process, e.g., the precipitation QPO at a fine time scale is often nested in a broad time scale QPO in Fig. 3.In addition, the precipitation dynamic in the coastal region is often influenced by some unique climatic and environmental factors, such as the typhoon and the storm tide.These factors may lead some stochastic impacts on the periodic characteristics in the coastal precipitation and decrease the significance of fine-scale periodicities, which may aggravate the mode mixing in the CWT.
Since the improved HHT method, which is based on the ensemble empirical mode decomposition (EEMD) algorithm, has been popularly used in recent years in decomposing data into multi-scale IMF components with prominent advantage in eliminating the mode mixing drawback in the traditional Empirical Mode Decomposition (EMD), the following study tries to introduce the improved HHT method into eliminating the mode mixing problem in the CWT analysis of the studied precipitation signals multi-scale periodic characteristics.
For each summer precipitation signal, the statistical results of five decomposition components obtained from the EEMD are displayed in Table 5.The corresponding HHT marginal spectra based on EEMD are shown in Fig. 6a3, b3 and c3).In the EEMD process, an ensemble size of 100 is used, and the added white noise in each ensemble member has a standard deviation of 0.2.Among the five decomposition components of each summer precipitation signal, IMF1-IMF4 are intrinsic mode function components at multiple frequencies and R5 is a trend component.For each summer precipitation signal, the fist decomposed IMF1, which denotes an average period at approximately 2 years, has the highest variance contribution ratios and the most significant correlation coefficients with the original data.For the corresponding HHT marginal spectrum of each summer precipitation signal, the most significant amplitude (or energy) accumulation occurs at the period of approximately 2 years.Further more, the EEMD-based HHT marginal spectra of the annual and the other seasonal precipitation signals (Fig. 6) show similar and significant fine-scale periodicities as the summer precipitation signal.Comparing the CWT method, the EEMD-based HHT method well improves the capabilities in eliminating the fine-scale mode mixing problem and in revealing the precipitation signals fine-scale periodic characteristics.As for the broad-scale periodicities analysis, the HHT results, that the annual and seasonal precipitation signals have approximate 9-and 25-year periodicities, match the CWT results.However, the significances in the broad-scale periodicities revealed by the HHT are slightly lower than that revealed by the CWT.
The above discussion proves that it is reasonable to accurately analyze the precipitation signals multi-scale periodicities by combining the CWT and HHT methods, i.e., using the EEMD-based HHT method to eliminate the fine-scale mode mixing drawbacks in the CWT and to improve the signal's multi-scale periodic characteristics.In view of this, the precipitation trend forecast based on the CWT should not ignore the precipitation signals fine-scale periodicities.In the case study, the three stations annual and seasonal precipitation trends at fine scales can be forecasted by deducing the future fluctuating of the wavelet transform coefficient curves at the corresponding periods of 2-3 years.Results show an overall decreasing trend of the precipitation in the next future 1-2 years in the studied stations, because all the finescale wavelet coefficients in 2011 locate at the end of each high-precipitation half period.

Applying research results of precipitation multi-scale time-frequency characteristics to soil hydrological process study
Before evaluating the impacts of coastal-region multi-scale precipitation trends on the soil hydrological process, the relationship between the precipitation dynamics and soil and groundwater salinization dynamics needs to be quantitatively analyzed first, which is shown in the following study by taking the relatively short monitoring data analysis in a newly cultivated farmland at the typical reclamation area Dongtai as an example.
The comparison results between the dynamics of monthly precipitation and monthly average soil electric conductiv-ity (EC), standing for soil salinity, at different soil layers are shown in Fig. 7a.Results show that the surface soil EC (10 cm layer) decreases significantly with the precipitation increase and increases with the precipitation decrease, with a significant correlation coefficient of 0.76 under the 0.01 confidence level (P = 0.004) by the logarithmic function fitting (Fig. 7b).On the contrary, see Fig. 7a, the 50 and 100 cm layers soil EC dynamics are in keeping with the precipitation dynamics.The positive correlation coefficient of 0.63 between the 50 cm layer soil EC and precipitation is significant under the 0.05 confidence level (P = 0.029) by the logarithmic function fitting, while the positive correlation coefficient of 0.44 between the 100 cm layer soil EC and precipitation is not significant under the 0.05 confidence level (P = 0.173) by the logarithmic function fitting.It seems that the soil salt is leached from the surface part to lower part when the precipitation increases.Further analysis shows that in average every 10 mm monthly precipitation increase in the monitoring site causes correspondingly 0.1-0.2dS m −1 decrease of the surface soil EC (Fig. 7b).
Because of the shallow groundwater table, usually changing from 80 to 150 cm layer over the whole year at the studied coastal reclamation region, the groundwater salinity dynamics has very close relationship with the precipitation change.The daily precipitation and groundwater EC dynamics based on the original monitoring data covering 1 whole year are shown in Fig. 8a.Results show that the groundwater EC usually increases significantly when the precipitation increases and leaches soil salt downwards.While, when the precipitation is low and evaporation is high, the salt in groundwater moves upwards into soil and the groundwater EC decreases.Based on the daily precipitation and groundwater EC data series, the cross wavelet power spectrum from cross wavelet transform (XWT) and the corresponding wavelet transform coherence (WTC) results are shown in Fig. 8b and c, respectively.Since the r(1) of daily groundwater EC and precipitation series are 0.98 and 0.30, respectively, the red noise test in XWT is chosen.Prominent resonant periods of the two signals around the 5-, 15-, 25-and 70-day bands stand out as being significant at the 0.05 significance level.The average squared wavelet coherence coefficient between the two signals around the 70-day band is the highest and equals to 0.88, which proves the most significant oscillation coherence between the two series at the periods around 70-day band.In view of this, the XWT phase angles within the 0.05 significant regions around the 70-day band are reasonable to be calculated.The mean phase −34.1 ± 5.3 • (where ± designates the circular standard deviation) at the periods around 70-day band is in-phase, which indicates about 5-7-day lag of the impact of precipitation dynamics on the groundwater EC dynamics.
The above discussion indicates that the precipitation increase at the coastal reclamation area can be very helpful for the surface soil salt leaching, while may increase the salt accumulation in the deeper soil profile and groundwater, which asks for reasonable drainage measures especially in rainy seasons to prevent the coastal groundwater salinization.The revealed future decreasing trend of precipitation, especially in spring and summer, in the studied coastal region has negative impact on the surface soil salt leaching, thus may aggravate the salinization in the salt-affected coastal land and probably lead a new secondary salinization process in the newly desalted and cultivated land in future years.This suggests reasonable salt leaching and evaporation suppression measures to prevent soil salinization, especially the root-zone salinization, in future years.
It should be noted that the precipitation trend forecast based on CWT and HHT in this study is qualitative.Thus, the further quantitative evaluation of precipitation trend and its impacts on soil and groundwater salinization in future research may need more advanced mathematical models and algorithms, such as discrete wavelet transform and neurofuzzy models (Yang et al. 2013), and have to pay more attention to the impacts of stochastic factors.In addition, because the soil and groundwater salinities monitoring site in this study is located at a newly cultivated farmland, the longer time series of soil and groundwater salinities monitoring data need to be collected in the future research for the improvement of quantitative analysis between the precipitation and salinization dynamics.

Conclusions
This study conducts a detailed analysis of multi-scale periodic characteristics involved in the annual and seasonal precipitation signals at different meteorological stations along the coastal reclamation region by using and comparing two types of popular continuous wavelet transforms (CWTs) and their wavelet power spectrums.Furthermore the future annual and seasonal precipitation trends in the studied coastal stations are deduced and forecasted by selecting a suitable wavelet to reveal the dominant precipitation periods.The critical fine-scale mode mixing problem in the CWT process is discussed and innovatively improved by combining the EEMD-based Hilbert-Huang transform (HHT) method.In addition, by analyzing the quantitative relationship between the precipitation dynamics and soil and groundwater salinities dynamics, some novel attempts are made to evaluate the impact of the revealed precipitation trend on the future soil hydrological process in a newly cultivated coastal farmland.Some interesting conclusions can be drawn as follows.pecially in increasing the confidence level of wavelet power spectrum test.
2. Comparing the traditional method of judging the precipitation trend by the climate tendency rate, which shows relatively poor significances, the CWT-based forecast approach by deducing the future wavelet coefficients dynamics at the dominant precipitation period with the highest significance level is proved more reasonable.By means of this approach, a general decreasing trend in the coastal annual and seasonal precipitation in future years is demonstrated.
3. Although the superiorities of CWT method in revealing the precipitation signals broad-scale periodicities and forecasting the precipitation trend are proved, the critical fine-scale mode mixing problem in CWT causes the poor significances of fine-scale periodicities, which cannot well match the previous studies on the whole valleys precipitation periodicities and on the climate backgrounds.Thus, the EEMD-based HHT method is innovatively introduced into the effective elimination of fine-scale mode mixing drawback in CWT.It is demonstrated reasonable to accurately analyze the precipitation signals multi-scale periodicities by combining the CWT and HHT methods.
4. The demonstrated prominent correlations between the precipitation dynamics and soil and groundwater salinities dynamics indicate that the precipitation increase can effectively leach the surface soil salt downwards, but may aggravate the salt accumulation in the deeper soil layers and groundwater, which asks for reasonable drainage measures in rainy seasons to prevent the possible groundwater salinization.By innovatively using the cross wavelet transform (XWT) and wavelet transform coherence (WTC) approaches, the significant coherence around the 70-day resonant period between the precipitation dynamics and groundwater salinity dynamics is revealed, which indicates the approximate 5-7day lag of the impact of precipitation dynamics on the groundwater salinity.Thus, the revealed future decreasing trend of precipitation, especially in spring and summer, may aggravate the soil salinization at the studied coastal reclamation region, which asks for reasonable salt leaching and evaporation suppression measures to prevent the future soil secondary salinization process.
In addition, some ideas about the further improvements of quantitative evaluation of precipitation trend and its impacts on the soil hydrological process and of longer monitoring data collection in the new cultivated land are proposed.

Figure 1 .
Figure 1.Map showing the locations of the studied meteorological stations and monitoring site.
a, b), where * denotes complex conjugation.The cross wavelet power spectrum is defined as W xy (a, b) 2 = W xy (a, b)W * xy (a, b) .The theoretical distribution of the cross wavelet power of two signals with background power spectrum P x k and P y k is given by Torrence and Compo (1998) as D W xy (a, b)

Figure 3 .
Figure 3. Real parts of the Morlet (a) and Mexican hat (b) wavelet transform coefficients of the annual precipitation at Dongtai during 1953-2011.

Figure 4 .
Figure 4. Global wavelet power spectra and their confidence tests using the white noise model for the Morlet (a) and Mexican hat (b) wavelet transforms of the annual precipitation at Dongtai during 1953-2011.

Figure 5 .
Figure 5. Morlet wavelet transform coefficients for the Dongtai annual and summer precipitation signals at the 24-year (a, c) and 16year (b, d) dominant periods, for the spring precipitation at the 16-year dominant period (e), for the autumn precipitation at the 23-year dominant period (f), and for the winter precipitation at the 17-year dominant period (g).

Figure 7 .
Figure 7. Monthly dynamics of the precipitation and soil EC of different layers in the new reclaimed farmland(a), and their correlation (b).

Figure 8 .
Figure 8. Original precipitation and groundwater EC data series (a), and their cross wavelet transform (b) and the corresponding average squared wavelet coherence coefficients (c).The thick black contour in (b) denotes the 0.05 significance level against the red noise test, the cone of influence where the edge effects might distort the picture is shown as a lighter shade, and the relative phase relationship is shown as the arrows (with in-phase pointing right, anti-phase pointing left).

1.
By combining the CWT and the corresponding wavelet power spectrum test, the multi-scale periodic characteristics, mainly the broad-scale periodicities, of annual and seasonal precipitation signals in the studied coastal reclamation region are revealed.Comparing the performances of two popular CWTs, the Morlet wavelet transform is proved more suitable than the Mexican hat in the studied area, because the Morlet shows more advantage in capturing the oscillatory behaviors and revealing the local feature of coastal precipitation periodicity, es-

valley and is about 5 km from the Yellow Sea. The Shanghai station (31.16 • N, 121.43 • E) is located at the south part of the lower reach of Yangtze River and is about 27 km from the East China Sea.
• N, 119.12 • E) is located at the north part of Huaihe River

Table 1 .
Fitting equations of the annual and seasonal precipitation series at the Ganyu, Dongtai and Shanghai stations based on the original data and 10-year running means data (y denotes the annual or seasonal precipitation (mm), x denotes the time series (1, 2, 3,...), r and P denote the correlation coefficient and P value in the regression model.

Table 2 .
Annual and seasonal precipitation periods and their significances at different time scales based on the Morlet and Mexican hat wavelet transforms at the Ganyu station.

Table 3 .
Annual and seasonal precipitation periods and their significances at different time scales based on the Morlet and Mexican hat wavelet transforms at the Dongtai station.

Table 4 .
Annual and seasonal precipitation periods and their significances at different time scales based on the Morlet and Mexican hat wavelet transforms at the Shanghai station.

Table 5 .
Statistical results of the multi-scale time-frequency characteristics from the EEMD of summer precipitation signals.