1 引言
海平面变化是气候变化的重要参数之一,受大气圈、水圈、岩石圈以及人类活动等多要素的综合作用,时空特性鲜明,变化机理复杂[1, 2, 3],海平面变化预测是海平面变化研究乃至全球研究中的热点也是难点问题[4, 5, 6, 7]。目前国际上关于器测时期海平面变化预测可分为数理统计预测和动力模型预测两种方法。对于大尺度如全球性海平面变化,可利用动力模型根据影响海平面变化的驱动因子(如全球平均气温)及作用机制进行过程清晰的长期预测[8]。然而对于区域影响要素复杂多样、物理过程及作用机制均不够清晰的区域性海平面变化研究,动力模型预测精度有待提高。数理统计预测是基于海平面过去时间变化序列进行建模和外推,具有资料可靠,精度高的优点,是现阶段区域海平面变化预测的有效途径[9]。
早期统计预测通常假定海平面过去变化特征在未来一段时期保持一致,线性拟合其趋势加以外推,是最简单可行的方法[10],但预测结果也相对较粗略。随着傅利叶变换、小波分析等时频分析方法被逐渐应用于海平面变化研究,验潮站序列被分解为确定性趋势项、周期项以及剩余随机序列[11, 12]。预测时将趋势项及周期项外推,并对剩余随机序列进行线性/非线性拟合预测[13]。近年,针对傅里叶变化基于固定时频窗以及小波分析的基函数选择难度大等问题,Huang等[14]提出新的时频分析方法——经验模态分解(Empirical Mode Decomposition,简称EMD),依据各周期信号振幅频率特征不同将原始数据分解为一系列具有不同特征尺度的序列,即本征模函数(Intrinsic Mode Function,简称IMF)。EMD方法提取的信号具有很高的信噪比,目前已被广泛应用于地学、水文、气象等领域[15, 16, 17]。王卫强等[18]利用EMD方法对中国沿岸海平面变化进行研究,并与灰色系统方法和随机动态预测模型比较,认为EMD方法不但能够得到较精确的年变率,而且能够分离出各种天文周期,较好地解决了资料序列的非平稳和非线性问题。集合经验模态分解(Ensemble Empirical Mode Decomposition,简称EEMD)是基于EMD方法的改进,利用噪声辅助减弱模态混叠现象的影响从而改善分解效果[19]。然而傅利叶变换和EEMD等时频解析方法始终是对现状海平面变化的分解研究,以此直接外推生成未来时间变化序列其精度有待提高。
另一方面,随着时间序列预测模型日益发展成熟,神经网络作为一门新兴边缘学科,因能较好解决非平稳序列间的非线性映射被广泛应用于地学、海洋、大气等领域[20, 21, 22, 23, 24]。BP(Back Propagation)神经网络是其中应用最为广泛、发展最为成熟的模型之一[25, 26, 27, 28]。但神经网络基于纯数据驱动预测,建模过程中无法考虑时间序列本身可能存在的物理意义及规律,无法解析海平面变化的过程及特征。因此,建立一种新的思路,打破海平面时间序列的时频分析方法及纯数据驱动预测方法的局限性,在科学解析过去海平面变化过程及特征的基础上,准确预估其变化规律是当前国际关于区域海平面变化研究需要解决的重要科学问题[29, 30]。
珠江三角洲位于中国广东沿海,区域内人口、财富密集,受海平面变化威胁严重[31, 32, 33]。三角洲内海平面在全球及区域尺度要素等共同影响下表现出相当的复杂性和不确定性,国内外已开展大量研究并取得一定成果[34, 35, 36, 37]。黄镇国等[10]基于香港及澳门站验潮序列预测2030年珠江三角洲海平面将上升30cm; 冯伟等[38]利用卫星测高技术研究发现1993~2009 年南海海平面变化为 5.5±0.7mm/a,明显高于同期全球平均海平面变化趋势(3.3±0.4mm/a)。然而以往关于海平面变化预测的研究通常利用现状时间序列进行线性拟合外推或者基于时频分析方法进行周期外推,研究精度有待进一步提高[39]。因此为明晰海平面变化过程并科学预测其趋势,本研究拟采用EEMD结合BP神经网络建立EEMDBP预测模型,以研究区内标准海洋站闸坡站验潮数据为例,对月均海平面高度序列进行拟合预测。 一方面分析海平面变化的规律,另一方面探讨提高海平面变化预测精度的思路,为区域海平面长期变化预测研究提供参考借鉴。
2 数据方法介绍 2.1 资料来源处理本研究收集全球海平面监测网络(Global Sea Level Observing System,简称GLOSS)标准站闸坡站1959~2011年海平面月均变化序列资料作为实证研究对象,对海平面长期变化进行预测研究。闸坡站位于中国南部广东沿海(21°35′N,111°50′E; 图1),是中国南海验潮站中历史最长、记录最完整的验潮站序列[40]。本文采用三次样条插值方法补齐序列中的缺测数据(1~2个月)[41]。
![]() |
图1 闸坡站位置图 Fig.1 Location of Zhapo tide gauge |
采用EEMD-BP模型进行海平面变化趋势预测主要包括时间序列的分解和预测两部分。 首先对原始海平面时间序列进行EEMD分解,提取海平面变化过程具有物理含义的信号并降低噪音影响,得到不同尺度周期或趋势(IMFs),实现序列的平稳化; 继而采用BP神经网络对分解后的各IMFs分别建立预测模型研究其未来变化趋势; 最后将各IMFs的预测结果重构,得到原始序列的预测值。
2.2.1 EEMD方法采用EEMD提取海平面变化序列中的周期项和趋势项,并经显著性检验,确定海平面变化隐含的周期及最终趋势,具体方法如下[42, 43]: 首先在原始序列中K次加入均值为0,幅值标准差为常数的高斯白噪声ni(i=1,2,…,K)得到时间序列fi(t)。然后对fi(t)进行EMD分解: 找出序列中所有极大值和极小值点,分别用三次样条函数合成上下包络线,并得到平均包络线mi1,将fi(t)减去mi1得到去掉低频趋势的hi1,多次重复上述过程,使平均包络线趋于零从而得到第一个IMF分量ci1(t),代表原始序列中最高频的分量,有:
![](PIC/20150212-M1.jpg)
对剩余序列ri1(t)重复上述分解过程,直到j次分解后所得剩余分量满足给定的结束标准[14],分解结束。
则原始序列可表示为:
![](PIC/20150212-M2.jpg)
根据不相关随机序列的统计值为0的原理[44],将上述步骤对应的IMF进行总体平均运算,以消除多次加入高斯白噪声对真实IMF的影响。
2.2.2 BP神经网络方法BP神经网络的预测思路是认为未来值与现状值之间隐含某种函数关系,通过建立神经网络来拟合这种关系并推导出未来值[45]。典型的BP网络包括输入层、隐含层(中间层)和输出层(见 图2),建模过程由信息的正向传播和误差的反向传播两部分组成。输入层将接收到的输入信息传递到隐含层,将信息处理变换后传递到输出层输出处理结果。当实际输出与预期输出不符时,进入误差的反向传播阶段。误差通过输出层,按误差梯度下降的方式修正各层权值,向隐含层、输入层逐层反传。神经网络通过反复学习训练,进行多次信息正向传播和误差反向传播,不断调整各层权值,直至输出误差减小到预设值,或者达到预设的学习次数。有关BP神经网络的基本原理已在众多文献中介绍[46, 47],本文不再赘述。
![]() |
图2 神经网络模型示意图 Fig.2 Structure of BP neural network |
闸坡站验潮资料经过EEMD分解生成9个IMFs,如下 图3所示,IMF1~IMF9依次对应8个周期项和1个长期趋势项。周期项分别包括3个月、7个月、13个月、26个月、49个月、101个月、220个月和314个月周期,基本对应常用理论天文周期[48]。其中,前3个周期分别对应季节,半年和1年的变化周期。26个月变化周期为准两年大气振荡,49个月周期与3~7年ENSO周期变化相符。101个月周期为准9年的半交点周期,与黄镇国等[49]关于广东沿岸珠江三角洲验潮记录的周期分析一致; 而220个月周期可能表现的是潮汐变化周期,与天文学上黄白交点运行的18.61年周期[50]相近。
![]() |
图3 闸坡站验潮资料的EEMD分解结果 Fig.3 Results of the EEMD of Zhapo tide gauge records |
为检验IMFs的信息的有效性,对分解结果进行白噪音检验,如下 图4所示,IMF1~3、IMF5~7均通过了置信度为95 % 的显著性检验,说明这6个IMFs具有统计学意义,它们分别对应的季节变化、半年变化、年变化、ENSO事件、准9年半交点周期及18.61天文年变化等6个周期项变化显著。未通过检验的其余2个IMFs则被认为是噪音序列,在实际预测中应予以剔除。从各IMF的方差贡献率来看,季节、半年和年周期分别占海平面变化的30.71 % ,42.88 % 以及20.05 % ,说明年内变化对中国南海平面变化影响较大,与前人[51]研究关于中国海平面变化以年内周期为主一致。同时根据趋势项IMF9的线性拟合显示,1959~2011中国南海近53年年变化率为1.6mm/a,与陈特固等[52]关于中国南海海域珠江三角洲海平面变化研究较一致。结果显示EEMD能较完整分解出海平面变化的常规周期项,提取的长期变化趋势合理,分解效果较好。
![]() |
图4 EEMD分解结果显著性检验 Fig.4 Significance test of IMFs from the EEMD of Zhapo tide gauge records |
考虑到海平面变化具有强烈的季节变化特征,本研究采用对应月时间序列组织数据建模,并通过数据序列滑动平移构建样本。原始资料01/1959~12/2011共636个月海平面高度资料被分为两部分,其中前576个资料(01/1959~12/2006)作为训练样本,后60个资料(01/2007~12/2011)作为检验样本[53]。经过反复调试,建立隐含层为1,输入层神经元为30,输出层为1(如输入层为1959~1988年每年1月海平面高度,输出层则为第31年1989年的1月海平面高度)的神经网络预测模型[54]。训练过程中逐步调整隐含层的节点数目,得到最佳隐含层节点数为9。在模型的隐含层中选用S型正切函数(tansig)作为输入层至隐含层的传递函数; 选用线性函数(purelin)作为隐含层与输出层间的传递函数; learngdm函数作为学习函数; trainlm函数作为训练函数以及mes函数作为性能函数。研究将学习率确定为0.02,并采用0.003作为误差边界值,即当误差值小于0.003时,认为学习完成,停止计算,输出结果,其他参数采用系数默认值。关于BP神经网络的建模预测过程在Matlab(R 2012a)中完成。
EEMD-BP建模预测即将分解所得IMF1~9分别按上述步骤进行建模预测。IMF1~9预测值与样本资料对比如下 图5所示。从 图5中可以看出,高频分量IMF1( 图5a)的拟合效果最差,但序列整体振幅及走向与样本变化较一致。由IMF2至IMF9,随着各IMF变化趋于平稳,预测效果逐渐变好。说明神经网络虽然适用于非平稳非线性序列预测,但当序列趋向平稳时,预测精度也随之上升。考虑到研究海平面长期变化趋势时,通常应先去除季节等年内高频振荡信号[55],因此,IMF1拟合效果不甚理想等对于长期如年代际海平面变化预测研究影响不大。
![]() |
图5 EEMD分解各IMF预测值与样本对比 Fig.5 Comparison between the prediction and samples of IMFs from EEMD |
为全面评价EEMD-BP模型预测精度,本文同时采用BP神经网络直接建模方法进行对比预测,即不经过EEMD分解,直接利用BP神经网络对闸坡站原始序列进行预测。为避免因降噪处理等造成与原始观测数据比较时产生偏差,进行预测精度对比时将全部IMF1~9预测结果叠加重构,形成EEMD-BP模型预测的最终结果。BP神经网络与EEMD-BP模型的预测结果见 图6,其中蓝色曲线代表观测值,红色曲线代表EEMD-BP模型预测结果,绿色曲线代表BP神经网络预测结果。 图6中可以看出3条曲线的整体趋势较一致,两种预测方法的拟合精度都较高。其中绿色曲线更易出现极值偏差,如在2007年9月,BP神经网络的预测值明显高于观测样本,而EEMD-BP预测结果则与观测样本较接近。另外在2009年11月及2010年11月,绿色曲线都明显偏离红蓝曲线。一方面,红色与蓝色曲线的拟合较绿色与蓝色曲线拟合的整体效果更好,即采用EEMD-BP模型预测比单独使用BP神经网络预测的效果好; 另一方面,BP神经网络预测过程中与原始观测值偏离较大处,即预测误差偏大的点较易出现在时间序列变化的极值点上,如图6中2008年8~10月及2011年8~10月。由此可知,尽管建模过程已通过学习过去变化规律对未来趋势进行模拟预测,但因海平面高度复杂多变,存在大量突变及极值点,相对于长期变化的趋势预测研究,对极值的预测处理相对要差一些,可另开展极值频率、重现期等预测研究。
![]() |
图6 BP神经网络预测及EEMD-BP模型预测结果比较 Fig.6 Comparison of the predictions from BP and EEMD-BP models |
为进一步分析两种方法的拟合精度及预测效果,本文采用相关系数(R)、平均误差(Mean Error,简称ME)和均方根误差(Root Mean Square Error,简称RMSE)对预测结果进行评价,见 表1。可以看出,两组预测值与检验样本均有较高的相关系数(R>0.75),较低的均方根误差(RMSE<37mm)及平均误差(ME<4mm)。Cheng等[56]利用验潮站数据及卫星资料预测英国沿海海平面变化,预测值与观测值最小RMSE为43.9mm,本文整体预测效果较理想。EEMD-BP模型预测结果的相关系数高于BP神经网络模型,而均方根误差及平均误差绝对值则低于BP神经网络模型,总体表现优于BP神经网络。因此,采用EEMD-BP模型预测,首先对数据进行平稳化及降噪处理,依据信号的不同变化特征进行模拟预测优于直接使用BP神经网络对原始数据进行预测。
![]() |
表1 BP神经网络预测与EEMD-BP神经网络预测误差统计比较 Table 1 Error statistics of the predictions from BP and EEMD-BP models |
考虑数据的有效性,本研究采用EEMD-BP模型,选用IMF2~3、IMF5~7及趋势项IMF9对珠江三角洲未来2050年海平面变化进行预测,结果如 图7所示。其中,蓝色为闸坡站1959~2011年海平面高度观测值,线性拟合结果显示该时期珠江三角洲海平面变化呈现上升趋势,速率约为2.22mm/a。 图7中红色曲线为预测值,线性拟合显示至2050年,上升速率将达2.72mm/a。
![]() |
图7 基于EEMD-BP模型的珠江三角洲1959~2050年海平面变化预测 Fig.7 Past sea level observation and sea level prediction from 1959 to 2050 in the Pearl River delta by EEMD-BP |
根据《2010年中国海平面公报》[57]显示,近30年广东沿海海平面以2.5mm/a速率上升,并预测未来30年将上升75~145mm。本文以闸坡站为代表站点研究珠江三角洲海平面变化,与以上研究相近,预测结果可靠。
4 结论本文利用珠江三角洲闸坡站1959~2011年月均海平面高度观测资料,结合EEMD及BP神经网络建立EEMD-BP预测模型进行海平面长期变化趋势预测研究,并将预测结果与直接使用BP神经网络预测对比,分析研究模型预测精度及适用性,在此基础上进行了珠江三角洲2050年海平面变化预测。研究结果表明:
(1)BP神经网络技术因其无需先验知识等优势,可以从纯数据驱动的出发,较好地对海平面高度时间序列进行学习训练、预测,适应海平面变化这类多因素、非平稳和非线性时间序列预测研究。但对于极值模拟预测能力有限,可另行开展极值、重现期等研究。
(2)EEMD-BP模型首先对原序列进行EEMD分解,得到具有不同时频特征的周期项及趋势项之后,再利用BP神经网络对各IMFs建立预测模型,预测精度高于直接使用BP神经网络预测。对各IMFs预测结果可看出,序列越趋于平稳,BP神经网络预测效果越好。EEMD-BP在对数据进行平稳化及降噪处理的基础上完成预测,模型构建合理可行,模型精度显著提高。
(3)闸坡站1959~2011年海平面高度时间序列经EEMD分解得到季节变化、半年变化、年变化、ENSO事件、准9年半交点周期及18.61天文年变化等6项周期项。其中,前3个周期的方差贡献率总和达93.64 % ,说明珠江三角洲海平面变化以季节、半年及年等年内周期变化为主。
(4)闸坡站1959~2011年海平面高度资料线性拟合结果显示该时段珠江三角洲以2.22mm/a速率上升,EEMD-BP模型预测至2050海平面上升速率增大至2.72mm/a,有加速上升趋势。
1 | 俞肇元, 袁林旺, 闾国年等. 西北太平洋边缘海区海面变化多尺度解析及空间分异. 地理研究, 2009, 28 (6):1644-1655 Yu Zhaoyuan, Yuan Linwang, Lü Guonian et al. Multi-scale analysis of sea-level change and its spatial characteristics in northwest Pacific Ocean marginal seas. Geographical Research, 2009, 28 (6):1644~1655 |
2 | 肖笃宁, 韩慕康, 李晓文等. 环渤海海平面上升与三角洲湿地保护. 第四纪研究, 2003, 23 (3):237-246 Xiao Duning, Han Mukang, Li Xiaowen et al. Sea level rising abound Bo Hai Sea and deltaic wetlands protection. Quaternary Sciences, 2003, 23 (3):237~246 |
3 | 谢志仁, 袁林旺. 略论全新世海面变化的波动性及其环境意义. 第四纪研究, 2012, 32 (6):1065-1077 Xie Zhiren, Yuan Linwang. Fluctuation characteristics of Holocene sea-level change and its environmental implications. Quaternary Sciences, 2012, 32 (6):1065~1077 |
4 | Barnett T P. The estimation of global sea-level change——A problem of uniqueness. Journal of Geophysical Research: Oceans, 1984, 89 (5):7980~7988 |
5 | Tilburg C E, Garvine R W. A simple model for coastal sea level prediction. Weather and Forecasting, 2004, 19 (3):511~519 |
6 | Suzuki T, Hasumi H, Sakamoto T T et al. Projection of future sea level and its variability in a high-resolution climate model:Ocean processes and Greenland and Antarctic ice-melt contributions. Geophysical Research Letters, 2005, 32 (19):1~5 |
7 | Slangen A B A, Katsman C A, van de Wal R S W et al. Towards regional projections of twenty-first century sea-level change based on IPCC SRES scenarios. Climate Dynamics, 2012, 38 (5-6):1191~1209 |
8 | Wild M, Calanca P, Scherrer S C et al. Effects of polar ice sheets on global sea level in high-resolution greenhouse scenarios. Journal of Geophysical Research:Atmospheres, 2003, 108 (D5):1~10 |
9 | Ayyub B M, Braileanu H G, Qureshi N. Prediction and impact of sea level rise on properties and infrastructure of Washington, DC. Risk Analysis, 2012, 32 (11):1901~1918 |
10 | 黄镇国, 张伟强, 吴厚水等. 珠江三角洲2030年海平面上升幅度预测及防御方略. 中国科学(D辑), 2000, 30 (2):202-208 Huang Zhenguo, Zhang Weiqiang, Wu Houshui et al. Prediction of sea level rising amplitude in 2030 and defensive countermeasures in the Pearl River delta. Science in China(Series D), 2000, 30 (2):202~208 |
11 | 欧素英, 陈子燊. 小波变换在相对海平面变化研究中的应用. 地理科学, 2004, 24 (3):358-364 Ou Suying, Chen Zishen. Application of wavelet transform to relative sea level research. Scientia Geographica Sinica, 2004, 24 (3):358~364 |
12 | 段华琛. 长江口附近验潮站海平面资料的周期图分析和周期信号检验. 海洋科学, 1982,(3):6-10 Duan Huachen. Analysing periodgrams and inspecting periodic signal of sea level deta supplied by the tide gauge stations near the mouth of the Yangtze River. Marine Sicences, 1982,(3):6~10 |
13 | 夏华永, 李树华. 带周期项的海平面变化灰色分析模型及广西海平面变化分析. 海洋学报, 1999, 21 (2):9-17 Xia Huayong, Li Shuhua. A grey model with periodic term for sea level analysis and its application to Guangxi coast. Acta Oceanologica Sinica, 1999, 21 (2):9~17 |
14 | Huang N E, Shen Z, Long S R et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings of the Royal Society A-Mathematical Physical and Engineering Sciences, 1998, 454 (1971):903~995 |
15 | 冯 颖. 经验模分解方法在卫星测高数据中的应用. 青岛: 国家海洋局第一海洋研究所硕士论文, 2011. 2-3 Feng Ying. Analysis of Satellite Altimeter Data Using Empirical Mode Decomposition and Hilbert-Huang Transformation. Qingdao:Master Thesis of the First Institute of Oceanography, State Oceanic Administration, 2011. 2~3 |
16 | 张衍广, 林振山, 李茂玲等. 基于EMD的中国生态足迹与生态承载力的动力学预测. 生态学报, 2008, 28 (10):5027-5032 Zhang Yanguang, Lin Zhenshan, Li Maoling et al. Prediction on the dynamics of ecological footprint and ecological capacity based on EMD in China. Acta Ecologica Sinica, 2008, 28 (10):5027~5032 |
17 | 杨 周, 林振山, 俞鸣同. 东亚冬季风演变和亚洲内陆干旱化信号的多尺度分析. 第四纪研究, 2011, 31 (1):73-80 Yang Zhou, Lin Zhenshan, Yu Mingtong. Multi-scale analysis of East Asian winter monsoon evolution and Asian inland drying force. Quaternary Sciences, 2011, 31 (1):73~80 |
18 | 王卫强, 陈宗镛, 左军成. 经验模态法在中国沿岸海平面变化中的应用研究. 海洋学报, 1999, 21 (6):102-109 Wang Weiqiang, Chen Zongyong, Zuo Juncheng. Empirical mode decomposition method for study of the mean sea level variation along the China coast. Acta Oceanologica Sinica, 1999, 21 (6):102~109 |
19 | Wu Z, Huang N E. Ensemble empirical mode decomposition:A noise-assisted data analysis method. Advances in Adaptive Data Analysis, 2009, 1 (1):1~41 |
20 | 薛联青, 崔广柏, 陈凯麒. 非平稳时间序列的动态水位神经网络预报模型. 湖泊科学, 2002, 14 (1):19-24 Xue Lianqing, Cui Guangbai, Chen Kaiqi. Dynamic water level neural network forecast model on non-stationary time series. Journal of Lake Sciences, 2002, 14 (1):19~24 |
21 | Malmgren B A, Kucera M, Nyberg J et al. Comparison of statistical and artificial neural network techniques for estimating past sea surface temperatures from planktonic foraminifer census data. Paleoceanography,2001, 16 (5):520~530 |
22 | 李双成, 郑 度. 人工神经网络模型在地学研究中的应用进展. 地球科学进展, 2003, 18 (1):68-76 Li Shuangcheng, Zheng Du. Applications of artificial neural natworks to geosciences:Review and prospect. Advance in Earth Sciences, 2003, 18 (1):68~76 |
23 | 郭庆春, 李 力, 张 冉等. 我国气温变化的神经网络预测模型. 热带气象学报, 2009, 25 (4):483-487 Guo Qingchun, Li Li, Zhang Ran et al. A neural-network-based forecasting model for temperature change in China. Journal of Tropical Meteorology, 2009, 25 (4):483~487 |
24 | 田光进, 刘纪远, 庄大方等. 基于遥感与GIS的20世纪90年代中国城镇用地时空特征. 第四纪研究, 2003, 23 (4):421-427 Tian Guangjin, Liu Jiyuan, Zhuang Dafang et al. The tmeporal-spatial characteristics of urban land in China in 1990s by remote sensing and GIS. Quaternary Sciences, 2003, 23 (4):421~427 |
25 | Wang R J, Feng X X, Xia Y J et al. A Back Propagation Artificial Neural Network Prediction Model of the gate freeze time for Injection Molded Polypropylenes. Journal of Macromolecular Science Part B-Physics, 2013, 52 (10):1414~1426 |
26 | 曾 文, 张训械, 胡 雄等. 利用人工神经网络建立我国南海地区电离层模式. 地球物理学报, 1999, 42 (5):581-589 Zeng Wen, Zhang Xunxie, Hu Xiong et al. An ionospheric model in the South China Sea region made by artificial neural network. Chinese Journal of Geophysics, 1999, 42 (5):581~589 |
27 | 匡文慧, 刘纪远, 邵全琴等. 区域尺度城市增长时空动态模型及其应用. 地理学报, 2011, 66 (2):178-188 Kuang Wenhui, Liu Jiyuan, Shao Quanqin et al. Dynamic urban growth model at regional scale and its application. Acta Geographica Sinica, 2011, 66 (2):178~188 |
28 | 张 赛, 廖顺宝. 多年平均气温空间化BP神经网络模型的模拟分析. 地球信息科学学报, 2011, 13 (4):534-538 Zhang Sai, Liao Shunbao. Simulation and analysis of spatialization of mean annual air temperature based on BP neural network. Journal of Geo-information Science, 2011, 13 (4):534~538 |
29 | Obeysekera J, Park J, Irizarry-Ortiz M et al. Probabilistic projection of mean sea level and coastal extremes. Journal of Waterway Port Coastal and Ocean Engineering, 2013, 139 (2):135~141 |
30 | Jevrejeva S, Moore J C, Grinsted A. Sea level projections to AD2500 with a new generation of climate change scenarios. Global and Planetary Change, 2012, 80-81:14~20 |
31 | 李 阔, 李国胜. 广东沿海地区风暴潮易损性评估. 热带地理, 2011, 31 (2):153-158,177Li Kuo, Li Guosheng. A storm surge vulnerability assessment for coastal Guangdong. Tropical Geography, 2011, 31 (2):153-158,177 |
32 | 曾昭璇, 丘世钧. 珠江口海平面上升对三角洲的影响. 人民珠江, 1994,(1):6-9 Zeng Zhaoxuan, Qu Shijun. Influence of rising sea level of Pearl River estuary on the delta. Pearl River, 1994,(1):6~9 |
33 | 董艺辛, 刘春莲, 陈木宏等. 珠江三角洲中部大鳌平原晚第四纪古生物记录及环境演化. 第四纪研究, 2012, 32 (6):1183-1198 Dong Yixin, Liu Chunlian, Chen Muhong et al. Late Quaternary paleontology and environmental changes in Da'ao Plain, middle Pearl River delta, Southern China. Quaternary Sciences, 2012, 32 (6):1183~1198 |
34 | 陈特固, 时小军, 余克服. 近50年全球气候变暖对珠江口海平面变化趋势的影响. 广东气象, 2008, 4 (2):1-3 Chen Tegu, Shi Xiaojun, Yu Kefu, The impacts of climate warming on sea-level trends at Pearl River estuary during 1957-2006. Guangdong Meteorology, 2008, 4 (2):1~3 |
35 | 董德化, 何焯霞, 黄希敏等. 珠江三角洲局部水位壅高及河口治理的几个问题. 人民珠江, 1999,(5):15-17 Dong Dehua, He Zhuoxia, Huang Ximin et al. On several problems in local rise of water level in the Pearl River dleta and improvement of Estuary. Pearl River, 1999,(5):15~17 |
36 | Huang Z G, Zong Y Q, zhang W Q. Coastal inundation due to sea level rise in the Pearl River delta, China. Natural Hazards, 2004, 33 (2):247~264 |
37 | 邵 磊, 赵 梦, 乔培军等. 南海北部沉积物特征及其对珠江演变的响应. 第四纪研究, 2013, 33 (4):760-770 Shao Lei, Zhao Meng, Qiao Peijun et al. The characteristics of the sediment in northern South China Sea and its response to the evolution of the Pearl River. Quaternary Sciences, 2013, 33 (4):760~770 |
38 | 冯 伟, 钟 敏, 许厚泽. 联合卫星重力、 卫星测高和海洋资料研究中国南海海平面变化. 中国科学: 地球科学, 2012, 42 (3):313-319 Feng Wei, Zhong Min, Xu Houze. Sea level variations in the South China Sea inferred from satellite gravity, altimetry, and oceanographic data. Science China:Earth Science, 2012, 42 (3):313~319 |
39 | 黄立人. 海面变化趋势的动态预测. 海洋通报, 1991, 10 (1):1-6 Huang Liren. Dynamical prediction of sea level change trend. Marine Science Bulletin, 1991, 10 (1):1~6 |
40 | 何 蕾, 李国胜, 李 阔等. 1959 年来珠江三角洲地区的海平面变化与趋势. 地理研究, 2014, 33 (5):988-1000 He Lei, Li Guosheng, Li Kuo et al. Changes and trends of sea level in the Pearl River delta in the last 50 years. Geographical Research, 2014, 33 (5):988~1000 |
41 | Church J A, White N J, Coleman R et al. Estimates of the regional distribution of sea level rise over the 1950-2000 period. Journal of Climate. 2004, 17 (13):2609~2625 |
42 | 冯 雯, 李可军, 翟 倩等. 中国北方地区(33°-41°N, 108°-115°E) 过去500年夏季降水量的周期成分分析. 第四纪研究, 2013, 33 (1):126-135 Feng Wen, Li Kejun, Zhai Qian et al. Multi-scale analysis of the precipitation index in north-central China(33°-41°N, 108°-115°E) over the past 500 years. Quaternary Sciences, 2013, 33 (1):126~135 |
43 | 郭志永, 翟秋敏. 北大西洋IODP U1313站深海沉积中环境敏感粒度组分研究. 第四纪研究, 2012, 32 (6):1107-1120 Guo Zhiyong, Zhai Qiumin. A study on environmentally sensitive grain-size groups of deep-sea sediment at IODP site U1313, North Atlantic. Quaternary Sciences, 2012, 32 (6):1107~1120 |
44 | 陈仁祥, 汤宝平, 马婧华. 基于EEMD的振动信号自适应降噪方法. 振动与冲击, 2012, 31 (15):82-86 Chen Renxiang, Tang Baoping, Ma Jinghua. Adaptive de-noising method based on ensemble empirical mode decomposition for vibration signal. Journal of Vibration and Shock, 2012, 31 (15):82~86 |
45 | 刘 莉, 叶 文. 基于BP神经网络时间序列模型的降水量预测. 水资源与水工程学报, 2010, 21 (5):156-159 Liu Li, Ye Wen. Precipitation prediction of time series model based on BP artificial neural network. Journal of Water Resources & Water Engineering, 2010, 21 (5):156~159 |
46 | Pradhan B, Lee S. Regional landslide susceptibility analysis using Back-Propagation Neural Network model at Cameron Highland, Malaysia. Landslides, 2010, 7 (1):13~30 |
47 | Filippo A, Torres A R, Kjerfve B et al. Application of Artificial Neural Network(ANN)to improve forecasting of sea level. Ocean & Coastal Management, 2012, 55:101~110 |
48 | 田素珍, 马继瑞. 太平洋月平均水位变化中的周期成分. 海洋学报, 1996, 18 (4):11-18 Tian Suzhen, Ma Jirui. Period elements of month average tide change in the Pacific Ocean. Acta Oceanologica Sinica, 1996, 18 (4):11~18 |
49 | 黄镇国, 张伟强, 陈特固. 香港/澳门近72年海平面变化曲线. 地理科学, 1999, 19 (3):91-94 Huang Zhengguo, Zhang Weiqiang, Chen Tegu. Sea level change curve in recent 72 years in Hongkong and Macau. Scientia Geographica Sinica, 1999, 19 (3):91~94 |
50 | Haigh I D, Eliot M, Pattiaratchi C. Global influences of the 18.61 year nodal cycle and 8.85 year cycle of lunar perigee on high tidal levels. Journal of Geophysical Research:Oceans, 2011, 116 (C6):1~12 |
51 | 郑文振. 我国海平面年速率的分布和长周期分潮的变化. 海洋通报, 1999, 18 (4):1-10 Zheng Wenzhen. Distribution of annual rates of sea level and variation of long period constituents in China. Marine Science Bulletin, 1999, 18 (4):1~10 |
52 | 陈特固, 杨清书, 徐锡祯. 广东沿海相对海平面变化特点. 热带海洋, 1997, 16 (1):95-100 Chen Tegu, Yang Qingshu, Xu Xizhen. The characteristics of sea level change along the coast of Guangdong Province. Tropic Oceanology, 1997, 16 (1):95~100 |
53 | 赵 青, 梁 娟. 基于遗传BP神经网络的非平稳时间序列预测. 四川理工学院学报(自然科学版), 2008, 21 (3):82-84 Zhao Qing, Liang Juan. Forecasting of nonstationary time series based on BP neural network of Genetic algorithms. Journal of Sichuan University of Science & Engineering(Natural Science Edition), 2008, 21 (3):82~84 |
54 | Hornik K, Stinchcombe M, White H. Multilayer feedforward networks are universal approximators. Neural Networks, 1989, 2 (5):359~366 |
55 | Asmar B N, Ergenzinger P. Long-term prediction of the water level and salinity in the Dead Sea. Hydrological Processes, 2002, 16 (14):2819~2831 |
56 | Cheng Y C, Andersen O B, Knudsen P. Integrating non-tidal sea level data from altimetry and tide gauges for coastal sea level prediction. Advances in Space Research, 2012, 50 (8):1099~1106 |
57 | 国家海洋局, 2010年中国海平面公报. 网址:http://www .soa .gov .cnzwgkhygb/zghpmgb/2010nzghpmgb/201212/t20121217_22808 .html 2011 State Ocean Administration, PRC. 2010 report on China's sea level. WWW Pages: http://www .soa .gov .cnzwgkhygb/zghpmgb/ 2010 nzghpmgb/201212/t20121217_22808.html2011 |
Abstract
The statistical model is one of the primary methods for sea level prediction. However, the conventional methods are still considered insufficient due to the complexity of the unstable and nonlinear sea level time series, and the prediction accuracy needs to be improved. Zhapo tide gauge (21°35'N,111°50'E)is a standard tide gauge with long records in the Pearl River delta, China. In this study the monthly records (01/1959~12/2011) of Zhapo, which was collected from Global Sea Level Observing System (GLOSS), was taken as an example, a hybrid method combining the data analysis method ensemble empirical mode decomposition (EEMD) and the artificial neural network (ANN) with back propagation algorithm (BP) was proposed to improve the prediction accuracy of sea level. EEMD extracted the frequency contributions called intrinsic mode functions (IMF) from the targeted time series oscillatory signal according to their unique fluctuation characteristic. 9 IMFs were obtained from the series of Zhapo tide gauge by using EEMD, and 6 of them passed the significance test on the level of 95%, suggesting that these 6 IMFs have statistical meaning. And the corresponding cycles such as seasonal cycle, semi-annual cycle, annual cycle, ENSO event cycle, semi nodical period and nutation cycle and so on are the significant characteristics in sea level change around Zhapo. The 6 IMFs were taken as the input factors of the BP neural network used for prediction. A BP neural network with three layers, which was optimized by setting the essential parameters such as the number of the hidden layers and the number of neurons in each layer, the learning rates, the momentum factor, and the number of training iterations and so on, was adopted in this study for modeling the sea level prediction. By reconstructing the outputs of the BP neural network, the prediction of the original time series was obtained. To evaluate the accuracy of the hybrid method, two predictions produced by the hybrid method and BP neural network directly were conducted in this paper. The results show that EEMD can extract the multiple time scales from the original time series, and the BP neural network is applicable for future sea level prediction. Comparing with the prediction by using BP neural network directly (R=0.76, RMSE=36.74mm, ME=-3.46), the EEMD-BP hybrid method performs better (R=0.89, RMSE=28.16mm, ME=2.31). The results suggest that EEMD-BP model significantly improves the prediction accuracy, and the preprocessing such as smoothing and denoising is strongly recommended for the prediction of unstable and nonlinear time series. The results of this study also reveals that the sea level has increased at a rate of 2.22mm/a over the period 1959~2011 in the Pearl River delta, and the prediction shows that it is expected to increase to 2.75mm/a by 2050.