长江河口段平均海面数值模拟研究

第３６卷第３期　水道港　口　Ｖ０１．３６　ＮＯ．３　２０１５年６月　Ｊｏｕｒｎａｌ　ｏｆ　Ｗａｔｅｒｗａｙ　ａｎｄ　Ｈａｒｂｏｒ　Ｊｕｎ．２０１５　长江河口段平均海面数值模拟研究　杨锋　，谭　亚。，王志伟　（１．９１６５０部队，广州５１０３２０；２．河海大学　港口海岸与近海工程学院，南京２１００９８）　摘要：为更好研究长江河口段平均海面的特性，建立了大通至长江口感潮河段的二维水动力数学模　型，通过模拟径流与潮汐共同作用下的水动力过程，并设置上下游不同条件的对照组，分析了河口段平　均海面的影响因子及响应规律。结果表明，长江河口段年平均海面Ａ。值由上游径流和外海潮波共同作　用产生，外海潮波给予了河道内平均水位一个沿程定值，而径流使得其产生了沿程衰减的趋势，这是河　口段Ａ。特征不同于外海的原因，也说明径流是其主要影响因子。Ａ。只受径流和外海潮位的年平均值影　响，而不受两者的年内变化影响。径流量的变化会导致上游Ａ０值比下游发生更加显著的响应；外海　值　的变化会使沿程发生大致等量的变幅。未来全球海平面上升后，河口段Ａｎ也将上升相近幅度。　关键词：数值模拟；平均海面；海平面上升；长江河口　中图分类号：Ｔｖ　１４２；Ｏ　２４２．１　文献标识码：Ａ　文章编号：１００５—８４４３（２０１５）０３—０２０４—０６　长江河口段西起江阴鹅鼻嘴，东至河口拦门沙，两岸地势较低，且人口稠密、经济发达，是黄金水道的人　海通道。因受到上游径流和外海潮波共同影响，长江口水动力条件复杂。平均海面是陆地高程和海图深度　基准面的起算面，长江口平均海面的变化规律如何，直接影响着当地的航道利用、盐水入侵、市区取水等经　济生活问题，并且对于当地提高潮汐预报精度、防灾减灾都具有十分重要意义。　沈焕庭…认为平均海面是人海河流河口区的侵蚀基准面，其升降会使河口纵比降发生变化，改变水流下　切能力，影响侵蚀与堆积。沈健、王宝灿　探讨了长江河口区的月平均海面时空变化规律，认为其月平均海　面变化主要受径流及其他海洋水文气象因子的作用，而且具有明显的季节性变化。陈宗镛等　通过计算长　江口附近４个验潮站过去２０　ａ的海面记录分析出海面变化总趋势是上升的，其变化受径流量大小和外海水　位共同影响，并且厄尔尼诺现象能使长江口年平均海面发生异常。陈西庆　通过分析表明，长江每年的人海　径流对该年的海面有重要影响，幅度为流量每变化１　０００　ｍ　／ｓ，年平均海面变化８．５　ｍｍ，但流量对长江口海面　的趋势性变化无显著影响。杨锋　通过对长江口实测潮位进行调和分析，发现受河床高度向上游抬升影响，　平均海面向上游呈递增趋势。　平均海面存在着日变化、月变化、年变化和多年变化，就长江河口而言，天文潮、径流、气压、气温和季风　等都是导致海面变化的基本因素，而且各因素的相互作用是相当复杂的。本文通过对江阴至共青圩站一整　年的逐时潮位数据进行调和分析　从而得到年平均海面Ａ。，并分析其在河口段的影响因子及响应程度。　收稿日期：２０１４—１０—１１；修回日期：２０１４—１卜ｏ３　基金项目：国家重点基础研究发展计划（９７３计划）课题：长江口海域水动力过程与生态系统演变机制（２０１０ＣＢ４２９００１）；国　家科技支撑计划（水沙变异条件下荆江与长江Ｉ＝Ｉ北支河道治理关键技术研究）课题：长江Ｉ＝ｔ北支建闸技术研究　（２０１３ＢＡＢ１２Ｂ０５）；江苏省水利科技项目：风暴潮耦合影响下长江镇扬河段水位预报技术研究（２０１２０３６）　作者简介：杨锋（１９８６一），男，河南省新乡人，助理工程师，主要从事验潮及水动力模拟分析工作。　Ｂｉｏｇｒａｐｈｙ：ＹＡＮＧ　Ｆｅｎｇ（１９８６一），ｍａｌｅ，ａｓｓｉｓｔａｎｔ　ｅｎｇｉｎｅｅｒ．　水道港　口　第３６卷第３期　中本文将平均海面　１００　徐六泾　口实测＿模拟　３００　徐六泾　口实测＿模拟　Ｏ　Ａ０作为一个角速率　８ｏ　为０的特殊分潮与　￣．ｒ　２５０　２００　其他分潮一起进行　．　１５０　了验证。　１００　由图２～图６可　５０　知，模拟值与实测值　Ｏ　Ｍ２　ｓ２　Ｎ２　Ｋ２　Ｋ１　Ｏ１　Ｐ１　Ｑ１　Ｍ４　ＭＳ４　Ｓａ　ＡＯ　吻合良好，逐时和Ｅｌ　分潮　均潮位偏差在１０　ｃｍ　图６分潮调和常数验证　以内，高低潮位时间　Ｆｉｇ．６　Ｖｅｒｉｉｆｃａｔｉｏｎ　ｏｆ　ｔｉｄａｌ　ｈａｒｍｏｎｉｃ　ｃｏｎｓｔａｎｔｓ　偏差小于３０　ａｒｉｎ；除Ｑ　等振幅量值较小的分潮外，其他分潮振幅偏差都在１０％Ｉ２内，迟角偏差小于１０。　。因　此，模型结果很好地复演了长江感潮河段的涨落潮过程。　２动力变化对Ａ。的影响．　大通站作为长江口的潮区界，一直被用作研究长江感潮河段的流量控制节点。连兴港位于长江口北支　末端，其水文特性为潮汐作用强于径流，因此本文模型视口门外的海洋潮汐特征与连兴港相同，并将其潮汐　特征赋予下游边界。本文在只考虑径流和潮流动力的情况下（图１），通过改变上游大通流量值和下游连兴　港潮位值来设置不同的上下游动力边界组　表１不同边界组合条件　合，每个组合计算时间长度为ｌ　ａ。Ｔａｂ．１　Ｃｏｍｂｉｎａｔｉｏｎ　ｃｏｎｄｉｔｉｏｎ　ｏｆｄｉｆｆｅｒｅｎｔｕｐｓｔｒｅａｍ　ａｎｄｄｏｗｎｓｔｒｅａｍｂｏｕｎｄａｒｉｅｓ　为了更直观地观察和分析年平均海面　垒　二　二　二二二二二二二二：：二　Ａ。在长江河口段的沿程变化趋势，现以各　，　大通２０００年流量过程　连＇￣２ｏｏｏ年逐时潮位　（日流量有年内变化）　（日均海面有年内变化）　潮位站距离上游大通站的距离，即河长长　＆，　大通２０００年流量过程　连兴港２０００年逐时潮位一２５　ｈ平均值　度作为横坐标，将江阴、天生港、徐六泾、杨　‘　（１３流量有年内变化）　（日均海面无年内变化）　林、吴淞和共青圩共６个测站的结果着重　，　恒定为０　／ｓ　连兴港２０００年逐时潮位　（日流量无年内变化）　（日均海面有年内变化）　标示出来，这六站距离大通河长分别约为：　＆，　，　大通２０００年流量过程连兴港２０００年逐时潮位一２５　ｈ平均值＋２９．３　ｃｌｎ　４１２　ｋｍ、４６２　ｋｍ、４９７　ｋｍ、５３４　ｋｍ、５６８　ｋｍ和　…　（日流量有年内变化）　（日均海面无年内变化）　ｆ　，．　恒定为２９　３５１　ｍ　／ｓ　连兴港２０００年逐时潮位　５９０ｋｍ。　～　（日流量无年内变化）　（日均海面有年内变化）　２．１动力变化对Ａ０的定性影响　为了确认河口段Ａ０是否是受上游径流下泄和外海潮波　向口门内上溯的共同影响，现以２０００年数据为基础来设置　上下游边界条件。　说明：逐时潮位减去２５　ｈ平均值可以过滤掉大部分长　周期低频潮波的影响，条件１．２中水位即在年平均海面２９－３　ｃｍ上下振动。２０００年大通平均流量为２９　３５１　ｍ。／ｓ；连兴港　站平均海面高度为２９．３　ｃｍ。　条件１．２和１．３是条件１的两种极端情况，目的是研究上　距大通距离（ｋｍ）　下游边界日均值的年内变化对平均海面是否有影响；条件　图７平均海面Ａ。沿程分布　１、２、３的对比则是为了探明河口段Ａ。是受径流和外海潮波　Ｆｉｇ．７　Ｄｉｓｔｒｉｂｕｔｉｏｎ　ａｌｏｎｇ　ｔｈｅ　ｉｒｖｅｒ　ｏｆ　Ａｏ　的单一作用还是共同影响。　表２主要站点Ａ０高度　由图７可知，条件１、１．２、１．３所计算出的Ａ０趋势线几乎　Ｔａｂ。２　Ｔｈｅ　ｖａｌｕｅｓ　ｏｆ　Ａｏ　ｉｎ　ｍａｊｏｒ　ｓｔａｔｉｏｎｓ　ｃｍ　重合在一起，说明长江河口段的沿程Ａ０值只取决于　上游径流量和外海潮位的平均值，而径流量值和潮　条件１　位值在年内是否出现长周期波动对其并无影响。　条件２　条件３　２０１５年６月　杨锋，等长江河口段平均海面数值模拟研究　２０７　对比条件１与２，当下游平均海面为０　ｃｍ时，径流单独产生的Ａ０值向下游依然保持着相同的衰减趋势和衰减　幅度，但各站数值都比原来减小了２５　ａｍ左右。当上游径流年均流量为０　ｍ　／ｓ时，条件３中的外海潮波向口　内上溯，但沿程Ａ０数值却基本保持不变，向上游呈现的趋势并未衰减，反而略微增大，这与河床高度向上游　抬升有关。从表２中的数据可知，各个站点条件１的Ａ。值都低于条件２与条件３之和１５　ｃｍ左右，且上游站点　中径流产生的水位值所占比重大，下游站点中外海潮波产生的水位值占总水位的比重更大。　这说明河口段平均海面Ａ。值由上游径流和外海潮波共同作用产生，径流对上游站点水位影响更大、外　海潮波对下游站点影响更大。在趋势上，外海的潮　表３不同边界组合条件　波给予了河道内平均水位一个沿程定值，而径流使Ｔａｂ．３　Ｃｏｍｂｉｎａｔｉｏｎ　ｃｏｎｄｉｔｉｏｎ　ｏｆｄｉｆｆｅｒｅｎｔｕｐｓｔｒｅａｍ　ａｎｄ　得其产生了沿程衰减的趋势，这是河口段Ａ。特征不　ｄｏｗｎｓｔｒｅａｍ　ｂｏｕｎｄａｒｉｅｓ　同于外海的原因，也说明径流是其主要的影响因子。　２．２动力变化对Ａ。的定量影响　为了研究长江河口段Ａ０高度对径流变化和外海　Ａ。变化响应的具体数值，现设置两组对照组与条件１　Ｒｎ　进行对比分析。　７０　说明：据１９９７～２００３年实测数据显示，大通年　ｕｕ　平均流量最小值为２００１年的２６　１５５　ｍ　／ｓ，最大值４０　为１９９８年的３９　３２１　ｍ　／ｓ。２０００年大通平均流量为　２９　３５１　ｍ　／ｓ放大１＿３倍后值为３８　１５６　ｍ　／ｓ，仍在往年　１０　范围内，符合实际（图８）。　据１９８９～２００９年实测数据显示，连兴港站平均海　面Ａ０最小值为１９９５年的１５．４８　ｃｍ，最大值为２００７年　图８大通不同条件的流量过程线　的３７．２３　ｃｍ。２００Ｏ年连兴港站Ａ０高度为２９．３　ｃｍ，放　Ｆｉｇ・　Ｐｒｏ　。。。　ｉｎｅ。　ｄｉ雎ｒｅｎ　ｄｉ　。ｈ　ｒｇ。　ａ　Ｄ咖ｎｇ　ｓ　砒ｉ。ｎ　大　口　多　小　游　站　数值均增大了相近的数量值，即在外海Ａｎ值增大了７．３　ｃｍ后，沿　距大通距离（ｋｍ）　程各站点Ａｎ值均增大了９　ｃｍ左右，且两者数值相近，这说明河口　图９平均海面　沿程分布　段站点Ａ０值对外海平均海面变化响应的敏感度几乎相同。　Ｆｉｇ．９　Ｄｉｓｔｒｉｂｕｔｉｏｎ　ａｌｏｎｇ　ｔｈｅ　ｒｉｖｅｒ　ｏｆ　Ａ０　２＿３全球海平面上升对Ａ。的影响　表４主要站点Ａ０高度（ａｍ）及变化量　由于人类活动对地球环境的影响迅速扩大，全　Ｔａｂ．４　Ｔｈｅ　ｖａｌｕｅ　ａｎｄ　ｖａｒｉａｔｉｏｎ　ｏｆ　Ａ０　ｉｎ　ｍａｊｏｒ　ｓｔａｔｉｏｎｓ　球变暖以及海平面上升是一全球性系列重大的环　境问题之一。在过去的１００年里，全球的海平面平　均上升速率为１－，－２　ｍｍ／ａ　。随着全球变暖，海平　面上升将加剧，大多数研究者认为，下一世纪全球　海平面上升速率可能达到过去１００　ａ平均上升速率　的３￣６倍ｎ　。据权威的ＩＰＣＣ报告称，到２１００年，全球海平面上升预测值的最佳估计值为６６　ｃｍ，最高估计值　为１１０　ｃｍ　…。　长江三角洲及其邻近地区是我国现代地壳沉降运动速率最大区域之一，也是人类开采地下水等造成地　面沉降较为严重的区域，近些年来，全区地面平均沉降速率为２～５　ｍｍ／ａ，其中上海市的最大沉降漏斗区为７　２０８　水道港　口　第３６卷第３期　４５　ｏｏｏ　４０ｏｏｏ　３５　ｏｏｏ　３００００　２５　ｏｏ０　暑　２Ｏ０００　１５　ｏｏｏ　惶　１０　ｏｏｏ　５　ｏｏｏ　０　距大通距离（ｋｉｎ）　图１０大通年平均流量与连兴港Ａ０多年变化　图１　１全球海面上升后的河口段Ａ。分布　Ｆｉｇ．１０　Ｔｅｍｐｏｒａｌ　ｖａｒｉａｔｉｏｎ　ｏｆ　ａｎｎｕａｌ　ｍｅａｎ　ｄｉｓｃｈａｒｇｅ　ａｔ　Ｄａｔｏｎｇ　ａｎｄ　Ａ０　ａｔ　Ｆｉｇ．１　１　Ｄｉｓｔｉｒｂｕｔｉｏｎ　ａｌｏｎｇ　ｔｈｅ　ｉｒｖｅｒ　ｏｆ　Ａｏ　ａｆｔｅｒ　ｔｈｅ　Ｌｉａｎｘｉｎｇｇａｎｇ　ｇｌｏｂａｌ　ｓｅａ　ｌｅｖｅｌ　ｒｉｓｅｓ　ｍｍ／ａ＿１　。因而，在２１世纪全球海平面的绝对上升　表５主要站点Ａ０高度（ｃｍ）及变化量　加上本地区的沉降所得到的海平面相对上升幅度　Ｔａｂ．５　Ｔｈｅ　ｖａｌｕｅ　ａｎｄ　ｖａｒｉａｔｉｏｎ　ｏｆ　Ａｏ　ｉｎ　ｍａｊｏｒ　ｓｔａｔｉｏｎｓ　将远超过全球平均值。据此，本文以２０００年数据为　条件　江阴　天生港徐六泾　杨林　吴淞　共青圩　基础，利用模型分别模拟长江口海域海平面上升　２０００年　１３５．９　１０２．７　８１．３　６８．５　５５．６　５０．２　５０　ａｍ和１００　ａｍ后的河口段Ａ０变化。　Ａ。抬升５０　ａｍ　１８３．５　１５２．０　１３１．１　１１８．４　１０４．９　９９．２　６　１７６．０　１６３．６　１４９．８　１４４＿３　１９９３　２００９年共１７　ａ的实测资料在图１０中显　∞　∞　∞　∞　加　Ａ０抬升１ｏｏ　ｃｍ　２２７．２　１９６．示，除了１９９８年等极致年份外，大通年平均流量与连兴港年平均海面Ａ０的变化趋势并不一致，且两者总的趋　势线倾斜方向刚好相反，即大通年平均流量呈逐年下降趋势，而连兴港Ａ。以９　ｍｍ／ａ的速率逐年上升。这与　全球海平面上升的大趋势相符，同时也说明连兴港站的Ａ。受上游径流影响极小、可以代表外海潮汐特征。　从图１１和表５可以看出，在外海海平面上升５０　ｃｍ后，长江河口段各站点Ａ０均增大了４９　ｃｍ左右；在外　海海平面上升１００　ａｍ后，长江河口段各站点Ａ。均增大了９３　Ｃｌｌｌ左右。这说明在外海海平面上升后，河口段　Ａ０均上升相近的高度，这就给该区域的人们提前敲响了警钟，应减小城市地面下沉并及早研究应对措施，以　防止重大海洋灾害的发生。　３结论　通过设置上下游不同条件的对照组，长江河口段平均海面Ａ０值由上游径流和外海潮波共同作用产生，　径流对上游站点水位影响更大、外海潮波对下游站点影响更大。外海的潮波给予了河道内平均水位一个沿　程定值，而径流使得其产生了沿程衰减的趋势，这是河口段Ａ０特征不同于外海的原因，也说明径流是其主要　的影响因子。径流量的变化会导致沿程Ａ。值发生响应变化，且上游比下游站点响应更敏感；外海Ａ０值的变　幅会使沿程振幅发生大致等量的变幅。未来全球海平面上升后，河口段Ａ０也上升相近幅度，这将提醒人们　及早采取措施预防海洋灾害。　参考文献：　［１］沈焕庭．平均海面研究的意义和计算方法［Ｊ］．海洋科学，１９８２（３）：３９—４３．　ＳＨＥＮ　Ｈ　Ｔ．Ｔｈｅ　ｃａｌｃｕｌａｔｉｏｎ　ｍｅｔｈｏｄ　ａｎｄ　ｒｅｓｅａｒｃｈ　ｓｉｇｎｉｆｉｃａｎｃｅ　ｏｆｍｅａｎ　ｓｅａ．１ｅｖｅｌ［ＪＪ．Ｍａｒｉｎｅ　Ｓｃｉｅｎｃｅｓ，１９８２（３）：３９－４３．　［２］沈健，王宝灿．长江河口区平均海面季节性变化的分析［Ｊ］．地理学报，１９９０（４）：４４卜４５０．　ＳＨＥＮ　Ｊ，ＷＡＮＧ　Ｂ　Ｃ．Ａｎａｌｙｓｉｓ　ｏｆ　ｓｅａｓｏｎａｌ　ｃｈａｎｇｅ　ｏｆｍｅａｎ　ｓｅａ—ｌｅｖｅｌ　ｉｎ　ｔｈｅ　Ｃｈａｎｇｊｉａｎｇ　ｅｓｔｕａｒｉｎｅ　ａｒｅａｌＪ　Ｊ．Ａｃｔａ　Ｇｅｏｇｒａｐｈｉｃａ　Ｓｉｎｉ—　ｃａ，１９９０（４）：４４１－４５０．　［３］陈宗镛，黄蕴和，周天华，等．长江口平均海面的初步研究［Ｊ］．海洋与湖沼，１９９１（４）：３１５—３２０．　ＣＨＥＮ　Ｚ　Ｙ，ＨＵＡＮＧ　Ｙ　Ｈ，ＺＨＯＵ　Ｔ　Ｈ，ｅｔ　ａ１．Ａ　ｐｒｅｌｉｍｉｎａｒｙ　ｓｔｕｄｙ　ｏｎ　ｍｅａｎ　ｓｅａ　ｌｅｖｅｌ　ｏｆ　ｔｈｅ　Ｃｈａｎｇｊｉａｎｇ　ｉｒｖｅｒ　ｅｓｔｕａｒｙ　ｌＪ　Ｊ．Ｏｃｅａｎｏｌｏ—　ｇｉａ　Ｅｔ　Ｌｉｍｎｏｌｏｇｉａ　Ｓｉｎｉｃａ，１９９１（４）：３１５－３２０．　［４］陈西庆．近七十年长江口海面变化研究及其意义［Ｊ］．上海水利，１９９２（２）：２１－３０．　ＣＨＥＮ　Ｘ　Ｑ．Ｔｈｅ　ｓｉｇｎｉｆｉｃａｎｃｅ　ａｎｄ　ｓｔｕｄｙ　ｏｎ　ｔｈｅ　ｃｈａｎｇｅ　ｏｆ　ｓｅａ　ｌｅｖｅｌ　ａｔ　ｔｈｅ　Ｃｈａｎｇｊｉａｎｇ　ｉｒｖｅｒ　ｅｓｔｕａｒｙ　ｉｎ　ｒｅｃｅｎｔ　７０　ｙｅａｒｓｌＪ　Ｊ．Ｓｈａｎｇｈａｉ　２０１５年６月　杨锋，等长江河口段平均海面数值模拟研究　２０９　Ｗａｔｅｒ，１９９２（２）：２１－３０．　［５］杨锋．长江河口段潮汐特性研究［Ｄ］．南京：河海大学，２０１４．　［６］赵有皓，王祥玲，张君伦．天文潮分析及预报实用系统［Ｊ］．河海大学学报：自然科学版，１９９９（４）：７３—７７．　ＺＨＡＯ　Ｙ　Ｈ，ＷＡＮＧ　Ｘ　Ｌ，ＺＨＡＮＧ　Ｊ　Ｌ．Ａ　ｐｒａｃｔｉｃａｌ　ｓｙｓｔｅｍ　ｏｆａｓｔｒｏｎｏｍｉｃａｌ　ｔｉｄａｌ　ａｎａｌｙｓｉｓ　ａｎｄ　ｐｒｅｄｉｃｔｉｏｎ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｈｏｈａｉ　Ｕｎｉ－　ｖｅｒｓｉｔｙ，１９９９（４）：７３—７７．　［７］Ｂｒｙｅ　Ｂ，Ｓｃｈｅｌｌｅｎ　Ｓ，Ｓａｓｓｉ　Ｍ，ｅｔ　ａ１．Ｐｒｅｌｉｍｉｎａｒｙ　ｒｅｓｕｌｔｓ　ｏｆ　ａ　ｆｉｎｉｔｅ．ｅｌｅｍｅｎｔ，ｍｕｌｔｉ．ｓｃａｌｅ　ｍｏｄｅｌ　ｏｆ　ｔｈｅ　Ｍａｈａｋａｍ　Ｄｅｌｔａ（Ｉｎｄｏｎｅｓｉａ）　［Ｊ］．Ｏｃｅａｎ　Ｄｙｎａｍｉｃｓ，２０１１，６１（８）：１　１０７－１　１２０．　［８］郑文振，吴乃华，金承钟，等．世界和中国的海平面变化［Ｊ］．海洋通报，１９９３（４）：９５—９９．　ＺＨＥＮＧ　Ｗ　Ｚ，ＷＵ　Ｎ　Ｈ，ＪＩＮＧ　Ｃ　Ｚ，ｅｔ　ａ１．Ｔｈｅ　ｓｅａ　ｌｅｖｅｌ　ｖａｒｉａｔｉｏｎ　ｉｎ　ｇｌｏｂａｌ　ａｎｄ　Ｃｈｉｎａ　ｓｅａ［Ｊ］．Ｍａｒｉｎｅ　Ｓｃｉｅｎｃｅ　Ｂｕｌｌｅｔｉｎ，１９９３（４）：９５—　９９．　［９］黄立人，马青．近几十年来的全球海面变化［Ｊ］．海洋学报：中文版，１９９３（６）：７６－８２．　ＨＵＡＮＧ　Ｌ　Ｒ，ＭＡ　Ｑ．Ｔｈｅ　ｓｅａ　ｌｅｖｅｌ　ｖａｒｉａｔｉｏｎ　ｉｎ　ｇｌｏｂａｌ　ｉｎ　ｒｅｃｅｎｔ　ｄｅＩｚａｄｅｓ［Ｊ］．Ａｃｔａ　Ｏｃｅａｎｏｌｏｇｉｃａ　Ｓｉｎｉｃａ，１９９３（６）：７６—８２．　［１０］Ｗａｒｒｉｃｋ　Ｒ　Ａ，Ｏｅｒｌｅｍａｎｓ　Ｈ．Ｓｅａ　ｌｅｖｅｌ　ｒｉｓｅ．Ｃｌｉｍａｔｅ　Ｃｈａｎｇｅ—Ｔｈｅ　ＩＰＣＣ　Ｓｃｉｅｎｉｉｆｃ　Ａｓｓｅｓｓｍｅｎｔ［ｚ］．Ｃａｍｂｒｉｄｇｅ　Ｕｎｉｖｅｒｓｉｔｙ　Ｐｒｅｓｓ，　１９９０：２５７－２８２．　［１　１］ＩＰＣＣ　Ｒｅｓｐｏｎｓｅ　Ｓｔｒａｔｅｇｉｅｓ　Ｗｏｒｋｉｎｇ　Ｇｒｏｕｐ．Ｓｔｒａｔｅｇｉｅｓ　ｆｏｒ　ａｄａｐｔａｔｉｏｎ　ｔｏ　ｓｅａ　ｌｅｖｅｌ　ｒｉｓｅ［Ｒ］．Ｃｏｍｂｒｉｄｇｅ：Ｃａｍｂｒｉｄｇｅ　Ｕｎｉｖｅｒｓｉｔｙ　Ｐｒｅｓｓ，　１９９０：１４３．　［１２］胡惠民，黄立人，杨国华．长江三角洲及其邻近地区的现代地壳垂直运动［Ｊ］．地理学报，１９９２（１）：２２—３０．　ＨＵ　Ｈ　Ｍ，ＨＵＡＮＧ　Ｌ　Ｒ，ＹＡＮＧ　Ｇ　Ｈ．Ｒｅｃｅｎｔ　ｃｒｕｓｔａｌ　ｖｅｒｔｉｃａｌ　ｍｏｗ￣ｍｅｎｔ　ｉｎ　ｔｈｅ　Ｃｈａｎｇｊｉａｎｇ　ｒｉｖｅｒ　ｄｅｌｔａ　ａｎｄ　ｉｔｓ　ａｄｊａｃｅｎｔ　ａｒｅａ［Ｊ］．Ａｃ—　ｔａ　Ｇｅｏｇｒａｐｈｉｃａ　Ｓｉｎｉｃａ，１９９２（１）：２２－３０．　Ｎｕｍｅｒｉｃａｌ　ｓｉｍｕｌａｔｉｏｎ　ｓｔｕｄｙ　ｏｎ　ｍｅａｎ　ｓｅａ　ｌｅｖｅｌ　ａｔ　Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　ｅｓｔｕａｒｙ　ｒｅａｃｈ　ＹＡＮＧ　Ｆｅｎｇ　，ＴＡＮ　Ｙａ　，ＷＡＮＧ　Ｚｈｉ－ｗｅｉ　（１．Ｐ　Ｎａｖｙ　９１６５０　Ｔｒｏｏｐｓ，Ｇｕａｎｇｚｈｏｕ　５１０３２０，Ｃｈｉｎａ；２．Ｃｏｌｌｅｇｅ　ｏｆＨａｒｂｏｒ，Ｃｏａｓｔａｌ　ａｎｄ　Ｏｆｆｓｈｏｒｅ　Ｅｎｇｉｎｅｅｒｉｎｇ，　Ｈｏｈａｉ　Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ　２１００９８，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｉｎ　ｏｒｄｅｒ　ｔｏ　ｓｔｕｄｙ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏ｛ｍｅａｎ　ｓｅａ　ｌｅｖｅｌ　ａｔ　Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　ｅｓｔｕａｒｙ　ｒｅａｃｈ．ａ　ｔｗｏ－ｄｉｍｅｎ－　ｓｉｏｎａｌ　ｈｙｄｒｏｄｙｎａｍｉｃ　ｍｏｄｅｌ　ｆｏｒ　ｔｈｅ　ｔｉｄａｌ　ｒｅａｃｈ　ｏｆ　ｔｈｅ　Ｄａｔｏｎｇ—Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　ｅｓｔｕａｒｙ　ｗａｓ　ｅｓｔａｂｌｉｓｈｅｄ．Ｂｙ　ｓｉｍｕｌａｔｉｎｇ　ｔｈｅ　ｈｙｄｒｏｄｙｎａｍｉｃｓ　ｉｎ　ｒｅｓｐｏｎｓｅ　ｔｏ　ｒｕｎｏｆｆ　ａｎｄ　ｔｉｄｅｓ　ａｎｄ　ｓｅｔｔｉｎｇ　ｄｉｆｆｅｒｅｎｔ　ｕｐｓｔｒｅａｍ　ａｎｄ　ｄｏｗｎｓｔｒｅａｍ　ｂｏｕｎｄａｒｉｅｓ，ｔｈｅ　ｉｎ—　ｌｆｕｅｎｃｅ　ｆａｃｔｏｒｓ　ａｎｄ　ｒｅｓｐｏｎｓｅｓ　ｏｆ　Ａ０　ａｔ　Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　ｅｓｔｕａｒｙ　ｒｅａｃｈ　ｗｅｒｅ　ａｎａｌｙｚｅｄ．Ｔｈｅ　ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ，ｔｈｅ　Ａｏ　ａｔ　Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　ｅｓｔｕａｒｙ　ｒｅａｃｈ　ｉｓ　ｃｏ—ｐｒｏｄｕｃｅｄ　ｂｙ　ｒｕｎｏｆｆ　ａｎｄ　ｔｉｄａｌ　ａｃｔｉｏｎ．Ｔｈｅ　ｏｆｆｓｈｏｒｅ　ｌｅｖｅｌ　ｇｉｖｅｓ　ａ　ｃｏｎｓｔａｎｔ　ｖａｌｕｅ　ｔｏ　ｔｈｅ　Ａ０　ａｌｏｎｇ　ｔｈｅ　ｅｓｔｕａｒｙ　ｒｅａｃｈ；ｍｅａｎｗｈｉｌｅ，ｔｈｅ　ｄｉｓｃｈａｒｇｅ　ｏｆ　ｕｒｎｏｆｆ　ｐｒｏｄｕｃｅｓ　ｔｈｅ　ｄｅｃａｙ　ｔｒｅｎｄ　ａｌｏｎｇ　ｔｈｅ　ｗａｙ　ｔｏ　ｉｔ，ｉｔ　ｓ　ｔｈｅ　ｒｅａｓｏｎ　ｗｈｙ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｔｈｅ　Ａ０　ａｌｏｎｇ　ｔｈｅ　ｅｓｔｕａｒｙ　ｒｅａｃｈ　ｉｓ　ｄｉｆｆｅｒｅｎｔ　ｆｒｏｍ　ｉｎ　ｏｆｆｓｈｏｒｅ　ａｒｅａ，ａｎｄ　ｉｔ　ｐｒｏｖｅｓ　ｔｈａｔ　ｔｈｅ　ｒｕｎｏｆｆ　ｉｓ　ｔｈｅ　ｍａｉｎ　ｉｎｆｌｕｅｎｃｅ　ｆａｃｔｏｒ．Ｔｈｅ　Ａｏ　ｉｎ　ｅｓｔｕａｒｙ　ｉｓ　ｏｎｌｙ　ａｆｆｅｃｔｅｄ　ｂｙ　ｔｈｅ　ａｖｅｒａｇｅ　ｌｅｖｅｌ　ｉｎ　ｏｆｆｓｈｏｒｅ　ａｎｄ　ｔｈｅ　ａｖｅｒａｇｅ　ｄｉｓｃｈａｒｇｅ　ｏｆ　ｕｒｎｏｆｆ，ｂｕｔ　ｎｏｔ　ｂｙ　ｔｈｅ　ｃｈａｎｇｅ　ｏｆ　ｒｕｎｏｆｆ　ａｎｄ　ｏｆｆｓｈｏｒｅ　ｔｉｄｅ　ｄｕｒｉｎｇ　ａ　ｙｅａｒ．Ｔｈｅ　ａｖｅｒａｇｅ　ｄｉｓｃｈａｒｇｅ　ｏｆ　ｒｕｎｏｆｆ　ｄｕｒｉｎｇ　ａ　ｙｅａｒ　ｗｉｌｌ　ａｆｆｅｃｔ　ｔｈｅ　Ａ０，ａｎｄ　ｔｈｅ　ｕｐｓｔｒｅａｍ　ｓｉｔｅｓ　ａｒｅ　ｍｏｒｅ　ｓｅｎｓｉｔｉｖｅ　ｔｈａｎ　ｄｏｗｎｓｔｒｅａｍ；ｔｈｅ　ｏｆｆｓｈｏｒｅ　Ａ０　ｗｉｌｌ　ｐｒｏｄｕｃｅ　ａ　ｓｉｍｉｌａｒ　ａｍｏｕｎｔ　ｏｆ　ｃｈａｎｇｅ　ｔｏ　ｉｔ．Ｗｉｔｈ　ｔｈｅ　ｇｌｏｂａｌ　ｓｅａ　ｌｅｖｅｌ　ｒｉｓｅ，ｔｈｅ　Ａ０　ｉｎ　ｅｓｔｕａｒｙ　ｗｉｌｌ　ｒｉｓｅ　ｂｙ　ａ　ｓｉｍｉｌａｒ　ｔｏ　ｏｆｆｓｂｏｒｅ．　Ｋｅｙ　ｗｏｒｄｓ：ｎｕｍｅｒｉｃａｌ　ｓｉｍｕｌａｔｉｏｎ；ｍｅａｎ　ｓｅａ　ｌｅｖｅｌ　Ａ０；ｓｅａ　ｌｅｖｅｌ　ｒｉｓｅ；Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　ｅｓｔｕａｒｙ　ｒｅａｃｈ