不同植被類型植物物種多樣性

袁王俊, 盧訓(xùn)令, 張維瑞, 王智慧, 張靈敏, 丁圣彥

1 河南大學(xué)環(huán)境與規(guī)劃學(xué)院, 開封 475004 2 河南大學(xué)中藥研究所, 開封 475004

不同植被類型植物物種多樣性

袁王俊1,2, 盧訓(xùn)令1, 張維瑞2, 王智慧2, 張靈敏2, 丁圣彥1，2，*

1 河南大學(xué)環(huán)境與規(guī)劃學(xué)院, 開封 475004 2 河南大學(xué)中藥研究所, 開封 475004

為了評(píng)價(jià)大面積人工種植杜仲對(duì)當(dāng)?shù)刂参锒鄻有缘挠绊?，以河南省汝?yáng)縣不同植被類型為研究對(duì)象，通過(guò)群落學(xué)調(diào)查，運(yùn)用重要值、Shannon-Wiener物種多樣性指數(shù)(H)、Simpson物種多樣性指數(shù)(D)和均勻度指數(shù)(JH′)等指標(biāo)，統(tǒng)計(jì)分析杜仲人工林、溫帶落葉闊葉林(以下簡(jiǎn)稱次生林)和撂荒地3 種不同植被類型的植物物種多樣性，探討杜仲種植對(duì)植物多樣性的影響。調(diào)查發(fā)現(xiàn)，杜仲林樣地中出現(xiàn)植物 82種，隸屬39科63屬，草本層為最發(fā)達(dá)的一層；次生林樣地中出現(xiàn)植物70 種，隸屬32科62 屬，喬木層為最發(fā)達(dá)層；撂荒地樣地中出現(xiàn)植物84種，隸屬35科69屬，無(wú)喬木層。杜仲林物種豐富度和多度均不亞于次生林和撂荒地，且杜仲林草本層物種豐富度和植株總數(shù)均高于次生林。統(tǒng)計(jì)分析顯示，杜仲林喬木層、灌木層和草本層多樣性指數(shù)H和D值、均勻度指數(shù)JH′值與其它兩種植被類型相比無(wú)顯著性差異(P>0.05)。因此，種植杜仲過(guò)程中采用合理密度，適當(dāng)管理，不僅能提供葉、花、果等資源，而且能夠豐富草本植物的種類和數(shù)量，增加植物物種多樣性。

杜仲林; 次生林; 撂荒地; 物種多樣性; 群落結(jié)構(gòu)

杜仲(EucommiaulmoidesOliv.)為杜仲科(Eucommiaceae)杜仲屬植物，是我國(guó)特有的第三紀(jì)孑遺植物[1]。杜仲皮和葉具有補(bǔ)肝腎、強(qiáng)筋骨功效[2-3]，種仁油中α-亞麻酸含量居植物之首[4-5]。杜仲膠具有獨(dú)特的橡膠和塑料二重性。2011年3月中國(guó)工業(yè)協(xié)會(huì)牽頭組織相關(guān)企業(yè)、科研單位組建了杜仲產(chǎn)業(yè)技術(shù)創(chuàng)新戰(zhàn)略聯(lián)盟[6]；2013年9月我國(guó)出版了杜仲產(chǎn)業(yè)綠皮書《中國(guó)杜仲橡膠資源與產(chǎn)業(yè)發(fā)展報(bào)告》，以推動(dòng)杜仲產(chǎn)業(yè)良性發(fā)展[1]。我國(guó)杜仲種植面積占全球90 %以上，栽培區(qū)域南到廣西、北達(dá)吉林，冬自上海，西抵新疆喀什[1]，為了滿足醫(yī)藥和橡膠產(chǎn)業(yè)對(duì)杜仲資源的需求，各地仍在大量栽培杜仲。

生態(tài)系統(tǒng)的物種多樣性和結(jié)構(gòu)特征是生態(tài)學(xué)研究中的熱點(diǎn)問(wèn)題之一[7-8]。物種多樣性是生態(tài)系統(tǒng)功能的重要基礎(chǔ)，生態(tài)系統(tǒng)中物種多樣性豐富，不僅有利于提高系統(tǒng)本身抵御外界不確定性因素的干擾和干擾后系統(tǒng)恢復(fù)能力，而且還有助于生態(tài)系統(tǒng)中營(yíng)養(yǎng)元素和肥力的可持續(xù)利用，生產(chǎn)力水平隨著物種多樣性的增加而提高[9-12]。人工栽培林木可為當(dāng)?shù)貛?lái)經(jīng)濟(jì)和社會(huì)效益。但在種植過(guò)程中需要進(jìn)行的林地清理、基礎(chǔ)設(shè)施建設(shè)、林地管理等，這些措施對(duì)當(dāng)?shù)厣锒鄻有缘挠绊懷芯肯鄬?duì)較少，且主要研究對(duì)象為桉樹，一些研究認(rèn)為桉樹人工林會(huì)降低當(dāng)?shù)厣锒鄻有?、土壤肥力和水分的作用[13-15]，另有研究認(rèn)為人工林能夠改善林地小氣候、提高土壤肥力、保護(hù)其林下生物多樣性[16-18]。

目前，杜仲人工林的研究主要集中在造林技術(shù)，種植技術(shù)，環(huán)境因子對(duì)生物量的影響等方面[19-20]，而人工栽培杜仲對(duì)當(dāng)?shù)刂参镂锓N多樣性的影響尚缺乏研究。本文擬以‘全國(guó)杜仲之鄉(xiāng)’汝陽(yáng)多年生杜仲林、次生林和撂荒地為研究對(duì)象，通過(guò)樣地調(diào)查和比較，分析人工杜仲林的對(duì)當(dāng)?shù)刂参锒鄻有缘挠绊?，為今后杜仲林的發(fā)展提供理論依據(jù)。

1 研究區(qū)概況

汝陽(yáng)縣位于河南省西部，地處東經(jīng)112°8′—112°38′，北緯33°49′—34°21′。最高山峰雞冠山海拔1602 m，最低點(diǎn)杜康河底海拔220 m，平均海拔543 m。屬暖溫帶大陸性季風(fēng)氣候，光照充足，氣候溫和，四季分明，年平均日照時(shí)數(shù)2177.3 h，日照百分率達(dá)49 %，年平均氣溫14 ℃，年均降雨量690 mm。森林覆蓋率為44.7 %。1987年被林業(yè)部批準(zhǔn)為杜仲基地縣。2000年被中國(guó)經(jīng)濟(jì)林協(xié)會(huì)命名為“全國(guó)杜仲之鄉(xiāng)”。目前，全縣已發(fā)展杜仲1.5萬(wàn) hm2，占全縣經(jīng)濟(jì)林總面積的58.4 %，杜仲皮、葉及系列產(chǎn)品年產(chǎn)量達(dá)1560萬(wàn)kg，年總產(chǎn)值1.2億元。

2 研究方法

2.1 樣地設(shè)置

2013年7月分別在種植杜仲較多的3個(gè)鄉(xiāng)鎮(zhèn)(王坪鄉(xiāng)、靳村鄉(xiāng)和三屯鄉(xiāng))各設(shè)置1個(gè)樣地(均于20世紀(jì)90年代初種植杜仲)，每個(gè)樣地杜仲林、次生林及撂荒地中各設(shè)置樣方3 個(gè)，共調(diào)查9個(gè)樣地，27個(gè)大樣方。喬木層樣方10 m × 10 m，在喬木樣方內(nèi)沿著樣地對(duì)角線分別設(shè)置5 個(gè)3 m × 3 m 灌木樣方，再在灌木樣方中央設(shè)置1 m × 1 m 草本樣方，共設(shè)置灌木和草本樣方各135個(gè)。

2.2 群落調(diào)查

調(diào)查指標(biāo)包括：?jiǎn)棠緦拥姆N類、數(shù)量、樹高、胸徑、基徑、冠幅；灌木層的種類、數(shù)量、高度、基徑、冠幅；草本層的種類、數(shù)量、高度、叢徑、蓋度等指標(biāo)；同時(shí)，記錄每個(gè)樣方地理坐標(biāo)、海拔、坡向、坡度、坡位等立地因子，樣方環(huán)境條件見表 1(僅列出杜仲林植被樣方點(diǎn)位信息，其他點(diǎn)位情況與其基本相同)。

表1 杜仲林樣地基本情況

2.3 數(shù)據(jù)分析

(1)物種豐富度

此次野外調(diào)查中各樣方面積是一致的，使用樣方內(nèi)物種數(shù)目S表示物種豐富度。

(2)重要值

喬木層重要值IV=(相對(duì)密度+相對(duì)顯著度+相對(duì)頻度)/3

灌木層重要值IV=(相對(duì)密度+相對(duì)高度+相對(duì)頻度)/3

草本層重要值IV=(相對(duì)密度+相對(duì)蓋度+相對(duì)頻度)/3

(3)Shannon-Wiener多樣性指數(shù) (H)

式中,Pi為第i物種的重要值，即Pi=ni/N，S為物種數(shù)。

(4)Simpson多樣性指數(shù)(D)

式中，Ni為第i物種的個(gè)體數(shù)，N為所有物種的個(gè)體總數(shù)，S為物種數(shù)。

(5)均勻度指數(shù)

本文選用基于Shannon -Wiener多樣性指數(shù)的均勻度指數(shù)進(jìn)行研究：

使用SPSS11.5軟件包中的單因素方差分析(α=0.05)檢驗(yàn)群落物種多樣性及其關(guān)系，所有圖表均在Excel 2007中完成。

3 結(jié)果與分析

3.1 物種豐富度

匯總3 個(gè)樣地3 種植被類型的喬木層、灌木層和草本層的科、屬及物種豐富度，結(jié)果如表2所示。杜仲林樣方內(nèi)共有植物82 種，隸屬于39 科63 屬；次生林樣方內(nèi)出現(xiàn)植物70 種，隸屬32 科62 屬；撂荒地樣方內(nèi)共有植物84 種，隸屬于35 科69 屬，次生林物種豐富度最小。杜仲林草本層58 種，2657 株，科、種和株數(shù)均大于喬木層和灌木層，表明杜仲林中草本層為最發(fā)達(dá)層。次生林喬木層有5個(gè)種，分別來(lái)自5 科5 屬，最主要的種是栓皮櫟，草本層50 種，960 株，其科、屬、種類仍然最大。撂荒地?zé)o喬木層，灌木層植物種類和數(shù)量都較少，而草本層極為發(fā)達(dá)，76 種，2469 株，其科、屬、種數(shù)都在3種樣地中是最大的。總體而言，在整個(gè)群落層面上，杜仲林的各層結(jié)構(gòu)與次生林相比并不簡(jiǎn)單，相反，草本層物種豐富度和植株總數(shù)都高于次生林。

表2 3種植被類型內(nèi)植物的科屬種數(shù)

3.2 群落物種重要值

喬木層共15 種，杜仲181 株，其重要值為70.73，其余種類如油松Pinustabulaeformis、栓皮櫟Quercusvariabilis、豆梨Pyruscalleryana等零星存在，重要值介于1.56—3.75之間。次生林喬木層共5 種，其中栓皮櫟130 株，其植株高大，最高者達(dá)23 m，為建群種，其重要值為81.01。杜仲林灌木層中鹽膚木Rhuschinensis35 株，重要值為44.96，牡荊Vitexnegundo重要值為11.86。重要值最小的綠葉胡枝子Lespedezabuergeri僅2.64。次生林灌木層中有黃櫨Cotinuscoggygria47 株，重要值為26.54，其次為牡荊和連翹Forsythiasuspensa，重要值為18.78和13.1。撂荒地灌木層有酸棗ZiziphusjujubaMill. Var.spinosa16 株，重要值為36.77，其次為牡荊，重要值為19.4。杜仲林草本層中野艾蒿Artemisialavandulaefolia530 株，重要值為17.63，細(xì)弱隱子草Cleistogenesgracilis344 株，重要值為8.53，鵝觀草Roegneriakamoji325 株，重要值為7.34。次生林草本層中細(xì)葉苔草Carexduriuscula144 株，重要值18.76，野艾蒿109 株，重要值10.36。撂荒地草本層中野艾蒿318 株，其重要值為12.42，其次為狗尾草Setariaviridis，重要值為7.91(表3)。

表3 3種植被類型主要物種組成及其重要值

3.3 物種多樣性和均勻度

本研究分別計(jì)算了不同植被類型喬木層、灌木層和草本層的Shannon-Wiener物種多樣性指數(shù)(H)和Simpson多樣性指數(shù)(D)。由表4可見，杜仲林喬木層多樣性指數(shù)H為1.361次生林多樣性指數(shù)H為0.705，二者無(wú)顯著性差異。二者多樣性指數(shù)D為 0.192和0.192，同樣不存在顯著性差異。灌木層中次生林的多樣性指數(shù)H為2.24，杜仲林和次生林分別為1.858和1.789，同樣沒(méi)有顯著性差異。多樣性指數(shù)D次生林最為0.847，次生林其次為0.75，杜仲林為0.549。草本層多樣性指數(shù)H撂荒地為3.635，杜仲林為3.263，次生林最低，為3.151，三者無(wú)顯著差異。多樣性指數(shù)D撂荒地、次生林、杜仲林依次為0.945、0.931和0.912，三者也無(wú)顯著差異。

由表4還可以看到杜仲林喬木層的均勻度指數(shù)(JH′)為0.503，稍大于次生林的0.438，杜仲林灌木層灌木層均勻度指數(shù)為0.807，次生林灌木層均勻度為0.808，二者基本一致，杜仲林草本層均勻度指數(shù)為0.804，次生林草本層均勻度為0.805，二者無(wú)差異。撂荒地灌木層和草本層均勻度指數(shù)分別為0.86和0.839，與杜仲林和次生林也無(wú)顯著性差異。

表4 3種植被類型物種多樣性指數(shù)及均勻度指數(shù)

4 討論

本研究9個(gè)杜仲林樣方共有植物82 種，杜仲林喬木層有15 種植物，并非單一樹種，原因是王坪鄉(xiāng)和靳村鄉(xiāng)是在原有山地基礎(chǔ)上栽培杜仲，林地清理不徹底，總的來(lái)說(shuō)林相整齊，喬木層一般高度在3—8 m, 和灌木層高度差別明顯。灌木層重要值較高的依次為鹽膚木、牡荊、黃櫨和連翹，草本層重要值較高是野艾蒿、隱子草和鵝觀草，草本層植物種類和株數(shù)都較多，蓋度大。次生林喬木層中栓皮櫟為建群種，植株蓋度大，林相整齊，灌木層16 種，草本層共50 種植物，僅960 株。撂荒地灌木層僅8 種，共35 株，高度整體較低，灌木層和草本層分界不是很明顯，而草本層極為發(fā)達(dá)，有76 種，蓋度高。總的來(lái)說(shuō)，杜仲林和撂荒地草本層種類和數(shù)量都較大。Chen等認(rèn)為草本植物多屬于逆境耐受型，尤其是一年生草本植物，主要靠種子繁殖，其特點(diǎn)是生活史短、種子多而繁殖力強(qiáng)，最能使用持久和強(qiáng)烈干擾[21]。因此，草本植物具有更廣泛的適應(yīng)性，對(duì)整個(gè)生態(tài)系統(tǒng)而言，它起到一種基礎(chǔ)保障作用。Grime認(rèn)為草本植物是生態(tài)系統(tǒng)的先鋒植被，也是退化生態(tài)系統(tǒng)的保留植被，適應(yīng)力強(qiáng)，有較寬的生態(tài)位。由于其生活史短暫，且能及時(shí)調(diào)整以應(yīng)對(duì)不同生活環(huán)境的變化，因而群落物種組成具有不穩(wěn)定性和多變性[22],本文結(jié)果與他們的觀點(diǎn)一致。次生林草本層種類和數(shù)量較少，糾其原因是栓皮櫟蓋度高，造成林分郁密度較高，使得光照條件不足，不利于種子的萌發(fā)和生長(zhǎng)，抑制了種子庫(kù)中非耐陰性草本先鋒物種種子的入侵和萌發(fā)。

分析Shannon-Wiener物種多樣性指數(shù)(H)和Simpson多樣性指數(shù)(D)，杜仲林的喬木層、灌木層和草本層與次生林、撂荒地均無(wú)顯著差異。在相同立地條件下，生境的差異是引起群落植物多樣性不同的主要原因之一。所處的生境條件存在較大差異時(shí)，群落就向不同的方向演替，繼而結(jié)構(gòu)和功能發(fā)生變化[23-24]。本研究海拔變化不大，對(duì)植物多樣性無(wú)影響，主要差異是種植了杜仲，而種植杜仲后沒(méi)有降低當(dāng)?shù)厣锒鄻有裕炊捎诙胖偕w度較小，增加了植物的種類和數(shù)量，尤其是草本植物，更有利于降低水土流失。同時(shí)，在調(diào)查中發(fā)現(xiàn)有些杜仲林由于精耕細(xì)作，杜仲生長(zhǎng)量大，林分郁密度極高，地面幾乎寸草不生，嚴(yán)重影響了當(dāng)?shù)刂参镂锓N多樣性。

綜上所述，比較杜仲林、次生林和撂荒地3種不同的植被類型，合理的種植密度下杜仲林植物物種多樣性并沒(méi)有減少，而且草本植物種類和數(shù)量相對(duì)于次生林有所增加，杜仲林不僅能夠提供葉、花、果等資源，并能夠增加植物物種多樣性。

[1] 杜紅巖, 胡文臻, 俞銳. 杜仲產(chǎn)業(yè)綠皮書—中國(guó)杜仲橡膠資源與產(chǎn)業(yè)發(fā)展報(bào)告. 北京: 社會(huì)科學(xué)文獻(xiàn)出版社, 2013: 1-5.

[2] 國(guó)家藥典委員會(huì). 中國(guó)藥典(一部). 北京: 中國(guó)醫(yī)藥科技出版社, 2010: 154-155.

[3] Zhang Q, Su Y Q, Zhang J F. Seasonal difference in antioxidant capacity and active compounds contents ofEucommiaulmoidesoliver leaf. Molecules, 2013, 18(2): 1857-1868.

[4] 杜紅巖. 杜仲活性成分與藥理研究的新進(jìn)展. 經(jīng)濟(jì)林研究, 2003, 21(2): 58-61, 82-82.

[5] 李欽, 杜紅巖, 杜蘭英, 赫錦錦, 段經(jīng)華. HPLC法測(cè)定杜仲雄花和杜仲雄花茶中京尼平苷酸、綠原酸和京尼平苷. 中草藥, 2009, 40(1): 71-72.

[6] 范仁德. 堅(jiān)持不懈, 繼續(xù)大力推進(jìn)杜仲膠產(chǎn)業(yè)化發(fā)展——在“2012中國(guó)杜仲產(chǎn)業(yè)發(fā)展高峰論壇”上的講話. 中國(guó)橡膠, 2012, 28(19): 4-6.

[7] Tilman D, Reich P B, Knops J M H. Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature, 2006, 441(7093): 629-632.

[8] Ives A R, Carpenter S R. Stability and diversity of ecosystems. Science, 2007, 317(5834): 58-62.

[9] Naeem S, Thompson L J, Lawler S P, Lawton J H, Woodfin R M. Declining biodiversity can alter the performance of ecosystems. Nature, 1994, 368(6473): 734-737.

[10] Tilman D, May R M, Lehman C L, Nowak M A. Habitat destruction and the extinction debt. Nature, 1994, 371(6492): 65-66.

[11] Kennedy T A, Naeem S, Howe K M, Knops J M H, Tilman D, Reich P. Biodiversity as a barrier to ecological invasion. Nature, 2002, 417(6889): 636-638.

[12] Wang C T, Long R J, Ding L M, Lai D Z, Li Y F. Species diversity, community stability and ecosystem function -extension of the continuous views. Pratacultural Science, 2005, 22(6): 1-7.

[13] Smith O H, Petersen G W, Needelman B A. Environmental indicators of agroecosystems. Advances in Agronomy, 2000, 69: 75-97.

[14] Wall A, Heiskanen J. Water-retention characteristics and related physical properties of soil on afforested agricultural land in Finland. Forest Ecology and Management, 2003, 186(1-3): 21-32.

[15] Sicardi M, García-Préchac F, Frioni L. Soil microbial indicators sensitive to land use conversion from pastures to commercialEucalyptusgrandis(Hill ex Maiden) plantations in Uruguay. Applied Soil Ecology, 2004, 27(2): 125-133.

[16] Keenan R, Lamb D, Woldring O, Irvine T, Jensen R. Restoration of plant biodiversity beneath tropical tree plantations in Northern Australia. Forest Ecology and Management, 1997, 99(1-2): 117-131.

[17] Manning A D, Fischer J, Lindenmayer D B. Scattered trees are keystone structures-implications for conservation. Biological Conservation, 2006, 132(3): 311-321.

[18] Zhang D J, Zhang J, Yang W Q, Wu F Z. Effects of afforestation withEucalyptusgrandison soil physicochemical and microbiological properties. Soil Research, 2012, 50(2): 167-176.

[19] 劉淑明, 梁宗鎖, 董娟娥. 不同水分條件下皮葉兩用杜仲林的生長(zhǎng)效應(yīng). 中南林業(yè)科技大學(xué)學(xué)報(bào), 2007, 27(5): 49-53.

[20] 聶朝俊, 彭智堅(jiān). 喀斯特山地人工杜仲林枯落物和土壤持水特性初步研究. 林業(yè)資源管理, 2003, (6): 27-31.

[21] 陳杰, 郭屹立, 盧訓(xùn)令, 丁圣彥, 蘇思, 郭靜靜, 李乾璽. 伊洛河流域草本植物群落物種多樣性. 生態(tài)學(xué)報(bào), 2012, 32(10): 302

[22] Grime J P. Plant Strategies, Vegetation Processes, and Ecosystem Properties. 2nd ed. USA: Wiley, 2006.

[23] Wang D P, Ji S Y, Chen F P. A review on the species diversity of plant communityal superiority. Chinese Journal of Ecology, 2001, 20(4): 55-60.

[24] 劉平, 秦晶, 劉建昌, 王華鋒, 王效科. 桉樹人工林物種多樣性變化特征. 生態(tài)學(xué)報(bào), 2011, 31(8): 2227-2235.

Plant′s diversity of different vegetation types

YUAN Wangjun1，2, LU Xunling1, ZHANG Weirui2, WANG Zhihui2, ZHANG Lingmin2, DING Shengyan1，2,*

1CollegeofEnvironmentandPlanning,HenanUniversity,Kaifeng475004,China2InstituteofChineseMateriaMedica,HenanUniversity,Kaifeng475004,China

Human-driven ecosystem simplification has highlighted questions about how the number of species in an ecosystem influences its functioning. Biodiversity is now known to affect ecosystem productivity and stability of the ecosystem. In this study, different vegetation types in Ruyang County, Henan Province, China, were examined to evaluate the impact of large-scale planting ofEucommiaulmoideson local species diversity. The species diversity was analyzed using Shannon-Wiener diversity (H), Simpson diversity (D), and evenness (JH′) indices in three different vegetation types:Eucommiaulmoidesplantation, temperate fallen broad-leaved forest (referred to as secondary forest), and abandoned land (old field). The results indicated that there were 82 species belonging to 39 families and 63 genera in theEucommiaulmoidesforest, in which the herbaceous layer was the most developed. The secondary forest contained 70 species belonging to 32 families and 62 genera, and the tree layer was more developed than more developed than the herbaceous and shrub layers. The old field included 84 species of 35 families and 69 genera but lacked a tree layer. The richness and abundance of species in theEucommiaulmoidesforest were similar to those in the secondary forest and old field, whereas the richness and abundance of species in the herbaceous layer of theEucommiaulmoidesforest were higher than those in the secondary forest. Our results showed that theH,D, andJHindices of the three layers of theEucommiaulmoidesforest were not significantly different from those of the secondary forest and old field (P> 0.05). Collectively, these results illustrate that, with a reasonable planting density and proper management,Eucommiaulmoidesforests can not only provide us with leaves and fruit but can also enrich the amount and diversity of herbaceous plants, thereby increasing biological diversity.

Eucommiaulmoidesforest; secondary forest; old field; species diversity; community structure

中國(guó)博士后基金項(xiàng)目(2013M531668); 河南省教育廳自然科學(xué)重點(diǎn)項(xiàng)目(13A180057)

2013-12-03;

2014-09-09

10.5846/stxb201312032870

*通訊作者Corresponding author.E-mail: syding@henu.edu.cn

袁王俊, 盧訓(xùn)令, 張維瑞, 王智慧, 張靈敏, 丁圣彥.不同植被類型植物物種多樣性.生態(tài)學(xué)報(bào),2015,35(14):4651-4657.

Yuan W J, Lu X L, Zhang W R, Wang Z H, Zhang L M, Ding S Y.Plant′s diversity of different vegetation types.Acta Ecologica Sinica,2015,35(14):4651-4657.