长白山北坡3个森林群落主要树种种间联结性
摘要:长白山是东亚大陆在我国境内唯一具有高山冻原的山地,分析长白山北坡3个群落乔木种的种间关联性尚未见报道。调查了长白山北坡海拔700-1700m的森林群落,根据调查资料,进行了TWINSPAN分类和种间联结性分析。结果表明:(1)35个样方被分为3个群落类型,即原始阔叶红松林、天然次生林和云冷杉林。34个树种被分为五类,第一类和第二类主要为阴性或中性偏阴性的硬阔叶树种,第三类和第四类主要为阳性或中性偏阳的软阔叶树种,第五类主要为针叶树种。(2)原始阔叶红松林建群种总体关联性呈不显著负相关,天然次生林和云冷杉林呈不显著正相关。(3)原始阔叶红松林、天然次生林和云冷杉林分别有89%、76%、47%的种对间关联程度紧密,植物群落稳定性与种间关系并非趋于正相关。
Analysis of interspecific associations among major tree species in three forest communities on the north slope of Changbai Mountain
ZHANG Yue1, 2, GUO Liping2, 3, YI Xuemei1, 2, CAO Wei1, WANG Yuanxia1, 2, WU Peili4, JI Lanzhu1
Abstract:Changbai Mountain has the only mountainous tundra habitat found in China. Interspecific associations among major tree species in three forestry communities on the north slope of Changbai Mountain have not yet been analyzed. In this study, we surveyed 35 quadrats in the three forest communities from 700 m to 1700 m of altitude, classified the samples using TWINSPAN, and analyzed for possible interspecific associations. The results show that (1) the 35 quadrats were divided into 3 different forest communities, a natural broad-leaved Korean pine forest, a natural secondary forest, and a spruce-fir forest. The 34 species were divided into 5 groups. The major trees species of the first and second groups were negative or neutral partial negative hard hardwood species. The major species of the third and fourth groups were positive or neutral partial Yang soft hardwood species. The fifth group was mainly composed of coniferous species. (3) Associations among major tree species in the natural broad-leaved Korean pine forest did not show any significant negative correlations, and no significant positive correlations were observed in the natural secondary forest and spruce-fir forest. (3) The percentage of closely correlated species pairs was 89%, 76%, and 47% for the natural broad-leaved Korean pine forest, natural secondary forest, and spruce-fir forest, respectively. Plant community stability and interspecies association did not tend to be positively correlated.
Key words:forest community TWINSPAN elevation gradient interspecific associations correlation coefficient
种间联结性是对一定时期内植物群落组成物种之间相互关系的静态描述[1],是以二元数据为基础,检验两个物种是否存在关联,表达的是一种定性关系[2, 3, 4]。种间联结通常是由群落演替过程中物种对环境的趋同性及其相互作用而引起的,它不仅是群落形成、演替的基础,也是划分群落的依据及物种间相互关系的表现形式之一[5]。不同物种间关联程度的测定,对于研究各物种之间的协作关系及其对环境的适应性和在群落中的分布、种群的进化以及群落的演替规律和分类都有重要意义[6, 7, 8]。
长白山是东亚大陆在我国境内唯一具有高山冻原的山地,随海拔高度的变化,呈现出垂直分布植被带[9, 10, 11],北坡的群落组成、结构沿海拔梯度的分布差异最为明显和完整[12]。长白山的重要用材树种主要有红松(Pinus koraiensis)、水曲柳(Fraxinus mandshurica)、紫椴(Tilia amurensis)、胡桃楸(Jugland mandshurica)、红皮云杉(Picea koraiensis)、长白鱼鳞云杉(Picea jezoensis)、臭冷杉(Abies nephrolepis)等,起到保持水土、涵养水源、调节小气候等作用[13, 14]。然而,20个世纪以来长白山原始天然林经历了3个阶段的过度采伐,残存的森林大部分已演变为天然次生林[15, 16, 17]。分析现存的原始阔叶红松林、云冷杉林和天然次生林的乔木种分类特征,了解树种间的关联性至关重要。以往的研究工作曾对长白山北坡锻树阔叶红松林群落、长白山北坡暗针叶林、长白山云冷杉针阔混交林的种间关联性、群落组成结构的特点及其沿海拔梯度的变化规律、两个演替阶段乔木的种间联结性等展开[7, 18, 19]。然而,同时分析长白山北坡3个植物群落,结合TWINSPAN分类、种间联结系数研究不同群落内乔木种的种间关联性尚未见报道。通过数量分类,可以认识长白山北坡海拔700 m到1700 m的3个植物群落的分布特征,对各群落中的主要乔木种进行种间联结性分析,可以反映出群落生境差异以及物种生态位的重叠与分离,为长白山地区植被的经营管理、植被恢复和生物多样性保护提供依据。
1 研究方法1.1 野外调查与样方设置
野外调查于2010、2011年进行,研究地点在长白山北坡安图县,在海拔700—1700m之间,共设置35个大小为30m×30m的样方,每个样方分为36个5m×5m的小样方,在每个小样方中对树高≥1.3m的树木,测定胸径和树高;对树高<1.3m的乔木,记录种名及株数。同时测定各样方的经纬度、海拔、坡度、坡向和坡位。统计各样方中出现的树种数,树种出现的样方数以及种对间的出现关系,排除偶见种,选取各森林群落中重要值较高的树种进行种间关联分析。
1.2 数据分析
种群指标包括各样方的树种数、个体总数、相对频度、多度与显著度;并以重要值作为反映树种在群落中的地位和作用的指标[20]:
1.3 群落类型划分及树种分类
本研究在长白山海拔700—1700m调查了35个样方(表 1),采用双向指示种分析法划分群落类型。将35个样方中重要值<0. 05 的物种剔除,筛选出乔木种共计34种(表 2),以此构成34×35的物种重要值数据矩阵,利用采用VESPAN 软件包中Hill 1979 年设计的TWINSPAN[21]划分群落类型和树种分类。
表1 35个样方信息Table 1 The 35 samples information
样方编号Sample No.海拔Elevation/m样方编号Sample No.海拔Elevation/m样方编号Sample No.海拔Elevation/m sp1700sp13898sp251056 sp2713sp14898sp261118 sp3736sp15900sp271120 sp4749sp16901sp281124 sp5769sp17935sp291126 sp6771sp18937sp301200 sp7793sp19995sp311300 sp8800sp201000sp321400 sp9803sp211012sp331500 sp10805sp221040sp341600 sp11821sp231047sp351700 sp12826sp241049——表2 34种乔木种Table 2 The 34 tree species
编号No.乔木种名Tree species编号No.乔木种名Tree species t1白桦Betula platyphyllat18糠椴Tilia mandshurica t2东北槭Acer mandshuricumt19裂叶榆Ulmus laciniata t3茶条槭Acer ginnala t20黄花落叶松Larix olgensis t4长白鱼鳞云杉Picea jezoensis var. komaroviit21蒙古栎Quercus mongolica t5稠李Prunus padust22三花槭Acer triflorum t6臭冷杉Abies nephrolepist23千金榆Carpinus cordata t7春榆Ulmus japonicat24青楷槭Acer tegmentosum t8大青杨Populus ussuriensist25色木槭Acer mono t9枫桦Betula costatat26沙松冷杉Abies holophylla t10黑桦Betula dahuricat27櫰槐Maackia amurensis t11红皮云杉Picea koraiensis t28山荆子Malus baccata t12红松Pinus koraiensist29山梨Pyrus ussuriensis t13胡桃楸Jugland mandshuricat30山杨Populus davidiana t14花楷槭Acer ukurunduenset31水曲柳Fraxinus mandshurica t15花楸Sorbus pohuashanensist32小楷槭Acer komarovii t16黄檗Phellodendron amurenset33髭脉槭Acer barbinerve t17紫花槭Acer pseudo-sieboldianumt34紫椴Tilia amurensis1.4 种间联结测定方法1.4.1 总体相关性测定
种间总体关联性反映的是在某地出现的多物种种间是否存在显著关联性,其测定采用[22]提出的方差比率法(VR),并采用统计量W 检验VR 值偏离1 的显著程度[23]。先作零假设,即研究树种间无显著关联,选用下列公式计算检验统计量[1]:
式中,S为总的树种数,N为总样方数;Tj为样方j内出现的目标树种总数,ni树种i出现的样方数;t为样方中种的平均数。在独立假设条件下期望值为1,VR>1表示树种间表现出正关联;VR<1表示树种间存在负的净关联。使用统计量W=N×(VR) 来检验VR值偏离1的显著程度。如果树种不显著关联,W则落入由 χ2 分布给出的界限的概率为90%; χ0.952(N) <W< χ0.052(N) 。
1.4.2 种间联结性分析
为了求出种间的联结系数和测定其联结的显著性,用2×2联列表统计各种对的具体值,即a、b、c、d 的取值。联结系数AC用来说明种间联结程度的相关系数。
若ad≥bc,则:
ad<bc,且d≥a,则:
ad<bc,且d<a,则:
式中,a、b、c、d分别表示两个树种同时出现,A出现B不出现,A不出现B出现以及两树种都不出现的样方数。AC的值域为[-1, 1],越接近于1,表明树种的正联结性越强;相反,AC值越接近于-1,则表明二者的负联结性越强;AC为0,树种间完全独立。
为了克服点联结系数AC受到d值影响大而造成偏差,本文选用测定两树种间关联度较好的Ochiai指数计算种间联结度:
OI指数[1]表示种对相伴出现机率和联结性程度,当a=0时,取值为0,表示种间完全相异,不同时出现在同一样方中;当a=N(总样方数)时,取值为1,表示同时出现在样方中。
2 结果分析2.1 TWINSPAN分类2.1.1 样方分类
对长白山北坡调查的35个样方用TWINSPAN进行等级分类(图 1),结合调查区域的实际生态意义,采用第三级的划分结果,将海拔700—1700m的森林群落划分为3个群落类型,群落的命名主要依据优势种的类型以及采伐历史作为群落划分结果的命称。

(1)群落Ⅰ 天然次生林群落
该群落包含1,5—6,8—15,17—18,21—23,26—27共计18个样方,除了样方9、10外,其余16个样方均分布在海拔700—1118m。主要优势树种有白桦(Betula platyphylla)、红松、山杨(Populus davidiana)、黄花落叶松(Larix olgensis)、蒙古栎(Quercus mongolica)和青楷槭(Acer tegmentosum)等,大多数为20世纪30、40年代日本盘踞东北时皆伐后形成的白桦天然次生林。
(2)群落Ⅱ 原始阔叶红松林群落
该群落包含2—4,7,16,19—20,24—25,28—29共计11个样方,分布在海拔713—1124m的长白山原始阔叶红松林内,主要优势树种有红松、紫椴、水曲柳、蒙古栎和槭树类树种。20世纪以来群落未受到采伐干扰,保存较为完整。
(3)群落Ⅲ 云冷杉林群落
该群落包含30—35 共计6个样方,分布在海拔1200—1700m。主要优势种为长白鱼鳞云杉、臭冷杉、黄花落叶松等针叶树种。
2.1.2 树种分类
长白山北坡34个树种的TWINSPAN分类结果见表 3。
表3 长白山北坡主要乔木种TWINSPAN分类矩阵Table 3 Quantitative classification matrix of the main tree species in the north slope of Changbai mountain
111111 1122222 1122222 333333 15690123478581236747690458923012345 左边数字为树种编号(表 2) ,最上方为样方编号,最下方是样方分组结果,最右边是树种分类结果,中间数字采用“伪种”方法将各指标数值分级,“-”表示植物在该样方不存在,“1”表示一次分类结果,“0”表示另一组 T3-------------2-------------3-------000 T5-1-----21---24------1-----131------000 T9--------------------1--22--42------000 T11-----------------------1-----------000 T134433545423334324444425443453322----000 T142--------------------------21------000 T20--------------------1-3-2--32------000 T24---------------------------31------000 T27----------------------2---1--------000 T3242--3122422-2-1----43-3--3-42------000 T33-1-------1--1--1--1--4-22-1--------000 T34-------------1----2354-4445-4------000 T2111----32--------------------------001 T12------1--2--3311441----------------001 T1711--11-2-1-----1--------11--2------001 T183345512331521---------2------------001 T19323-----1---1-23-1--1----2-1-------001 T222343333233342-----423-1--4--1------001 T23122--12-2------------------13------001 T265444433454333232--13314134134------001 T28-12-2121242131----331--------------001 T29-111-222----1--1-1-----------------001 T30-----11--1-------------------------001 T31-44--34--22-22322---------1--------001 T1-44--544453-5454552---1-1-1--34----01 T7232-222221---111---24--1-1--21-2---01 T83--3-2--11--1--3-1------4-----2----01 T25-21-2-1421112133231133234222-32----01 T16-------------1----1--------1--21--210 T21---------2244-2423-----------53442410 T4-----11---3--243--3-23444-2--55555511 T6----------432244553-3423--54-55555411 T10--------------2--1---1--3---11121--11 T15----------14--32-1-3-1422--1-33452111 00000000000000000000000000000111111 00000000000000000011111111111001111 00000000001111111100000000011本研究采用5级分类水平将长白山北坡主要乔木种分为5个植被类型,以生活型和生境为依据,根据外业调查数据,将34个乔木种划分为3个类群。
(1)类群Ⅰ
包括000和001内的所有树种,主要为阴性或中性偏阴性的树种,如稠李(Prunus padus)、黑桦(Betula dahurica)、红松、胡桃楸、裂叶榆(Ulmus laciniata)、千金榆(Carpinus cordata)、沙松冷杉(Abies holophylla)、水曲柳、紫椴和槭树类树种。
(2)类群Ⅱ
包括01内的所有树种,如白桦、春榆、大青杨、青楷槭等,该组树种主要为阳性或中性偏阳的软阔叶树种。
(3)类群Ⅲ
包括10和11内的所有树种,如花楸、落叶松、长白鱼鳞云杉、臭冷杉、枫桦、花楷槭等,主要为高海拔地区的针叶树种。
根据《中国植被》,植被型是一定特殊生境的产物[24],从上述样方和树种的分类结果可知,原始阔叶红松林包括水曲柳、胡桃楸、紫椴、黄檗等硬阔叶树种,主要以红松等阴性或中性偏阴性的树种为主;而天然次生林以白桦、大青杨等阳性或中性偏阳的软阔叶树种为主。虽然原始阔叶红松林与天然次生林的海拔相似,但是受历史人为采伐因素的影响,群落内乔木种组成差异较大。云冷杉林以分布在高海拔地区的落叶松、长白鱼鳞云杉、臭冷杉等针叶树种为主。以上分类结果说明TWINSPAN 数量分类结果与实际植被空间格局分布较为一致,是一种行之有效的群落分类方法。
2.2 种间联结性2.2.1 原始阔叶红松林
重要值是综合反映群落中不同树种的相对重要性,为了客观地选出能概括各森林群落特征的主要树种,选取最小值3 作为阈值,将阔叶红松林内重要值大于3的树种(表 4)进行种间关联分析(图 2)。红松的重要值是蒙古栎的4.24倍,原始阔叶红松林内树种组成复杂、种类丰富,各树种所占资源差异不大,建群种优势度差异不显著。
表4 阔叶红松林内重要值排在前10位的树种Table 4 The top ten tree species in the broad-leaved korean pine forest
排序Rank树种Species重要值Ⅵ排序Rank树种Species重要值Ⅵ 1红松14.026色木槭7.11 2髭脉槭11.557长白鱼鳞云杉4.94 3紫花槭9.528青楷槭3.92 4紫椴8.639水曲柳3.82 5臭冷杉7.3910蒙古栎3.31
根据10 个优势种群的存在与不存在矩阵,计算VR 值( VR = ST2/ δ2T = 0.6033< 1),说明10个种群间表现出负联结。对以上测定结果,采用统计量W 检测VR偏离1的显著性, W=N×(VR) =6.64;查表得相应的 χ2 值,χ0.052(11)=19.68 ,χ0.952(11)=4.57 ,因此 W 落入 χ0.952(11)与χ20.05(11) 之间,说明阔叶红松林内主要树种总体上表现出不显著的负关联。反映出阔叶红松林群落主要乔木种对环境的适应性相同或相似,但又存在资源利用上的竞争,总体间不显著的负关联,则是利用种间一定程度的分离,以降低竞争强度。
联结系数(AC)度量种间联结性的强度,可作为对种间联结性的一种说明。阔叶红松林内AC值范围在<-0.5、[-0.5,-0.1]、0、[0.1,0.5]、>0.5、不存在的种对分别占27%、6%、18%、4%、11%、33%,正负关联比为0.45,说明51%的种对间呈独立分布。阔叶红松林内种间联结度OI指数半矩阵图(图 4),表明89%的种对(40对)的OI值在[0.6,0.1],反映出原始阔叶红松林内主要树种89%的种对间关联程度紧密。
2.2.2 天然次生林
天然次生林内选取重要值大于3的树种(表 5)进行种间关联分析(图 3)。白桦的重要值是青楷槭的4.84倍,天然次生林内树种组成简单,建群种优势度差异显著。
表5 天然次生林内重要值排在前10位的树种Table 5 The top ten tree species in the natural secondary forest
排序Rank树种Tree species重要值VI排序Rank树种Tree species重要值VI 1白桦15.556臭冷杉5.92 2红松11.767蒙古栎5.42 3色木槭10.728山杨3.98 4紫花槭7.89黄花落叶松3.66 5紫椴7.6910青楷槭3.21根据12个优势种群的存在与不存在矩阵,计算VR 值( VR = S2T / δ2T =1.35>1),说明种群间表现出负联结。对以上测定结果,采用统计量W 检测VR偏离1的显著性, W=N×(VR) =24.3;查表得相应的 χ2 值, χ20.05(18)=28.87 , χ0.952(18)=9.39 ,因此 W 落入 χ0.952(18) 与 χ20.05(18) 之间,说明天然次生林主要树种总体上表现出不显著的正关联,反映出各物种对群落环境有较一致的要求。
天然次生林内联结系数(AC)的范围在< - 0.5、[-0.5,-0.1]、[0.1,0.5]、>0.5、不存在的种对分别占18%、17%、18%、14%、33%,正负关联比为0.91,33%的种对间呈独立分布。种间联结度OI指数半矩阵图(图 3),表明76%的种对(50对)的OI值在[0.6,1],说明主要乔木树种间76%的种对关联程度紧密。

2.3.3 云冷杉林
云冷杉林内重要值大于等于3的树种有5种(表 6),高海拔地区林下树种组成简单,长白鱼鳞云杉的重要值是白桦的12.9倍,群落内建群种优势度差异很大,优势树种极为明显。
表6 云冷杉林内重要值排在前5位的树种Table 6 The top five tree species in the spruce-fir forest
排序Rank树种Tree species重要值VI 1长白鱼鳞云杉38.38 2臭冷杉23.80 3黄花落叶松12.62 4花楷槭8.92 5白桦2.98根据5个优势种群的存在与不存在矩阵,计算VR 值(VR = S2T / δ2T = 2> 1),说明11个种群间表现出正联结。对以上测定结果,采用统计量W 检测VR偏离1的显著性, W=N×(VR) =12;查表得相应的 χ2 值, χ20.05(6)=12.59 , χ0.952(6)=1.64 ,因此 W 落入 χ0.952(6) 与 χ20.05(6) 之间,说明云冷杉林内主要树种总体上表现出不显著的正关联。
为了客观的分析树种间关系,将群落内所有树种进行种间关联性分析(图 4)。云冷杉林内AC值范围在<-0.5、0、[0.1,0.5]、>0.5、不存在的种对分别占2%、31%、13%、16%、38%,正负关联比为14.5,69%的种对间呈独立分布。种间联结度OI指数半矩阵图(图 4),表明26种对(占47%)的OI值在[0.6,1],说明11个树种间47%的种对关联程度紧密。

3 结论与讨论3.1 TWINSPAN分类
本文采用双向指示种分析法对调查样地和树种进行分类,35个样地被分为3个群落,即海拔700—1120m为原始阔叶红松林和天然次生林,因为20世纪30、40年代日本盘踞东北时部分阔叶红松林被皆伐后形成的杨桦、白桦天然次生林;1120—1700m为云冷杉林,划分结果与前人的结论大致相同[10, 25, 26, 27]。
闫琰等[28]的研究表明长白山温带次生杨桦林和次生针阔混交林中,小尺度(10m×10m)中性过程和生态位过程均能解释群落物种多度分布,但生态位过程的解释能力相对较大,而在中等尺度(30m×30m)上中性过程开始逐渐取代生态位过程程成为维持物种多度分布的主要机制。本研究的调查采用中等尺度(30m×30m),原始阔叶红林内主要树种中第一、二类的树种占74%,第三、四类树种占20%;天然次生林内主要树种中第一、二类的树种占50%、第三、四类的树种占15%。历史采伐干扰改变了林分中的资源和空间状况[29],原始林与天然次生林在资源上优势树种宽度值不同,物种多度分布也不同。云冷杉林内主要树种中第三、五类的树种占72%,是由于臭冷杉、长白鱼鳞云杉具有最大的生态位宽度[30],是该区森林群落的主要建群种。生态位宽度较大物种对林内资源环境的适应力较强,占有的空间范围较广,在群落内的重要值也较大。应用生态位理论可以较好的解释中等尺度研究的树种分类结果。
3.2 种间关联性分析
从群落中乔木种的重要值看,随着海拔高度的升高,建群种数量显著减少,云冷杉林建群种数量是原始林和次生林的一半。原始林和次生林建群种数量虽然相同,但受历史采伐的影响,次生林内白桦的重要值明显高于红松,建群种之间的优势度差异较大。从群落内种间总体关联性结果看,在较大的尺度下(样方单位面积30 m×30 m),原始阔叶红松林的建群种总体呈不显著负关联,天然次生林建群种总体呈不显著正关联。一般认为,物种总体联结随着演替的发展变化反映出次生演替各阶段群落的稳定性,特别是演替发展至后期总体上往往表现出显著正相关[31]。本研究中原始阔叶红松林未受采伐等人为影响,是处于演替后期的顶级群落;天然次生林伐后恢复近半个世纪,是处于演替中期的群落,群落内建群种总体关联性的研究结果不支持上述观点。林木间竞争是森林生态系统中的普遍现象[32],建群种是对群落结构和群落环境的形成有明显控制作用的物种,原始阔叶红松林经过长期的演替,建群种间的竞争使得该生态系统达到平衡状态。
从种间联结性结果看,3个群落内种对之间多呈正联结关系,树种间正联结是生态学特性及对环境要求的差异和互补性,导致各种之间相互兼容,互相促进生长[32],常常相伴出现。原始阔叶红松林内,黄檗与色木槭、水曲柳都为主要伴生树种,在群落中占有重要的优势度,表现出显著的正联结关系。天然次生林中,臭冷杉与落叶松都喜湿生生境,为主要伴生种。云冷杉林中长白鱼鳞云杉与臭冷杉、红松的生态习性也较为一致,都为阴性、喜冷湿环境及酸性土壤的针叶树种。杜志等[13]研究发现从云冷杉针阔混交林到云冷杉针叶混交林和近原始云冷杉林,云杉、冷杉和红松两两间的空间关联性,表现为负相关的现象增多。本研究中的云冷杉林分布在海拔1200—1700m范围内,研究范围较广,云杉与冷杉相伴出现的几率较大;随着海拔升高,红松的数量减少,因而云杉和冷杉分别于红松相伴出现的几率降低。3个群落中,绝大多数种对的联结强度存在明显的差异,因为树种对微环境适应性不同,导致它们在不同群落中存在几率不同[3, 18, 33]。只有臭冷杉与青楷槭在阔叶红松林与天然次生林中,白桦与青楷槭在云冷杉林和天然次生林中联结性均较高。了解不同群落不同演替阶段的种间联结性和联结强度,可以为森林演替不同阶段的经营管理提供理论依据。
相关知识
九龙山南方铁杉群落物种多样性及乔木种种间联结性
长白山北坡林线对气候变化的响应过程研究
白花鬼针草入侵对本土植物群落主要物种种间联结和稳定性的影响
长白山北坡珍稀濒危植物及其自然保护
长白山 北半球幸存的温带森林
长白山哈泥泥炭地七种苔藓植物种间联结和生态位研究
长白山阔叶红松林主要树种种群结构及其林木径向生长对气候响应
长白山北坡访花昆虫研究(Ⅲ)——访花昆虫多样性
冀北山区混交林植被生态位和种间联结分析.pdf文档全文免费阅读、在线看
长白山北坡春季物候特征及其对气候变化的响应
网址: 长白山北坡3个森林群落主要树种种间联结性 https://www.huajiangbk.com/newsview766832.html
上一篇: 昆明植物所在群落水平植物花部性状 |
下一篇: 宁夏荒漠草原典型群落的植被特征及 |
推荐分享

- 1君子兰什么品种最名贵 十大名 4012
- 2世界上最名贵的10种兰花图片 3364
- 3花圈挽联怎么写? 3286
- 4迷信说家里不能放假花 家里摆 1878
- 5香山红叶什么时候红 1493
- 6花的意思,花的解释,花的拼音 1210
- 7教师节送什么花最合适 1167
- 8勿忘我花图片 1103
- 9橄榄枝的象征意义 1093
- 10洛阳的市花 1039