黃河三角洲地區(qū)沿海養(yǎng)殖池塘空間演變階段性劃分及特征分析

王紫陽(yáng)　吳曉青　豆曉　王曉杰

摘要基于1984—2020年遙感影像數(shù)據(jù)，獲取連續(xù)長(zhǎng)時(shí)間序列的養(yǎng)殖池塘空間分布數(shù)據(jù)，并采用時(shí)空熱點(diǎn)分析和PLR_FP方法對(duì)黃河三角洲地區(qū)養(yǎng)殖池塘空間演變的階段性進(jìn)行劃分。結(jié)果表明，黃河三角洲地區(qū)養(yǎng)殖池塘面積呈現(xiàn)波動(dòng)增長(zhǎng)態(tài)勢(shì)，年均增加29.3 km²，時(shí)空增長(zhǎng)模式多樣化，且具有明顯的空間分異性；增長(zhǎng)熱點(diǎn)分布在濱州沿海、東營(yíng)河口區(qū)西部、墾利黃河口南部，成為目前黃河三角洲地區(qū)池塘養(yǎng)殖重點(diǎn)分布區(qū)，而萊州灣沿岸池塘養(yǎng)殖規(guī)模逐漸萎縮。養(yǎng)殖池塘的空間演變過(guò)程可劃分為起步期（1984—1992年）、平穩(wěn)期（1993—1998年）、擴(kuò)張期（1999—2013年）和縮減期（2014—2020年）4個(gè)階段，不同階段養(yǎng)殖池塘?xí)r空增長(zhǎng)模式和養(yǎng)殖業(yè)發(fā)展特征差異明顯。

關(guān)鍵詞養(yǎng)殖池塘；空間演變；階段劃分；變化特征；黃河三角洲

中圖分類號(hào)S967.4文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)0517-6611（2023）11-0050-06

doi：10.3969/j.issn.0517-6611.2023.11.012開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Stage Division and Characteristics Analysis of the Spatial Evolution of Coastal Aquaculture Ponds in Yellow River Delta

WANG Zi-yang WU Xiao-qing DOU Xiao4 et al（1. Yantai Institute of Coastal Zone Research， Chinese Academy of Sciences， Yantai，Shandong 264003;2. University of Chinese Academy of Sciences， Beijing 100049;3. Key Laboratory of Coastal Environmental Processes and Ecological Remediation， Chinese Academy of Sciences， Yantai，Shandong 264003;4. College of Public Administration， Shandong Agricultural University，Taian，Shandong 271000）

AbstractBased on the remote sensing image data from 1984 to 2020， this study obtained the spatial distribution data of aquaculture ponds in a long-term continuous sequence，and the spatio-temporal hotspot analysis and PLR_FP method were used to divide the stages of the spatial evolution of aquaculture ponds in Yellow River Delta.The results showed that the area of aquaculture ponds in the Yellow River Delta showed a fluctuating growth trend， with an average annual increase of 29.3 km2. The spatiotemporal growth patterns were diverse， with obvious spatial differentiation. Growth hotspots were distributed along the coast of Binzhou， the west of Dongying Estuary and the south of Kenli Yellow River Estuary， which had become the key distribution areas for pond aquaculture in the Yellow River Delta， while the scale of pond aquaculture along the Laizhou Bay was gradually shrinking. The spatial evolution process of aquaculture ponds could be divided into four stages： the initial stage （1984-1992）， the stable stage （1993-1998）， the expansion stage （1999-2013） and the shrinking stage （2014-2020）， and there were significant differences in the spatiotemporal growth patterns and aquaculture characteristics of aquaculture ponds at different stages.

Key wordsAquaculture ponds;Spatial evolution;Staged division;Change characteristics;Yellow River Delta

池塘養(yǎng)殖是海水養(yǎng)殖的一種重要模式，作為“優(yōu)質(zhì)糧食”的生產(chǎn)場(chǎng)所，其生態(tài)災(zāi)害風(fēng)險(xiǎn)較低，經(jīng)濟(jì)收益較高，對(duì)我國(guó)“藍(lán)色糧倉(cāng)”建設(shè)具有重要意義［1］?，F(xiàn)代黃河三角洲及萊州灣沿岸為典型的三角洲海岸和淤泥質(zhì)海岸，灘涂資源豐富、池塘養(yǎng)殖歷史悠久、養(yǎng)殖規(guī)模大，是山東省重要的養(yǎng)殖基地［2］。近幾十年來(lái)，隨著社會(huì)經(jīng)濟(jì)發(fā)展水平的提高和養(yǎng)殖技術(shù)的進(jìn)步，黃河三角洲地區(qū)池塘養(yǎng)殖不斷向渤海擴(kuò)張，占用了大量自然灘涂濕地，造成濱海濕地萎縮、生物多樣性下降、土壤次生鹽漬化、碳封存環(huán)境破壞等問(wèn)題［3-7］。與此同時(shí)，養(yǎng)殖池塘也是重要的人工濕地，具有補(bǔ)充地下水和調(diào)節(jié)氣候等生態(tài)功能［8-9］，科學(xué)合理地進(jìn)行池塘養(yǎng)殖還能夠成為踐行碳中和戰(zhàn)略的重要負(fù)排放途徑［10-11］。在推動(dòng)黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展戰(zhàn)略背景下，實(shí)現(xiàn)池塘養(yǎng)殖功能優(yōu)化布局，成為化解濱海濕地生態(tài)保護(hù)與高效生態(tài)化養(yǎng)殖之間矛盾的重要管理手段之一。而研究沿海池塘養(yǎng)殖時(shí)空演變特征作為其優(yōu)化布局的一項(xiàng)重要基礎(chǔ)性工作，有助于揭示該地區(qū)海水養(yǎng)殖業(yè)發(fā)展歷程，獲取地區(qū)池塘養(yǎng)殖的獨(dú)特性，為黃河三角洲地區(qū)灘涂濕地池塘養(yǎng)殖空間優(yōu)化布局和生態(tài)保護(hù)修復(fù)工程規(guī)劃建設(shè)提供數(shù)據(jù)和理論支持。

當(dāng)前學(xué)者多從土地利用和景觀格局變化的角度揭示黃河三角洲地區(qū)養(yǎng)殖池塘的時(shí)空演變特征及其對(duì)濱海濕地的生態(tài)影響［12-15］，但是由于其時(shí)空分布數(shù)據(jù)不具備連續(xù)性，或僅使用經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)，造成對(duì)養(yǎng)殖池塘?xí)r空演變特征和變化趨勢(shì)的把握不夠準(zhǔn)確，不能對(duì)海水養(yǎng)殖業(yè)發(fā)展過(guò)程作出有效判斷。該研究基于遙感、相關(guān)專項(xiàng)調(diào)查資料等數(shù)據(jù)，重點(diǎn)揭示黃河三角洲地區(qū)1984—2020年池塘海水養(yǎng)殖時(shí)空演變過(guò)程，劃分其空間演變階段，探討不同階段空間增長(zhǎng)特征及其驅(qū)動(dòng)因素，以期為池塘養(yǎng)殖可持續(xù)管理與開(kāi)發(fā)和環(huán)境保護(hù)提供技術(shù)支持。

1資料與方法

1.1研究區(qū)概況研究區(qū)位于山東半島北部，西起魯冀交界的漳衛(wèi)新河河口，東至煙臺(tái)、濰坊交界的膠萊河口，地理位置為117°76′～119°61′E、36°92′～38°28′N，包括黃河三角洲、萊州灣自海岸線向陸延伸15 km的海岸帶地區(qū)，行政區(qū)涉及濱州、東營(yíng)、濰坊3個(gè)地級(jí)市。此區(qū)域?yàn)榉凵坝倌噘|(zhì)平原海岸，地勢(shì)低平，灘涂范圍寬廣，有黃河、小清河、濰河、膠萊河等十余條淡水河流入海，近岸海域營(yíng)養(yǎng)鹽、有機(jī)質(zhì)多，餌料充盈，是山東省重要的海水養(yǎng)殖基地。

1.2數(shù)據(jù)源與預(yù)處理遙感數(shù)據(jù)源主要為L(zhǎng)andsat系列衛(wèi)星，時(shí)間跨度為1984—2020年，通過(guò)地理空間數(shù)據(jù)云（http：//www.gscloud.cn/）平臺(tái)下載獲取，共計(jì)74景影像。在ENVI中進(jìn)行預(yù)處理，包括輻射校正、大氣校正、全色波段和多光譜的融合等。輔助數(shù)據(jù)包括地形圖、歷史年份的高清影像和海岸帶專項(xiàng)調(diào)查數(shù)據(jù)等。海水養(yǎng)殖相關(guān)統(tǒng)計(jì)數(shù)據(jù)來(lái)源于省、市統(tǒng)計(jì)局歷年統(tǒng)計(jì)年鑒，部分年份缺失數(shù)據(jù)采用移動(dòng)平均法進(jìn)行補(bǔ)充。

1.3研究方法

1.3.1遙感信息提取。養(yǎng)殖池塘的提取首先借助水體指數(shù)并參照明顯的解譯標(biāo)志［16］識(shí)別其潛在區(qū)域；然后，結(jié)合地物遙感解譯判別標(biāo)志、專家經(jīng)驗(yàn)知識(shí)、海岸帶土地利用調(diào)查資料以及實(shí)地調(diào)研，采用人工目視解譯方法對(duì)養(yǎng)殖池塘和鹽田、水庫(kù)坑塘進(jìn)行區(qū)分，處理獲得1984—2020年黃河三角洲地區(qū)連續(xù)37年的養(yǎng)殖池塘分布數(shù)據(jù)。利用同時(shí)期高清谷歌影像、高分一號(hào)等高分辨率遙感數(shù)據(jù)及研究團(tuán)隊(duì)近幾年實(shí)地調(diào)查點(diǎn)位進(jìn)行解譯精度驗(yàn)證，結(jié)果顯示37年養(yǎng)殖池塘總體解譯精度均在79%以上，滿足研究需要。

1.3.2階段性劃分方法。該研究基于養(yǎng)殖池塘空間分布連續(xù)多年的遙感觀測(cè)數(shù)據(jù)，采用基于特征點(diǎn)的分段線性表示法（piecewise linear representation based feature points， PLR_FP）提取時(shí)間序列數(shù)據(jù)分段特征點(diǎn)［17］，進(jìn)行池塘養(yǎng)殖空間擴(kuò)展階段性劃分。PLR_FP方法可以概化時(shí)間序列數(shù)據(jù)的主要狀態(tài)特征，進(jìn)行有效的數(shù)據(jù)壓縮和局部噪聲過(guò)濾。具體做法：首先，識(shí)別時(shí)間序列數(shù)據(jù)的局部極值點(diǎn)；然后，采用閾值R限定前后極值點(diǎn)數(shù)據(jù)的垂直差比，將符合要求的點(diǎn)識(shí)別為不同尺度下的分段特征點(diǎn)［17-18］；最后，借助雙樣本W(wǎng)ilcoxon秩和檢驗(yàn)將不具有顯著性的特征點(diǎn)進(jìn)行排除［19］，依據(jù)最終分段線性表示的擬合系數(shù)與常數(shù)0之間的關(guān)系判斷每階段的具體趨勢(shì)。

其中，秩和檢驗(yàn)用于檢驗(yàn)2個(gè)獨(dú)立的養(yǎng)殖序列樣本是否來(lái)自相同的總體。從2個(gè)獨(dú)立總體抽取樣本并將所有觀測(cè)值由大到小排列后，每一個(gè)觀測(cè)值在混合排列中所占位置的名次稱為它的秩，當(dāng)其中某樣本秩和（statistic）較小，混合樣本為一個(gè)隨機(jī)混合樣本，秩和檢驗(yàn)中P值小于0.05，可以認(rèn)為該分段點(diǎn)前后時(shí)間序列數(shù)據(jù)具有明顯差異，可對(duì)前后2個(gè)時(shí)段內(nèi)養(yǎng)殖池塘演變特征相同這一假設(shè)提出質(zhì)疑。

1.3.3空間演變模式判別方法。對(duì)于連續(xù)性面狀空間分布數(shù)據(jù)，采用時(shí)空熱點(diǎn)分析方法量化分析近37年養(yǎng)殖池塘的空間格局類型及其出現(xiàn)頻率和時(shí)空位置［20］，揭示池塘養(yǎng)殖空間演變分異規(guī)律。

該方法首先基于ArcGIS熱點(diǎn)分析方法計(jì)算各時(shí)間截面Getis-Ord G*統(tǒng)計(jì)量并進(jìn)行顯著性檢驗(yàn)；空間權(quán)重矩陣采用鄰接標(biāo)準(zhǔn)中的“Queen”，即面要素有公共邊或結(jié)點(diǎn)即為相鄰，采用“CONTIGUITY_EDGES_CORNERS”工具計(jì)算，鄰域步長(zhǎng)設(shè)置為1年；然后，分析各區(qū)域截面Getis-Ord G*統(tǒng)計(jì)量的Mann-Kendall趨勢(shì)；其中，區(qū)域截面的構(gòu)造采用覆蓋研究區(qū)范圍的1 km格網(wǎng)；最后，結(jié)合養(yǎng)殖池塘面積遙感觀測(cè)數(shù)據(jù)，按照模式分類規(guī)則（表1）對(duì)比歸類為不同時(shí)空模式［20］。

2結(jié)果與分析

2.1養(yǎng)殖池塘?xí)r空演變總體特征從圖1可以看出，1984—2020年黃河三角洲地區(qū)沿海養(yǎng)殖池塘規(guī)模呈現(xiàn)波動(dòng)增長(zhǎng)的態(tài)勢(shì)，從期初的66.7 km²增長(zhǎng)至1 150.7 km²，增長(zhǎng)約16.3倍，年均增速為29.3 km²。其中，1986—1989、2000—2005、2010—2013年養(yǎng)殖池塘面積增速較快，并于2013年達(dá)到最高峰（1 359.9 km²）；2016—2018年養(yǎng)殖池塘規(guī)模出現(xiàn)縮減，3年共減少231.1 km²；而1990—1999年養(yǎng)殖池塘面積保持在相對(duì)穩(wěn)定的水平，可見(jiàn)養(yǎng)殖池塘規(guī)模長(zhǎng)時(shí)間序列年度增長(zhǎng)呈現(xiàn)明顯的階段性變化特點(diǎn)。

黃河三角洲地區(qū)沿海養(yǎng)殖池塘空間演變也呈現(xiàn)出明顯的時(shí)空分異特征，不同時(shí)間段養(yǎng)殖池塘增長(zhǎng)的冷熱點(diǎn)發(fā)生遷移變化，時(shí)空演變模式呈現(xiàn)多樣化（圖2）。研究期初，濱州的無(wú)棣、沾化和濰坊市的昌邑、寒亭等縣（市、區(qū)）的Gi*統(tǒng)計(jì)量較高，為區(qū)域內(nèi)沿海池塘養(yǎng)殖的起步區(qū)；隨后一段時(shí)間，各縣（市、區(qū)）養(yǎng)殖池塘增長(zhǎng)變化的Gi*統(tǒng)計(jì)量分布相對(duì)均勻，無(wú)明顯的熱點(diǎn)區(qū)，僅在壽光、河口、利津等縣（市、區(qū)）保持較低水平。但是在2000年之后，養(yǎng)殖池塘區(qū)域增長(zhǎng)的冷熱點(diǎn)分布逐漸變得兩級(jí)分明，沾化、無(wú)棣的Gi*統(tǒng)計(jì)量驟增，出現(xiàn)了新一輪的池塘養(yǎng)殖高潮，并帶動(dòng)與其毗鄰的東營(yíng)河口區(qū)池塘養(yǎng)殖的擴(kuò)張，成為“逐漸增強(qiáng)的熱點(diǎn)”，也是目前池塘養(yǎng)殖的重點(diǎn)分布區(qū)；而萊州灣西南部沿岸縣（市、區(qū)）Gi*統(tǒng)計(jì)量明顯降低，歸為“逐漸減弱的熱點(diǎn)”，因鹽化工、臨港產(chǎn)業(yè)發(fā)展和城鎮(zhèn)化建設(shè)，該區(qū)域養(yǎng)殖池塘規(guī)模大面積萎縮。養(yǎng)殖池塘空間增長(zhǎng)的冷點(diǎn)分布較為分散，其中“逐漸增強(qiáng)的冷點(diǎn)”和“分散的冷點(diǎn)”位于黃河三角洲國(guó)家級(jí)自然保護(hù)區(qū)，說(shuō)明自然保護(hù)地建設(shè)管理對(duì)養(yǎng)殖池塘空間擴(kuò)展具有明顯的約束作用。

2.2養(yǎng)殖池塘規(guī)模變化階段性劃分為準(zhǔn)確劃分養(yǎng)殖池塘規(guī)模變化的階段，分析不同階段空間演變模式及特征，該研究采用PLR_FP方法對(duì)養(yǎng)殖池塘面積連續(xù)遙感觀測(cè)數(shù)據(jù)進(jìn)行特征點(diǎn)識(shí)別（圖1），通過(guò)調(diào)整常數(shù)R閾值（0≤R≤20.39）的大小，將符合要求的特征點(diǎn)識(shí)別為不同分割等級(jí)下的分界點(diǎn)（表2）。由于從第五級(jí)分割等級(jí)開(kāi)始，不顯著點(diǎn)多于50%，故保留前4個(gè)等級(jí)的特征點(diǎn)，進(jìn)行秩和檢驗(yàn)后，最終確定養(yǎng)殖池塘面積變化分段特征顯著的特征點(diǎn)為1992、1998、2005和2013年。

以上述分段特征點(diǎn)進(jìn)行養(yǎng)殖池塘面積變化的分段線性擬合，結(jié)果顯示1984—1992、1999—2005、2006—2013年3個(gè)時(shí)段的擬合系數(shù)大于0，說(shuō)明這些時(shí)段池塘養(yǎng)殖規(guī)模呈現(xiàn)增長(zhǎng)趨勢(shì)；1993—1998年的線性擬合斜率近乎為0，池塘養(yǎng)殖規(guī)?；緵](méi)發(fā)生變化；2014—2020年的線性擬合斜率為負(fù)值，說(shuō)明此階段池塘養(yǎng)殖規(guī)模出現(xiàn)縮減。由此將黃河三角洲地區(qū)沿海養(yǎng)殖池塘空間擴(kuò)展過(guò)程劃分為4個(gè)時(shí)段，即池塘養(yǎng)殖起步期（1984—1992年）、平穩(wěn)期（1993—1998年）、擴(kuò)張期（1999—2013年）、縮減期（2014—2020年）。

2.3黃河三角洲地區(qū)養(yǎng)殖池塘空間演變階段性特征依據(jù)養(yǎng)殖池塘演變劃分的4個(gè)階段，進(jìn)一步分析不同階段黃河三角洲地區(qū)沿海池塘養(yǎng)殖空間演變特征、養(yǎng)殖業(yè)發(fā)展特點(diǎn)及驅(qū)動(dòng)因素。

（1）池塘養(yǎng)殖起步期（1984—1992年）。養(yǎng)殖池塘空間增長(zhǎng)總體上自濱州沿海地區(qū)開(kāi)始起步，隨后在濰坊沿岸多點(diǎn)興起，最終在濱州和萊州灣西南部海岸帶地區(qū)呈現(xiàn)大幅擴(kuò)張（圖3a）；至期末，養(yǎng)殖規(guī)模已達(dá)到559.33 km²，年均增速為54.74 km²。山東黃河三角洲國(guó)家級(jí)自然保護(hù)區(qū)的建設(shè)，限制了東營(yíng)市河口區(qū)與墾利區(qū)的圍海養(yǎng)殖活動(dòng)。這一時(shí)期，受益于對(duì)蝦、扇貝養(yǎng)殖技術(shù)的突破和改革開(kāi)放的良好發(fā)展環(huán)境，沿海池塘對(duì)蝦、貝類養(yǎng)殖獲得快速發(fā)展，全省漁業(yè)總產(chǎn)值達(dá)到期初的3.83倍。但是，養(yǎng)殖品種單一，養(yǎng)殖方式粗放，且由于過(guò)度追求產(chǎn)量和規(guī)模，高密度養(yǎng)殖和盲目用藥造成水環(huán)境惡化，給附近灘涂及海域造成巨大生態(tài)壓力［22］。

（2）池塘養(yǎng)殖平穩(wěn)期（1993—1998年）。養(yǎng)殖池塘空間增長(zhǎng)熱點(diǎn)規(guī)模小，且分布相對(duì)分散（圖3b），主要位于東營(yíng)和濱州沿海，養(yǎng)殖規(guī)模穩(wěn)定在547 km²左右。這一時(shí)期，在“海上山東”和黃河三角洲綜合開(kāi)發(fā)戰(zhàn)略的推動(dòng)下，經(jīng)過(guò)養(yǎng)殖結(jié)構(gòu)調(diào)整，池塘養(yǎng)殖種類逐漸多樣化，實(shí)現(xiàn)了從追求產(chǎn)量和規(guī)模向注重質(zhì)量和效益的重大轉(zhuǎn)變［23］。此階段后期受1993和1996年對(duì)蝦、櫛孔扇貝大規(guī)模死亡事件的影響，黃河三角洲地區(qū)開(kāi)始積極探索各種魚(yú)蝦貝立體混養(yǎng)、輪養(yǎng)的模式，依靠科技進(jìn)步開(kāi)展對(duì)蝦池的綜合利用，不斷挖掘生態(tài)效益［24］。

（3）池塘養(yǎng)殖擴(kuò)張期（1999—2013年）。養(yǎng)殖池塘呈現(xiàn)大面積快速擴(kuò)張態(tài)勢(shì)，養(yǎng)殖規(guī)模由1999年的542.18 km²迅速增長(zhǎng)至2013年的1 359.87 km²，年均增速54.51 km²，遠(yuǎn)高于其他時(shí)段；擴(kuò)張熱點(diǎn)區(qū)變得更為集中，濱州東部、墾利沿岸成為新的增長(zhǎng)核心區(qū)域，萊州灣南岸池塘養(yǎng)殖規(guī)?？s小，最終形成了濱州北部的廣闊灘涂以及現(xiàn)今黃河入?？谝阅系臇|營(yíng)市沿岸兩大重點(diǎn)養(yǎng)殖區(qū)（圖3c）。這一時(shí)期，多種池塘混養(yǎng)模式技術(shù)探索獲得成功及應(yīng)用，養(yǎng)殖種類日益多樣化，三疣梭子蟹、中國(guó)明對(duì)蝦、日本囊對(duì)蝦、凡納濱對(duì)蝦、刺參和鮑都具有一定規(guī)模，養(yǎng)殖生態(tài)意識(shí)進(jìn)一步增強(qiáng)，可持續(xù)、集約化發(fā)展成為各地海水養(yǎng)殖追求的新發(fā)展方向［25-26］。

（4）池塘養(yǎng)殖縮減期（2014—2020年）。2014—2016年養(yǎng)殖池塘面積縮減230.14 km²，之后基本穩(wěn)定在1 113.93 km²；養(yǎng)殖池塘空間增長(zhǎng)受到城市化、工業(yè)化進(jìn)程不斷推進(jìn)和黃河三角洲灘涂濕地生態(tài)修復(fù)工程建設(shè)的影響，壽光、寒亭、廣饒以及利津、墾利北端的池塘規(guī)?？s減（圖3d），陸域空間發(fā)展受到擠壓，而海域開(kāi)放式養(yǎng)殖獲得快速發(fā)展［27-28］。這一時(shí)期，“東參西養(yǎng)”以及各種生態(tài)混養(yǎng)模式技術(shù)獲得成功，并且在現(xiàn)代生態(tài)高效農(nóng)業(yè)示范區(qū)建設(shè)帶動(dòng)下，各地注重引導(dǎo)發(fā)展生態(tài)與經(jīng)濟(jì)利益相協(xié)調(diào)的海水養(yǎng)殖業(yè)［24，29］，黃河三角洲地區(qū)海參、凡納濱對(duì)蝦養(yǎng)殖以及青蛤、菲律賓蛤仔、褐牙鲆、縊蟶、鹵蟲(chóng)、海蜇養(yǎng)殖業(yè)蒸蒸日上，一定程度上實(shí)現(xiàn)了經(jīng)濟(jì)效益和生態(tài)效益的雙贏［30］。

3討論

該研究基于連續(xù)長(zhǎng)時(shí)間序列的養(yǎng)殖池塘面積遙感觀測(cè)數(shù)據(jù)和PLR_FP方法對(duì)黃河三角洲地區(qū)沿海養(yǎng)殖池塘空間演變過(guò)程進(jìn)行階段性劃分，識(shí)別出4個(gè)明顯的變化階段，并分析得出不同時(shí)段養(yǎng)殖池塘?xí)r空增長(zhǎng)模式和養(yǎng)殖特征具有明顯差異，較好地揭示了區(qū)域池塘養(yǎng)殖業(yè)發(fā)展歷程。階段性劃分結(jié)果比較客觀，與已有研究也存在一致性。如宋媛等［31］研究揭示1994—1999年萊州灣、渤海灣等海灣養(yǎng)殖池塘面積保持均衡變動(dòng)，2004年起因萊州灣鹽業(yè)迅速發(fā)展，養(yǎng)殖池向鹽田發(fā)生轉(zhuǎn)換，以及2009—2014年渤海灣養(yǎng)殖業(yè)重大改革帶來(lái)的快速發(fā)展等，都符合該研究區(qū)養(yǎng)殖池塘規(guī)模變化的重要時(shí)間節(jié)點(diǎn)。

漁業(yè)經(jīng)濟(jì)統(tǒng)計(jì)指標(biāo)的變化也可以在一定程度上反映地區(qū)漁業(yè)發(fā)展歷程，如許學(xué)工［32］、周鑫等［33］進(jìn)行黃河三角洲地區(qū)漁業(yè)綜合效益評(píng)估時(shí)發(fā)現(xiàn)經(jīng)濟(jì)效益在1992、1998和2006年等時(shí)間節(jié)點(diǎn)同樣呈現(xiàn)出明顯的階段性差異，此結(jié)果與該研究基于遙感觀測(cè)數(shù)據(jù)的時(shí)間分段點(diǎn)有所偏差，主要原因是漁業(yè)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)涵蓋范圍廣，而且研究區(qū)內(nèi)海水養(yǎng)殖模式多種多樣，除了池塘養(yǎng)殖，還包括工廠化養(yǎng)殖、海域底播養(yǎng)殖、海域筏式養(yǎng)殖等類型，池塘養(yǎng)殖以外的這些海水養(yǎng)殖模式所帶來(lái)的養(yǎng)殖產(chǎn)量和養(yǎng)殖效益在海水養(yǎng)殖業(yè)中占據(jù)較大比重，采用遙感技術(shù)手段無(wú)法精確獲取連續(xù)時(shí)間段高時(shí)空精度的各類養(yǎng)殖模式空間分布信息，而養(yǎng)殖產(chǎn)量和漁業(yè)總產(chǎn)值統(tǒng)計(jì)指標(biāo)更是難以剝離出池塘養(yǎng)殖所帶來(lái)的收益。

4結(jié)論

該研究基于遙感、時(shí)空熱點(diǎn)分析和PLR_FP方法，揭示了1984—2020年黃河三角洲地區(qū)沿海養(yǎng)殖池塘?xí)r空演變過(guò)程及階段性變化特征，主要結(jié)論如下：

（1）黃河三角洲地區(qū)養(yǎng)殖池塘規(guī)模呈現(xiàn)波動(dòng)中增長(zhǎng)態(tài)勢(shì)，近37年間增長(zhǎng)了16.3倍，年均增速29.3 km²，形成1984—1992年和1999—2013年兩大增長(zhǎng)高峰段，2013年養(yǎng)殖面積達(dá)到最高峰（1 359.9 km²）。養(yǎng)殖池塘空間增長(zhǎng)冷熱點(diǎn)變化具有明顯的時(shí)空分異性，空間增長(zhǎng)模式也呈現(xiàn)多樣化。養(yǎng)殖池塘空間增長(zhǎng)自濱州、濰坊沿海興起，隨后增長(zhǎng)熱點(diǎn)在北部地區(qū)由濱州自西向東轉(zhuǎn)移至河口區(qū)西部，而南部地區(qū)則由濰坊沿岸向西、向北遷移到東營(yíng)黃河口南部，最終形成了濱州沿海、東營(yíng)河口區(qū)西部、墾利黃河口南部等池塘養(yǎng)殖重點(diǎn)區(qū)，而萊州灣南岸池塘養(yǎng)殖規(guī)模逐漸萎縮。

（2）黃河三角洲地區(qū)沿海養(yǎng)殖池塘的空間演變過(guò)程可劃分為4個(gè)階段，即起步期（1984—1992年）、平穩(wěn)期（1993—1998年）、擴(kuò)張期（1999—2013年）和衰減期（2014—2020年）；由于海水養(yǎng)殖技術(shù)的進(jìn)步、漁業(yè)經(jīng)濟(jì)促進(jìn)政策和濕地生態(tài)保護(hù)措施的實(shí)施，不同階段池塘養(yǎng)殖空間演變特征、主要驅(qū)動(dòng)因素和漁業(yè)生產(chǎn)經(jīng)營(yíng)狀況有較大差異。

參考文獻(xiàn)

［1］ 楊紅生.我國(guó)藍(lán)色糧倉(cāng)科技創(chuàng)新的發(fā)展思路與實(shí)施途徑［J］.水產(chǎn)學(xué)報(bào)，2019，43（1）：97-104.

［2］ 徐巖，朱玉貴.山東省海水貝類養(yǎng)殖產(chǎn)業(yè)變化特征分析［J］.現(xiàn)代農(nóng)業(yè)科技，2020（12）：215-216，220.

［3］ 曹瓊，黃佳芳，羅敏，等.濱海沼澤濕地轉(zhuǎn)化為養(yǎng)殖塘對(duì)其碳儲(chǔ)量的影響［J］.中國(guó)環(huán)境科學(xué)，2022，42（3）：1335-1345.

［4］ 賈軍，崔正國(guó)，唐小雙，等.不同水力條件下海水人工濕地凈化效果及其影響分析［J］.漁業(yè)現(xiàn)代化，2021，48（5）：26-35.

［5］ CHI Y，SHI H H，ZHENG W，et al.Spatiotemporal characteristics and ecological effects of the human interference index of the Yellow River Delta in the last 30 years［J］.Ecological indicators，2018，89：880-892.

［6］ WANG C D，LI X，YU H J，et al.Tracing the spatial variation and value change of ecosystem services in Yellow River Delta，China［J］.Ecological indicators，2019，96：270-277.

［7］ ZHANG X J，WANG G Q，XUE B L，et al.Dynamic landscapes and the driving forces in the Yellow River Delta wetland region in the past four decades［J/OL］.Science of the total environment，2021，787［2022-03-25］.https：//doi.org/10.1016/j.scitotenv.2021.147644.

［8］ 吳志峰，曹崢，宋松，等.粵港澳大灣區(qū)濕地遙感監(jiān)測(cè)與評(píng)估：現(xiàn)狀、挑戰(zhàn)及展望［J］.生態(tài)學(xué)報(bào)，2020，40（23）：8440-8450.

［9］ 張緒良，谷東起，豐愛(ài)平，等.萊州灣南岸濱海濕地的氮、磷循環(huán)過(guò)程及調(diào)控對(duì)策［J］.中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2008，16（5）：1127-1133.

［10］ 焦念志，劉紀(jì)化，石拓，等.實(shí)施海洋負(fù)排放踐行碳中和戰(zhàn)略［J］.中國(guó)科學(xué)：地球科學(xué)，2021，51（4）：632-643.

［11］ 張安琪.海水養(yǎng)殖助力“藍(lán)碳計(jì)劃”的途徑［J］.安徽農(nóng)業(yè)科學(xué)，2022，50（2）：89-91.

［12］ 孔梅，曹惠明，高興國(guó)，等.黃河三角洲自然保護(hù)區(qū)國(guó)土空間變化遙感監(jiān)測(cè)與評(píng)價(jià)［J］.海洋環(huán)境科學(xué)，2021，40（2）：272-276.

［13］ OTTINGER M，KUENZER C，LIU G H，et al.Monitoring land cover dynamics in the Yellow River Delta from 1995 to 2010 based on Landsat 5 TM［J］.Applied geography，2013，44：53-68.

［14］ 仇志強(qiáng)，毛德華，相恒星，等.5個(gè)時(shí)期黃河流域濕地格局及其變化研究［J］.濕地科學(xué)，2021，19（4）：518-526.

［15］ 黃建濤，鄭偉，萬(wàn)年新，等.近30年來(lái)萊州灣濱海濕地景觀格局變化特征研究［J］.海洋科學(xué)，2021，45（2）：76-90.

［16］ 任玲玲，栗云召，于淼，等.1984—2015年黃河三角洲人工濕地的演變及主要驅(qū)動(dòng)因子［J］.農(nóng)業(yè)資源與環(huán)境學(xué)報(bào)，2020，37（4）：493-502.

［17］ 喻高瞻，彭宏，胡勁松，等.時(shí)間序列數(shù)據(jù)的分段線性表示［J］.計(jì)算機(jī)應(yīng)用與軟件，2007，24（12）：17-18.

［18］ HUANG L，CHEN X，ZHANG Y，et al.Identification of topic evolution：Network analytics with piecewise linear representation and word embedding［J］.Scientometrics，2022，127（9）：5353-5383.

［19］ 代穩(wěn)，呂殿青，李景保，等.1951—2014年洞庭湖水沙階段性演變特征及驅(qū)動(dòng)因素分析［J］.水土保持學(xué)報(bào)，2017，31（2）：142-150，176.

［20］ 單寶艷，張智璇，陳艷秋，等.時(shí)空格局分析方法及其實(shí)證應(yīng)用：以山東省制造業(yè)為例［J］.測(cè)繪科學(xué)技術(shù)學(xué)報(bào)，2021，38（6）：624-630，638.

［21］ GUO M，LI J，HE H S，et al.Detecting global vegetation changes using Mann-Kendal（MK）trend test for 1982-2015 time period［J］.Chinese geographical science，2018，28（6）：907-919.

［22］ 趙廣苗.海水池塘養(yǎng)殖模式的演變與發(fā)展趨勢(shì)［J］.齊魯漁業(yè)，2006，23（7）：1-2.

［23］ 李大海，潘克厚，陳玲玲.改革開(kāi)放以來(lái)我國(guó)海水養(yǎng)殖政策的演變與發(fā)展［J］.中國(guó)漁業(yè)經(jīng)濟(jì)，2008，26（3）：57-61.

［24］ 李靜.海水養(yǎng)殖污染與生態(tài)修復(fù)對(duì)策［J］.農(nóng)家參謀，2020（15）：114，116.

［25］ 董雙林.多維視角下的新時(shí)代水產(chǎn)養(yǎng)殖業(yè)發(fā)展［J］.水產(chǎn)學(xué)報(bào)，2019，43（1）：105-115.

［26］ 方建光，李鐘杰，蔣增杰，等.水產(chǎn)生態(tài)養(yǎng)殖與新養(yǎng)殖模式發(fā)展戰(zhàn)略研究［J］.中國(guó)工程科學(xué)，2016，18（3）：22-28.

［27］ 李大海，韓立民.陸海統(tǒng)籌構(gòu)建糧食安全保障新體系研究［J］.社會(huì)科學(xué)輯刊，2019（6）：109-117.

［28］ 麥康森，徐皓，薛長(zhǎng)湖，等.開(kāi)拓我國(guó)深遠(yuǎn)海養(yǎng)殖新空間的戰(zhàn)略研究［J］.中國(guó)工程科學(xué)，2016，18（3）：90-95.

［29］ XU J，HAN L M，YIN W.Research on the ecologicalization efficiency of mariculture industry in China and its influencing factors［J］.Marine policy，2022，137［2022-03-25］.https：//doi.org/10.1016/j.marpol.2021.104935.

［30］ LIANG Y X，CHENG X W，ZHU H，et al.Historical evolution of mariculture in China during past 40 years and its impacts on eco-environment［J］.Chinese geographical science，2018，28（3）：363-373.

［31］ 宋媛，李國(guó)慶，吳晶，等.1989—2019年環(huán)渤海地區(qū)養(yǎng)殖池的空間變化特征［J］.海洋通報(bào)，2021，40（1）：92-100.

［32］ 許學(xué)工，梁澤，周鑫.黃河三角洲陸海統(tǒng)籌可持續(xù)發(fā)展探討［J］.資源科學(xué)，2020，42（3）：424-432.

［33］ 周鑫，許學(xué)工.黃河三角洲（東營(yíng)市）高效生態(tài)漁業(yè)綜合效益評(píng)估［J］.北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，51（3）：518-524.