蚯蚓是土壤生态系统主要动物类群之一，具有重要的生态学与环境指示功能，能够促进植物根系营养成分吸收，改善土壤结构等(丁亦男，王帅，2010)。蚯蚓隶属于环节动物门Annelida寡毛纲Oligochaeta后孔寡毛目Opisthopora，是一种重要的大型土壤动物。截止到2011年，我国共记录陆栖蚯蚓303种，隶属9科31属(徐芹，肖能文，2011)。

近年来，蚯蚓分类学研究主要侧重于通过形态学或辅以分子生物学方法对研究地区的蚯蚓物种多样性进行系统调查与分析，总结出该地区最新的蚯蚓名录和新种的详细记录与描述，如：王海军(2005)利用形态学方法对采集于湖北、贵州和广东等部分地区的蚯蚓样本进行物种鉴定，共记录6个新种，更新了研究地区的蚯蚓物种名录。同时，应用对应因子分析(FAC)综合分析文献数据，初步阐明我国环毛类蚯蚓的系统演化关系。 孙静(2009)通过对2006年在海南地区采集的大量蚯蚓样本进行形态学鉴定，发现蚯蚓新种10种，同时将数值分类学原理成功应用到蚯蚓分类学。孙倩(2010)对在胶东与辽东地区采集的蚯蚓样本进行鉴定，编制了最新的蚯蚓物种多样性名录。Blakemore(2013)通过调查蒙古地区陆栖蚯蚓多样性，共记录到诺登爱胜蚓Eisenia nordenskioldi的2个新亚种，分别是E. nordenskioldi mongol和E. nordenskioldi onon ssp. n.，并对其进行了详细的形态学描述与DNA条形码分析。

然而，鲜有利用形态学方法鉴定并分析蚯蚓优势物种的相关研究，尤其是河北地区的相关研究目前还未见详细报道，故本文对河北廊坊地区蚯蚓优势种形态多样性进行了初步分析。

河北廊坊地区位于华北平原北部，属大陆性季风气候，年平均降雨量540 mm，年平均气温11.5 ℃，土壤类型为潮壤土(杜希勋，吴志刚，1986)(表 1)。

本次调查在每种生境类型中随机选取3个样方，每个样方面积50 cm×50 cm，土层深度为0~15 cm。蚯蚓分布生境主要包括次生林地、菜园、农田、潮湿的草垛和枯枝落叶层等。挖掘法采集蚯蚓个体。野外采集到的蚯蚓置于潮湿土壤带回实验室，在采集当日用70%酒精处死，以95%酒精或无水乙醇保存，全部个体现保存于廊坊师范学院生命科学学院浸制动物标本室。

利用游标卡尺测量蚯蚓个体的体长、体宽，体视显微镜观察其他外部特征(体节数、体色、背孔始位、环带形状及所占体节数、雌孔、雄孔、受精囊孔位置与数目)和内部解剖特征(受精囊形状、大小与数量等)，依据《中国陆栖蚯蚓》(徐芹，肖能文，2011)介绍的蚯蚓形态学特征为鉴定标准。对蚯蚓优势种不同采集群体的形态结构(体长、体宽、体节数)进行比较，量化形态多样性指标。

运用Statistic 8.0和Excel 2007进行聚类、相关性和方差分析。

本次调查共记录蚯蚓3科6属9种(附录)，其中优势种为天锡杜拉蚓Drawida gisti (n=55)，日本杜拉蚓D. japonica(n=78)，梯形流蚓Aporrectodea trapezoides(n=37)。

天锡杜拉蚓 隶属链胃蚓科Moniligastridae杜拉蚓属Drawida。口前叶为前叶式。背孔自2/3节间始。环带位于Ⅹ~节；雄孔1对，位于10/11节间，阴茎高而尖，尖端朝向体前；雌孔1对，位于11/12腹面节间，较雄孔不明显。受精囊孔1对，位于7/8腹面节间沟。受精囊1对，导管在7/8隔膜后多圈盘绕。身体暗褐色或灰色，环带灰白(95%酒精浸泡)。

日本杜拉蚓 隶属链胃蚓科杜拉蚓属。无背孔。口前叶为前叶式。环带位于Ⅹ~节，雄孔1对，位于第10节腹面后缘。砂囊2~3个，位于Ⅶ~Ⅵ Ⅹ节。受精囊较天锡杜拉蚓小而圆，位于7/8隔膜后方。身体灰褐色，环带白色(95%酒精浸泡)。

梯形流蚓 隶属正蚓科Lumbricidae流蚓属Aporrectodea。口前叶为上叶式。 背孔自12/13节间始。环带位于Ⅹ Ⅹ Ⅶ~Ⅹ Ⅹ Ⅹ Ⅳ节。性隆脊位于Ⅹ Ⅹ Ⅺ~Ⅹ Ⅹ Ⅹ Ⅲ节。雄孔位于第15节腹面。受精囊孔2对，位于9/10与10/11腹面节间。身体后部背腹扁平。身体红褐色，头部偏暗(95%酒精浸泡)。

经统计分析，研究数据符合正态分布。天锡杜拉蚓的体节数为51~207，体长为59.80 mm±19.15 mm，体宽为4.30 mm±1.21 mm；日本杜拉蚓体节数为54~174，体长为42.00 mm±10.90 mm，体宽为3.00 mm±0.69 mm；梯形流蚓体节数为67~191，体长为63.40 mm±20.72 mm，体宽为4.50 mm±0.92 mm(表 2)。

3个优势种的个体体长与体宽之间的显著正相关反映了蚯蚓外部结构的发育特点，即体长随体宽的增大而增加(图 1~图 3)，体长与体节数之间极显著的回归关系亦表明体长随体节数的增加而增长(图 4~图 6)。

图 1 天锡杜拉蚓体长与体宽的相关性(n=55) Fig. 1 Correlation between the length and width of Drawida gisti (n=55)

图选项

图 2 日本杜拉蚓体长与体宽的相关性(n=78) Fig. 2 Correlation between the length and width of Drawida japonica (n=78)

图选项

图 3 梯形流蚓体长与体宽的相关性(n=37) Fig. 3 Correlation between the length and width of Aporrectodea trapezoides (n=37)

图选项

图 4 天锡杜拉蚓体长与体节数的相关性(n=55) Fig. 4 orrelation between the length and segment number of Drawida gisti (n=55)

图选项

图 5 日本杜拉蚓体长与体节数的相关性(n=78) Fig. 5 Correlation between the length and segment number of Drawida japonica (n=78)

图选项

图 6 梯形流蚓体长与体节数的相关性(n=37) Fig. 6 Correlation between the length and segment number of Aporrectodea trapezoides (n=37)

图选项

3种蚯蚓的体长聚类分析符合种属差异，即同属的天锡杜拉蚓与日本杜拉蚓 身体长度较正蚓科的梯形流蚓更为相似(图 8)，然而天锡杜拉蚓较粗，日本杜拉蚓与梯形流蚓的体宽较为接近(图 7)。

图 8 3种蚯蚓体节数的聚类分析(欧式距离) Fig. 8 Cluster analysis of the length of three dominant species (Euclidean distance)

图选项

图 7 3种蚯蚓体宽的聚类分析(欧式距离) Fig. 7 Cluster analysis of the width of three dominant species (Euclidean distance)

图选项

徐芹和肖能文(2011)指出，蚯蚓个体的背孔始位、环带形状及所占体节数、雌孔、雄孔、受精囊孔位置与数目、受精囊形状、大小与数量等是蚯蚓物种分类的基本特征，即同一种蚯蚓的上述结构相同或相似，没有显著差异。而蚯蚓的体长、体宽和体节数在不同的蚯蚓物种或同一物种的不同地理群体可能变化较大。故本研究选取上述3个指标以量化河北廊坊地区陆栖蚯蚓优势种的形态多样性。一般来说，蚯蚓的体宽并不和体长成正比。但本研究的3种蚯蚓个体体长与体宽存在一定的正比关系，并且回归关系显著(P＜0.05)。其中以天锡杜拉蚓体长与体宽的拟合度最高(R²=0.671)，日本杜拉蚓的拟合度最低(R²=0.335)，表明天锡杜拉蚓体长与体宽的相关性最高，而日本杜拉蚓的相关度最低。而天锡杜拉蚓的体宽较大，未与同属的日本杜拉蚓聚在一起的原因可能与其在土壤环境中的生态型有关(张卫信等，2007)。同时，鉴于文献记录资料相对有限，并且记录的蚯蚓种类较少、形态分类特征不完善、分布范围较窄，未涉及本研究地区的细致比较也可能是结论出现差异的原因。天锡杜拉蚓、日本杜拉蚓与梯形流蚓个体的体长与体节数之间的拟合度均较高(R²＞0.85)，并且呈极显著正相关(P＜0.01)，这一结论与文献记录结果(徐芹，肖能文，2011)相一致，并且反映了环节动物身体的长度随体节数的增多而变大的基本发育规律。此外，3种蚯蚓在不同采集点的样本数量均存在差异，故其体长与体宽的回归曲线拟合度(R²)也不尽相同。在未来的研究中应注意进一步拓展采集区域，扩大样本数量。

我国的蚯蚓分类学研究，如王海军(2005)、孙静(2009)和孙倩(2010)均侧重于通过个体的形态结构特征完成陆栖蚯蚓的物种鉴定，进而更新研究地区的蚯蚓新物种或新记录物种的地理分布名录。而本次研究首次在系统调查廊坊地区蚯蚓物种多样性的基础之上，重点分析了该地区的优势物种个体之间形态多样性，探讨其外部表观结构差异的显著性。Donnelly等(2014)结合形态学、生物统计学及DNA条形码等方法分析并比较了广布于欧洲地区的红正蚓Lumbricus rubellus的形态结构差异与地理群体间高度的遗传多样性。尽管该项研究的研究区域、研究对象与研究方法和本次研究均有所不同，但是其形态学数据及其分析方法为本研究的顺利开展具有重要的指导意义。此外，本文填补了河北廊坊地区大型土壤动物研究的空白，并为该地区陆栖蚯蚓生态学的深入研究提供了基础数据。