任廣波，劉艷芬，馬 毅，張 杰

(1．國家海洋局第一海洋研究所，山東 青島266061;2．東營市海洋與漁業(yè)局，山東 東營257091)

1964年，米草屬植物首次由南京大學(xué)從英國引種到我國的江蘇海岸帶［1］，用于護堤和改良土壤，之后的30年里，包括大米草和互花米草在內(nèi)的米草植物不斷在我國沿海省份推廣、引種，到目前已遍布我國除海南和臺灣的所有沿海省份［2］。2003年，國家環(huán)境保護總局正式將米草屬的互花米草納入我國第一批外來入侵物種名單［3］?，F(xiàn)代黃河三角洲自1990年前后開始在孤東采油區(qū)北側(cè)五號樁附近引種互花米草［4］，隨后便迅速生長蔓延，至今已遍布現(xiàn)代黃河三角洲自然保護區(qū)潮間帶區(qū)域。黃河三角洲自然保護區(qū)是以保護珍稀水禽丹頂鶴和黑嘴鷗等為主的國家級自然保護區(qū)，而上述珍禽以蘆葦、堿蓬和檉柳為主要棲息地，若互花米草繼續(xù)擴張最終將影響到該區(qū)域生態(tài)格局，可能使該區(qū)域陷入嚴(yán)重的生態(tài)危機［5，6］。有研究已經(jīng)表明，互花米草的生長已對濱海濕地保護區(qū)內(nèi)保護鳥類的種群數(shù)量造成了影響［7-11］。

目前，互花米草在現(xiàn)代黃河三角洲的瘋狂漲勢和蔓延已經(jīng)引起了國內(nèi)外很多學(xué)者以及當(dāng)?shù)卣块T的高度重視，有針對性地開展了一些生態(tài)學(xué)方面的監(jiān)測和防治對策研究［5，6，12-16］，于祥等［14］對2005至2007年該區(qū)域互花米草的分布與面積進行過分析，而對于黃河三角洲自然保護區(qū)所在的現(xiàn)代黃河三角洲自引種至今互花米草的分布和變遷情況，卻未見研究報道。

本文基于遙感和GIS技術(shù)，收集了覆蓋黃河三角洲區(qū)域1989年至2013年24年間的24個時相的Landsat5 TM、Landsat7 ETM+和Landsat8 OLI遙感影像，結(jié)合2013年的現(xiàn)場調(diào)查數(shù)據(jù)，應(yīng)用纓帽變換和SVM的互花米草識別技術(shù)，通過分類后處理，完成了24年間逐年的互花米草分布信息提取;根據(jù)互花米草的面積隨時間的變化特點以及其分布隨河口區(qū)海岸變化的特征，分別從時間上和空間上分為4個時段和4個區(qū)域，開展了對互花米草在現(xiàn)代黃河三角洲的發(fā)展和變遷分析。最后，結(jié)合黃河每年的水沙量測量數(shù)據(jù)和氣象數(shù)據(jù)，分析了其變化原因。

1 區(qū)域和數(shù)據(jù)

1．1 研究區(qū)域

現(xiàn)代黃河三角洲，在地理上指1934年黃河分流點由墾利縣寧海下移至漁洼時開始建造的、以漁洼為頂點、北起挑河口南至宋春榮溝的扇形區(qū)域，完全包括黃河三角洲國家級自然保護區(qū)的區(qū)域范圍，同時又是勝利油田的主產(chǎn)區(qū)。近年來，由于黃河來水來沙量的銳減［12，17］和地面沉降［18］，致使大面積的現(xiàn)代黃河三角洲區(qū)域發(fā)生不同程度的侵蝕。同時，勝利油田的開發(fā)生產(chǎn)和東營市的城市建設(shè)都要求要有一片穩(wěn)固的海岸，于是相關(guān)部門于1990年前后在孤東采油區(qū)以北的五號樁附近引進互花米草進行繁育種植，并獲得成功。至2007年互花米草面積已達571．6 hm2［14］。在這一時間段中，互花米草的面積和分布經(jīng)歷了一個動態(tài)變化的過程［4］。田家怡等［6］的研究表明，在五號樁區(qū)域，互花米草入侵后，以此為棲息地的傳統(tǒng)鳥類的種樹減少了約75%。若不加控制，將會對保護區(qū)的生態(tài)安全造成極大的威脅［5，19，20］。

1．2 數(shù)據(jù)

(1)遙感影像數(shù)據(jù)

收集的24景遙感影像和遙感影像傳感器參數(shù)詳細信息分別如表1和表2所示。

表1 遙感影像類型及成像時間Tab．1 Landsat remote sensing images and their acquired dates

圖1 研究區(qū)域位置圖Fig．1 Location of the study area

表2 Landsat5 TM，Landsat7 ETM+和Landsat 8 OLI傳感器參數(shù)［21］Tab．2 Sensor parameters of Landsat 5 TM，Landsat7 ETM+and Landsat 8 OLI

由表1和表2，所使用的衛(wèi)星遙感數(shù)據(jù)完整覆蓋了互花米草在現(xiàn)代黃河三角洲引種至今的24年，除1990年(6月中旬)、1994年(10月上旬)和1999年(6月下旬)外，有21個年份的影像都是7～9月互花米草生長旺盛季節(jié)獲取的，利于對互花米草分布范圍和變遷情況的準(zhǔn)確監(jiān)測。21景影像中9月份的數(shù)據(jù)占2/3，原因是在黃河三角洲區(qū)域9月份處于傳統(tǒng)意義上的秋季，少雨少云，易獲得高質(zhì)量的遙感影像數(shù)據(jù)。本研究中，為使從不同時相影像中獲得的互花米草的遙感識別結(jié)果之間具有可比性，使用上述各年份影像中具有相同空間分辨率(30 m)的多光譜數(shù)據(jù)進行分析。由于2005年之后因Landsat5衛(wèi)星的超期服役，數(shù)據(jù)漸少，進而以Landsat7 ETM+數(shù)據(jù)作為補充，但其2003年之后的影像因SLC停止工作而出現(xiàn)數(shù)據(jù)條帶的缺失，故需在互花米草自動分類提取后，對結(jié)果進行后處理。

(2)現(xiàn)場數(shù)據(jù)

分別于2013年5月、9月和2014年5月三次對現(xiàn)代黃河三角洲區(qū)域進行了現(xiàn)場踏勘，應(yīng)用多種手段實地獲取了互花米草生長區(qū)域的現(xiàn)場數(shù)據(jù):應(yīng)用便攜式地物光譜儀和成像光譜儀獲取了互花米草光譜數(shù)據(jù)，作為互花米草衛(wèi)星遙感信息提取的依據(jù);應(yīng)用2014年無人機獲取的互花米草分布區(qū)域位置等信息，配合現(xiàn)場抵近照片，作為互花米草分布遙感信息提取結(jié)果的驗證數(shù)據(jù)。

2 方法

2．1 互花米草的遙感解譯

互花米草具有耐鹽、耐淹、耐掩埋的特點，可生長于海岸帶的潮間帶的高潮帶上部至中潮帶下部［22-25］，具有很強的入侵性［26］，屬于先鋒物種，我國本土特別是北方不存在與之在區(qū)域上競爭的植被物種。在我國多分布在淤泥質(zhì)潮灘上，受潮水周期性沖刷的潮間帶和潮下帶。故在遙感影像中，互花米草存在的區(qū)域多濕度較大，且會有存在于海水中或被水淹沒的情況?；诨セ撞莸纳鲜錾L特點，首先對遙感影像進行特征波段選取和組合，根據(jù)互花米草植被特殊的生長環(huán)境，對其進行纓帽線性變換，生成亮度、綠度和濕度三個分量;同時，選擇對葉綠素含量(即植被蓋度)和水分敏感(常處于淺水或濕度較大的潮間帶)的第4和第5波段，對于OLI圖像，則采用對應(yīng)的第5和第6波段，作為特征波段。應(yīng)用以上5個波段的組合，結(jié)合2013年的現(xiàn)場踏勘的數(shù)據(jù)資料選擇訓(xùn)練樣本，以SVM為分類器，進行地面覆蓋類型分類和互花米草信息提取。利用無人機獲取的互花米草分布和范圍數(shù)據(jù)、現(xiàn)場照片數(shù)據(jù)和互花米草光譜數(shù)據(jù)，對互花米草遙感信息提取結(jié)果進行驗證和修正。

由于少部分影像受云霧等影響，特別是對于因SLC設(shè)備停止工作而導(dǎo)致的條帶缺失的Landsat 7 ETM+影像，需要在SVM自動分類之后對分類結(jié)果進行后處理。

2．2 互花米草變遷分析

現(xiàn)代黃河三角洲海岸潮間帶隨時間蝕淤變化較為劇烈［27］，且因處于黃河入?？陂T處，又受勝利油田石油開發(fā)的影響，不同區(qū)域的變化情況也不一。故本研究擬分別按時間和空間維度分析互花米草變遷的特征。

在時間上，擬根據(jù)互花米草面積變化特點，劃分時間段開展分析工作;在空間上，擬根據(jù)互花米草不同生長區(qū)域的蝕淤變遷特點，劃分子區(qū)域，如圖1中右圖所示。

在互花米草分布變化分析的基礎(chǔ)上，結(jié)合1990年以來的黃河來水來沙量、河口尾閭改道蝕淤、風(fēng)暴潮和氣象數(shù)據(jù)等歷史資料，分析其變化原因。

3 結(jié)果與分析

3．1 互花米草面積和分布范圍監(jiān)測結(jié)果

圖2 互花米草分布范圍監(jiān)測結(jié)果示意圖，底圖為2013年7月15日Landsat OLI綠波段圖像Fig．2 The Spartina alterniflora distribution monitoring results from 1990 to 2013

圖3 研究區(qū)互花米草面積監(jiān)測結(jié)果Fig．3 Spartina alterniflora areas monitoring results from 1990-2013

總體上，自1990年至2013年的24年間，現(xiàn)代黃河三角洲互花米草面積和分布范圍變化劇烈。首先，在分布范圍上逐年擴大，如圖2所示，1990年僅在最先引種的五號樁和孤東采油區(qū)東南側(cè)有分布，而從2007年開始，分布范圍擴大到了現(xiàn)行黃河入海流路東南側(cè)。其次，互花米草分布總面積隨時間波動較大(圖3)。在時間上，可分為4個不同的變化階段，分別是1990引種到1994年的第一次減少階段，1994年到2001年的短暫恢復(fù)后第二次減少階段，2001年到2006年的短暫恢復(fù)后第三次減少階段和2007年至2013年的擴張階段。

在空間上，根據(jù)現(xiàn)代黃河三角洲24年來的互花米草分布變化特點，可將互花米草分布區(qū)劃分成4個子區(qū)域(圖2)。其中區(qū)域A為最先引種互花米草的區(qū)域，在整個研究時段內(nèi)位置最為穩(wěn)定，但面積卻在不斷的減小;區(qū)域B為1990年之后海岸不斷受侵蝕后退的區(qū)域;區(qū)域C為由于黃河入?？诳陂T不斷發(fā)生變化而極不穩(wěn)定的區(qū)域;區(qū)域D為近年來不斷向海淤長的區(qū)域。4個區(qū)域不同的蝕淤變化特點會給互花米草24年來的發(fā)展變化帶來截然不同的影響。

3．2 互花米草面積長時間變化的原因分析

圖4 20世紀(jì)90年代黃河下游斷流天數(shù)Fig．4 The Yellow River drying up days in one yearfrom 1990 to 1999

圖5 1990-2010年黃河汛期(7月～10月份)徑流量和輸沙量(利津站統(tǒng)計值)Fig．5 Runoff and sediment measurements in Yellow River flood seasons from 1990 to 2010(July to October)at Lijin station

表3 20世紀(jì)90年代以來黃河三角洲的風(fēng)暴潮災(zāi)害［28］Tab．3 Yellow River Delta storm surges in 1990s

1990年互花米草在黃河三角洲引種成功后，到1992年短時間內(nèi)達到了整個研究時段內(nèi)面積的極值，約1 137 hm2。而此前這一時段內(nèi)，研究區(qū)黃河入海徑流量較為充沛(圖5)，斷流天數(shù)較少(圖4)，同時，這一時段內(nèi)也無明確記載的海洋災(zāi)害(如風(fēng)暴潮等)，給互花米草的生長營造了一個相對理想的環(huán)境。1992年，黃河三角洲區(qū)域遭受了一次嚴(yán)重的風(fēng)暴潮(表3)，導(dǎo)致1993年監(jiān)測到互花米草分布面積有所下降。而1993年和1994年黃河的斷流日期處在互花米草生長旺盛的7月到10月份，進一步造成了互花米草面積的連續(xù)下降，到1994年降到了這一時期的低值點，如圖3所示。

經(jīng)過1995年和1996年兩年的恢復(fù)，到1997年，研究區(qū)又遭受了一次較嚴(yán)重的風(fēng)暴潮過程，同時，該年度黃河斷流天數(shù)達到了歷史極值的226天，這導(dǎo)致了從該年度開始，至2001年，互花米草面積在研究區(qū)不斷減少，以致在2000年和2001年的遙感影像中，幾乎無法識別明顯呈規(guī)模的互花米草。2000年前后，黃河水量達到了歷史極低值點，入海泥沙的不足直接導(dǎo)致了黃河口北岸的海岸侵蝕，這是造成這一時期互花米草退化的原因之一。

2002年之后，黃河開始大規(guī)模的在小浪底以下流域?qū)嵤┱{(diào)水調(diào)沙，此后雖黃河未再斷流，但水沙量依舊不高，加上2003年的又一次嚴(yán)重風(fēng)暴潮，自2004年，互花米草面積又開始了新一輪的減少。

到2007年，在沒有較大災(zāi)害影響的情況下，互花米草面積正不斷得到恢復(fù)，并呈爆發(fā)性增長的態(tài)勢。

綜合分析互花米草面積的變化，其和這一時段內(nèi)黃河河口段汛期徑流量和汛期輸沙量具有較大的相關(guān)性，經(jīng)線性回歸分析，互花米草面積與黃河汛期(7月～10月份)總徑流量和輸沙量之間的相關(guān)系數(shù)R2為0．56(P＜0．01)。

圖6 不同區(qū)域互花米草面積演變Fig．6 The area changes of Spartina alterniflora in 4 sub-regions

3．3 不同區(qū)域互花米草面積的對比分析

區(qū)域A位于勝利油田孤東采油區(qū)五號樁附近，通過遙感監(jiān)測發(fā)現(xiàn)，互花米草分布面積于1992年達到極大值，之后處于不斷的侵蝕過程中。根據(jù)張帆等［4］的研究，該區(qū)域互花米草的退化歸因于貝殼堤的不斷發(fā)育，灘涂中貝殼密度大的區(qū)域，互花米草生長情況較差。最終形成了如圖2區(qū)域A中2010年的新月形分布，且自2006年以來，面積和分布變化均較小。同時，在1994年、1997年和2000年區(qū)域A中互花米草近于消失(圖6)。對于1994年，概因所用圖像為10月17日獲取，處深秋季節(jié)，互花米草在遙感影像中已不明顯;對于1997年，該年度8月19日發(fā)生了一次特大風(fēng)暴潮，導(dǎo)致了互花米草被掩埋而在遙感影像中未被發(fā)現(xiàn);進入2000年后黃河下游實現(xiàn)了常年不斷流，但在2000年到2002年，黃河利津水文站觀測的年流量和6月份、7月份流量均達到了不斷流年份歷史同期的最低值［29］，同時，根據(jù)氣象觀測資料［30］，2000年到2002年三角洲地區(qū)的全年降水量和6月份、7月份降水量均為1961年來最低值，在淡水注入和降水都較小的情況下，該時間段內(nèi)互花米草的生存環(huán)境鹽度升高，抑制了該年份互花米草的生長。

在區(qū)域B中(圖2)，互花米草的面積和分布表現(xiàn)出了較大的不穩(wěn)定性。2000年之前的面積明顯大于之后，并于1997年達到了極值，在1994年和1999年有兩次面積的急劇下降。于2000年達到低谷之后，從2007年開始分布面積逐年增大。該區(qū)域互花米草分布變化的原因，除上述的黃河水沙量變化和降水問題外，還包括該區(qū)域常年遭受的嚴(yán)重侵蝕。僅從互花米草的位置來看(圖2)，24年來其向海一側(cè)的最外緣向陸地方向(西南)移動了約2．6 km。同時這期間還經(jīng)歷了幾次較大的風(fēng)暴潮，綜合起來導(dǎo)致了區(qū)域B中互花米草分布面積的起伏變化。

區(qū)域C在4個區(qū)域中受黃河尾閭影響最大。1996年，黃河由清8汊流路改道，形成了最新的向東北方向發(fā)展的新三角洲陸地區(qū)域，這一改變使原本在1995年較為茂盛的互花米草因河流沖刷而消失殆盡。而自2002年始，該區(qū)域互花米草分布面積得到了恢復(fù)，由黃河汛期徑流量監(jiān)測結(jié)果可見(圖5)，自2002年，黃河開始了調(diào)水調(diào)沙，使得該年度之后的每年6月底開始，在互花米草生長季之初，口門區(qū)域得到了大量淡水的集中補充。淡水的注入既緩解了2000年以來的高鹽度，又帶來了豐富的陸源營養(yǎng)物質(zhì)。河口口門區(qū)域的互花米草到2004年開始了逐年的衰退，從之后的遙感影像中發(fā)現(xiàn)，原因是該區(qū)域該時段口門區(qū)域沖淤變化劇烈，每年都有不同的流路分叉發(fā)育，導(dǎo)致區(qū)域C在這一時期內(nèi)得不到一個穩(wěn)定的生長環(huán)境，而逐漸消失。

對于區(qū)域D，2004年之前，該區(qū)域一直是海域，2004年之后，經(jīng)過3年的發(fā)展，于2007年的圖像中發(fā)現(xiàn)了第一片成規(guī)模的互花米草。之后至今，黃河口門流路不斷向北發(fā)展，且調(diào)水調(diào)沙使該區(qū)域泥沙量補充豐富，為互花米草的生長、發(fā)展和擴張營造了一個相對穩(wěn)定的環(huán)境，致使區(qū)域D中的互花米草位置穩(wěn)定，面積不斷擴大。

綜合上述子區(qū)域的變遷特點，除區(qū)域D由于近口門可得到大量的泥沙補充而暫未出現(xiàn)侵蝕外，其余3個子區(qū)域的互花米草生長均未有效遏制海岸的侵蝕，表現(xiàn)為區(qū)域A和區(qū)域B的互花米草面積在24年里不斷的縮小和后退，而區(qū)域C雖因河口改道而消失(在生長變遷初期也發(fā)生了侵蝕)。

4 討論

對于引起互花米草面積隨時間變化的原因分析，我們應(yīng)用了黃河水沙量、黃河改道、風(fēng)暴潮和氣象條件等4種佐證資料。通過分析發(fā)現(xiàn)，互花米草的變化并非僅由其中1種因素影響而發(fā)生，但由于后3種原因難以量化，故未作統(tǒng)計意義上的與面積變化的相關(guān)性分析。

本文所用數(shù)據(jù)少部分存在云霧覆蓋，還有部分Landsat 7 ETM+影像存在條帶缺失，對于這部分影像，不可避免的會造成統(tǒng)計誤差。應(yīng)用遙感技術(shù)開展互花米草的監(jiān)測工作，有以下2方面的優(yōu)勢:其一，能高效、準(zhǔn)確、客觀的獲取不同歷史時期的互花米草分布和變化狀況;其二，可有效的完成難以進入?yún)^(qū)域的互花米草監(jiān)測，特別是對于黃河三角洲這樣潮灘廣闊、難以實施現(xiàn)場測量的區(qū)域。同時也存在2方面的技術(shù)限制:第一，受傳感器重訪周期和云霧影響，難以在一個長時間序列里獲取到互花米草同一物候期的影像，在一定程度上給植被的識別和面積判定帶來影響;第二，基于歷史影像的信息提取結(jié)果往往缺乏驗證數(shù)據(jù)，導(dǎo)致信息提取結(jié)果存在一定的不確定性。

5 結(jié)論

基于自互花米草在現(xiàn)代黃河三角洲引種的1990年到2013年共24景Landsat衛(wèi)星遙感影像，結(jié)合2013年和2014年的現(xiàn)場踏勘數(shù)據(jù)，完成了互花米草在現(xiàn)代黃河三角洲的分布和變遷監(jiān)測;在此基礎(chǔ)上，根據(jù)互花米草隨時間變化的特點和其在不同分布區(qū)域的變化特征，在時間上分4個階段和在空間上分4個區(qū)域?qū)セ撞莸淖冞w進行了分析，并結(jié)合黃河水沙量、河口改道、風(fēng)暴潮和氣象等數(shù)據(jù)資料，對變遷原因進行了分析。主要結(jié)論包括:

24年來互花米草在現(xiàn)代黃河三角洲分布范圍逐年擴大，但總面積隨時間出現(xiàn)明顯波動，最大面積出現(xiàn)在1992年，約1 137 hm2，而在部分年份如1994年、2000年和2006年互花米草出現(xiàn)了較大規(guī)模的漸少，2007年至今呈爆發(fā)性增長。

在區(qū)域分布上，首先引種互花米草的五號樁附近區(qū)域A的互花米草一直處于退化狀態(tài)，現(xiàn)在面積僅為1992年峰值時的1/10;現(xiàn)行黃河入?？跂|南側(cè)的區(qū)域D的互花米草自2007年發(fā)現(xiàn)至今呈爆發(fā)性生長態(tài)勢;位于孤東采油區(qū)東南側(cè)的區(qū)域B和位于黃河現(xiàn)行口門處的區(qū)域C，同時受到海岸蝕淤和河口改道的影響，隨時間呈不規(guī)律變化狀態(tài)。整體上，互花米草的引種并未達到有效的固岸或者促淤的作用。

互花米草在現(xiàn)代黃河三角洲的發(fā)展和變遷與黃河水沙輸送量有較大相關(guān)性，決定系數(shù)R2達到了0．56(P＜0．01)，口門改道、海岸侵蝕和風(fēng)暴潮都是造成互花米草退化的直接原因。

［1］ CHUNG C H．Forty years of ecological engineering withSpartinaplantations in China［J］．Ecological Engineering，2006，27(1):49-57．

［2］ 左平，劉長安，趙書河，等．大米草屬植物在中國海岸帶的分布現(xiàn)狀［J］．海洋學(xué)報，2009，31(5):101-111．ZUO Ping，LIU Changan，ZHAO Shuhe，et al．Distribution ofSpartinaplants in China coastal zone［J］．Acta Oceanologica Sinica，2009，31(5):101-111．

［3］ 中華人民共和國環(huán)境保護部．關(guān)于發(fā)布中國第一批外來入侵物種名單的通知［EB/OL］．2003-01-10．Ministry of Environmental Protection of PRC．The notification of China’s first batch invasive alien species list［EB/OL］．2003-01-10．

［4］ 張帆，劉長安，姜洋．灘涂鹽沼濕地退化機制研究［J］．海洋開發(fā)與管理，2008，8:99-101．ZHANG Fan，LIU Changan，JIANG Yang．Study of the tidal salt marsh wetland degradation mechanisms［J］．Ocean Development and Management，2008，8:99-101．

［5］ 楊光，張錫義，宋志文．黃河三角洲地區(qū)大米草入侵與防治對策［J］．青島建筑工程學(xué)院學(xué)報，2005，26(2):57-59．YANG Guang，ZHANG Xiyi，SONG Zhiwen．Biological invasion countermeasure ofSpartina sppin the Yellow River Delta Coast［J］．Journal of Qingdao Institute of Architecture and Engineering，2005，26(2):57-59．

［6］ 田佳怡，于祥，申保忠，等．黃河三角洲外來入侵物種大米草對灘涂鳥類的影響［J］．中國環(huán)境管理干部學(xué)院學(xué)報，2008，18(3):87-90．TIAN Jiayi，YU Xiang，SHEN Baozhong，et al．Effect of an alien invasive speciesSpartina anglicaon birds in Shoal in Yellow River Delta［J］．Journal of EMCC，2008，18(3):87-90．

［7］ GRAY A J，MARSHALL D F，RAYBOULD A F．A century of evolution inSpartina anglica［J］．Advanced in Ecological Research，2004，21:1-62．

［8］ FOSS S．Spartina:threat to Washington’s saltwater habitat［R］．Washington State Department of Agriculture Bulletin．Olympia．1992．

［9］ DAEHLER C C，STRONG D R．Status，prediction and prevention of introduce cord grassSpartina spp．invaisions in pacific estuaries，USA［J］．Biological Conservation，1996，78:51-58．

［10］ GOSS-CUSTARD J D，CLARKE R T，DURELL S E，et al．Population consequences of winter habitat loss in a migratory shore bird II Model prediction［J］．Journal of Applied Ecology，1995，32(2):337-351．

［11］ 張東，陳小勇．應(yīng)全面評估大米草在我國沿海地區(qū)擴展繁殖帶來的利弊［J］．上海建設(shè)科技，2005，1:35-37．ZHANG Dong，CHEN Xiaoyong．Evaluation on advantages and disadvantages of spartina expansion along China coast［J］．Journal of Shanghai Construction Science and Technology，2005，1:35-37．

［12］ 孫衛(wèi)東，彭子成．治理黃河三角洲海岸退蝕的生物措施-米草生態(tài)防護工程［J］．中國地質(zhì)災(zāi)害與防治學(xué)報，1996，(3):45-48．SUN Weidong，PENG Zicheng．An useful biological management on coastal protection in yellow river delta:ecological engineering ofSpartina［J］．The Chinese Journal of Geological Hazard and Control，1996，7(3):45-48．

［13］ 石東里，田家怡，陳印萍．黃河三角洲外來入侵物種大米草的生物學(xué)與生態(tài)學(xué)特性［J］．濱州學(xué)院學(xué)報，2009，25(3):27-32．SHI Dongli，TIAN Jiayi，CHEN Yinping．Biological and ecological characteristics of an invasive alien speciesSpartinain Yellow River Delta［J］．Journal of Binzhou University，2009，25(3):27-32．

［14］ 于祥，田家怡，李建慶，等．黃河三角洲外來入侵物種大米草的分布面積與擴展速度［J］．海洋環(huán)境科學(xué)，2009，28(6):684-686．YU Xiang，TIAN Jiayi，LI Jianqing，et al．Distribution area and expanding speed of alien invasive speciesSpartina sppin Yellow River Delta［J］．Marine Environmental Science，2009，28(6):684-686．

［15］ 呂培茹，田家怡，申保忠．外來生物米草入侵對黃河三角洲鹽沼生物的影響［J］．海洋湖沼通報，2010，3:157-165．LV Peiru，TIAN Jiayi，SHEN Baozhong．Effect ofSpartina sppinvasion on saltmarsh community of yellow river delta［J］．Transactions of Oceanology and Limnology，2010，3:157-165．

［16］ 朱士文，潘秀蓮，李秀啟，等．外來物種米草對黃河三角洲生態(tài)環(huán)境的影響［J］．山東農(nóng)業(yè)科學(xué)，2012，44(3):73-75．ZHU Shiwen，PAN Xiulian，LI Xiuqi，et al．Effects of exotic speciesSpartina anglicaon ecological environment of the Yellow River Delta［J］．Shandong Agricultural Science，2012，44(3):73-75．

［17］ 張曉龍，李培英．現(xiàn)代黃河三角洲的海岸侵蝕及其環(huán)境影響［J］．海洋環(huán)境科學(xué)，2008，27(5):475-479．ZHANG Xiaolong，LIPeiying．The coastal erosion and its environmental impacts of the modern Yellow River delta［J］．Ocean Environmental Sciense，2008，27(5):475-479．

［18］ 趙全基．我國海岸侵蝕災(zāi)害及對策［J］．災(zāi)害學(xué)，1992，7(4):39-42．ZHAO Quanji．China coastal erosion hazards and countermeasures［J］．Journal of Catastrophology，1992，7(4):39-42．

［19］ CHEN Z Y，Li B，ZHONG Y，et al．Local competitive effects of introducedSpatina alternifloraonScirpusmariqueterat Dongtan of Chongming island，the Yangtze River estuary and their potential ecological consequences［J］．Hydrobiologia，2004，528(1-3):99-106．

［20］ GAN X，ZHANG K，MA Z，et al．The effect of invasions of the grassSpatina alternifloraon wintering birds on Chongming Island，Dongtan Reserve，China［J］．Journal of Ornithology，2006，147(25):169．

［21］ Landsat Mission．http://landsat．usgs．gov/about_project_descriptions．php．2013．

［22］ 徐國萬，卓榮宗．我國引種互花米草的初步研究［J］．南京大學(xué)學(xué)報，1985，6(3):212-225．XU Guowan，ZHUO Rongzong．The preliminary study of the introduction ofSpartina［J］．Journal of Nanjing University．1985，6(3):212-225．

［23］ HOWESB L，DACEY J W H，GOEHRINGER D D．Factors controlling the growth form ofSpartina alterniflora:feedbacks between above-ground production，sediment oxidation，nitrogen and salinity［J］．Journal of Ecology，1986，74:881-898．

［24］ 陳吉余，時鐘．中國濱海沼澤的初步研究［J］．地理科學(xué)，1999，10(1):58-68．CHEN Jiyu，SHI Zhong．Preliminary study on Chinese coastal marshes［J］．Scientia Geographica Sinica，1999，10(1):58-68．

［25］ 劉永學(xué)，陳君，張忍順，等．江蘇海安鹽沼植被演替的遙感圖像分析［J］．農(nóng)村生態(tài)環(huán)境，2001，17(3):39-41．LIU Yongxue，CHEN Jun，ZHANG Renshun，et al．Analysis of salt marsh vegetation succession by remote sensing images in Jiangsu Hai’an［J］．Journal of Ecology and Rural Environment，2001，17(3):39-41．

［26］ 陳中義，李博，陳家寬．米草屬植被入侵的生態(tài)后果及管理對策［J］．生物多樣性，2004，12(2):280-289．CHEN Wenyi，LI Bo，CHEN Jiakuan．The ecological consequences and management countermeasures of invasionSpartinavegetation［J］． Biodiversity Science，2004，12(2):280-289．

［27］ MENG X M，JIA Y G，SHAN H X，et al．An experimental study on erodibility of intertidal sediments in the Yellow River delta［J］．International Journal of Sediment Reasearch，2012，27(2):240-249．

［28］ 張曉龍，李培英，劉月良．黃河三角洲風(fēng)暴潮災(zāi)害及其對濱海濕地的影響［J］．自然災(zāi)害學(xué)報，2006，15(2):10-13．ZHANG Xiaolong，LI Pei ying，LIU Yueliang．Storm surge disaster and its impact on coastal wetlands in Yellow River Delta［J］．Journal of Natural Disasters，2006，15(2):10-13．

［29］ 王學(xué)金，陳立強，宋玉敏，等．三角洲地區(qū)黃河水資源利用現(xiàn)狀及對策［J］．水利規(guī)劃與設(shè)計，2013，1:18-19．WANG Xuejin，CHEN Liqiang，SONG Yumin，et al．Yellow River Delta water resources utilization and countermeasures［J］．Water Resources Planning and Design，2013，1:18-19．

［30］ 林琳，劉健，陳學(xué)群，等．黃河三角洲1961-2000年水資源時空變化特征［J］．水資源保護，2012，28(1):29-37．LIN Lin，LIU Jian，CHEN Xuequn，et al．The changing features of the water resources of the Yellow River delta in 1961-2000［J］．Water Resources Production，2012，28(1):29-37．