摘要：【目的】明確云南省山區(qū)旱地稻田蠐蠟種類組成、種群空間分布型及抽樣方法，為云南山區(qū)旱地稻田蠐蠟類害蟲的預(yù)測預(yù)報及綜合防控提供參考依據(jù)?！痉椒ā窟x擇云南省普洱市瀾滄縣竹塘鄉(xiāng)蒿枝壩村、哈果馬村和諾國村3個主要旱地稻種植區(qū)為試驗地點，于2022和2023年采用五點取樣法對旱地稻田中的蠐嘈種類進(jìn)行系統(tǒng)調(diào)查及鑒定；通過聚集度指標(biāo)法、Iwao回歸分析、Taylor冪函數(shù)法及種群聚集均數(shù)分析不同試驗點旱地稻田蠐蠟空間分布型和聚集原因，并基于空間分布結(jié)果，擬合建立理論抽樣和序貫抽樣模型。【結(jié)果】云南省普洱市瀾滄縣3個試驗地點的旱地稻田共發(fā)現(xiàn)9種蠐蠟，隸屬于3科6屬，其中黑阿鰓金龜（Apogonia cupreoviridis）和東方絹金龜（Maladera orientalis）是主要的優(yōu)勢種，相對豐度分別為30.91%和24.18%。聚集度指標(biāo)法、Iwao回歸分析及Taylor冪法結(jié)果顯示，2022和2023年蒿枝壩村、哈國馬村和諾國村的蠐蠟空間分布型均表現(xiàn)為聚集分布（alt;0，βgt;1）;聚集均數(shù)λ值均大于2，表明昆蟲本身的生活習(xí)性和環(huán)境氣候、土壤濕度、植株生長狀況等因素是影響旱地稻田蠐蠟個體聚集的主要因素。通過建立理論抽樣和序貫抽樣模型可知，當(dāng)允許誤差D為0.3，蒿枝壩村、哈國馬村和諾國村的蟲口密度均為2頭/m2時，2022年所需的理論抽樣數(shù)分別為353、393和155m2，2023年所需抽樣數(shù)為76、44和49m2。【結(jié)論】云南省普洱市瀾滄縣竹塘鄉(xiāng)3個村的旱地稻種植田蠐蠟共有9種，主要優(yōu)勢種為黑阿鰓金龜和東方絹金龜。蠐增種群在田間呈聚集分布，引起聚集的原因是昆蟲本身的生活習(xí)性和環(huán)境氣候、土壤濕度、植株生長狀況等因素。

關(guān)鍵詞：旱地稻；蠐蠟；地下害蟲；空間分布；序貫抽樣；云南

中圖分類號：S431.116文獻(xiàn)標(biāo)志碼：A文章編號：2095-1191（2024）01-0128-09

Composition structure and spatial distribution of white grubs in upland rice field in Yunnan

LI Zheng-fei1，YANG Bao-yun1，HE Yan-yan12，YANG Hu-quan'，YAN Chao-jian-ping'，CAO Hu-ying'，PENG Yue-jin1，CHEN Bin1，DU Guang-zu1

（'College of Plant Protection，Yunnan Agricultural University/State Key Laboratory for Conservation and Utilization of Bio-resources in Yunnan，Kunming，Yunnan650201，China;2Institute of Resource Plants，Yunnan University，Kunming，Yunnan650500，China）

Abstract：[Objective】To identify the species composition，population spatial distribution pattern and sampling method of white grubs in upland rice field of mountain area in Yunnan，which provided areference for theprediction and comprehensive control of white grubs in upland rice field.【Method]The white grub population was systematically investi gated and identified by five-point sampling method in2022and2023in three main upland rice growing areas，Haozhiba Village，Haguoma Village and Nuoguo Village of Zhutang Township，Lancang County，Pu'er City，Yunnan Province.Aggregation index method，Iwao regression analysis，Taylor power function method and population aggregation mean were used to analyze the spatial distribution patterns and aggregation reasons of white grubs in different experimental sites.Based on the spatial distribution results，theoretical sampling and sequential sampling models were formed.[Result]A total of 9 species of white grubs belonging to6genera and3families were found in3experimental sites in Lancang County，Pu'er City，Yunnan Province.Apogonia cupreoviridis and Maladera orientalis were the dominant species.The relative abundances were30.91%and24.18%，respectively.The aggregation index method，Iwao regression analysis and Taylor power analysis showed that the spatial distributions of white grub in Haozhiba Village，Haguoma Village and Nuoguo Village in2022and2023were all aggregation distributions（alt;0，βgt;1）.And the aggregation mean λvalues were greater than2，the results showed that the living habits of insects，environmental climate，soil moisture and plant growth were the main factors affecting the aggregation of grubs in upland rice field.Through the establishment of theoretical sampling and sequential sampling model，when the allowable error Dwas0.3and the white grub population density in Haozhiba Village，Haguoma Village and Nuoguo Village was2insects/m2，the required theoretical sampling numbers were353，393and155m2in2022，and76，44and49m2in2023respectively.【Conclusion]Atotal of9white grub species are found in upland rice fields in3villages of Zhutang Township，Lancang County，Pu'er City，Yunnan Province.The main dominant species are A.cupreoviridis and M.orientalis，and the grub populations are all aggregated in the field in indi-vidual groups，which are caused by the living habits of insects，environmental climate，soil moisture and plant growth conditions.

Keywords：upland rice;white grubs;undergroundpests;spatial distribution;sequential sampling;Yunnan

Foundation items：Yunnan Major Science and Technology Project（202202AE090036）;Open Project of State Key Laboratory of BiologicalResources Protection and Utilization Jointly Built by Yunnan（gzkf2022004）

0引言

【研究意義】蠐嘈是鞘翅目（Coleoptera）金龜總科（Scarabaeoidea）幼蟲的統(tǒng)稱，是國內(nèi)外公認(rèn)的難預(yù)測及防治的主要地下害蟲。世界上已記載的蠐蠟種類有3萬種以上，我國已記載的種類1800余種，其中為害農(nóng)、林和牧草的蠐蠟有110種以上（張美翠等，2014）。蠐蠟類害蟲食性較雜，很多種類可取食多種作物的種子、塊根、塊莖和幼苗等地下部位，給農(nóng)業(yè)生產(chǎn)造成巨大損失（楊廣東等，2010;吳鶴敏等，2020;Katumanyane et al.，2023）。蠐蠟類害蟲在田間呈不均勻分布，且多個種群混合發(fā)生，生活為害習(xí)性各異，發(fā)生規(guī)律錯綜復(fù)雜等，導(dǎo)致對蠐蠟的田間取樣及調(diào)查比較困難（王文娟等，2012;Jordan et al.，2012）。水稻旱作是一種“雨養(yǎng)旱種”的雜交稻栽培技術(shù)，可讓旱地稻與水田稻一樣分蘗。與傳統(tǒng)水稻相比，旱地稻土層氧氣足，為蠐蠟的發(fā)生提供了良好的環(huán)境。蠐蠟棲于土壤中危害旱地稻的根系，在水稻生長期，蠐曹啃食根和莖，輕則造成植株枯黃，重則造成整株枯死，嚴(yán)重影響旱地稻的產(chǎn)量和品質(zhì)，給種植戶造成嚴(yán)重的經(jīng)濟損失。因此，研究云南山區(qū)旱地稻田蠐蠟組成結(jié)構(gòu)及空間分布特征型，對旱地稻田蠐蠟類害蟲的預(yù)測預(yù)報及綜合防控，保障云南部分山區(qū)和半山區(qū)群眾口糧安全，助推鄉(xiāng)村振興等方面具有重大意義?！厩叭搜芯窟M(jìn)展】目前對蠐蠟的空間分布格局研究主要采用回歸模型分析方法，馬艷芳等（2009）用回歸分析研究不同深度土層蠐蠟的空間分布類型，發(fā)現(xiàn)大栗鰓金龜（Melolontha hippo-castani）3齡幼蟲在0～10cm土壤中呈聚集分布。王文娟等（2012）通過回歸分析法研究蠐蠟在裸大麥田間的空間分布，結(jié)果表明蠐蠟在裸大麥田間屬聚集分布，且蠐蠟個體間相互吸引。劉思雨等（2020）采用回歸分析法研究煙田蠐蠟種類和混合種群的空間分布型，共發(fā)現(xiàn)5種蠐蠟，蠐蠟個體間相互排斥，主要呈聚集分布，且聚集度依賴于密度。陸俊姣等（2021）也采用回歸分析法對小麥田地下害蟲的空間分布進(jìn)行研究，發(fā)現(xiàn)小麥田地下蠐嘈空間分布型為聚集分布。通過對金龜子成蟲的空間分布型進(jìn)行回歸分析，發(fā)現(xiàn)白星花金龜（Protaetia brevitarsis）成蟲在玉米田間和荒漠草原的分布型呈聚集分布（劉政等，2010;賈彥霞等，2011）;麗腹弓角鰓金龜（Toxo-spathius auriventris）成蟲在西藏林芝農(nóng)牧林區(qū)也呈聚集分布（臧建成等，2018）。近年來，對于蠐蠟的研究主要集中于花生田、玉米田和金銀花地等，對這些類型地塊蠐蠟的種類組成、生物學(xué)特性及發(fā)生動態(tài)均有報道，如劉富順等（2022）采用田間調(diào)查和室內(nèi)試驗對河北花生田蠐蠟進(jìn)行研究，明確了該地區(qū)花生田共有9種蠐蠟，暗黑鰓金龜（Holotrichia paral-lela）為主要優(yōu)勢種。李霞等（2023）調(diào)查金銀花田蠐曹的種類及危害程度，結(jié)果表明主要的蠐蠟種類為華北大黑鰓金龜（Holotrichia oblita）、東方絹金龜（Maladera orientalis）和福婆鰓金龜（Brahmina fal-dermanni），危害高峰期為7—9月。【本研究切入點】目前對于蠐蠟的組成結(jié)構(gòu)、空間分布型和抽樣技術(shù)的研究主要集中于煙田、花生田及麥田等，對于旱地稻田中蠐蠟的研究尚未見文獻(xiàn)報道?！緮M解決的關(guān)鍵問題】以云南省普洱市瀾滄縣3個主要旱地稻種植區(qū)為試驗地點，于2022和2023年采用五點取樣法對旱地稻田中的蠐蠟進(jìn)行系統(tǒng)性調(diào)查，明確蠐嘈的種類組成及其在田間的空間分布結(jié)構(gòu)，并模擬形成旱地稻田蠐蠟種群的抽樣技術(shù)，以期為云南山區(qū)旱地稻田蠐蠟類害蟲的預(yù)測預(yù)報及綜合防控提供科學(xué)參考。

1材料與方法

1.1試驗材料及試驗地概況

水稻品種為滇禾優(yōu)615。周圍農(nóng)田種植玉米、甘蔗。試驗地點位于云南省普洱市瀾滄縣竹塘鄉(xiāng)蒿枝壩村（22°44'5\"N，99°47'59\"E，海拔1490m）、哈果馬村（22°46'52\"N，99°47'13\"E，海拔1530m）和諾國村（22°40'41\"N，99°48'4\"E，海拔1370m）旱地稻田。

1.2調(diào)查采集方法

試驗于2022和2023年8月旱地稻拔節(jié)期進(jìn)行。其中2022年每試驗點調(diào)查15塊樣地，2023年每試驗點調(diào)查20塊樣地，每塊樣地面積不小于1ha。調(diào)查方法采用五點取樣法，每塊旱地稻田選取5個樣點，每點調(diào)查1m2，取樣深度20～30cm，分別統(tǒng)計記錄每m2旱地稻根部蠐蠟的種類和數(shù)量。

1.3蠐蠟種類調(diào)查及鑒定

將采集到的蠐蠟按照不同的大小和特征分開收集，每種抽取一部分用沸水煮沸片刻后浸入75%酒精中保存。通過體視鏡對其形態(tài)進(jìn)行觀察和鑒定，鑒定參考《中國經(jīng)濟昆蟲志》（第二十八冊）（張芝利，1984）和《中國常見地下害蟲圖鑒》（曹雅忠和李克斌，2017）等。對于無法鑒定的種類，將活蟲帶回實驗室，并接入牙簽盒內(nèi)單頭飼養(yǎng)。盒內(nèi)約盛4/5容積的濕土，相對濕度控制在18%，飼喂新鮮馬鈴薯塊，飼養(yǎng)至成蟲，進(jìn)行種類鑒定（周洪旭等，2009;許浩等，2012）。

1.4聚集度指標(biāo)計算

將2022—2023年不同試驗點旱地稻田蠐蠟的調(diào)查數(shù)據(jù)整理分析，以每個試驗點所選樣地的不同田塊為單位，分別計算各塊樣地內(nèi)蠐蠟的平均蟲口密度（m）和樣本方差（S2），并在此基礎(chǔ)上，測定平均擁擠度（m*）和蠐蠟種群空間分布的各項指標(biāo)，最終確定蠐在不同樣地的空間分布格局（張孝羲，2002;戈峰，2008;劉思雨等，2020）。計算平均擁擠度（m*）、擴散系數(shù)（C）、叢生指數(shù)（I）、負(fù)二項分布值（K）、久野指數(shù)（CA）和聚塊指數(shù)（C*）5個聚集度指標(biāo)確定旱地稻田蠐蠟混合種群的空間分布型（沈佐銳，2009）。計算公式如下：

m*=m+S2/m-1

C=S2/m

當(dāng)Clt;1、C=0、Cgt;1時，種群分別為隨機分布、均勻分布和聚集分布。

I=（S2-m）/m

當(dāng)K0、I=0、Igt;0時，種群分別為均勻分布、隨機分布和聚集分布。

K=m2/（S-m）

當(dāng)Klt;0、Kgt;0、K→00時，種群分別為均勻分布、聚集分布和隨機分布。

CA=1/K=（S2-m）/m2

當(dāng)CAlt;1、CA=1、CAgt;1時，種群分別為均勻分布、隨機分布和聚集分布。

C*=m*/m

當(dāng)C*lt;1、C*=1、C*gt;1時，種群分別為均勻分布、Poisson分布和聚集分布。

1.5回歸分析

采用回歸分析法擬合蠐蠟在旱地稻田間的空間分布型。

Iwao回歸分析法：m*=α+βm，當(dāng)m*和m呈線性關(guān)系時，用于研究m*與平均蟲口密度之間的相互關(guān)系。式中，α表示平均擁擠度，當(dāng)αgt;0時，表示分布的基本成分是個體群，個體之間相互吸引；α=0時，表示分布的基本成分是單個體；αlt;0時，表示分布的基本成分是單個體，且個體之間相互排斥。β表示相對聚集度，代表著個體群的空間分布型，當(dāng)βlt;1為均勻分布，β=1為隨機分布，βgt;1為聚集分布（陸俊姣等，2021;He et al.，2022）。

Taylor冪函數(shù)法：S2=amb，可轉(zhuǎn)化為回歸方程lgS2=1ga+blgm。式中，a、b為參數(shù)，當(dāng)1ga=0，b=1時，種群為隨機分布；當(dāng)lgagt;0，b=1時，種群為聚集分布，聚集度無密度依賴性；lgagt;0，bgt;1時，種群為聚集分布，聚集度有密度依賴性；Igalt;0，blt;1時，種群為均勻分布，且種群密度越高分布越均勻（He et al.，2022）。

1.6聚集原因分析

利用Blackkith的種群聚集均數(shù)（λ）檢驗蠐蠟在旱地稻田的聚集原因，計算公式：λ=m×r/2K。式中，m為各樣地蠐蠟密度的平均值；K為各樣地蠐蠟負(fù)二項分布的平均值；r為自由度等于2K時的c20s值（概率為0.05）。當(dāng)λlt;2時，表明個體聚集是由氣候、土壤濕度和植株生長狀況等諸多環(huán)境因素引起的；λ≥2時，表明個體聚集是由昆蟲本身的生活習(xí)性和上述環(huán)境因素綜合作用的結(jié)果（王利軍等，2011）。

1.7理論抽樣模型

基于線性回歸方程分析的結(jié)果，平均擁擠度m*和平均密度m呈顯著相關(guān)，可利用m*-m線性回歸聚集度參數(shù)α和β確定蠐蠟在不同蟲口密度下的最適抽樣數(shù)（王琨等，2022），理論最適抽樣數(shù)公式為：N=（t/D）2[（α+1）m+（β-1）]，其中N為最適理論抽樣數(shù)，t為95%置信度下的分布值（取1.96），D為允許誤差值（陸俊姣等，2021）。

1.8序貫抽樣模型

采用Iwao（1976）提出的復(fù)合序貫抽樣技術(shù)：Tou=nm±t√n（a+1）m?+（β-1）m2，擬合得到旱地稻田蠐嘈種群理論抽樣數(shù)上、下限。式中：累積蟲量上限T?取模型加號計算的數(shù)值，下限T?取模型減號計算的數(shù)值。n為田間抽樣數(shù)，m?為臨界防治指標(biāo)，t為相應(yīng)置信度下的分布值，α和β為Iwao回歸模型的聚集參數(shù)（孫小旭等，2019）。

1.9統(tǒng)計分析

試驗數(shù)據(jù)采用DPS14.0進(jìn)行統(tǒng)計分析。

2結(jié)果與分析

2.1蠐嘈種類組成及相對豐度結(jié)果

通過2023年的田間挖土調(diào)查并對旱地稻田的蠐蠟進(jìn)行種類鑒定，云南省普洱市瀾滄縣竹塘鄉(xiāng)3個試驗點的旱地稻種植田蠐蠟共有9種，隸屬于3科6屬（表1）。其中鰓金龜亞科（Melolonthidae）7種，分別為小黃鰓金龜（Metabolusflavescens）、闊脛瑪絹金龜（Maladera verticalis）、小闊脛瑪絹金龜（Maladera ovatula）、東方絹金龜（Maladera orienta-lis）、棕脊鰓金龜（Holotrichia castanea）、暗黑鰓金龜（Holotrichia parallela）、黑阿鰓金龜（Apogonia cupreo-viridis）;麗金龜亞科（Rutelidae）1種，為銅綠麗金龜（Anomala corpulenta）;花金龜亞科（Cetoniinae）1種，為白星花金龜（Potosia brevitarsis）。

蒿枝壩村共采集到蠐蠟標(biāo)本214頭，包含9個種；哈國馬村共采集到標(biāo)本233頭，包括5個種；諾國村共采集到標(biāo)本103頭，包括7個種。蒿枝壩村的優(yōu)勢種為黑阿鰓金龜、東方絹金龜和暗黑鰓金龜，相對豐度分別為35.05%、16.82%和14.49%;哈國馬村蠐嘈的優(yōu)勢種為東方絹金龜、黑阿鰓金龜和闊脛瑪絹金龜，相對豐度分別為35.62%、27.04%和22.75%;諾國村主要蠐蠟優(yōu)勢種為黑阿鰓金龜和闊脛瑪絹金龜，相對豐度分別為31.06%和20.39%。

從3個試驗點采集的情況來看，蠐蠟的物種數(shù)量為：蒿枝壩村gt;諾國村gt;哈國馬村，蠐蠟個體數(shù)量為：哈國馬村gt;蒿枝壩村gt;諾國村。3個試驗點共采集到550頭蠐蠟標(biāo)本，其中，黑阿鰓金龜170頭，相對豐度為30.91%;東方絹金龜133頭，相對豐度為24.18%;闊脛瑪絹金龜93頭，相對豐度為16.91%;暗黑鰓金龜42頭，相對豐度7.64%;這4個種類是瀾滄旱地稻田主要的蠐蠟種類。

2.2蠐蠟的空間分布型—聚集度指標(biāo)分析結(jié)果

從2022和2023年各試驗點旱地稻田蠐蠟種群的各項聚集度指標(biāo)來看（表2），不同年度各試驗點蠐嘈種群的Cgt;1、Igt;0、Kgt;0、C*gt;1，均符合聚集分布的指標(biāo)要求，表明不同年度各試驗點旱地稻田蠐蠟種群在土壤中的分布類型均屬于聚集分布。

2.3蠐蠟空間分布型一回歸分析結(jié)果

Iwao回歸分析法表明，在2022年調(diào)查的15塊樣地和2023年調(diào)查的20塊樣地中，不同試驗點旱地稻田蠐嘈種群m*-m回歸方程均表現(xiàn)為：αlt;0，βgt;1，表明不同年度各試驗點旱地稻田的蠐嘈種群的空間分布類型為聚集分布，且個體間相互排斥，分布的基本成分是單個體（表3）。

Taylor冪函數(shù)法表明，各年份不同試驗點樣地調(diào)查的參數(shù)均為bgt;1，說明旱地稻田中的蠐蠟種群均為聚集分布，與用聚集度方法測定的結(jié)果一致（表4）。

2.4聚集原因分析結(jié)果

從2022和2023年各試驗點旱地稻田蠐蠟種群聚集數(shù)可知（表5）。2023年各試驗點旱地稻田的蟲口密度均高于2022年，諾國試驗點蟲口密度最低。不同年度3個試驗點的旱地稻田λ值均gt;2，表明昆蟲本身的生活習(xí)性和環(huán)境氣候、土壤濕度、植株生長狀況等因素是影響旱地稻田蠐蠟個體聚集的主要因素。

2.5蠐蠟理論抽樣模型分析結(jié)果

2022年不同試驗點旱地稻田蠐的理論抽樣模型分別為：蒿枝壩村N=（t/D）2（0.118/m+8.222）、哈國馬村N=（t/D）2（-0.184/m+9.308）、諾國村N=（t/D）2（0.679/m+3.282）;2023年的理論抽樣模型分別為：蒿枝壩村N=（t/D）2（0.623/m+1.459）、哈國馬村N=（t/D）2（0.460/m+0.810）、諾國村N=（t/D）2（0.670/m+0.805）（表6）。當(dāng)蟲口密度一定時，隨著允許誤差增大，所需的抽樣數(shù)減少。當(dāng)允許誤差和蟲口密度一定的情況下，各年度不同試驗點旱地稻田所需的抽樣數(shù)不同，表現(xiàn)為2022gt;2023年。不同試驗點旱地稻田所需最適抽樣數(shù)也不同，在2023年，所需最適抽樣數(shù)表現(xiàn)為蒿枝壩村gt;諾國村gt;哈國馬村；2022年，所需最適抽樣數(shù)表現(xiàn)為哈國馬村gt;蒿枝壩村gt;諾國村。例如，當(dāng)允許誤差為0.3時，蠐蠟的蟲口密度為2頭/m2時，在2022年，蒿枝壩村、哈國馬村、諾國村所需的理論抽樣數(shù)分別為353、393和155m2。而2023年，蒿枝壩村、哈國馬村、諾國村所需的理論抽樣數(shù)分別為76、44和49m2。因此，在不同試驗點旱地稻田的調(diào)查中，可根據(jù)允許誤差的范圍及田間蠐蠟蟲口密度，參照表6確定田間調(diào)查時的抽樣數(shù)。

2.6蠐蠟序貫抽樣模型分析結(jié)果

2022年不同試驗點旱地稻田蠐蠟的序貫抽樣模型分別為：蒿枝壩村Tou=4n±2.36√n、哈國馬村Tou）=4n±2.21√n、諾國村To（1）=4n±1.14√n;2023年的序貫抽樣模型分別為：蒿枝壩村Toc=4n±5.08√n、哈果馬村Tou=4n±3.85√n、諾國村Tau=4n±12.80√n（表7）。各年份不同試驗點的旱地稻田蠐蠟的防治指標(biāo)也有所差異。如在2023年，當(dāng)調(diào)查蒿枝壩村2m2內(nèi)蠐蠟混合種群時，如果累積蟲量大于13頭，則需要進(jìn)行防治；若小于3頭，則不需要進(jìn)行防治；累積蟲量在二者之間，則繼續(xù)調(diào)查。直至達(dá)到理論抽樣數(shù)80m2，此時累積蟲量如果大于354頭，則需要進(jìn)行防治。當(dāng)調(diào)查哈國馬村和諾國村蠐蠟混合種群密度時，如果累積蟲量大于抽樣上限，則需要進(jìn)行防治；若小于抽樣下限，則不需要進(jìn)行防治；如果累積蟲量介于二者之間，則需繼續(xù)調(diào)查。

3討論

昆蟲的空間分布型是其種群個體在自然空間的散布狀態(tài)，其形成不僅與每種昆蟲自身習(xí)性有密切的關(guān)系，也受環(huán)境因子的影響（Leite et al.，2006;朱金良等，2007;Voss et al.，2009）。研究昆蟲空間分布型不僅能揭示昆蟲在空間的分布特性，還能為昆蟲在田間的抽樣和防治提供科學(xué)依據(jù)（Martins et al.，2018;李召波等，2019;李少華等，2022）。研究蠐嘈在土壤中的分布格局，對制定調(diào)查方法、預(yù)測蟲情和制定防治措施均具有重要意義。本研究結(jié)果表明2022和2023年云南省普洱市瀾滄縣旱地稻種植區(qū)生境下蠐蠟的主要優(yōu)勢種為黑阿鰓金龜和東方絹金龜。與劉思雨等（2020）調(diào)查發(fā)現(xiàn)云南大理煙田蠐蠟的主要優(yōu)勢種為小黃鰓金龜和銅綠麗金龜；李霞等（2023）發(fā)現(xiàn)的山東平邑金銀花田主要蠐嘈種類為華北大黑鰓金龜和東方絹金龜?shù)南撓瀮?yōu)勢種間存在差異。猜測與地域特點、氣候條件、土壤濕度、土壤養(yǎng)分、土壤理化性質(zhì)和土壤微生物群落結(jié)構(gòu)（劉福順等，2022）以及當(dāng)?shù)氐姆N植模式、前茬作物類型和寄主植物種類的差異有關(guān)（Cherman et al.，2014;曲春娟等，2023;Castelo and Crespo，2023）。

本研究發(fā)現(xiàn)，2022和2023年云南省普洱市瀾滄縣3個村旱地稻上的蠐蠟種群均呈聚集分布，與云南大理煙田、山西小麥田和西藏裸大麥田蠐蠟種群呈聚集性分布的結(jié)果（王文娟等，2012;劉思雨等，2020;陸俊姣等；2021）一致。猜測蠐蠟在我國不同地理位置、不同寄主植物上的空間分布型多表現(xiàn)為聚集性分布，可能與蠐蠟成蟲金龜子的產(chǎn)卵偏好有關(guān)。本研究通過對聚集原因的分析，發(fā)現(xiàn)不同試驗點旱地稻田的蠐蠟聚集因素λgt;2，表明樣地蠐蠟個體的聚集是由昆蟲本身的生活習(xí)性和環(huán)境氣候、土壤濕度、植株生長狀況等因素引起，與陸俊姣等

（2021）對小麥田地下害蟲混合種群的研究結(jié)果一致。而與劉思雨等（2020）調(diào)查發(fā)現(xiàn)煙田蠐蠟聚集是由環(huán)境導(dǎo)致的結(jié)論有所不同，可能由于調(diào)查生境下蠐蠟種類的差異造成；不同金龜子產(chǎn)卵習(xí)性存在很大差異，如粗狹肋鰓金龜（Holotrichiascarobiculata）交配后，卵多數(shù)散產(chǎn)于土中，而一些金龜子則是將卵聚產(chǎn)在植株根部（王興紅等，2022）。此外，蠐曹生活史上的差異也可能對聚集原因造成影響，如黑阿鰓金龜為2年1代，而小黃鰓金龜則為1年1代。

作物的不同會對金龜子產(chǎn)卵行為產(chǎn)生影響，曲春娟等（2023）研究表明暗黑鰓金龜產(chǎn)卵選擇明顯偏好花生，而銅綠麗金龜產(chǎn)卵更偏向于選擇玉米。劉福順等（2022）研究表明土壤質(zhì)地會影響蠐蠟的發(fā)生水平，壤土中的蠐蠟發(fā)生量遠(yuǎn)高于黏土和砂土；秸稈還田有利于蠐蠟種群生存的環(huán)境條件，會增加花生田蠐蠟的數(shù)量。宿主寄生模式會對蠐蠟的空間分布產(chǎn)生影響，寄生強度會造成局部蠐蠟密度的變化（Castelo and Crespo，2023）。調(diào)查不同生境下蠐蠟的空間分布，還應(yīng)考慮周圍的外界環(huán)境條件影響。山區(qū)旱地稻田周圍生態(tài)環(huán)境較好，植被茂盛，下一步調(diào)查時還應(yīng)具體分析周邊寄主植物以及土壤質(zhì)地等因素與蠐蠟空間分布格局之間的關(guān)系。

蠐嘈生活在土壤中，調(diào)查監(jiān)測和防治難度大（仵均祥，2016）。因此，確定科學(xué)的理論抽樣量，可保證調(diào)查時的準(zhǔn)確性和時效性并能節(jié)省人力（沈佐銳，2009;鄭伯平等，2011）。蠐蠟類害蟲在田間常以多個種群混合發(fā)生危害，不同種類蠐蠟造成的為害損失不同（Echeverri-Molina and Govender，2016）。本研究調(diào)查的旱地稻田蠐蠟種類中，白星花金龜以成蟲危害，而其幼蟲為腐食性，不造成危害（楊誠，2015）。東方絹金龜成蟲是林業(yè)上的大害蟲，而其幼蟲對作物的危害不大（仵均祥，2016）。黑阿鰓金龜體型較小，對根系造成的危害也較小。暗黑鰓金龜則是造成嚴(yán)重危害的種類，成蟲和幼蟲均嚴(yán)重危害植物，會造成農(nóng)林業(yè)的重大損失（朱長建，2023）。因此，對旱地稻田的蠐曹進(jìn)行調(diào)查時，還應(yīng)結(jié)合本研究旱地稻田蠐蠟的種類來綜合確定蠐曹的抽樣數(shù)量和防治標(biāo)準(zhǔn)，如蒿枝壩村旱地稻田相比其他2個村暗黑鰓金龜?shù)臄?shù)量更多，受害可能性較大，防治時應(yīng)重點關(guān)注；而哈國馬村主要為危害較小的種類，應(yīng)適當(dāng)放寬防治的標(biāo)準(zhǔn)，減少化學(xué)農(nóng)藥的施用。在蠐蠟的防治中，還應(yīng)結(jié)合不同田塊優(yōu)勢種的差異，制定最佳防治策略，比如，根據(jù)不同種類蠐蠟對不同農(nóng)藥的抗藥性和敏感性的差異，選擇合適的農(nóng)藥混搭和配比方案；同時結(jié)合生物防治，考慮蠐蠟個體差異篩選對不同田塊蠐蠟優(yōu)勢種高致病性的白僵菌、綠僵菌菌株。

本研究未詳細(xì)統(tǒng)計旱地稻田不同蠐蠟種類的分布情況，且未對危害較嚴(yán)重的種類進(jìn)行系統(tǒng)分析，之后的研究中，將進(jìn)一步探索蠐蠟在田間的分布與土壤微生物群落結(jié)構(gòu)的相關(guān)性；研究造成不同種類危害差異的機理；詳細(xì)分析主要危害類群在田間的空間分布結(jié)構(gòu)；模擬更準(zhǔn)確的旱地稻田蠐蠟種群的抽樣技術(shù)及建立蠐蠟主要危害種類在不同蟲口密度下的序貫抽樣模型和理論最適抽樣數(shù)，為該試驗點旱地稻田蠐蠟類害蟲的預(yù)測預(yù)報及田間防控提供更精準(zhǔn)的理論指導(dǎo)。

4結(jié)論

云南省普洱市瀾滄縣竹塘鄉(xiāng)3個村的旱地稻種植田蠐蠟共有9種，分別為小黃鰓金龜、闊脛瑪絹金龜、小闊脛瑪絹金龜、東方絹金龜、棕脊鰓金龜、暗黑鰓金龜、黑阿鰓金龜、銅綠麗金龜和白星花金龜。其中，黑阿鰓金龜和東方絹金龜為主要優(yōu)勢種。3個村的旱地稻中的蠐蠟種群均以個體群在田間呈聚集分布，引起聚集的原因是昆蟲本身的生活習(xí)性、環(huán)境氣候、土壤濕度和植株生長狀況等因素。

參考文獻(xiàn)：

戈峰.2008.昆蟲生態(tài)學(xué)原理與方法[M].北京：高等教育出版社.[Ge F.2008.Principle and methods of insect eco-logy[M].Beijing：Higher Education Press.]

賈彥霞，賀奇，王新譜.2011.寧夏荒漠草原白星花金龜成蟲空間分布型的初步研究[J].草地學(xué)報，19（1）：177-180.[Jia YX，He Q，Wang XP.2011.A preliminary study on the spatial distribution pattern of Potosia brevitaris adults[J].Acta Agrestia Sinica，19（1）：177-180.]doi：10.3969/j.issn.1007-0435.2011.01.030.

李少華，王云鵬，王榮成，尹萍，李向東，鄭方強.2022.玉米田桃蛀螟幼蟲的空間分布型與抽樣技術(shù)[J].中國農(nóng)業(yè)科學(xué)，55（10）：1961-1970.[LiSH，Wang YP，Wang RC，Yin P，LiX D，Zheng FQ.2022.Spatial distribution pat-tern and sampling technique of Conogethes punctiferalis larvae inmaize fields[J].Scientia Agricultura Sinica，55（10）：1961-1970.]

李霞，于毅，宋振巧，門興元，崔洪瑩，李麗莉，宋瑩瑩，王建華，郭文秀.2023.山東平邑地區(qū)金銀花田蠐蠟的種類調(diào)查[J].山東農(nóng)業(yè)科學(xué)，55（3）：142-147.[LiX，YuY，Song"Z Q，MenXY，Cui HY，Li LL，Song YY，Wang JH，Guo WX.2023.Species investigation of grubs in honey suckle fields in Pingyi，Shandong Province[J].Shandong Agricultural Sciences，55（3）：142-147.]doi：10.14083j.issn.1001-4942.2023.03.021.

李召波，夏濤，王應(yīng)清，曾千春，馬春花，馬麗宣，李強，高熹，田洋.2019.辣木瑙螟幼蟲的空間分布型及抽樣技術(shù)研究[J].云南農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)），34（4）：597-602.[LiZ B，Xia T，Wang YQ，Zeng QC，Ma CH，Ma LX，LiQ，Gao X，Tian Y.2019.Research on the spatial distri-bution pattern and sampling technique of Noorda blitealis walk larva[J].Journal of Yunnan Agricultural University（Natural Science），34（4）：597-602.]doi：10.12101j.issn.1004-390X（n）.201804020.

劉福順，馮曉潔，劉春琴，吳娛，席國成，王慶雷.2022.河北滄州花生田蠐蠟發(fā)生動態(tài)及影響因素分析[J].中國植保導(dǎo)刊，42（4）：33-37.[Liu FS，F(xiàn)eng XJ，Liu CQ，Wu Y，Xi GC，Wang QL.2022.Analysis on the occurrence dyna-mics and influencing factors of grubs in peanut fields in Cangzhou，Hebei Province[J].China Plant Protection，42（4）：33-37.]doi：10.3969/j.issn.1672-6820.2022.04.006.

劉思雨，戶艷霞，范志勇，杜廣祖，孫軍偉，陳斌.2020.煙草田間蠐蠟種類及其空間分布與抽樣技術(shù)研究[J].南方農(nóng)業(yè)學(xué)報，51（11）：2722-2730.[LiuSY，HuYX，F(xiàn)anZY，Du"GZ，SunJW，Chen B.2020.Species composition，spatial distribution pattern and sampling technique of white grubs in tobacco fields[J].Journal of Southern Agriculture，51（11）：2722-2730.]doi：10.3969/j.issn.2095-1191.2020.11.015.

劉政，孫艷，王少山，李國英.2010.白星花金龜在玉米田的空間分布和抽樣技術(shù)[J].植物保護(hù)，36（6）：125-127.[Liu Z，Sun Y，Wang SS，LiGY.2010.Spatial distribution pat-tern and sampling model for the adults of Liocola brevitar-sis incornfields[J].Plant Protection，36（6）：125-127.]

陸俊姣，任美鳳，李大琪，李霞，史高川，董晉明.2021.小麥田地下害蟲幼蟲的空間分布型與抽樣技術(shù)[J].植物保護(hù)學(xué)報，48（2）：340-347.[Lu JJ，Ren MF，LiD Q，LiX，Shi GC，Dong JM.2021.Spatial distribution pattern and sam-pling technique of the larvae of subterranean pest insects in wheat fields[J].Journal of Plant Protection，48（2）：340-347.]doi：10.13802/j.cnki.zwbhxb.2021.2020073.

馬艷芳，謝宗謀，張永強，郭彥林.2009.大栗鰓金龜幼蟲空間分布型研究[J].植物保護(hù)，35（4）：140-142.[Ma YF，Xie ZM，Zhang YQ，Guo YL.2009.Spatial distribution pat-terns of Melolontha hippocastani larvae[J].PlantProtec-tion，35（4）：140-142.]doi：10.3969/j.issn.0529-1542.2009.04.033.

曲春娟，姜曉靜，萬書波，張正，杜龍，鞠倩，李曉，曲明靜.2023.花生/玉米間作對暗黑鰓金龜和銅綠麗金龜發(fā)生的影響[J/OL].中國油料作物學(xué)報.https：/link.cnki.net urlid/42.1429.S.20230810.1550.002.[Qu CJ，Jiang XJ，Wan SB，Zhang Z，Du L，Ju Q，Li X，Qu MJ.2023.Effects of peanut/maize intercropping system on theoccur-rence of Holotrichia parallela and Anomala corpulenta"[J/OL].Chinese Journalof Oil CropSciences.https：//link.cnki.net/urlid/42.1429.S.20230810.1550.002.]doi：10.19802/j.issn.1007-9084.2023039.

沈佐銳.2009.昆蟲生態(tài)學(xué)及害蟲防治的生態(tài)學(xué)原理[M].北京：中國農(nóng)業(yè)大學(xué)出版社.[Shen ZR.2009.The ecologi-cal principle of insect ecology and pest control[M].Bei-ing：China Agricultural University Press.]

孫小旭，趙勝園，靳明輝，趙慧媛，李國平，張浩文，姜玉英，楊現(xiàn)明，吳孔明.2019.玉米田草地貪夜蛾幼蟲的空間分布型與抽樣技術(shù)[J].植物保護(hù)，45（2）：13-18.[Sun XX，Zhao SY，Jin MH，Zhao HY，LiGP，Zhang HW，Jiang YY，Yang XM，Wu KM.2019.Larval spatial distribution pattern and sampling technique of the fall armyworm"Spodoptera frugiperdain maize fields[J].Plant Protection，45（2）：13-18.]doi：10.16688/j.zwbh.2019115.

王錕，張麗蘭，黎斌，杜廣祖，陳斌.2022.山地玉米植株及田間草地貪夜蛾幼蟲的空間分布特征[J].南方農(nóng)業(yè)學(xué)報，53（10）：2894-2903.[Wang K，Zhang LL，LiB，Du GZ，Chen B.2022.Spatial distribution of larvae Spodoptera frugiperda（J.E.Smith）in mountainous corn plants and field[J].Journal of Southern Agriculture，53（10）：2894-2903.]doi：10.3969/j.issn.2095-1191.2022.10.021.

王利軍，郭文超，徐建軍，肖利鳳，韓振忠，張建萍.2011.馬鈴薯甲蟲空間分布型及抽樣技術(shù)研究[J].環(huán)境昆蟲學(xué)報，33（2）：147-153.[Wang LJ，Guo WC，XuJJ，Xiao LF Han ZZ，Zhang JP.2011.The spatial distribution pattern and sampling techniques of the colorado potato beetle，Leptinotarsa decemlineata（Say）[J].Journal of Environ-mental Entomology，33（2）：147-153.]doi：10.3969/j.issn.1674-0858.2011.02.002.

王文娟，唐小琴，孟長軍.2012.裸大麥田蠐蠟的空間分布型及抽樣技術(shù)研究[J].江蘇農(nóng)業(yè)科學(xué)，40（8）：105-107.[Wang WJ，Tang XQ，Meng CJ.2012.Spatial distribu-tion and sampling technology of grubs innaked barley field[J].Jiangsu Agricultural Sciences，40（8）：105-107.]doi10.15889/j.issn.1002-1302.2012.08.077.

王興紅，劉曼，黃振興，楊廣明，羅麗林，張歆逢.2022.粗狹肋鰓金龜?shù)纳飳W(xué)特性研究[J].中國南方果樹，51（3）：174-178.[Wang XH，Liu M，Huang ZX，Yang GM，Luo LL，Zhang XF.2022.Biological characteristics ofHolo-trichiascrobiculata Brenske[J].South China Fruits，51（3）：174-178.]doi：10.13938/j.issn.1007-1431.20210616.

吳鶴敏，陳琦，蔣月麗，范志業(yè)，沈海龍，劉迪，侯艷紅，陳莉，李雷雷，王文豪，李世民.2020.河南漯河地區(qū)金龜優(yōu)勢種及種群發(fā)生動態(tài)[J].中國植保導(dǎo)刊，40（7）：38-43.[Wu HM，Chen Q，Jiang YL，F(xiàn)an ZY，Shen HL，Liu D，Hou YH，Chen L，LiLL，Wang WH，LiS M.2020.Dominant species and population dynamics of scarabaeoidea in Luohe area of Henan Province[J].China Plant Protection40（7）：38-43.]doi：10.3969/j.issn.1672-6820.2020.07.008.

仵均祥.2016.農(nóng)業(yè)昆蟲學(xué)[M].北京：中國農(nóng)業(yè)出版社.[Wu JX.2016.Agricultural entomology[M].Beijing：China Agri-culture Press.]

許浩，文禮章，陳坤.2012.湖南煙田常見金龜子種類及其識別調(diào)查研究[J].中國農(nóng)學(xué)通報，28（21）：221-228.[XuH，Wen LZ，Chen K.2012.Identification of common scara-baeid occurrence in tobacco fields in Hunan Province，China（Coleoptera，Scarabaeidae）[J].Chinese Agricultural Science Bulletin，28（21）：221-228.]doi：10.3969/j.issn1000-6850.2012.21.040.

楊誠.2015.白星花金龜生物學(xué)及其對玉米秸稈取食習(xí)性的研究[D].泰安：山東農(nóng)業(yè)大學(xué).[Yang C.2015.Studies on the biology of Liocola brevitasis（Lewis）and its fee-ding habit on maize straw[D].Tai'an：Shandong Agricul-tural University.]

楊廣東，熊桂蓮，張文輝.2010.三種常見地下害蟲的發(fā)生規(guī)律及綜合防治方法[J].現(xiàn)代園藝，（4）：32-33.[Yang GD，Xiong GL，Zhang WH.2010.The occurrence regula-rity and comprehensive control measures of three common underground pests[J].Contemporary Horticulture，（4）：32-33.]doi：10.14051/j.cnki.xdyy.2010.04.006

曹雅忠，李克斌.2017.中國常見地下害蟲圖鑒[M].北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社.[YuYZ，LiK B.2017.A Guide tocommon underground pests in China[M].Beijing：China Agricultural Science and Technology Press.]

臧建成，張帥，相棟.2018.麗腹弓角鰓金龜成蟲在柳樹上的空間分布型及抽樣技術(shù)[J].高原農(nóng)業(yè)，2（5）：507-511.[Zang JC，Zhang S，Xiang D.2018.Spatial distribution patternsand sampling techniques of Toxospathiusauriven-tris bates adults on willows[J].Journal of Plateau Agricul-ture，2（5）：507-511.]doi：10.19707/j.cnki.jpa.2018.05.008.

張美翠，尹姣，李克斌，曹雅忠，仵均祥.2014.地下害蟲蠐蠟的發(fā)生與防治研究進(jìn)展[J].中國植保導(dǎo)刊，34（10）：20-28.[Zhang MC，Yin J，LiK B，Cao YZ，Wu JX.2014.Research progress on theoccurrences of white grub and its control[J].China Plant Protection，34（10）：20-28.]doi：10.3969/j.issn.1672-6820.2014.10.005.

張孝羲.2002.昆蟲生態(tài)及預(yù)測預(yù)報第三版[M].北京：中國農(nóng)業(yè)出版社.[Zhang XX.2002.Insect ecology and fore-cast（The3\"edition）[M].Beijing：China Agriculture Press.]

張芝利.1984.中國經(jīng)濟昆蟲志：鞘翅目金龜總科幼蟲[M].北京：科學(xué)出版社.[Zhang ZL.1984.Economic entomo-logy of China：Larvae of Coleoptera Scarabaeoidea[M].Beijing：Science Press.]

鄭伯平，鄭長英，顧松東，王俊平.2011.西花薊馬（Franklini-ella occidental）在白三葉草上的空間分布[J].青島農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)版），28（1）：47-49.[Zheng BP，

Zheng CY，Gu SD，Wang JP.2011.The spacial distribu-tion pattern of Franklimiella occidentalis on Trifolium""repens[J].Journal of Qingdao Agricultural University（Natu"ral Science），28（1）：47-49.]doi：10.3969/J.ISSN.1674-148X.2011.01.013.

周洪旭，譚秀梅，李長友，鄭桂玲，李國勛.2009.營養(yǎng)和濕度對華北大黑鰓金龜生長發(fā)育和生殖的影響[J].華北農(nóng)學(xué)報，24（4）：201-204.[Zhou HX，Tan XM，LiCY，Zheng"G L，LiGX.2009.Effect of nourishment and humidity on"the development and reproduction of Holotrichia oblita[J].Acta Agriculturae Boreali-Sinica，24（4）：201-204.]

朱金良，王華弟，陳躍，徐云，吳降星，蔡國梁.2007.麥田灰飛虱種群空間分布型及抽樣技術(shù)探討[J].中國農(nóng)學(xué)通報，23（8）：382-386.[Zhu JL，Wang HD，Chen Y，XuY，Wu JX，Cai GL.2007.Spatial distribution of small brown"planthopper，Laodelphaxstriatellus Fallen，population and"its sampling technique in wheat and barley fields[J].Chi-nese Agricultural Science Bulletin，23（8）：382-386.]doi：10.3969/j.issn.1000-6850.2007.08.083.

朱長建，馬雨夢，李斌，梁佳惠，曹付，楊洋，劉蘇，李茂業(yè).2023.對暗黑鰓金龜幼蟲高毒力的球孢白僵菌篩選及其模擬田間防效評價[J].昆蟲學(xué)報，66（8）：1009-1019.[Zhu CJ，MaYM，Li B，Liang JH，Cao F，Yang Y，Liu S，Li MY.2023.Screening and simulated field efficacy evaluation of Beauveria bassianawih high virulence to""Holotrichia parallela（Coleoptera：Scarabaeidae）larvae[J].Acta Entomologica Sinica，66（8）：1009-1019.]doi：10.16380/j.kcxb.2023.08.002.

CasteloM K，Crespo JE.2023.Habitats and parasitoid abun-dance influence spatial density dependence patterns，ren-dering an asilid fly as apotential biological controller of white grubs[J].Frontiers in Agronomy，5：1029232.doi：10.3389/fagro.2023.1029232.

Cherman MA，Morón MA，Dal Prá E，Valmorbida I，Guedes JV C.2014.Ecological characterization of white grubs（Coleoptera：Melolonthidae）community in cultivated and noncultivated fields[J].Neotropical Entomology，43：282-288.doi：10.1007/s13744-014-0214-0.

Echeverri-Molina D，GovenderP.2016.Community structure and morphospecies composition of whitegrubs（Coleop-tera：Scarabaeidae）attacking plantation Acacia mearnsii seedlings inKwaZulu-Natal，South Africa[J].African Ento-mology，24（1）：170-179.

He Y，Wang K，Du G，Zhang Q，Li B，Zhao L，He P，Chen B.2022.Temporal and spatial distributionpatterns ofSpodop tera frugiperda in mountain maize fields inChina[J].In-sects，13（10）：938.doi：10.3390/insects13100938.

Jordan TA，Youngman RR，LaubC L，Tiwari S，Kuhar TP，Balderson TK，Moore DM，Saphir M.2012.Fall soi sampling method for predicting spring infestation of white grubs（Coleoptera：Scarabaeidae）in corn and the benefits of clothianidin seed treatment inVirginia[J].Crop protec tion，39：57-62.doi：10.1016/j.cropro.2012.04.006.

Katumanyane A，Slippers B，Wondafrash M，Malan AP，Hur-ley BP.2023.Susceptibility of white grubs from forestry and sugarcane plantations in South Africa to entomopatho-genic nematodes[J].Biocontrol，68（2）：155-167.doi：org/10.1007/s10526-023-10185-7.

Leite GL D，Picanco M，Zanuncio JC，Ecole CC.2006.Fac-tors affecting herbivory of Tarips palmi（Thysanoptera：Thripidae）and Aphis gossypii（Homoptera：Aphididae）on"the eggplant（Solanum melongena）[J].Brazilian Archives"of Biology and Technology，49（3）：361-369.doi：10.1590/s?516-89132006000400002.

Martins JC，Pican?o MC，Silva RS，Alfredo G，Tarcisio G，Guedes RN.2018.Assessing the spatial distribution of Tuta absoluta（Lepidoptera：Gelechiidae）eggs in open-field tomato cultivation through geostatistical analysis[J].Pest Management Science，74，30-36.doi：10.1002/ps.4664

VossS C，Spafford H，Dadour IR.2009.Annual and seasonal patterns of insect succession on decomposing remains at two locations in Western Australia[J].Forensic Science International，193（1-3）：26-36.doi：10.1016/j.forsciint.2009.08.014.

（責(zé)任編輯 李洪艷）