海南省木蘭灣海岸沙地輸沙量分布特征研究

引言研究背景輸沙量是研究風(fēng)蝕地貌、沙漠變化以及防治沙漠化問題的關(guān)鍵指標(biāo)。目前對低海拔地區(qū)沙漠、戈壁和沙漠化地區(qū)的輸沙量進行了大量研究[1]。但對海岸沙地的輸沙量研究仍然不足。海南省木蘭灣作為中國的海岸線之一，具有月牙形的砂質(zhì)海灘，其獨特的地理位置和氣候特征對木蘭灣海岸沙地的風(fēng)沙活動具有主要影響，是典型的海岸沙地。然而，隨著人類活動的加劇和自然環(huán)境的變化，木蘭灣海岸沙地的輸沙量分布特征也發(fā)生了變化，這給海岸帶的生態(tài)環(huán)境和經(jīng)濟發(fā)展帶來了重要影響。目前，對于海南省木蘭灣海岸沙地的研究還處在初期階段，為了推進對海南省木蘭灣海岸沙地的相關(guān)研究，填補海南省木蘭灣海岸沙地輸沙量分布特征的空缺。本項目擬選取海南省木蘭灣作為典型研究區(qū)，明晰該地區(qū)海岸沙地輸沙量分布特征對闡明木蘭灣的風(fēng)沙活動規(guī)律，深入了解海岸線的演變規(guī)律和保護海岸生態(tài)環(huán)境具有重要意義，對海南省制定合理的海岸帶開發(fā)和海岸風(fēng)沙資源利用與防治具有重要的現(xiàn)實意義。國內(nèi)外研究進展國內(nèi)研究進展對于風(fēng)沙活動的研究，國內(nèi)最早于20世紀(jì)中葉，而且研究地區(qū)主要集中在沙漠半沙漠地區(qū)和江河流域地區(qū)，相比之下國內(nèi)對于海岸沙地的風(fēng)沙活動研究較少。其中，已有的研究表明，沙漠地區(qū)關(guān)于地表風(fēng)沙活動特征的研究多集中于風(fēng)況、輸沙勢、輸沙方向及輸沙量的研究[2]。如杜鶴強等學(xué)者通過野外實測塔克拉瑪干沙漠腹地新月形沙丘近地表100cm高度風(fēng)沙流的輸沙量，發(fā)現(xiàn)在近地表100cm高度內(nèi)，除了背風(fēng)坡外，新月形沙丘各部位風(fēng)沙流的輸沙量隨著高度的增加呈現(xiàn)出不同程度的減小[3]。張偉民等學(xué)者通過對莫高窟頂戈壁輸沙量的系統(tǒng)監(jiān)測應(yīng)用輸沙勢理論定量研究復(fù)雜風(fēng)向（速）下輸沙量的變化規(guī)律[4]。孫小雲(yún)等學(xué)者基于塔克拉瑪干沙漠22個氣象站2005—2007年的逐時風(fēng)記錄數(shù)據(jù)和輸沙勢定義計算了2007年各測站的16個方位的輸沙勢，并分析了輸沙勢的空間分布特征與季節(jié)變化[5]。張昊等人利用烏珠穆沁沙地的風(fēng)況數(shù)據(jù)，研究沙地起沙風(fēng)況和輸沙勢變化特征，同時結(jié)合沉積物粒徑數(shù)據(jù)，揭示該地區(qū)風(fēng)沙環(huán)境的空間分異規(guī)律[6]。汪海嬌，田麗慧等學(xué)者以青海湖東岸的克土沙區(qū)為研究區(qū)域，利用風(fēng)速風(fēng)向數(shù)據(jù)及全方位集沙儀的集沙數(shù)據(jù)，比較分析本區(qū)的風(fēng)沙活動特征[7]。落桑曲加，張焱等學(xué)者利用旋轉(zhuǎn)多路方口集沙儀對雅魯藏布江流域中游河岸沙地、河漫灘沙地和山麓沙地3種地表類型、6個點2020年9-12月的輸沙量進行野外實地觀測，對比不同下墊面的輸沙量，分析地表特征對輸沙量的影響過程[1]。國內(nèi)對海岸沙地的研究方向主要集中在海岸風(fēng)沙地貌類型、發(fā)育模式、風(fēng)沙運動、沉積特征、發(fā)育演變及海岸風(fēng)沙防治等，但近20年來在海岸沙地起沙風(fēng)速、風(fēng)沙流結(jié)構(gòu)、沙丘移動與形態(tài)變化等方面開展了一定的觀測研究，來探尋海岸沙地形成發(fā)育的動力學(xué)過程[8]。近年來，隨著對海岸帶環(huán)境保護意識的提高，越來越多的學(xué)者開始關(guān)注海岸沙地的動態(tài)變化和輸沙機制。如董玉祥，馬駿研究了輸沙量對海岸沙丘表面風(fēng)沙流中不同粒徑組沙粒輸沙量的垂向分布影響的原因，主要應(yīng)與不同風(fēng)速氣流的攜沙極限以及不同粒徑組沙粒的運動方式密切相關(guān)[9]。在對中國河北省昌黎海岸的沙丘分布區(qū)域進行研究時，李志星等學(xué)者采用了鄰近的三個氣象站所收集的2008至2018年間風(fēng)力數(shù)據(jù)，分析了沙土搬運勢能的時空變動，他們還通過衛(wèi)星圖像，并結(jié)合多種影響因素，對海岸沙丘的分布模式及其形態(tài)演變背后的驅(qū)動因素進行了綜合探討[10]。木蘭灣作為海南省的重要海岸帶之一，也逐漸引起了學(xué)者們的關(guān)注。吳多譽，張航飛等學(xué)者利用海灘斷面高層數(shù)據(jù)，對木蘭灣海岸海灘斷面地形進行了研究，木蘭灣海灘斷面地形具有季節(jié)性的變化特征，取得了一定的成果[11]。國外研究進展相對于國內(nèi)研究，國外學(xué)者對于風(fēng)沙活動的研究較為深入和系統(tǒng)。對于風(fēng)沙現(xiàn)象及其攜帶沙粒運動的研究，國際上已經(jīng)在多個地區(qū)進行，其中包括尼羅河的河谷和三角洲、科威特沙漠以及埃及的El-Khanka地區(qū)[2]。如MohamedE.Hereher從五個氣象站獲得的數(shù)據(jù)和遙感圖像中提取沙丘形態(tài)和沙丘推進速率，以此用來評估和量化尼羅河河谷和三角洲西側(cè)流沙和沙丘活動的嚴(yán)重程度[12]。J.M.Al-Awadhi等人基于八個氣象站的風(fēng)況數(shù)據(jù)資料和Fryberger的輸沙勢定義，分析了科威特風(fēng)沙活動的輸沙勢分布特征[13]。NazariSamani等人在研究中將風(fēng)況數(shù)據(jù)和地貌測量相結(jié)合，對伊朗中部KashanErg的不同沙丘進行分離和分類[14]。國外學(xué)者對海岸沙丘的研究中結(jié)合了多種因素，如NoamLevin等人使用衛(wèi)星圖像和雷達數(shù)據(jù)對海岸沙丘活動進行量化，并結(jié)合風(fēng)況數(shù)據(jù)和波浪數(shù)據(jù)分析輸沙勢和風(fēng)向性的變化[15]。JasperKnight等人利用遙感和野外實測數(shù)據(jù)，提取地上沙丘和通風(fēng)口的地貌數(shù)據(jù)和結(jié)合風(fēng)力狀況，分析了海岸沙丘的遷移模式和動態(tài)變化[16]。研究的主要內(nèi)容通過觀測海南島東北部木蘭灣2023年全年不同高度、不同方向輸沙量，不同方向輸沙量占比圖和全年實測輸沙率玫瑰圖進行分析，分析木蘭灣輸沙量、總輸沙量、總輸沙率隨隨高度和方向的變化，以探究海岸沙地形成演化過程和機理。研究區(qū)概況本研究區(qū)位于海南省文昌市北部鋪前鎮(zhèn)的木蘭灣，實驗監(jiān)測地的經(jīng)緯度位置為（20°36′19″N，110°44′29″E）。本研究區(qū)域?qū)儆跓釒Ъ撅L(fēng)性氣候，全年高溫，年平均氣溫23℃，旱雨兩季分明，降水量大，雨季集中在夏季，日照充足，年平均日照時數(shù)大于2000小時，風(fēng)向隨著季節(jié)的更替有明顯轉(zhuǎn)變。臺風(fēng)是其主要的自然災(zāi)害。本研究區(qū)屬于濱海堆積地貌，由7個具有獨立性的月牙形沙灘組成，陸域地表地形平緩，沿海岸線伸展，形成平緩的沙灘，沙灘略向大海傾斜，后方地形逐漸升高，局部形成沙丘，沙丘發(fā)育良好，地表多為沙生植被覆蓋。圖1木蘭灣實驗地監(jiān)測位置數(shù)據(jù)來源與處理數(shù)據(jù)來源本研究數(shù)據(jù)來自于木蘭灣實驗地積沙儀2023年的監(jiān)測數(shù)據(jù)，數(shù)據(jù)包含研究區(qū)2023年積沙儀10、30、50、100cm四個不同高度的輸沙量，N、NE、E、SE、S、SW、W、NW八個不同方向的輸沙量以及2023年各個方向的總輸沙量。數(shù)據(jù)分析方法在研究區(qū)選擇典型砂質(zhì)海岸作為觀測區(qū)，面積250m×250m研究區(qū)部署全方位沙子收集裝置（該裝置包括進沙通道、儲沙容器、防塵頂蓋和一個圓筒狀的外殼等組件）以捕獲360°范圍內(nèi)的所有沙子傳輸過程。計算輸沙率的公式為：Q=W其中Q為輸沙率，單位為kg/m·a； W是集沙量，單位為kg；L為積沙儀進口寬度，單位為m；ΔT為時間。本次研究使用的全方位定點積沙儀可以監(jiān)測10、30、50、100cm等4個高度的風(fēng)沙流特征，并且積沙儀按N、NE、E、SE、S、SW、W、NW的8個方向擺放，可以收集不同方位的所有輸沙，測得的數(shù)據(jù)可與風(fēng)向資料建立關(guān)系。在安裝積沙儀時，確保其下端口的底面與地面保持水平，且積沙口朝向與風(fēng)向相對。啟動積沙過程前，須同時開啟所有積沙儀的口蓋，并在積沙過程結(jié)束后，統(tǒng)一關(guān)閉各積沙儀的口蓋。取樣前，需對自封袋進行稱重，然后分層次采集地表10、30、50、100cm的沙塵樣本，并在烘干后，使用精度為1/10000克的電子天平進行稱重。使用Excel軟件整理和分析輸沙率等數(shù)據(jù)，并用origin2022制作輸沙率玫瑰圖進行分析。結(jié)果與分析木蘭灣2023年不同高度輸沙量分布特征根據(jù)實驗數(shù)據(jù)以及不同高度輸沙量折線圖（圖2）可知：從不同高度的總輸沙量來看，在積沙儀10、30、50、100cm四個不同的高度中，積沙儀10cm高度在八個方向的總輸沙量為424.38g；積沙儀30cm高度在八個方向的總輸沙量為203.06g；積沙儀50cm高度在八個方向的總輸沙量為136.86；積沙儀100cm高度在八個方向的總輸沙量最為51.37g。其中，研究區(qū)10cm高度的總輸沙量最大，100cm高度的總輸沙量最小，積沙儀不同高度總輸沙量隨著積沙高度的增加而減少，兩者呈負相關(guān)關(guān)系。以此得出，研究區(qū)近地表10cm高度區(qū)域的風(fēng)沙活動強度最強，100cm高度區(qū)域的風(fēng)沙活動強度最弱，而且研究區(qū)風(fēng)沙活動的強弱隨著高度的增加而逐漸減弱。從不同高度的輸沙量變化幅度來看，積沙儀10cm高度的輸沙量變化幅度最大，其中最小輸沙量為15.64g，最大輸沙量為152.16g，最小輸沙量與最大輸沙量差值為136.52g，輸沙量波動最大；積沙儀30cm高度的輸沙量變化幅度僅次于10cm的變化幅度，其中最小輸沙量為3.02g，最大輸沙量為53.87g，最小輸沙量與最大輸沙量差值為50.85g，其輸沙量波動程度僅次于10cm的輸沙量波動程度；積沙儀50cm高度的輸沙量變化幅度僅次于10cm和30cm的變化幅度，其中最小輸沙量為2.03g，最大輸沙量為39.6g，差值為37.57g，其輸沙量波動程度次于10cm和30cm的輸沙量波動程度；積沙儀100cm高度的輸沙量變化幅度最小，最小輸沙量為0.89g，最大輸沙量為16.57g，差值為15.68g，輸沙量波動最小。積沙儀不同高度的輸沙量變化幅度隨著高度的增加而減小，不同高度的輸沙量波動程度也隨高度的增加而減小，呈負相關(guān)關(guān)系。從不同高度的輸沙量分布方位來看,積沙儀10cm高度在N、NE、E、SE、S、SW、W、NW八個方位的輸沙量分別為29.68g、52.61g、62.53g、130.28g、152.16g、96.33g、37.31g、15.64g，其中，10cm高度的最大輸沙量分布在S方向，最小輸沙量分布在NW方向；積沙儀30cm高度在N、NE、E、SE、S、SW、W、NW八個方位的輸沙量分別為6.7g、28.41g、25.04g、37.94g、53.87g、43.84g、4.24g、3.02g，其中30cm高度的最大輸沙量分布在S方向，最小輸沙量分布在NW方向；積沙儀50cm高度在N、NE、E、SE、S、SW、W、NW八個方位的輸沙量分別5.7g、39.6g、14.46g、24.73g、25.91g、19.94g、4.49g、2.03g，其中50cm高度的最大輸沙量分布在NE方向，最小輸沙量分布在NW方向；積沙儀100cm高度在N、NE、E、SE、S、SW、W、NW八個方位的輸沙量分別為4.7g、12.73g、16.57g、8.89g、3.87g、1.13g、2.59g、0.89g，其中100cm高度的最大輸沙量分布在E方向，最小輸沙量分布在NW方向。不同高度的輸沙量，均是以偏南方向和偏東方向為主要分布方向，西北方向的輸沙量分布最少；除了NE和W這兩個方向，其余六個方向隨著積沙儀高度的增加，同一方向上的輸沙量逐漸減少。總體來說，除極個別情況外，木蘭灣2023年輸沙量隨高度的變化大致體現(xiàn)為：高度越高，輸沙量越??；木蘭灣2023年風(fēng)沙活動強度隨高度的變化大致體現(xiàn)為:高度越高，風(fēng)沙活動強度越弱。圖2不同高度輸沙量折線圖木蘭灣2023年不同方向輸沙量分布特征結(jié)合圖3從不同方向的總輸沙量來看，積沙儀在N方向的總輸沙量為46.78g，占不同方向總輸沙量比例為5%；積沙儀在NE方向的總輸沙量為133.35g，占不同方向總輸沙量比例為14%；積沙儀在E方向的總輸沙量為118.6g，占不同方向總輸沙量比例為12%；積沙儀在SE方向的總輸沙量為201.84g，占不同方向總輸沙量比例為21%；積沙儀在S方向的總輸沙量為235.81g，占不同方向總輸沙量比例為24%；積沙儀在SW方向的總輸沙量為161.24g，占不同方向總輸沙量比例為17%；積沙儀在W方向的總輸沙量為48.63g，占不同方向總輸沙量比例為5%；積沙儀在NW方向的輸沙量為21.58g，占不同方向總輸沙量比例為2%。研究區(qū)S方向的總輸沙量最多，所占比例最大，NW方向的總輸沙量最少，所占比例最??；研究區(qū)2023年輸沙量主要以偏南（S、SW）方向和偏東（SE、NE、E）方向的輸沙為主,這兩組方向占2023年輸沙方向的88%，其中偏南方向以S方向所占比例最高（24%），偏東方向以SE方向所占比例最高（21%）。由此得出，研究區(qū)風(fēng)沙活動主要集中在偏南和偏東方向，S方向的風(fēng)沙活動強度最強，NW方向的風(fēng)沙活動強度最弱。圖3不同方向總輸沙量占比圖結(jié)合圖4從不同方向輸沙量的分布高度來看，積沙儀N方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為29.68g、6.7g、5.7g、4.7g，分別占N方向總輸沙量比例為64%、14%、12%、10%；積沙儀NE方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為52.61g、28.41g、39.6g、12.73g，分別占NE方向總輸沙量比例為39%、21%、30%、10%；積沙儀E方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為62.53g、25.04g、14.44g、16.57g，分別占E方向總輸沙量比例為53%、21%、12%、13%；積沙儀S方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為152.16g、53.87g、25.91g、3.87g，分別占S方向總輸沙量的比例為64%、23%、11%、2%；積沙儀SW方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為96.33g、43.84g、19.94g、1.13g，分別占SW方向總輸沙量比例為60%、27%、12%、1%；積沙儀W方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為37.31g、4.42g、4.49g、2.59g，分別占W方向總輸沙量的比例為77%、9%、9%、5%；積沙儀NW方向在10cm、30cm、50cm、100cm四個高度的輸沙量分別為15.64g、3.02g、2.03g、0.89g，分別占NW方向總輸沙量比例為73%、14%、9%、4%。由以上分析及圖表資料可知，在八個方向的輸沙量占比中，10cm高度輸沙量所占比例最多，在近地面10cm高度為八個方向輸沙量的主要集中分布高度，是主要的風(fēng)沙活動高度，風(fēng)沙活動強度較強；除E方向外，其余7個方向的輸沙量占比中，100cm高度輸沙量所占比例最少，風(fēng)沙活動強度較弱。在不同方向的輸沙量占比中，輸沙量的垂直分布存在差異。其中，在各個方向中，垂直輸沙量的最大差異如下：N方向為24.98g，NE方向為39.94g，E方向為48.07g，SE方向為121.39g，S方向為148.29g，SW方向為95.2g，W方向為34.72g，NW方向為14.75g。由此可得，導(dǎo)致垂直輸沙量差異的原因可能是不同季節(jié)，風(fēng)向不同，土壤含水量不同。推測其研究區(qū)的主風(fēng)向主要以南風(fēng)和東南風(fēng)為主，在濕季，土壤含水量大，土壤不易被風(fēng)力侵蝕，不易起風(fēng)沙活動，故輸沙量的高度越低，垂直方向上的輸沙量差異較大；在干季，土壤含水量少，土壤易被風(fēng)力侵蝕，易起風(fēng)沙活動，故輸沙量的高度越高，垂直方向上的輸沙量差異較小。圖4各個方向輸沙量占比圖木蘭灣2023年實測輸沙率分布特征結(jié)合實驗數(shù)據(jù)和圖5來看，2023年研究區(qū)年總輸沙率為48.39kg/m·a，其中在N、NE、E、SE、S、SW、W、NW八個方向的輸沙率分別為2.34kg/m·a、6.67kg/m·a、5.93kg/m·a、10.09kg/m·a、11.79kg/m·a、8.06kg/m·a、2.43kg/m·a、1.07kg/m·a，輸沙率介于0~7.5之間的方向有N、EN、E、W、NW，介于7.5~15之間的有S、SE、SW，輸沙率呈現(xiàn)出明顯的方向變化特征，其中S方向的輸沙率最大，NW方向的輸沙率最小。本研究區(qū)輸沙率以偏南方向（S、SW)和偏東方向(SE、E、NE）為主，其中偏南方向以S方向的輸沙率最大（11.79kg/m·a），偏東方向以SE方向的輸沙率最大（10.09kg/m·a）。這與前面分析的不同方向總輸沙量的分布特征基本一致。圖52023年木蘭灣輸沙率玫瑰圖結(jié)論木蘭灣海岸沙地輸沙量隨高度的變化顯著，輸沙量的分布在10cm高度最大，100cm高度最??；輸沙活動主要集中在近地表10cm高度，并且隨著高度的增加輸沙量遞減，研究區(qū)風(fēng)沙活動強度隨著高度的增加而減弱。木蘭灣海岸沙地輸沙量隨方向的變化顯著，輸沙量的分布在S方向最大，NW方向最小，研究區(qū)風(fēng)沙活動主要集中在偏南和偏東方向。不同方向的輸沙量垂直分布存在差異，原因是不同季節(jié)，風(fēng)向不同，土壤含水量不同，導(dǎo)致輸沙垂直分布的差異。木蘭灣海岸沙地2023年總輸沙率為48.39kg/m·a，研究區(qū)南部與東部輸沙率較大，其中最大輸沙率分布在S方向，北部與西部輸沙率較小，其中最小輸沙率分布在NW方向。2023年研究區(qū)總輸沙量在不同方向上的分布特征與總輸沙率分布特征基本一致。論文的不足與展望本文基于木蘭灣2023年全年不同高度、不同方向?qū)崪y輸沙量，系統(tǒng)地分析了輸沙量、總輸沙量、總輸沙率隨高度和方向的變化。通過對海南省木蘭灣海岸沙地輸沙量的分布特征進行研究，揭示了木蘭灣海岸沙地輸沙量分布規(guī)律，為后續(xù)推測木蘭灣海岸沙地演變機理的研究奠定了基礎(chǔ)。盡管對木蘭灣海岸沙地輸沙量分布特征的研究取得了一定的進展，但是由于作者本科階段學(xué)習(xí)的知識比較籠統(tǒng)，學(xué)識積累和自身水平有限，仍有許多問題有待進一步的探討和解決。研究過程中使用的方法比較簡單，目前的實驗結(jié)果分析可能還不夠透徹，需要學(xué)習(xí)更為嚴(yán)謹?shù)姆治龇椒?，更加深入、全面的分析實驗結(jié)果。海南省木蘭灣海岸屬于典型的熱帶海岸，其海岸沙地的形成是一個復(fù)雜的過程。研究過程中忽略了輸沙量的其他影響因素，如風(fēng)速的變化、沙源的分布情況和植被覆蓋率，需要更加深入地研究這些因素對輸沙量的影響機制和相互作用關(guān)系，這有助于我們更好地理解木蘭灣海岸沙地的動態(tài)變化過程，并為相關(guān)研究和應(yīng)用提供更加準(zhǔn)確地理論基礎(chǔ)。未來對海南省木蘭灣的研究可以從以下幾個方面展開：一是加強實地觀測的和采樣工作，深入了解木蘭灣海岸沙地的地形地貌、氣候氣象和水動力條件等特征；二是需要進一步探索新的研究方法的技術(shù)手段。如利用先進的遙感技術(shù)和GIS技術(shù)，建立木蘭灣海岸沙地的時空數(shù)據(jù)庫和數(shù)字高程模型（DEM），為研究輸沙量分布特征提供更加準(zhǔn)確和全面的數(shù)據(jù)支持；三是開展數(shù)值模擬研究，模擬木蘭灣海岸沙地的動力地貌過程，預(yù)測未來輸沙量分布趨勢。同時，還需要加強跨學(xué)科合作和交流，推動海南省木蘭灣相關(guān)領(lǐng)域的研究進展，同時加強研究成果的應(yīng)用和推廣。參考文獻落桑曲加，張焱，馬鵬飛等.雅魯藏布江中游不同地表輸沙量特征[J].中國沙漠，2022，42(02):6-13.羅鳳敏，高君亮，辛智鳴等.烏蘭布和沙漠東北緣起沙風(fēng)風(fēng)況及輸沙特征[J].農(nóng)業(yè)工程學(xué)報，2019，35(04):145-152.杜鶴強,韓致文，王濤，等.新月形沙丘表面風(fēng)速廓線與風(fēng)沙流結(jié)構(gòu)變異研究[J].中國沙漠,2012,32(01):9-16.張偉民，王濤，汪萬福等.復(fù)雜風(fēng)況條件下戈壁輸沙量變化規(guī)律的研究[J].中國沙漠，2011，31(03):543-549.孫小雲(yún)，房彥杰，趙景峰等.塔克拉瑪干沙漠輸沙勢時空分布特征[J].干旱區(qū)地理，2020，43(01):38-47.張昊，黨曉宏，蒙仲舉等.烏珠穆沁沙地主要風(fēng)沙環(huán)境特征及形成機制研究[J].干旱區(qū)研究,2023,40(10):1687-1697.DOI:10.13866/j.azr.2023.10.15.汪海嬌，田麗慧，張登山等.青海湖東岸沙地風(fēng)沙活動特征[J].中國沙漠，2020，40(01):49-56.董玉祥，李志忠.近40年中國海岸風(fēng)沙地貌研究回顧[J].中國沙漠,2022,42(01):12-22.董玉祥，馬駿.輸沙量對海岸沙丘表面風(fēng)沙流中不同粒徑沙粒垂向分布的影響[J].中山大學(xué)學(xué)報(自然科學(xué)版)，2009，48(03):102-108.李志星，