羅紫丹，孟憲偉，羅新正

（1.長(zhǎng)春師范大學(xué)，吉林長(zhǎng)春，130032；2.國(guó)家海洋局第一海洋研究所，山東青島，266061；3.煙臺(tái)大學(xué)，山東煙臺(tái)，264005）

百年內(nèi)全球海平面上升、地殼上升和潮灘沉積對(duì)廣西英羅灣紅樹林分布的影響

羅紫丹1，孟憲偉2，羅新正3

（1.長(zhǎng)春師范大學(xué)，吉林長(zhǎng)春，130032；2.國(guó)家海洋局第一海洋研究所，山東青島，266061；3.煙臺(tái)大學(xué)，山東煙臺(tái)，264005）

為有效保護(hù)和合理利用廣西紅樹林，探討了自2000年到2100年全球海平面上升、地殼上升和潮灘沉積對(duì)英羅灣紅樹林分布的影響。采用機(jī)理分析方法建立紅樹林邊界位置預(yù)測(cè)模型；利用預(yù)測(cè)模型計(jì)算在2000年數(shù)字高程模型上定位2100年紅樹林邊界的數(shù)據(jù)；利用Global Mapper軟件在2000年數(shù)字高程模型上確定紅樹林邊界的位置。利用Mapinfo professional地理信息系統(tǒng)軟件分析海平面上升引起的紅樹林邊界的移動(dòng)。通過研究得到在2000年數(shù)字高程模型上表達(dá)2100年紅樹林邊界位置的數(shù)據(jù)，以及自2000年到2100年紅樹林邊界移動(dòng)專題地圖。結(jié)果表明：自2000年到2100年，低、中模式海平面上升，英羅灣紅樹林邊界向外擴(kuò)展，面積增大；高模式海平面上升，紅樹林邊界向陸移動(dòng)，面積增大；極端模式海平面上升，紅樹林邊界向陸移動(dòng)，面積減小。

全球海平面上升；地殼上升；潮灘沉積；數(shù)字高程模型；紅樹林分布；英羅灣

Keywords：global sea level rise；crustal rise；tidal flat sedimentation；digital elevation model；mangrove distribution；Yingluo Bay

紅樹林具有重要的生態(tài)功能和經(jīng)濟(jì)價(jià)值。紅樹林能夠消能護(hù)岸（楊建民，2008），防御風(fēng)暴潮、海浪、海嘯的災(zāi)難性破壞；能夠促淤（譚曉林等，1997；Brodie et al，1996）、固碳（Ceron-Breton et al，2011），凈化水體和大氣；能夠?yàn)楹Ｑ笊锾峁⒌兀荒軌驗(yàn)槿祟愄峁┠静?、薪材、藥材和海產(chǎn)品（Faraco et al，2010）。

紅樹林分布在平均潮位（Mean Sea Level，縮寫為MSL）線與平均大潮高潮位（Mean High Water Springs，縮寫為MHWS）線之間的潮間帶（Matthew et al，2009），能夠隨平均潮位線和平均大潮高潮位線的移動(dòng)而移動(dòng)（Eric et al，2006）。因此，紅樹林向海邊界（Seaward Mangrove Margin，縮寫為SMM）與平均潮位線基本一致，紅樹林向陸邊界（Landward Mangrove Margin，縮寫為L(zhǎng)MM）與平均大潮高潮位線基本一致。

全球海平面上升（Global Sea Level Rise，縮寫為GSLR）抬升平均潮位和平均大潮高潮位，從而驅(qū)動(dòng)紅樹林邊界向陸地方向移動(dòng)。地殼上升（Crustal Rise，縮寫為CR）和潮灘沉積（Tidal Flat Sedimentation，縮寫為TFS）不改變平均潮位和平均大潮高潮位的高程，而是通過增高潮灘表面的高程，減小水深，從而驅(qū)動(dòng)紅樹林邊界向海洋方向移動(dòng)。地形和海堤既不改變潮位高程，也不改變水深，僅僅調(diào)控紅樹林邊界移動(dòng)的速率。

以英羅灣紅樹林為研究對(duì)象，預(yù)測(cè)自2000年到2100年全球海平面上升、地殼上升和潮灘沉積對(duì)紅樹林邊界移動(dòng)的影響。首先，建立紅樹林邊界位置預(yù)測(cè)模型；其次，利用預(yù)測(cè)模型計(jì)算在2000年數(shù)字高程模型上表達(dá)2100年紅樹林邊界位置的數(shù)據(jù)；第三，運(yùn)用Global Mapper和數(shù)字高程模型，確定紅樹林邊界的空間位置；最后，運(yùn)用Mapinfo professional地理信息系統(tǒng)軟件分析紅樹林邊界的移動(dòng)。

揭示百年內(nèi)全球海平面上升、地殼上升和潮灘沉積對(duì)紅樹林分布的影響，可為廣西紅樹林的有效保護(hù)和合理利用提供科學(xué)支持。

1 研究方法

1.1研究區(qū)概況

英羅灣位于南海西北部，廣西壯族自治區(qū)合浦縣山口鎮(zhèn)沿岸（圖1）。英羅灣的紅樹林主要有木欖（Bruguiera gymnorrhiza）、紅海欖（Rhizophora stylosa）、秋茄（Kandelia candel）、桐花樹（Aegiceras corniculatum）和白骨壤（Avicennia marina）。它們構(gòu)成了木欖群落、木欖+紅海欖群落、紅海欖群落、紅海欖+秋茄群落、秋茄群落、秋茄+桐花樹群落、秋茄+白骨壤群落、桐花樹群落，以及桐花樹+白骨壤群落（梁士楚，1996）。木欖種群面積為12.1 hm2；紅海欖種群面積為44.3hm2；秋茄為12.1hm2；桐花樹為8.1hm2；白骨壤為4.0 hm2（梁士楚，1999；范航清等，2005）。

圖1 研究區(qū)地理位置

紅樹林種群在英羅灣成帶狀分布。木欖主要分布在內(nèi)灘；紅海欖分布在內(nèi)灘和中灘；秋茄分布在內(nèi)灘、中灘和外灘；桐花樹分布在內(nèi)灘、中灘和外灘；白骨壤分布在外灘（溫遠(yuǎn)光等，2002）。內(nèi)灘位于平均高潮位線附近，物質(zhì)組成以淤泥質(zhì)為主，坡度為4.35%；外灘位于平均潮位線附近，物質(zhì)組成以泥沙質(zhì)為主，坡度為1.84%（梁士楚，1996）。

英羅灣潮位參照石頭埠驗(yàn)潮站的潮位數(shù)據(jù)。石頭埠驗(yàn)潮站是最接近英羅灣的驗(yàn)潮站，二者相距30 km。英羅灣位于博白-合浦-北海斷裂以東的云開大山斷?。üＯ椋?994）。新生代以來，英羅灣所處的地殼處于不斷地上升過程之中（尹克堅(jiān)，1991）。英羅灣紅樹林向陸邊界潮灘主要由風(fēng)化物質(zhì)組成。特別是向陸邊界不斷向陸推進(jìn)，邊界處幾乎無沉積積累。因此，潮灘沉積驅(qū)動(dòng)英羅灣紅樹林向海邊界向海方向移動(dòng)，幾乎不對(duì)向陸邊界的移動(dòng)產(chǎn)生影響。

解放塘海堤位于英羅灣紅樹林向陸邊界的中段，長(zhǎng)1 552 m，為圍海造田于上世紀(jì)50年代修建，因此而清除了76.47 hm2的紅樹林。

從1991年到2000年，因?yàn)閲鷫?，英羅灣紅樹林每年減少4.2 hm2；從2000年到2010年，因?yàn)樽匀粩U(kuò)展和人工造林，英羅灣紅樹林每年增加2.9 hm2（朱耀軍等，2013）。由于英羅灣紅樹林區(qū)已納入山口國(guó)家級(jí)紅樹林生態(tài)自然保護(hù)區(qū)，未來英羅灣紅樹林景觀的演變將主要受自然要素的影響。

1.2軟件與數(shù)據(jù)

1.2.1 軟件

研究中使用了3款軟件。Envi圖像處理軟件，美國(guó)ITT visual information solutions公司發(fā)行；Global mapper制圖軟件，美國(guó)Blue marble Geographics公司發(fā)行；Mapinfo professional地理信息系統(tǒng)軟件，美國(guó)Mapinfo Corporation發(fā)行。

1.2.2 數(shù)據(jù)

（1）遙感影像與數(shù)字高程模型

1景Landsat遙感影像數(shù)據(jù)，軌道號(hào)為p124r045，分辨率為30 m，采集于2000年10月30日。1景SRTM DEM數(shù)字高程模型數(shù)據(jù)，軌道號(hào)為n21e109，分辨率為90 m，采集于2000年2月，高程基準(zhǔn)為大地水準(zhǔn)面。上述數(shù)據(jù)均從http：//glcfapp.glcf.umd.edu：8080/esdi/下載。

（2）全球平均海平面上升

IPCC AR4指出，從1990-1999年段到2090-2099年段，全球海平面按照6種情景上升（IPCC， 2007）。選擇B1、A1B和A1FI情景對(duì)應(yīng)的全球平均海平面的上升值的中值作為全球平均海平面上升的低、中、高模式的值，將“A1FI+冰原動(dòng)力增量”作為全球平均海平面上升的極端模式值。這4種模式全球平均海平面分別上升28 cm、35 cm、43 cm和58 cm。

（3）潮位

根據(jù)從1969到1993年25年的實(shí)測(cè)潮位數(shù)據(jù)統(tǒng)計(jì)結(jié)果（莫永杰等，1996），以及1985高程基準(zhǔn)與大地水準(zhǔn)面的差值（焦文海等，2002），得到英羅灣的平均潮位為63 cm，平均大潮高潮位為248 cm。2000年平均潮位和平均大潮高潮位采用自1969年到1993年連續(xù)25年的統(tǒng)計(jì)結(jié)果。

（4）地殼上升速率

1954年至1989年，英羅灣所處地塊地殼上升速率為0.5~3.5 mm/a，平均上升速率為2.0 mm/a（莫永杰等，1996）。到2100年，地殼將上升20 cm。

（5）潮灘沉積速率

2010年5月，在英羅灣紅樹林向海邊界處采集兩個(gè)沉積巖芯，長(zhǎng)度分別為8.8 m和5.2 m，編號(hào)分別為YLW02和YLW04。在華東師范大學(xué)河口海岸國(guó)家重點(diǎn)實(shí)驗(yàn)室進(jìn)行210Pb測(cè)年，以137Cs為時(shí)標(biāo)，計(jì)算出YLW02和YLW04巖芯的平均沉積速率分別為0.21 cm/a（表1）和0.16 cm/a（表2），平均為0.185 cm/a。按此平均沉積速率計(jì)算，到2100年，英羅灣潮灘沉積厚度將增加18.5 cm。

1.3技術(shù)路線

1.3.1 構(gòu)建在2000年DEM上表達(dá)2100年紅樹林向海邊界（SMM）和向陸邊界（LMM）位置數(shù)據(jù)的預(yù)測(cè)模型

通過機(jī)理分析，建立預(yù)測(cè)模型（公式（1）、公式（2）），計(jì)算在2000年DEM上表達(dá)2100年紅樹林向海邊界、向陸邊界位置的數(shù)據(jù)。

公式（1）中，SMM2100是在2000年SRTM DEM上表達(dá)2100年紅樹林向海邊界位置的數(shù)據(jù)；MSL2000是2000年的平均潮位。SLR2100表示自2000年到2100年的海平面上升值。CR2100是自2000年到2100年地殼上升幅度。TFS2100是自2000年到2100年紅樹林向海邊界處沉積層厚度增量。

表1 210Pb測(cè)定的YLW02巖芯沉積速率剖面

表2 210Pb測(cè)定的YLW04巖芯沉積速率剖面

公式（2）中，LMM2100是在2000年SRTM DEM上表達(dá)2100年平均大潮高潮位線位置的數(shù)據(jù)；MHWS2000是2000年的平均大潮高潮位。SLR2100、CR2100的含義與公式（1）中的相同。

受初始條件限制，公式（1）、（2）僅適用于英羅灣及其附近海區(qū)。

1.3.2 計(jì)算在2000年DEM上表達(dá)2100年SMM和LMM位置的數(shù)據(jù)

利用公式（1）和相應(yīng)的MSL2000、SLR2100、CR2100、TFS2100值，可計(jì)算出的SMM2100值。利用公式（2），和相應(yīng)的MHWS2000、SLR2100、CR2100值，可計(jì)算出的LMM2100值。

1.3.3 生成紅樹林的向海邊界和向陸地邊界

應(yīng)用Global Mapper制圖軟件，輸入2000年2月的SRTM數(shù)字高程模型數(shù)據(jù)，并將投影校正為北半球50區(qū)UTM投影，以與假彩色遙感影像底圖的投影相匹配。輸入2000年平均潮位和表達(dá)2100年年紅樹林向海邊界位置的數(shù)據(jù)，生成紅樹林向海邊界。輸入2000年平均大潮高潮位和表達(dá)2100年紅樹林向陸邊界的數(shù)據(jù)，生成紅樹林向陸邊界。以矢量、Mapinfo的格式輸出、保存。

1.3.4 SMM和LMM移動(dòng)距離的計(jì)算模型

應(yīng)用Mapinfo v7.0地理信息分析軟件，輸入由ENVI軟件合成的假彩色遙感影像，提取海堤信息。輸入2000年紅樹林邊界和2100年紅樹林邊界數(shù)據(jù)。

參照2000年紅樹林的向海邊界，分析自2000年到2100年全球海平面上升、地殼上升和潮灘沉積對(duì)紅樹林向海邊界移動(dòng)的影響。提取2100年紅樹林向海邊界與2000年向海邊界之間區(qū)域的面積（As）和周長(zhǎng)（Ps），利用公式（3）計(jì)算紅樹林向海邊界的移動(dòng)距離（DS）。正值代表向海前進(jìn)，負(fù)值代表向陸后退，零值代表穩(wěn)定。

公式（3）中，sgn（）是符號(hào)函數(shù)。對(duì)于紅樹林向海邊界，當(dāng)CR2100+TFS2100-SLR2100＞0時(shí)，sgn（）取正，紅樹林向海邊界向海方向移動(dòng)，表現(xiàn)為紅樹林向海側(cè)的擴(kuò)張；當(dāng)CR2100+TFS2100-SLR2100=0時(shí)，sgn（）取0，紅樹林向海邊界保持穩(wěn)定；CR2100+TFS2100-SLR2100＜0時(shí)，sgn（）取負(fù)，紅樹林向海邊界向陸方向移動(dòng)，表現(xiàn)為紅樹林向海側(cè)的退縮。

將紅樹林向海邊界的進(jìn)退引起的紅樹林向海側(cè)的范圍變化設(shè)想為長(zhǎng)為L(zhǎng)s，寬為Ds的規(guī)則矩形，則有：

解上述方程，即可得到公式（3）中紅樹林邊界變動(dòng)范圍的平均寬度：

同理，參照2000年紅樹林的向陸邊界，分析自2000年到2100年海平面上升和地殼上升對(duì)紅樹林向陸邊界移動(dòng)的影響。提取2100年紅樹林向陸邊界與2000年向陸邊界之間區(qū)域的面積（Al）和周長(zhǎng)（Pl），利用公式（4）計(jì)算紅樹林向陸邊界的平均移動(dòng)距離（Dl）。正值代表向陸前進(jìn)，負(fù)值代表向海后退，零值代表穩(wěn)定。

2 研究結(jié)果

2.12100年平均潮位線（SMM2100）和平均大潮高潮位線（LMM2100）位置

根據(jù)4種海平面上升模式，根據(jù)公式（1）、公式（2）計(jì)算在2000年SRTM DEM上表達(dá)2100年平均潮位線位置、平均大潮高潮位線位置的數(shù)據(jù)（表3）。

表3 2100年平均潮位線（SMM2100）和平均大潮高潮位線（LMM2100）位置

2.2 SMM和LMM移動(dòng)專題圖

將2000年SRTM DEM數(shù)據(jù)，MSL2000、MHWS2000、SMM2100和LMM2100數(shù)據(jù)輸入Global Mapper制圖軟件，制作基于2000年SRTM DEM的紅樹林向海邊界SMM和向陸邊界LMM。將SMM和LMM輸入Mapinfo professional地理信息系統(tǒng)軟件，制作SMM和LMM移動(dòng)專題地圖（圖2）。

平均海平面上升28cm、35cm、43cm或58cm，紅樹林向陸邊界中段受海堤阻擋而保持穩(wěn)定，東段（EL）和西段（WL）均向陸方向移動(dòng)，移動(dòng)距離隨海平面上升而增大。平均海平面上升28 cm或者35 cm，紅樹林向海邊界向海方向移動(dòng)，移動(dòng)距離隨海平面上升而減小。平均海平面上升43 cm或58 cm，紅樹林向海邊界均向陸地方向移動(dòng)。

3 討論

利用Mapinfo professional地理信息系統(tǒng)提取2000年SMM與2100年SMM之間區(qū)域的面積和周長(zhǎng)，應(yīng)用公式（3）計(jì)算SMM移動(dòng)的平均距離。提取2000年LMM與2100年LMM之間區(qū)域的面積和周長(zhǎng)，應(yīng)用公式（4）計(jì)算LMM移動(dòng)的平均距離。

圖2 不同海平面上升模式下英羅灣紅樹林空間分布的變化

3.1自2000年到2100年，英羅灣紅樹林分布的變化

當(dāng)全球海平面以低模式上升時(shí)，即到2100年上升28 cm，英羅灣紅樹林向陸邊界的東段與西段分別向陸地方向移動(dòng)1.6 m和1.4 m，紅樹林分別增加0.36 hm2和0.11 hm2；向海邊界向海方向移動(dòng)5.7 m，導(dǎo)致紅樹林增加1.94 hm2。在全球海平面以低模式上升時(shí)，英羅灣紅樹林將凈增加2.41 hm2。

當(dāng)全球海平面以中模式上升時(shí)，即到2100年上升35 cm，英羅灣紅樹林向陸邊界的東段與西段分別向陸地方向移動(dòng)2.7 m和2.5 m，使紅樹林分別增加0.62 hm2和0.21 hm2；向海邊界向海方向移動(dòng)1.8 m，使紅樹林增加0.63 hm2。在全球海平面以中模式上升時(shí)，英羅灣紅樹林將凈增加1.46 hm2。

當(dāng)全球海平面以高模式上升時(shí)，即到2100年上升43cm，英羅灣紅樹林向陸邊界的東段與西段分別向陸地方向移動(dòng)4.1 m和3.8 m，使紅樹林分別增加0.93 hm2和0.31 hm2；向海邊界向陸方向移動(dòng)2.3 m，使紅樹林減少0.79 hm2。在全球海平面以高模式上升時(shí)，英羅灣紅樹林將凈增加0.45 hm2。

當(dāng)全球海平面以極端模式上升時(shí)，即到2100年上升58 cm，英羅灣紅樹林向陸邊界的東段與西段分別向陸地方向移動(dòng)6.6 m和7.8 m，使紅樹林分別增加1.5 hm2和0.64 hm2；向海邊界向陸方向移動(dòng)9.7 m，使紅樹林減少3.29 hm2。在全球海平面以極端模式上升時(shí)，英羅灣紅樹林將凈減少1.15 hm2。

當(dāng)全球海平面上升28 cm或35 cm時(shí)，紅樹林向海邊界向海移動(dòng)，移動(dòng)距離相應(yīng)減小；當(dāng)海平面上升43 cm或58 cm時(shí)，紅樹林向海邊界向陸移動(dòng)，移動(dòng)距離相應(yīng)增大。自2000年到2100年，當(dāng)海平面上升28cm、35cm、43cm或58cm時(shí)，紅樹林向陸邊界向陸移動(dòng)，移動(dòng)距離相應(yīng)增大（圖3）。

圖3 自2000年到2100年紅樹林邊界移動(dòng)的平均距離

自2000年到2100年，當(dāng)平均海平面上升28cm或35cm時(shí)，紅樹林向海、向陸兩側(cè)面積均增加，總面積增大；當(dāng)海平面上升43 cm時(shí)，紅樹林向陸側(cè)面積增大，向海側(cè)面積減小，總面積仍增大；當(dāng)海平面上升58 cm時(shí)，紅樹林向陸側(cè)面積增加，向海側(cè)面積減小，總面積減?。▓D4）。

圖4 自2000年到2100年紅樹林面積變化

3.22000年英羅灣紅樹林分布與潮位線的關(guān)系

利用研究區(qū)2000年的遙感高程數(shù)據(jù)和潮位數(shù)據(jù)，在Global Mapper制圖軟件中生成平均潮位線和平均大潮高潮位線；利用研究區(qū)2000年遙感影像數(shù)據(jù)的4、5、3波段在ENVI遙感數(shù)據(jù)處理軟件中疊加，生成假彩色遙感影像；將潮位線與假彩色遙感影像在Mapinfo地理信息系統(tǒng)軟件中疊加，結(jié)果表明，平均潮位線與紅樹林向海邊界基本吻合，平均大潮高潮位線與紅樹林向陸邊界也基本吻合，表現(xiàn)出平均潮位線與紅樹林向海邊界基本一致、平均大潮高潮位線與紅樹林向陸邊界基本一致的基本趨勢(shì)。

3.3公式（1）、（2）左端項(xiàng)SMM和LMM的物理意義

由于地殼上升，2000年和2100年具有不相同的數(shù)字高程模型。在2100年的數(shù)字高程模型上，必須使用全球平均海平面的絕對(duì)變化值來分析2100年的潮位線的位置；在2000年的數(shù)字高程模型上，必須使用全球平均海平面的相對(duì)變化值（考慮地殼上升的影響）分析2100年的潮位線的位置。

因?yàn)橹挥?000年的數(shù)字高程模型，就必須使用全球平均海平面的相對(duì)變化在2000年數(shù)字高程模型上分析2100年的潮位線的位置。利用公式（1）、（2）計(jì)算的數(shù)值，是潮位的相對(duì)變化值，不是在2100年數(shù)字高程模型能夠使用的高程數(shù)據(jù)，而是只能夠在2000年數(shù)字高程模型上使用的、表達(dá)2100年潮位線位置的數(shù)據(jù)。

表4 紅樹林變化面積的對(duì)比

3.4地形、地貌對(duì)紅樹林面積變化值的影響

利用遙感數(shù)據(jù)計(jì)算的紅樹林面積的變化，僅是紅樹林地表面積變化在水平面上的投影面積。利用紅樹林面積變化區(qū)的地形坡度，根據(jù)紅樹林面積變化的投影，可以近似得到紅樹林面積變化的實(shí)際值（地表面積）（表4）。

表中數(shù)據(jù)表明，地形、地貌對(duì)紅樹林面積變化值的影響比較小。地表面積比投影面積大2‰到18‰。

3.5紅樹林面積變化值的誤差

從2000年到2100年，全球平均海平面上升的低、中、高和極端模式值和地殼上升速率，均為其分布范圍內(nèi)的中值。全球海平面上升低模式范圍為28±10 cm；中模式為34.5±13.5 cm；高模式為42.5±16.5 cm；極端模式為57.5±21.5 cm。研究區(qū)地殼上升范圍為20±15 cm。

根據(jù)誤差疊加規(guī)則、潮位相對(duì)變化的閾值及其對(duì)應(yīng)的面積，計(jì)算出紅樹林在低、中、高和極端模式下面積變化的分布范圍。全球海平面上升低模式時(shí)，紅樹林增加2.41±6.09 hm2；中模式時(shí)增加1.46±6.12 hm2；高模式時(shí)增加0.45±7.96 hm2；極端模式時(shí)減少1.15±8.40 hm2。

全球海平面上升每種模式都有一個(gè)較大的變化范圍，地殼上升也有一個(gè)較大的變化范圍，使得相應(yīng)的紅樹林的面積變化也有一個(gè)比較大的變動(dòng)范圍。因此，文中所預(yù)測(cè)的紅樹林變化情況，僅是一個(gè)平均趨勢(shì)。

4 結(jié)論

通過海平面上升、地殼上升和潮灘沉積自2000年到2100年對(duì)英羅灣紅樹林分布影響的研究，得到以下結(jié)論。

（1）平均海平面以低模式（28 cm）或中模式（35 cm）上升時(shí)，廣西英羅灣紅樹林向海、向陸邊界均向外擴(kuò)展，導(dǎo)致紅樹林面積分別增加2.41 hm2和1.46 hm2。

（2）平均海平面以中模式（43 cm）上升時(shí)，廣西英羅灣紅樹林向海、向陸邊界均向陸地方向移動(dòng)，導(dǎo)致紅樹林面積增加0.45 hm2。

（3）平均海平面以極端模式（58 cm）上升時(shí)，廣西英羅灣紅樹林向海、向陸邊界均向陸地方向移動(dòng)，導(dǎo)致紅樹林面積減少1.15 hm2。

Brodie M Y,Harvey L E,1996.A Spatial Analysis of the Relationship between Mangrove（Avicennia marina var australasica）Physiognomy and Sediment Accretion in the Hauraki Plains,New Zealand. Estuarine,Coastal and Shelf Science,(42):231-246.

Ceron-Breton R M,Ceron-Breton J G,Sanchez-Junco R C,et al,2011. Evaluation of carbon sequestration potential in mangrove forest at three estuarine sites in Campeche,Mexico.International Journal of Energy and Environment,5(4):487-494.

Eric L G,Joanna E,Vainuupo J,et al,2006.Adapting to Pacific Island mangrove responses to sea level rise and climate change.CLIMATE RESEARCH,32:161-176.

Faraco L F D,Jose M A,Paulo C L,2010.A methodology for assessing the vulnerability of mangroves and fisherfolk to climate change.Pan-American Journal of Aquatic Sciences,5（2）：205-223.

IPCC,2007.Climate Change 2007-The Physical Science Basis.Cambridge University Press,Cambridge.

Matthew H,Heather C,Pablo H,et al,2009.the impacts of sealevel rise on the california coast.California Climate Change Center,http：// pacinst.org/publication/the-impacts-of-sea-level-rise-on-thecalifornia-coast/

范航清，陳光華，何斌源，等，2005.山口紅樹林濱海濕地與管理.北京：海洋出版社.

郭福祥，1994.廣西大地構(gòu)造單元.桂林冶金地質(zhì)學(xué)院學(xué)報(bào)，14（3）：233-243.

焦文海，魏子卿，馬欣，等，2002.1985國(guó)家高程基準(zhǔn)相對(duì)于大地水準(zhǔn)面的垂直偏差.測(cè)繪學(xué)報(bào)，31（3）：196-201.

梁士楚，1996.廣西英羅灣紅樹植物群落的研究.植物生態(tài)學(xué)報(bào)，20（4）：310-321.

梁士楚，1999.廣西的紅樹林資源及其可持續(xù)利用.海洋通報(bào)，18（6）：77-83.

莫永杰，李平日，方國(guó)祥，等，1996.海平面上升對(duì)廣西沿海的影響與對(duì)策.北京：科學(xué)出版社.

譚曉林，張喬民，1997.紅樹林潮灘沉積速率及海平面上升對(duì)我國(guó)紅樹林的影響.海洋通報(bào)，16（4）：29-35.

溫遠(yuǎn)光，劉世榮，元昌安，2002.廣西英羅灣紅樹植物種群的分布.生態(tài)學(xué)報(bào)，22（7）：1105-1110.

楊建民，2008.海岸帶邊坡防浪林消浪理論與實(shí)驗(yàn)研究.海洋通報(bào)，27（2）：16-21.

尹克堅(jiān)，1991.廣西近期地殼垂直運(yùn)動(dòng)初探.大地測(cè)量與地球動(dòng)力學(xué)，4：86-90.

朱耀軍，郭菊蘭，武高潔，等，2013.近20年來英羅灣紅樹林景觀過程及周邊土地利用/覆蓋變化.北京林業(yè)大學(xué)學(xué)報(bào)，35（2）：22-29.

（本文編輯:袁澤軼）

Impact of global sea level rise，crustal rise and tidal flat sedimentation for one hundred years on the mangrove distribution in the Yingluo Bay of Guangxi

LUO Zi-dan1,MENG Xian-wei2,LUO Xin-zheng3
（1.Changchun Normal University，Changchun130032,China;2.The First Instituteof Oceanography,SOA, Qingdao 266061,China;3.Yantai University,Yantai 264005,China）

In this paper,the impact of global sea level rise,crustal rise and tidal flat sedimentation on the mangrove distribution is explored to protect the mangrove in the Yingluo Bay of Guangxi.The tidal level-predicting models are built through the mechanism analysis.The model is used to calculate data placing seaward and landward mangrove margins of 2100 on the digital elevation model of 2000.The seaward and landward mangrove margins are determined on the digital elevation model of 2000 with global mapper software.The data are got placing the seaward and landward mangrove margins of 2100 on the digital elevation model of 2000 and the maps of mangrove margin migration from 2000 to 2100.The results suggest three conclusions.First,when global sea level rises at low and moderate mode,the mangrove margins will migrate outward,and mangrove area will increase;Second,when the global sea level rises at high moderate mode,the mangrove margins will migrate landward，and the mangrove area will increase;Finally,when the global sea level rises at the extreme mode,the mangrove margins will migrate landward,and the mangrove area will decrease.

P735

A

1001-6932（2017）02-0209-08

10.11840/j.issn.1001-6392.2017.02.012

2015-09-08；

2015-12-16

全球變化研究國(guó)家重大科學(xué)研究計(jì)劃（2010CB951203）；海洋沉積與環(huán)境地質(zhì)國(guó)家海洋局重點(diǎn)實(shí)驗(yàn)室開放基金（MASEG201203）；山東省高等學(xué)校教學(xué)改革項(xiàng)目（2012192）；山東省研究生教育創(chuàng)新計(jì)劃（SDYY13058）。

羅紫丹（1990-），計(jì)算機(jī)科學(xué)與技術(shù)專業(yè)。電子郵箱：luozidanlan@sina.com。

羅新正，博士，教授。電子郵箱：xzhluo@163.com。