基于“3S”技術(shù)的南流江采砂點時空分布特征

朱文軒，黎樹式,2，馮炳斌，楊夏玲

（1. 南寧師范大學(xué)地理科學(xué)與規(guī)劃學(xué)院，南寧 5300011；2.廣西北部灣海洋環(huán)境變化與災(zāi)害研究重點實驗室，廣西欽州 535011；）

世界城鎮(zhèn)化處在快速發(fā)展的過程中，到2050年，預(yù)計68%的世界人口將居住在城市（Hackney et al,2020），未來90%的城鎮(zhèn)化將發(fā)生在亞非等發(fā)展中國家。而砂石是城市建設(shè)必不可少的原材料，在過去的20年間全球范圍內(nèi)砂石的使用量幾乎翻了2 倍，全球砂石骨料年產(chǎn)量和消費量高達400 億t，其中，中國為200 億t（王瓊杰，2019），而砂石的自然增長速度遠落后于消耗量。不可持續(xù)的采砂活動會造成一系列的環(huán)境危害，例如：在河床形成散亂的采砂坑，改變河道形態(tài)，造成河岸的不穩(wěn)定（高耶等，2017；Koehnken et al.,2020）；還會影響河水的理化性質(zhì)，導(dǎo)致水質(zhì)下降（楊惠雯等，2020；袁婷等，2020）；此外，也可能導(dǎo)致河流水位降低，降低洪水發(fā)生的頻率（Park et al.,2020），到達河口地區(qū)的泥沙量也會因此減少，影響河口三角洲沖淤的動態(tài)平衡（黎兵，2020）。南流江流域在改革開放以來的近幾十年間河砂的開采量也急劇上升，給流域生態(tài)環(huán)境的穩(wěn)定帶來巨大的挑戰(zhàn)。為應(yīng)對采砂活動帶來的環(huán)境挑戰(zhàn)，需對流域采砂活動活動的過程變化、動力因素和環(huán)境影響等特征進行充分認識，以便為采砂活動的科學(xué)有效管理、流域生態(tài)治理提供科學(xué)依據(jù)。

遙感影像數(shù)據(jù)相比較傳統(tǒng)的觀測數(shù)據(jù)具有時間跨度大、空間覆蓋范圍廣、獲取便利等優(yōu)點，在采砂活動的分析方面已有較多的應(yīng)用。在湖泊采砂活動分析方面，較多學(xué)者利用Landsat 系列影像通過識別采砂船分析湖泊采砂活動時空變化（Wu et al.,2009;Li et al.,2014），還有學(xué)者利用夜間燈光遙感數(shù)據(jù)依據(jù)光源分析夜間湖泊的采砂活動時空變化（Duan et al.,2019）和運用MODIS等遙感影像數(shù)據(jù)對水體懸沙濃度進行反演以分析采砂活動時空變化（Barnes et al.,2015）。河流采砂活動分析方面，有學(xué)者利用高分辨的影像識別河流中的采砂船來估算河流的采砂量（李志威等，2019），還有學(xué)者通過圖像紋理特征結(jié)合實地驗證識別采砂場來分析某一時間河流的采砂現(xiàn)狀（陳冬勤等，2013）。通過對已有研究的回顧，發(fā)現(xiàn)由于湖泊地區(qū)水域開闊，采用Landsat 等中高分辨率的影像，結(jié)合“識別采砂船”和“采砂區(qū)域光源”的方式進行采砂活動分析且具有較好的成效；而河流水域較狹長，采砂的規(guī)模相對較小，采砂活動的部分解譯標(biāo)志，如“采砂船”“水色”“光源”在Landsat 等中高分辨率影像中不顯著，無法滿足識別需求。因此，以往河流采砂活動的研究多采用高分辨率的影像開展。但受限于高分辨率影像的歷史跨度較短、早期影像稀缺等因素，還未見有河流采砂活動長期時空變化研究。通過對Landsat 影像中河流采砂活動的影像特征進一步分析發(fā)現(xiàn)，采砂區(qū)域的部分關(guān)鍵指示特征，如河床紋理特征和砂石的光譜特征，及部分水色特征，在影像中具有較好的表現(xiàn)，可滿足對采砂活動識別的需求。

鑒于此，本文擬采用Landsat 系列影像，以南流江為研究對象，分析河流采砂活動的時空變化特征，在此基礎(chǔ)上，進一步探討南流江采砂活動時空變化的驅(qū)動因素及產(chǎn)生的環(huán)境效應(yīng)，以期為南流江的生態(tài)保護和綜合治理提供科學(xué)依據(jù)。本案例研究有利于拓展Landsat 系列影像的應(yīng)用場景，為河流采砂活動的長期時空變化研究提供新的方法。

1 研究區(qū)域與數(shù)據(jù)

1.1 研究區(qū)域概況

南 流 江（21°30′—23°00′N，109°00′—110°30′E）位于廣西壯族自治區(qū)南部，是中國西南地區(qū)最大的獨流入海河流，發(fā)源于玉林市大容山南側(cè)，流經(jīng)北流、玉林等6縣（市）區(qū)，于合浦縣注入北部灣。南流江河長287 km，流域面積8 635 km2；上游河寬30～70 m，河流落差較大；中游河寬80～100 m，河流落差較小；下游河寬250～400 m，河面寬廣，流速平緩。流域的劃分，上游中游大致以“龍門江-張黃江-石灣鎮(zhèn)與泉水鎮(zhèn)的分界線”一線為界，中游與下游大致以“龍麗江部分河道-博白縣城所在地-三步江部分河道”一線為界。流域地形地貌復(fù)雜，上中游以平原丘陵地形為主，下游地區(qū)以平原為主，河口地區(qū)形成面積約為500 km2的三角洲（王丹媛，2019）。流域2020 年人口總數(shù)為1 089.78萬人，流經(jīng)北部灣經(jīng)濟區(qū)（圖1），區(qū)域內(nèi)經(jīng)濟增長速度顯著加快，發(fā)展前景廣闊，地理位置十分重要。

圖1 研究區(qū)域Fig.1 Location of study area

1.2 基礎(chǔ)數(shù)據(jù)與輔助數(shù)據(jù)

基礎(chǔ)數(shù)據(jù)為13幅Landsat5 TM數(shù)據(jù)和5幅Landsat8 OLI_TIR 數(shù)據(jù)，來源于美國地質(zhì)勘探局官網(wǎng)①https://earthexplorer.usgs.gov/bulk，影像數(shù)據(jù)信息見表1所示。本研究還用到以下輔助數(shù)據(jù)：

表1 研究區(qū)詳細遙感數(shù)據(jù)Table 1 Detailed remote sensing data of the study area

1）全國1∶100 萬矢量數(shù)據(jù)，來源為全國地理信息資源目錄服務(wù)系統(tǒng)②https://www.webmap.cn/。2）南流江下游常樂站的年平均流量和年平均含沙量數(shù)據(jù)（1988-2018 年）來源于《中國水文年鑒：珠江流域水文資料》（水利部珠江水利委員會，1988—2018）；1990-2002年水質(zhì)數(shù)據(jù)來源于《廣西年鑒》（廣西壯族自治區(qū)地方志編纂委員會，1990-2002），2003-2018 年水質(zhì)數(shù)據(jù)來源于《廣西壯族自治區(qū)環(huán)境統(tǒng)計公報》（廣西壯族自治區(qū)環(huán)境保護廳，2003—2018）；潮位數(shù)據(jù)來源于相應(yīng)年份的《潮汐表》（國家海洋信息中 心， 1997、 1988、 2007、2019）。3）玉林市、欽州市、北海市1989-2018 年的人均GDP數(shù)據(jù)，來源于《廣西年鑒》（廣西壯族自治區(qū)地方志編纂委員會，1989—2018）。

2 研究方法

2.1 通過目視解譯識別采砂點

研究區(qū)域的遙感影像經(jīng)過預(yù)處理后，進行采砂點的目視解譯。通過影像中疑似采砂點與現(xiàn)場勘驗結(jié)果對比，確定采砂點的圖像紋理與光譜特征，進而確定采砂點的判斷依據(jù)：河床形態(tài)、河床植被覆蓋情況、水色情況；與前一時相影像對比有無河床缺失。為便于量化統(tǒng)計，采砂區(qū)域岸線長度每400 m算作1處，不足400 m也算作1處。如圖2-a所示，Landsat8真彩色影像中采砂點處無植被覆蓋，呈淺黃色，且岸線曲折，水體偏渾濁。在圖2-b中Landsa8影像的562假彩色合成影像中采砂點呈淺綠色；在圖2-c 的谷歌影像中，清晰地顯示多艘采砂床正進行采砂作業(yè)，采砂點處堆積著成片的河沙；圖2-d為2015年的實地驗證圖，圖中采砂船正進行采砂作業(yè)。通過圖2-c、2-d 對圖2-a、2-b 的驗證，可基本確定采砂點的影像特征。之后利用ArcGIS10.7 對采砂點進行矢量化，得到1988、1998、2008、2018年的采砂點數(shù)據(jù)。

圖2 采砂點示例Fig.2 Example of sand mining points

2.2 南流江流域建筑用地提取

依據(jù)南流江流域的土地利用特點，構(gòu)建土地利用分類體系。體系包含6 個土地利用類型：耕地、林草地、水體、建筑用地、人工濕地和其他用地。之后對南流江流域遙感影像進行監(jiān)督分類，分類算法采用隨機森林法，該算法在數(shù)據(jù)分類方面有較多優(yōu)點，諸如人工干涉少、運算周期短、數(shù)據(jù)刻畫全面多樣等（趙宇虹，2021）。最后使用ENVI5.3 拓展工具中的Random Forest Classification工具進行分類，提取各年份建筑用地，并分析面積變化。

2.3 河床面積缺失計算及河口三角洲岸線變化分析

集聚的采砂點往往構(gòu)成1個大的采砂區(qū)域，多位于河流凸岸。利用ArcGIS10.7將各年份影像的采砂區(qū)域矢量化，之后通過屬性計算統(tǒng)計各年份采砂區(qū)域的面積變化，進而得出由于采砂活動所導(dǎo)致的采砂區(qū)域面積缺失量。同樣利用該軟件將各年份（根據(jù)衛(wèi)星過境時三角洲地區(qū)的潮位信息，經(jīng)過仔細甄別，選取潮位相近的19891202、19981109、20081220、20201207四期影像）的河口三角洲岸線矢量化，岸線的判斷采用平均高潮線作為海陸分界指示線，分析各年份河口三角洲的岸線時空變化（圖3）。

2.4 相關(guān)分析法

社會經(jīng)濟發(fā)展是采砂活動開展的原動力，經(jīng)濟發(fā)展帶來的城市基礎(chǔ)建設(shè)對砂石產(chǎn)生直接需求，是采砂活動直接動力。因而分析經(jīng)濟發(fā)展和城市面積變化與采砂活動的過程聯(lián)系，對于砂石供需矛盾的調(diào)節(jié)、采砂活動的規(guī)劃管理有重要意義。同時，完善合理的監(jiān)管是采砂活動持續(xù)健康發(fā)展的必要保證，分析監(jiān)管措施對采砂活動的影響，有助于評估監(jiān)管成效，發(fā)展監(jiān)管問題，從而提升監(jiān)管效率。因此，分別從流域人均GDP 變化、建筑用地面積變化、采砂監(jiān)管3個方面分析其與南流江采砂點變化的關(guān)系（見圖3），探討采砂點時空變化的動力因素。同時關(guān)注采砂活動的環(huán)境影響，分別從水質(zhì)變化、水沙特征、河床面積變化、三角洲岸線變化4個方面分析采砂活動所產(chǎn)生的環(huán)境效應(yīng)（見圖3）。

圖3 研究路線Fig.3 Research roadmap

3 結(jié)果與分析

3.1 采砂點的時空分布特征

根據(jù)采砂點的判斷規(guī)則，對遙感影像進行采砂點提取，得到1988-2018 年南流江采砂點分布。圖4、5 顯示，南流江全河段的采砂點在1998-2008 年數(shù)量快速上升，2018年采砂點數(shù)量大幅下降，上中下游采砂點數(shù)量呈相同的變化趨勢，均在2008 年達到數(shù)量的最高點，分別為9、55、89個。各縣采砂點數(shù)量的時間變化趨勢與各流域分區(qū)一致（見圖4），合浦縣、博白縣、浦北縣、福綿區(qū)均呈先增加后減少的趨勢，同樣在2008 年采砂點的數(shù)量達到最多，分別為81、46、21、4個。各個時期采砂點空間的分布情況為：下游＞中游＞上游，合浦縣＞博白縣＞浦北縣＞福綿區(qū)（合浦縣大致位于下游，浦北縣和博白縣大致位于中游，福綿區(qū)大致位于上游），呈現(xiàn)此分布特征的主要原因是，下游地區(qū)地勢逐漸平緩，流速下降，河床中的泥沙沉積量逐漸增多，可供開采的砂石資源豐富。此外，下游地區(qū)交通條件便利（圖6），便于河沙外運，這也進一步刺激砂石的開采。

圖4 南流江采砂點統(tǒng)計（a.分區(qū)統(tǒng)計；b.縣統(tǒng)計）Fig.4 Statistics of sand mining points in Nanliu River(a.Regional statistics;b.County statistics)

圖5 1988—2018年南流江簡化采砂點分布Fig.5 Distribution of simplified sand mining points in Nanliu River from 1988 to 2018

圖6 2018年距南流江10 km內(nèi)主要公路密度與采砂點分布Fig.6 Density of main roads and distribution of sand mining points within 10 km away from Nanliu River in 2018

3.2 采砂典型區(qū)域

圖7為選取的下游地區(qū)的采砂典型區(qū)域，位于常樂鎮(zhèn)和泉水鎮(zhèn)的交界處，典型區(qū)域內(nèi)河道較為曲折，水流流速緩慢，河床泥沙儲量豐富?？擅黠@地看出，由于采砂活動，河道形態(tài)發(fā)生較大形變，河床逐漸破碎呈鋸齒狀，河道明顯拓寬，河流生態(tài)環(huán)境受到較為嚴重的破壞；并且相較之前的影像，2018年河道的沙洲部分基本消失，可推測，由于常年過度開采，部分河道砂石儲量日漸枯竭；在受破壞程度方面，下游地區(qū)的受破壞程度要大于中游地區(qū)和上游地區(qū)。采砂點的時空分布方面，所示區(qū)域的采砂點數(shù)量同樣在1988-2008 年呈上升狀態(tài)，2018 與2008 年數(shù)量基本持平，與下游分區(qū)采砂點的數(shù)量變化趨勢相一致。

圖7 下游采砂典型區(qū)域（a.1988-1998年；b.2008-2018年）Fig.7 Typical area of downstream sand mining(a.1988-1998;b.2008-2018)

4 討論

4.1 南流江采砂點時空分布的驅(qū)動力

4.1.1 社會經(jīng)濟發(fā)展 圖8 顯示，南流江流域3 市的人均地區(qū)生產(chǎn)總值的變化總體呈不斷上升的趨勢，由開始的緩慢上升到后來的急速上升。欽州、北海、玉林市3 市2018 年的人均GDP 相比1989 年分別增長了56、82 和35 倍，年均增長率分別為192%、283%和122%。經(jīng)濟的持續(xù)增長對基礎(chǔ)設(shè)施建設(shè)產(chǎn)生旺盛的需求，砂石作為基礎(chǔ)設(shè)施建設(shè)的主要原材料的需求量也隨之激增。統(tǒng)計期間，南流江經(jīng)濟水平與采砂點的數(shù)量整體上均呈快速上升趨勢。因此，可認為南流江流域經(jīng)濟持續(xù)增長是推動采砂點增加的主要原因。

圖8 南流江流域市級1989-2018年人均地區(qū)生產(chǎn)總值變化Fig.8 Change of per capita regional GDP of cities in Nanliu River Basin from 1989 to 2018

4.1.2 建筑用地面積變化 對南流江流域遙感影像數(shù)據(jù)進行監(jiān)督分類，得到南流江流域1988-2018年土地利用狀況和建筑用地變化（圖9），南流江1988、1998、2008、2018年的建筑用地面積分別為159.359、232.721、426.000、526.122 km2，年均增長率為12.2%。沙石是城市建設(shè)的主要原料，城市基礎(chǔ)建設(shè)直接促使南流江的河沙開采量的快速增長。在研究期間采砂點數(shù)量與建筑用地面積的變化呈階段性一致，1988-2008 年呈正相關(guān)，尤其1998-2008 年建筑用地面積增長量達到最大，為193.279 km2。在此期間增加的采砂點數(shù)量也最多，達到57個，可認為建筑行業(yè)的需求直接刺激了采砂點數(shù)量的不斷增加。

圖9 1988-2018年南流江流域土地利用類型Fig.9 Land use types of Nanliu River basin from 1988 to 2018

4.1.3 采砂行業(yè)的監(jiān)管 南流江采砂點的數(shù)量于2018年出現(xiàn)急劇下降，考慮在該期間流域內(nèi)經(jīng)濟增長迅速，砂石的需求旺盛，推測可能是由于采砂監(jiān)管加強所致。經(jīng)對相關(guān)資料查證，在此期間，有關(guān)部門分別于2011 和2017 年頒布實施了《廣西壯族自治區(qū)河道采砂管理辦法》（廣西壯族自治區(qū)人民政府，2011）和《廣西壯族自治區(qū)采砂管理條例》（廣西壯族自治區(qū)人民代表大會常務(wù)委員會，2016），監(jiān)管力度不斷加強。2016 年出臺《關(guān)于辦理非法采礦、破壞性采礦刑事案件適用法律若干問題的解釋》（最高人民法院最高人民檢察院，2016）規(guī)定非法采砂行為全面入刑，對非法采砂活動起到強烈的震懾作用。這一系列采砂法律法規(guī)的完善，顯著地加強和完善了對采砂活動的監(jiān)管。

從執(zhí)法打擊非法采砂的成效看，2016年開始中央環(huán)保督察組進駐廣西壯族自治區(qū)，被稱為“史上最嚴環(huán)保督察”，反饋了包括南流江采砂問題在內(nèi)的眾多環(huán)保意見。各級政府對反饋的各項采砂問題進行了嚴格整改。2016年，玉林市各級相關(guān)部門共查處九洲江、南流江流域環(huán)境違法行為242起，搗毀非法采砂船353艘，北海市于11月成立采砂秩序治理專項小組，對14家違規(guī)采石場依法立案查處；2017 年，欽州市對南流江40 余個非法采砂點，組織開展了4 次打擊行動，清理取締了32 個；2018年，玉林市打擊取締違法采砂場67個，違法洗沙場28 個（玉林市地方志編纂委員會，2019），欽州市取締南流江非法采砂場184個（欽州市地方志編纂委員會，2019）。在相關(guān)部門的多次高強度打壓下，效果立竿見影，南流江的非法采砂活動得到有效控制，采砂點的數(shù)量有明顯下降。

4.2 南流江采砂活動的環(huán)境效應(yīng)分析

4.2.1 采砂活動對南流江水沙特征的影響 在

1988-2018年南流江常樂站年平均流量呈弱增加趨勢，而年平均含沙量呈明顯減少趨勢，2007年平均含沙量達到最低（圖10-a），年平均流量與年平均含沙量存在明顯的負相關(guān)性，由此推測年平均含沙量下降可能是由于人為因素所致；在1988-2018年采砂點整體的變化呈快速上升態(tài)勢，而年平均含沙量呈明顯下降趨勢（圖10-b），年平均含沙量與采砂點數(shù)量同樣存在負相關(guān)性，據(jù)此推斷，強烈的采砂等人類活動可能是含沙量下降的重要原因之一。過度采砂致使河底可供侵蝕的泥沙量大幅減少，導(dǎo)致在年平均流量變化不大的情況下，河流的挾沙量出現(xiàn)明顯下降。

圖10 1988—2018年南流江常樂站年平均流量與年平均含沙量變化（a）和年平均含沙量與采砂點變化（b）Fig.10 Change of annual average flow and annual average sediment concentration of Changle Station on Nanliu River(a)and change of annual average sediment concentration and sand extraction point(b)during 1988-2018

4.2.2 南流江水質(zhì)變化 中國于1988、1999、2002 年分別頒布了GB3838-1988、GHZB1-1999和GB3838-2002（現(xiàn)行）3版《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》，不同的標(biāo)準(zhǔn)所檢測的項目以及各水質(zhì)類別的標(biāo)準(zhǔn)限值發(fā)生較大改變。為了數(shù)據(jù)的可比性，將已有水質(zhì)數(shù)據(jù)分2個部分進行統(tǒng)計，對于1988-2002年的水質(zhì)變化，從超標(biāo)物質(zhì)變化的角度進行分析；對于2003-2018年的水質(zhì)變化，從水質(zhì)類別變化的角度進行分析。由表2 可知，1992、1998、1999、2002年都出現(xiàn)懸浮物質(zhì)超標(biāo)，2000、2001年出現(xiàn)石油類物質(zhì)超標(biāo)。懸浮物質(zhì)超標(biāo)可能與采砂活動對水體的擾動相關(guān)，石油物質(zhì)超標(biāo)可能與采砂船和運砂船的增加相關(guān)。在1988-2002年懸浮物與石油類物質(zhì)超標(biāo)頻次明顯增加，同時采砂點的數(shù)量處于不斷上升的趨勢。由此推測，采砂活動可能對上述物質(zhì)超標(biāo)作出較大貢獻。表3顯示，2003-2014年南流江的水質(zhì)整體上優(yōu)于Ⅲ類標(biāo)準(zhǔn)，但在2015-2018年水質(zhì)下降為Ⅳ類水質(zhì)。采砂變化方面，2008-2018年由于采砂影響，采砂區(qū)域面積減少了3.55 km2，高于之前2個10年的采砂區(qū)域面積缺失量。且考慮到對采砂活動的監(jiān)管力度整體加強主要在2016年之后，由此推斷2008-2015年采砂強度可能仍處在不斷加劇的狀態(tài)，強烈的采砂活動可能是導(dǎo)致該期間水質(zhì)下降的一個重要的原因。

表2 南流江1990—2002年水質(zhì)超標(biāo)物質(zhì)情況Table 2 Substances exceeding the standard in water quality of Nanliu River from 1990 to 2002

表3 南流江2003—2018年水質(zhì)類別情況Table 3 Water quality categories of Nanliu River from 2003 to 2018

4.2.3 河岸河床侵蝕 為分析采砂活動對河床的破壞情況，利用南流江流域的遙感影像，對采砂區(qū)域進行矢量化，構(gòu)建矢量面（圖11），之后對采砂區(qū)域矢量面進行面積計算匯總。結(jié)果顯示，1988、1998、2008 和2018 年采砂區(qū)域面積分別為44.64、42.45、40.95 和37.40 km2。從圖11 可以明顯看到，由于采砂活動導(dǎo)致河床岸線呈鋸齒狀，部分區(qū)域有明顯殘缺，采砂區(qū)域面積由于采砂每10年分別減少2.19、1.50和3.55 km2，年均減少約為0.24 km2，可見，過度的采砂活動對河床破壞十分嚴重。

圖11 南流江采砂區(qū)域Fig.11 Sand mining area of Nanliu River

4.2.4 河口三角洲岸線變化 為分析采砂點數(shù)量的變化對南流江河口三角洲的影響，對南流江河口三角洲岸線進行矢量化（圖12）。結(jié)果顯示，近30年來南流江三角洲的岸線變化還是以向海延伸為主，淤積主要發(fā)生在淤泥質(zhì)海岸，主要是由于紅樹林快速生長導(dǎo)致的岸線向海擴張。在海岸線的蝕退方面，雖因采砂活動等人為因素引起輸沙量的減少，但目前并沒有發(fā)現(xiàn)三角洲地區(qū)岸線受到明顯侵蝕，主要有以下原因：1）河口三角洲地區(qū)的生態(tài)環(huán)境不斷轉(zhuǎn)好，紅樹林面積不斷擴大，對海岸線起有效防護作用；2）人工防浪堤、防潮堤的修建使得海岸線不易被侵蝕；3）南流江入?？趯儆趶姵焙０?，強潮為河口帶來大量泥沙（黎樹式，2017）。

圖12 河口三角洲岸線變化Fig.12 Shoreline change of estuarine delta

5 結(jié)論與建議

5.1 結(jié)論

利用1988、1998、2008、2018年共4個時相的遙感影像分析南流江采砂點的時空分布，探討其背后驅(qū)動因素及所產(chǎn)生的環(huán)境效應(yīng)，得出主要的結(jié)論為：

1）30 年來南流江采砂點時空變化明顯。采砂點在1988-2008 年呈不斷增長態(tài)勢，2008 年達到頂點；之后，2018年采砂點的數(shù)量快速下降，各流域分區(qū)以及各縣域采砂點的變化類似。各個時期均為下游地區(qū)采砂點最多，中游次之，上游地區(qū)最少，在縣域的分布為：合浦縣＞博白縣＞浦北縣＞福綿區(qū)，這與河流泥沙的沉積規(guī)律一致。

2）采砂點時空變化的驅(qū)動因素是多樣的，經(jīng)濟發(fā)展帶來的城市化快速推進，是造成采砂點數(shù)量增長的主要原因。2018 年的采砂點數(shù)量出現(xiàn)下降，主要由于是2015年以來南流江采砂監(jiān)管力度不斷加強所致。

3）南流江采砂活動對周圍環(huán)境造成較大影響，首先，采砂活動造成采砂區(qū)域面積不斷減小，進而床載泥沙量不斷下降，導(dǎo)致河流的挾沙量不斷下降；其次，采砂活動是導(dǎo)致南流江水質(zhì)下降的重要原因；此外，雖然因采砂活動等人為因素引起輸沙量的減少，但目前并沒有發(fā)現(xiàn)三角洲地區(qū)岸線受到明顯的侵蝕。

5.2 南流江綜合治理建議

采砂活動發(fā)生區(qū)域和所造成的影響波及到整個流域，各地區(qū)在采砂監(jiān)管和生態(tài)治理方面應(yīng)繼續(xù)加強合作，提升綜合治理能力。在生態(tài)恢復(fù)方面，應(yīng)繼續(xù)加強對生活工業(yè)養(yǎng)殖廢水的排放監(jiān)管，規(guī)范采砂作業(yè)流程，提升水環(huán)境質(zhì)量；在生態(tài)保護方面，南流江入海的泥沙量減少將會是一個長期的趨勢，因此，需要繼續(xù)增加紅樹林的面積，增強海岸的防護能力。在砂石的可持續(xù)開采方面，一方面有關(guān)部門或許可以通過支持南流江泥沙輸移沉積狀況研究，確定泥沙的開采區(qū)域或者開采量，以便實現(xiàn)自然供應(yīng)與開采量之間的平衡；另一方面，應(yīng)積極尋找替代性的資源，如大力發(fā)展和推廣應(yīng)用機制砂（國家發(fā)展和改革委員會，2020），通過多種途徑滿足經(jīng)濟發(fā)展所需的砂石量。最后，南流江的保護和治理關(guān)系到流域內(nèi)全體民眾的福祉，應(yīng)充分發(fā)動群眾，對非法采砂等破壞流域生態(tài)環(huán)境的活動進行監(jiān)督，保障民意反饋通道的通暢高效，形成重視生態(tài)、共同治理的良好社會風(fēng)氣。