• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    The influence of canopy-layer composition on understory plant diversity in southern temperate forests

    2017-07-08 02:27:31LucianaMestrenicaToroManrquezRosinaSolerAlejandroHuertasHerreraGuillermoMartnezPasturandMarVanessaLencinas
    Forest Ecosystems 2017年2期

    Luciana Mestre,Mónica Toro-Manríquez,Rosina Soler,Alejandro Huertas-Herrera,Guillermo Martínez-Pastur and María Vanessa Lencinas

    The influence of canopy-layer composition on understory plant diversity in southern temperate forests

    Luciana Mestre*,Mónica Toro-Manríquez,Rosina Soler,Alejandro Huertas-Herrera,Guillermo Martínez-Pastur and María Vanessa Lencinas

    Background:Understory plants represents the largest component of biodiversity in most forest ecosystems and plays a key role in forest functioning.Despite their importance,the influence of overstory-layer composition on understory plant diversity is relatively poorly understood within deciduous-evergreen broadleaved mixed forests. The aim of this work was to evaluate how tree overstory-layer composition influences on understory-layer diversity in three forest types(monospecific deciduous Nothofagus pumilio(Np),monospecific evergreen Nothofagus betuloides(Nb),and mixed N.pumilio-N.betuloides(M)forests),comparing also between two geographical locations (coast and mountain)to estimate differences at landscape level.

    Results:We recorded 46 plant species:4 ferns,12 monocots,and 30 dicots.Canopy-layer composition influences the herb-layer structure and diversity in two different ways:while mixed forests have greater similarity to evergreen forests in the understory structural features,deciduous and mixed were similar in terms of the specific composition of plant assemblage.Deciduous pure stands were the most diverse,meanwhile evergreen stands were least diverse. Lack of exclusive species of mixed forest could represent a transition where evergreen and deciduous communities meet and integrate.Moreover,landscape has a major influence on the structure,diversity and richness of understory vegetation of pure and mixed forests likely associated to the magnitude and frequency of natural disturbances,where mountain forest not only had highest herb-layer diversity but also more exclusive species.

    Conclusions:Our study suggests that mixed Nothofagus forest supports coexistence of both pure deciduous and pure evergreen understory plant species and different assemblages in coastal and mountain sites.Maintaining the mixture of canopy patch types within mixed stands will be important for conserving the natural patterns of understory plant composition in southern beech mixed forests.

    Landscape,Understory communities,Understory structure,Forest structure,Nothofagus,Conservation

    Background

    The understory is a critical feature of forest ecosystems, which affects energy flow and nutrient cycling,biodiversity and regeneration capacity(Gilliam 2007).Also,understory responds fast to both natural and anthropogenic disturbances(Lencinas et al.2011),e.g.avoiding erosion and generating suitable microenvironments for other species development(Simonson et al.2014).However,understory plants have heterogeneous composition,structure and distribution patterns,which depend on the overstory species,forest structure,micro-environmental and stand conditions(e.g.,Márialigeti et al.2016;Tinya and ódor 2016).Moreover,understory could provide habitat and resource specificity to sustain rare or exclusive species in the herbaceous layer,which have excellent potential for use as biodiversity indicators(Spyreas and Matthews 2006). Thus,several forest understory species are well-known indicators of site conditions,impact of forest management and effects of climate change(Gracia et al.2007;Chavez and Macdonald 2012),and could act as a useful indicatorof forest sustainability,ecosystem health and conservation status(Dale et al.2002;Lencinas et al.2008a;Chavez and Macdonald 2012).

    Temperate forest is a worldwide distributed biome that usually exhibit strong seasonality,favoring those plant and animal species that were able to adapt to short-term climatic variation(Dreiss and Volin 2013).In temperate deciduous forests,trees drop their leaves in autumn allowing high seasonal variation in light availability for the understory;while in evergreen broadleaved forests-more commonly found in the Southern Hemisphere-trees keep their leaves year round,therefore understory plants are adapted to low light availability in all seasons.In both broadleaved forest types,understory plants have adapted to maximize their light/energy capture in different ways.Kudo et al.(2008)found that some understory plants inhabiting deciduous stands complete most of their annual growth cycle in the few weeks between snowmelt and the closure of the tree canopy.Meanwhile shaded plants inhabiting evergreen broadleaved forests tend to allocate more growth to leaves,displaying larger,thinner leaves which take up more horizontal area for light capture(Damascos 1996; Promis et al.2008).Even the most shade-intolerant species maydisappearcompletelyin pureevergreen stands developing a forest typical climax shade-tolerant species. However,there is little information about particularities of understory composition,structure and distribution patterns in mixed deciduous-evergreen broadleaved forests (see Kuninaga et al.2015).

    Magellanic forests,one of the most pristine areas in the world(Mittermeier et al.2003),are conformed mainly by Nothofagus tree species,of which deciduous N.pumilio had a wide distribution range and dominates inner the main island of the Archipelago of Tierra del Fuego,while evergreen N.betuloides develops in softer environmental conditions,as in a middle altitude on the mountain valleys,or next to the Beagle Channel and big lake coasts (Allué et al.2010).In a narrow area of their distribution range in the archipelago,N.betuloides grows together with N.pumilio forming a mixed deciduous-evergreen forest(Tuhkanen 1992),in the two both predominant landscapes where N.betuloides occurred(coast and mountains).Mixed forests of N.pumilio and N.betuloides also occurred in mild climate zones of the continent,for which the distribution and dominance of vascular plant species in the understory has been barely described (Veblen et al.1978,1979;Cruz et al.2007).Contrary,in the southernmost island conditions,understory composition and structure of mixed deciduous-evergreen Nothofagus forests have been little explored(Promis et al. 2008)compared to pure Nothofagus forests.Many researches in N.pumilio were developed studying natural unmanaged and several anthropogenic impacted forests (e.g.,Martínez Pastur et al.2002;Lencinas et al.2008a, 2008b,2011),but pure N.betuloides forests have received fewer attention.Both pure forest types usually develop a sparse and short understory layer(less than 1 m tall), mainly conformed by herbaceous and moss species (Moore 1983;Matteri and Schiavone 2002).But differences and similarities in understory composition and structure among pure and mixed Nothofagus forests are still unknown.Therefore,in this work we evaluate the influence of canopy-layer dominance(deciduous,evergreen, and mixed deciduous-evergreen)from old-growth Nothofagus forests on understory structure(frequency,cover, sub-layers),and vascular plant composition(richness and pattern assemblage),at two typical landscape locations (coast and mountain)of their natural distribution on Tierra del Fuego Island.

    We hypothesize that:(i)the composition and structure of understory in mixed forests includes characteristics of both pure forest types,being structurally similar to any of these two,and reaching greater species richness in mixed compared to pure forests;(ii)differences and similarities among pure and mixed forest are maintained in different landscape locations(coast-mountain)although the specific composition of vascular plants could vary between such locations;(iii)pure forests had high number of exclusive plant species(e.g.,shade tolerant species)compared to mixed forest,which can be considered indicator species of each canopy composition(deciduous or evergreen)and landscape locations(coastal or mountain sites).Understanding the influence of tree canopy-layer composition on the distribution of understory plants is important for biological conservation in these southern-beech forests. Moreover,in a possible scenario of global climate change mixed-species forests have been put forward as important risk reduction and adaptation strategies for environmental stress conditions(Forrester et al.2016).

    Method

    Study site description

    Native forests in Tierra del Fuego Island occupy more than 30%surface from inner mountains to Beagle Channel coasts(Allué et al.2010).At this latitude,these forests grow naturally under extreme environmental conditions as the climate is strongly influenced by oceanic current, Antarctica proximity but also the Andes Mountains (Kreps et al.2012).Climate is characterized by short,cool summers and long,snowy winters with frequent occurrence of frosts.Only three months per year have daily temperature above 0°C,and the growing season extends approximately from November to March.Mean monthly temperatures vary from–7°C to 14°C(extremes vary from–17°C in July to 22°C in January)and precipitations reaches up to 600 mm?year?1.Annual wind speed outside forests is 8 km?h?1in average,reaching up to 100 km?h?1during storms.The wind storms and the avalanches are the mostcommon natural disturbances causing tree blown-down in small and large areas(Martínez Pastur et al.2011;Kreps et al.2012).

    The study was carried out in mature mixed Nothofagus forests(>250 years)located south-west of Tierra del Fuego Island,Argentina (Fig.1),withoutprevious timber management in the last 50 years.Three Nothofagus forests situations were considered according to their canopy-layer composition:(1)pure deciduous forests with>80%N. pumilio crown cover;(2)pure evergreen forests with>80% N.betuloides crown cover;and(3)mixed N.pumilio-N. betuloides forests,with similar proportions of both species in the canopy.Likewise,two landscape locations were selected where these three forest types could be found:coast close to the Beagle Channel,and mountain toward the inner island.The coast sites were placed within the Tierra del Fuego National Park where the altitude is 50-100 m a.s.l.;the mean annual temperature varies between 2.1°C and 10.9°C and accumulated rainfall is 710 mm?year?1with abundant snow fall(Kreps et al.2012).Meanwhile,mountain forest sites were located through Garibaldi Pass in last foothills of the Andes Mountain Range where the altitude is 300–350 m a.s.l.;mean annual temperature varies between–0.6°C and 7.2°C and accumulated rainfall is 545 mm?year?1(Kreps et al.2012).

    Forest structure characterization

    Forest structure was characterized by angle-count sampling plots(n=10 for each forest type and location,N=60) following Bitterlich(1984).To apply this methodology we used a Criterion RD-1000(Laser Technology,USA)with a variable basal area factor(K)between 6 or 7,which allowed to recording tree basal area(BA)and density(N).Additionally,dominant height(DH)was estimated using a TruPulse 200 laser clinometer and distance rangefinder(Laser Technology,USA)by averaging height of the three taller trees per plot.Also,species and diameter at breast height(MD) was determined for each tree included in each plot.These data were collected between October and December 2012.

    Above the center of each forest plot,hemispherical photos of the canopy were taken.We used an 8 mm fisheye lens(Sigma Ex-AF4,Japan)mounted on a 35 mm digital camera(Nikon D50,Japan)with a tripod and level and they were oriented to magnetic north.The camera was set 1 m above the ground.This height is enough to avoid registering understory or shrub cover,and only register tree canopy cover.The program Gap Light Analyzer v2.0 0(Frazer et al.2001)was used to calculate canopy cover(CC)as a percentage,direct(RDR)and diffuse(RDF)solar radiation transmitted through canopy, and total radiation ratio(RT)(direct plus diffuse radiation transmitted through canopy,related to total radiation incident on a horizontal surface located above forest canopy). This data was obtained by taking just one photo per plot in January-February 2013(summer in South Hemisphere) avoiding the influence of direct sun light.The program offers the total value of RDR,RDF and RT for each month along the whole year,which are based on latitude and longitude,as well as cloudy index of each specific location. For details of inputs and models see Martínez Pastur et al. (2011).In our work,we averaged the values of the growing season(September to April).

    Fig.1 Location of study area in south-west Argentinean Tierra del Fuego.Rectangles indicate sampling locations:on the coast of the Beagle Channel(54°50'29"S and 68°31'54"W)and on the Andean Mountains region(54°41'05"S and 67°48'10"W)

    Understory structure and composition

    To evaluate the understory structure and composition we established two linear transects(10 m length each)associated to each forest structure plot(n=10 for each forest type and location,N=60).Transects were located parallel and 5 m apart,one from each other,and then averaged.We used the point-intercept method(Brancaleoni et al.2003) completing 50 intercept points per transect(every 20 cm) to record vascular plant species richness(calculated at the transect level)and ground cover by species and types(bare soil,debris and vegetation,which was differentiated in lichens,bryophytes,dicots,monocots and ferns).Cover was estimated by sub-layers(top or dominant layer,down or suppressed layer),and total(dominant plus suppressed sub-layers).We then calculated the number of layers at each intercepting point.Vascular plants(dicots,monocots and ferns)were taxonomically determined and classified by their origin(native or exotic)and life form(trees,shrubs, prostrate herbs,erect herbs,rhizomatous grasses and caespitose grasses)following Moore(1983)and Correa (1969–1998).Occurrence frequency for each species by forest type is shown in Additional file 1.We could not differentiate species of bryophyte and lichens so both of them were considered at group level.Bare soil refers to forest soil(including rocks)usually cover by litter but without vegetation.Debris included branches and trunks more than 3 cm in diameter.Also,an understory“dominance index”was calculated for each cover type(bare soil,debris,lichens, bryophytes,ferns,monocots and dicots),dividing each dominant cover type by total cover at each transect,and then the results were expressed as percentage.Fieldwork was conducted during January-February 2013(summer in South Hemisphere).This period is considered the most appropriated to conduct vegetation surveys in Nothofagus forests at this latitude to record all plant species inhabiting the understory layer and to avoid neglect species as can occur before or after that date within the growing season (Lencinas et al.2008a).The short duration of the growing period and the hard climate conditions influence on growth and develop of understory plant species to better express them(vegetative and reproductively)in this period.Juvenile trees(seedlings and saplings of N.pumilio and N.betuloides less than 0.5 m height)were also included as understory components.

    Data analysis

    Forest and understory structure variables were analyzed using a two-way ANOVA design,in which the forest type (three levels:deciduous,mixed,evergreen)and landscape location(two levels:coast,mountain)were the main factors of analysis.Likewise,structure differences between overstory species were compared in mixed forests,using also two-way ANOVAs,where species(two levels:N.pumilio, N.betuloides)and landscape location(two levels:coast, mountain)were the main factors of analysis.Variables were ln(Y)or ln(Y+1)transformed for the analysis to accomplish statistical ANOVA assumptions,but not transformed mean values are shown in Tables.The comparisons of mean values were done using the Tukey honestly significant difference test(p<0.05).Statistically significant interactions were plotted to get a better interpretation of our results.

    Complementary,four multivariate methods were applied: (i)a non-metric multidimensional scaling(NMS)ordination,with 9999 unrestricted random permutations of the raw data,(ii)a multi-response permutation procedures (MRPP),both based on Bray-Curtis(S?rensen)distance; (iii)a detrended correspond analysis(DCA)with down weighting of rare species and axes rescaling;and(iv)an Indicator Species Analysis(IndVal),including a random reallocation procedure with 4999 permutations(Monte Carlo test)to evaluate the significance of the maximum IndVals (p<0.005).All the multivariate analysis were conducted with species abundance(intercepted species in transects).

    NMS and MRPP were jointly used to assess similarity or difference of understory composition among forest types and locations,and support the visual assessment of these patterns.DCA allows evaluating differential botanical composition of each forest type.IndVal was useful to determine whether any species could be associated with certain habitats or treatments(Dufrêne and Legendre 1997).All multivariate methods were performed using vascular plant species cover,but NMS and MRPP also included general ground cover types(bare soil,debris, lichens and bryophytes)in the matrix.

    In MRPP,within-group homogeneity is measured as the chance-corrected within-group agreement(A):when all items are identical within groups,then A=1(perfect association);when heterogeneity within groups equals expectation by chance,then A=0(no association between sample units and defined groups);when there is less agreement within groups than expected by chance,then A<0 (McCune and Grace 2002).On the other hand,IndVal varied between 0(no indication)and 100(perfect indication), being highest when all individuals of a species are found in a single habitat(high specificity)and the species occurs in all samples of that habitat(high fidelity).Species with IndVal greater than 50 and p values lower than 0.05 can be considered as“indicator species”(characteristic species of a particular site or condition),while species with IndVal between 25 and 50 and p values lower than 0.05 can be considered as“detector species”(provide information from more than one habitat,and indicate the change direction in the habitat quality)(Tejeda-Cruz et al.2008).IndVal tests were performed for each forest type and location,jointly and separately.We used the software PC-ORD(Bruce and Mefford 1999)to conduct all multivariate analyses.

    Results

    Forest structure

    Forest structure showed significant differences among forest types and locations(Table 1).Deciduous forests had significantly fewer tree density and higher dominant height and mean diameter than evergreen forests.Mixed forests did not differentiate from evergreen forests in tree density and dominant height,while they differed from both deciduous and evergreen forests in mean diameter. On the other hand,canopy cover and basal area were significantly higher in mountain than in coastal sites,while differences were not significant between locations for the other variables.There were significant interactions for the basal area and direct solar radiation(Fig.2).Interactions in basal area occurred due to in coastal sites all forest types had similar values while in mountain sites deciduous stands had significant lower basal area than mixed and evergreeen stands.At the same time,basal area of coastal and mountain locations differed only in mixed forests. Interactions for direct solar radiation occurred due to in mountain sites all forest types had similar values while in coastal sites mixed stands showed significantly higher valuesthan deciduousand evergreen stands.And simultaneously,significant differences between locations occurred only in deciduous forests(Fig.2).

    The tree-layer of deciduous forest was 100%N.pumilio trees,while evergreen forests had 92%BA of N.betuloides and 8%of N.pumilio.The tree-layer of mixed forests had 67%trees of N.pumilio(58%BA)and 33%of N.betuloides (42%BA)(Table 2).Basal area,dominant height and mean diameter were similar for both Nothofagus species in mixed forests,but there were significant differences for tree density and percentage of basal area(greater values for N. pumilio than for N.betuloides).Likewise,basal area,density and dominant height of mixed stands were similar between locations,but there were significant differences in mean diameter of trees(greater values in the coast than in the mountain).Also,interactions occurred due to both species had similar values in mountain sites but in coastal sites in N.pumilio showed greater values than N.betuloides, whitout significant differences between locations for each tree species.

    Understory richness,structure and cover

    There were observed 46 vascular plant species in the understory of deciduous,mixed and evergreen Nothofagus forests,including four ferns,12 monocots and 30 dicots (Additionalfile1).Onlytwo species(Chiliotrichum diffusum and Taraxacum officinale)were observed but not intercepted(recorded)on transects.For this such species were included in the Annex but not in the statistical analysis.The overall richness was related to the forest type, following the order deciduous>mixed>evergreen,and ranging from 39 to 21 species,of which 33%of the total richness(15 species)were found in all forest types.The occurrence frequency of the three-environment-shared species was related to the forest type,being the more frequent species N.betuloides,Rubus geoides and Pernettya pumila in evergreen forests(with more than 70%occurrence frequency),and N.pumilio,Codonorchis lessonii and Adenocaulon chilense in deciduous forests(with more than 80%occurrence frequency).Four of these species were also the most frequent in mixed forests(N.betuloides,N. pumilio,C.lessonii and R.geoides).The exclusive species from each forest type were seven in deciduous,three in evergreen and only one in mixed forests,while 17 species were shared between mixed and deciduos forests only,and three species between mixed and evergreen forests only. There were not species shared between pure forests only.

    Likewise,richness was very similar between coast and mountain locations(33 and 36 species,respectively),and 50%of the total richness(23 species)were shared between the both locations.The more frequent shared species were N.pumilio in both locations,C.lessonii and N.betuloides in coast(with more than 67%occurrence frequency),and R.geoides and P.pumila in mountain(with more than 77%occurrence frequency).The exclusive species from each location were 10 species in coast and 13 in mountain.

    Table 1 Forest structure in Nothofagus forest of Tierra del Fuego,according to forest type(deciduous,mixed,or evergreen)and landscape location(coast or mountain)

    Fig.2 Interactions corresponding to Tables 1,2,3,4 and 5.For BA(basal area),RDR(direct solar radiation),Richness,Lichens and Bryophyte porcentual covers,Bryophyte and Fern understory indexes,letters indicates significant differences by Tukey test(p<0.05),with small letters for comparisons among forest types(deciduous,mixed or evergreen),and capital letters for comparisons between locations(coast or mountain).For BA(%)(proportional basal area),letters indicates significant differences by Tukey test(p<0.05),with small letters for comparisons between species(N.pumilio and N.betuloides). Letters are shown when differences were statistically significant.

    Only four of the observed vascular plants were alien species:Cerastium fontanum,Poa pratensis,Stellaria media and Taraxacum officinale(less than 10%occurrence frequency in the whole sampling).These did not occurr in evergreen forests at all,but they did in deciduous and mixed forests(less than 20%occurrence frequency and 1%cover),being found T.officinale and Cerastium fontanun in deciduous and mixed forests,S.media only in mixed forests,and Poa pratensis only in deciduous forests.On the other hand,P.pratensis were only found in coasts,S.media and T.officinale occurred only in mountains,and C.fontanum was present in both locations.

    Table 2 Forest structure in mixed Nothofagus forests of Tierra del Fuego,according to species(N.pumilio and N.betuloides)and landscape location(coast or mountain)

    Understory structure was influenced by both forest type and landscape location(Table 3).A maximum of 18 plant species per plot and 5 vertical strata were registered,with significantly highestrichness and strata quantity in deciduous forests.Interaction in richness occurred by significantly higher values in deciduos forests than in mixed and evergreen forests in coast sites(Fig.2),while in mountain,richness significantly differed among all forest types following the gradient evergreen

    About cover by taxonomic groups(Table 4),fern,monocot and dicot covers were significantly higher in deciduous forests than in mixed and evergreen forests,while monocot and dicot covers were also higher in moutain than in coast locations.Interaction in lichens occurred by significantly greater cover in coast than in mountain in mixed and evergreen forests(Fig.2),but not in deciduous one,while there were not detected statistical differences among forest types for each location.For bryophytes,interactions took place due to significantly higher cover in evergreen than in deciduous and mixed forests in coast locations(Fig.2),but similar and higher bryophyte cover in evergreen and mixed than in deciduous forests in mountain location.Likewise, all forest types presented significantly higher bryophyte cover in mountain than in coast locations.

    In relation to the understory index(Table 5),that reflects the relative dominance of each taxonomic group,dicots were dominant for all forest types and locations(percentage index =30%in average),without significant differences among them.Bare soil was also important in mixed and evergreen forests(percentage index=30%in average),and significatly lower in deciduos forests(percentage index=13%),followed by debris with also greater values in mixed and evergreen(percentage index=11%in average)than in deciduos forests (percentage index=6%).Monocots had the lowest index in evergreen(1%)and the greast index in deciduous(4%),with intermediate values for mixed forests,while lichens did not differentiated among forest types(percentage index=0.5%in average).On the other hand,bare soil,debris,monocots and lichens presented significantly higher indexes in coast than in mountain.Interactions in bryophyte index occurred by significant differences among deciduos vs.mixed and evergreen forests in mountain(Fig.2),but only among deciduous and evergreen in coast locations,while fern index was significantly higher in deciduous coast forests than in other forest type and locations.

    Table 3 Understory structure(species richness,strata quantity)and general ground cover types(bare soil,debris and vegetation)in Nothofagus forests of Tierra del Fuego,according to forest type(deciduous,mixed,or evergreen)and landscape location(coast or mountain)

    Table 4 Understory ground cover by taxonomic groups in Nothofagus forests of Tierra del Fuego,according to forest type (deciduous,mixed,or evergreen)and landscape location(coast or mountain)

    Understory assemblage patterns

    Ordination graphic by NMS highlighted relations among forest types and locations(Fig.3),since Axis 1 split between coast and mountain plots,while Axis 2 showed a gradient among deciduous,mixed and evergreen forests,which was more evident in mountain than in coast location.Moreover, the deciduous forests had higher internal variability while mixed and evergreen stands showed less internal variability forming closer groups.MRPP confirmed significant statistical differences between locations(T=–14.390;p<0.001; A=0.087),and also,but greater,among forest types(T=–18.260;p<0.001;A=0.158),with significant differences among all forest type pairwise comparisons(T<–4.062;p <0.004;A>0.040),where greatest A(0.193)was observed in evergreen vs.deciduous comparison.Likewise,MRPP showed significant differences among combined“forest type-location”(T=–20.980;p<0.001;A=0.296),with significant differences among all pairwise comparisons (T<–3.503;p<0.009;A>0.080),where greatest A (0.328)was observed in coast mixed vs.evergreen mountain comparison.

    For DCA analysis,only two ordination axes were retained for interpretation(eigenvalues of 0.5989 for Axis 1 and 0.0267 for Axis 2).These analysis emphasized species association with forest types,with more richness related to deciduous forests than evergreen or mixed forests(Fig.4).Likewise,there were more species shared only between mixed and deciduous forests(13 species, which did not occur in evergreen forests),than between mixed and evergreen forests(four species,which did not occur in deciduous forests),meanwhile none species were shared only between pure evergreen and deciduous forests.Among the species shared by all forest types, Maytenus disticha appeared as commonly observed inany of them,while R.geoides is more related to evergreen forests,P.alopecurus to mixed forests,and N. pumilio,Senecio acanthifolius,Codonorchis lessoni, Gunneramagellanica,Blechnum penna-marinaand Macrachaenium gracile are more associated to deciduous forests.Different life forms were distributed in all forest types,but erect herbs were more abundant in deciduous forests.

    Table 5 Understory index by groups(bare soil,debris,lichens,bryophytes,ferns,monocots and dicots)in Nothofagus forests, considering forest type(deciduous,mixed or evergreen)and landscape location(coast or mountain)

    Fig.3 Non-metric multidimensional scaling(NMS)ordination based on ground covers,to assess similarity or difference of understory structure among forest types and locations

    Indicator species analysis showed seven indicator species for deciduous forests(Table 6),of which five species did not occurred exclusively in this forest type(Fig.4),complemented with six detector species.Meanwhile,evergreen forests presented only two indicator and one detector species,and there was not detector or indicator species for mixed forests. Likewise,one indicator and two detector species were found in coastal sites,and four indicator and four detector species were found at mountain locations(Table 6).

    Discussion

    The three forest types investigated in this study differed in overstory tree composition and the selected forest structure variables(Table 1).Regarding solar radiation variables,it is known that beech crowns transmit only a small proportion of daylight to the forest floor(Ellenberg 1996;Barbier et al. 2008)and that the effective amount of radiation reaching the forest floor is determined by forest structural characteristics and tree species(van Oijen et al.2005;Barbier et al. 2008;Promis et al.2008;Martínez Pastur et al.2012;Kuninaga et al.2015).According with this,we expected the highest level of solar radiation in deciduous forest,because the canopy of N.pumilio is quite sparse and their leaves are substantially thinner than N.betuloides.However,we were unable to detect differences in solar radiation between deciduous and evergreen forests(Table 1),probably due to we taken hemispherical photos during summer but assumed the same(uniform)amount of canopy cover during the entire growing season for all forest types.A month-bymonth analysis would give us a better idea of the monthly differences in solar radiation due to the dynamics of canopies throughout the year(Promis et al.2010)mainly in deciduous and mixed forests.

    Fig.4 Detrended correspondence analysis(DCA)based on vascular plant species cover and life forms in Nothofagus forests.Species codes are in Additional file 1

    Table 6 Indicator species analysis(IndVal and probability between brackets)of understory vascular plant species,for each Nothofagus forest type and landscape location

    Our study showed that canopy-layer composition influences the understory structure and composition(Tables 3 and 4;Figs.3 and 4).While mixed forests have greater similarity to evergreen forests in the understory structural features,deciduous and mixed were similar in terms of the specific composition of plant assemblage.This could show deciduous and evergreen species differentially influence the understory communities beyond the simple comparison of the identity of plant species.

    In other forest ecosistems,the authors conclude that mixedwood forest may facilitate higher diversity and/or richness by providing a heterogeneous understory environment allowing the co-existence of common shade tolerant species along with resource-demanding,shade intolerant species and generalist species in patches with a mixed tree composition(Chavez and Macdonald 2012;Yu and Sun 2013).However,and contrary to our hypothesis,mixed broadleaved forests did not increase plant diversity at understory level compared to pure forests.We found maximum richness and vegetation cover in pure deciduos forests,not in mixed stands(Tables 3 and 4;Fig.4).Mixing of deciduous and evergreen tree species generate an intermediate level of plant diversity(richness and cover)derived from the mixture of shared plants,many of these herbs (Macrachaenium gracile,Senecio acanthifolius)shrubs (Maytenus disticha,Rubus geoides)or Nothofagus seedlings that grow under a closed canopy regardless of the canopy composition.In a review work,Barbier et al.(2008)concluded that mixing of deciduous and coniferous tree species generally affects understory diversity,but in almost all cases maximum diversity is observed in one of the pure stands, not in mixed stands.It is important to note that most of the current evidence about overstory-understory interactions in mixed forest derived from mixed-wood forests(deciduous broadleaf and conifer)of northern hemisphere.On the opposite side,the coverage of bryophytes group on the forest floor was improved in evergreen stands likely by lower competition with vascular plants for water and light at understory level(Gallo et al.2013).The light-sensitivity of bryophytes is supported by other field studies as well (Tinya and ódor 2016).

    Our results confirm the presence of exclusive species in pure forests(deciduous or evergreen)and landscape locations,and this species could indicate particular characteristics of each pure forest type(Table 6).Some common understory species on evergreen forests are also perennial species(e.g.,Gaultheria mucronata,Maytenus disticha), which could have some independence from environmental fluctuations because they have better nutrient retention over time(e.g.,N.betuloides can retain seven years a leaf in the tree)and can develop in poor content nutrient habits(Donoso and Donoso 2006;Promis et al.2008). The understory of deciduous Nothofagus forest is dominated by herbaceous plants associated with frequent small-scalle disturbances(e.g.,gap formation)or a more dynamic than the evergreen forest canopy(Lencinas et al. 2008b,2011;Promis et al.2010).The connection between deciduous and mixed forests is given by the presence and abundance of herbaceous species(e.g.,Dysopsis glechomoides,Gavilea lutea,Viola magellanica)that naturally grown under closed canopies but benefiting from smallscale disturbances at canopy level.Microclimatic conditions could varies more widely in deciduous and mixed forest stands than in evergreen stands,in response to the annual gradient of light,influencing soil properties(Gazol and Ibá?ez,2010;Yu and Sun 2013;Márialigeti et al. 2016).This generates more diverse site conditions over short distances and thus more habitat opportunities for a wide variety of plants.

    Other authors highlight the importance of canopyspecies leaf litter as a key factor influencing soil acidity and thereby nutrient stocks,whereas the upper 10 cm of soil are most significantly influenced by tree-species effects(Nordén 1994;Finzi et al.1998;Augusto et al.2002, 2003).Additionally,thickness of litter layer varies according to tree species(Neirynck et al.2000;Augusto et al. 2002;van Oijen et al.2005).A close positive correlation between soil pH and herb-layer diversity(e.g.H?rdtle etal.2003;Ferretti et al.2006;Barbier et al.2008),as well as litter layer thinness(e.g.Borchsenius et al.2004;van Oijen et al.2005;Barbier et al.2008)was particularly detected in forests dominated by beech,a group of trees which produces less decomposable and acidifying litter(Ellenberg 1996).However,our study did not include nutrient composition of forest soils neither the litter input from different canopies,but major contribution of the deciduous Nothofagus species in mixed stands(i.e.,highest proportion of basal area and tree density of deciduous species) may evidence major input of litterfall from deciduos species and therefore a strong influence on plant species assemblage.

    Regarding landscape locations,we found that understory structural and compositional patterns between pure and mixed forest are maintained in coast and mountain,except in richness,lichen and bryophytes cover,where we found interactions(Fig.2).As mountain forests grown under more stressful conditions,(e.g.,extreme temperature values, shorter growing season)than coastal forests,at similar latitude greater diversity would be expected in the later one (Lencinas et al.2008b).Contrary to this idea,we found that landscape location have an important effect over forest type and understory composition,showing that mountain was more diverse(e.g.,higher richness,strata quantity and vegetation cover)than coast,which could be due to a lowest content of litter in the coast as was described by Frangi et al.(2005)in evergreen forests floors considered less fertile. Moreover,the steepest slopes of mountain sites determines more frequent small-scale disturbances than coastal sites with gentle slopes(Roig 1998)generating greater microenvironmental heterogeneity.Following the environmental heterogeneity hypothesis(Huston 1994),small-scale environmental differences result in increased herb-layer diversity.Related to this hypothesis,we found some evidence (pers.obs.)indicating a more dynamic stand structure in mountain than in costal forests.In example,we observed big trees fallen over more than 50 years ago and dense groups of Nothofagus saplings(advanced regeneration)in deciduos and mixed stands on mountains,meanwhile the forest structure on the coast seems to be more stable and with recent fallen trees.These recent disturbances determines less stability of forest soil(e.g.,removal of the organic layer,exposed rocks,etc.)generating a poorest and nonstable understory in coastal sites compared to more stable conditions after older disturbances in mountain forests. However,more detailed studies are necessary to understand the effect of landscape location on understory species assemblages including a wide range of geographical locations and other landscape metrics(e.g.,forest patch size).

    The understanding of the interaction of overstory and understory species was key under future climate change.At these high latitudes,tree growth and distribution is most sensitive to climatic conditions(Hertel et al.2008;Massaccesi et al.2008;Ivancich et al.2012),and therefore,a shift in plant communities forced by climate change would be more evident.In Southern Patagonia,the maximum temperature and annual precipitation are predicted to increase until the end of the twenty-first century,possibly leading to a higher net primary productivity(Kreps et al.2012).This could cause a change on vegetation composition between coast and mountain forests,since exclusive plants of mountain forests are adapted to colder and dryer conditions.On the other hand,higher net primary productivity would produce important changes on soil properties,driving a shift in the proportion of deciduous and evergreen and therefore on understory communities.

    Conclusions

    We conclude that canopy-layer composition influences the herb-layer structure and diversity in two different ways:the presence of the evergreen component in mixed stands determines structural features at understory level,while the deciduous component determines the specific composition of plant assemblage.The deciduos forests were the most diverse and had higer number of exclusives species compared to mixed and evergreen forests.Lack of indicator species exclusive of mixed forest support it as transitional habitat where both evergreen and deciduous communities meet and integrate.

    At the same time,our results indicates that landscape has a major influence on the structure,diversity and richness of understory vegetation of pure and mixed forests mainly associated to the magnitude and frequency of natural disturbances.Moreover,we want to stress that mountain forest not only had higher herb-layer diversty but also more species exclusive of such habitat.Our study suggests that mixed Nothofagus forest provide suitable habitat for both pure deciduous and pure evergreen understory plant species and different assemblages in coastal and mountain sites.To better understand the light effects on interactions overstory-understory,we consider that a month-by-month analysis is necessary in a next research phase.We proposed that maintaining the mixture of canopy patch types within mixed stands will be important for conserving the natural patterns of understory plant composition in southern beech mixed forests.

    Additional file

    Additional file 1:List of vascular plant species observed in the understory of deciduous,mixed and evergreen Nothofagus forests of Tierra del Fuego,with occurrence frequency(%)in each forest type. (DOCX 31 kb)

    Acknowledgements

    The authors gratefully thank the Administración de Parques Nacionales and the staff of Tierra del Fuego National Park.Project APN-1236 supported this research during 2013-2014.

    Authors’contributions

    MVL and GMP conceived the idea.LM,RS GMP and MVL conducted the field work and statistical analyzes.MTM and AHH made figures and analyzed the bibliography.All authors analyzed the bibliography and contributed to the discussion of the results and wrote the manuscript.All authors read and approved the final manuscript.

    Competing interests

    The authors declare that they have no competing interests.

    Field permits

    Research permissions to sample understory plants at Tierra del Fuego National Park and Garibaldi were provided by the Administración de Parques Nacionales(Project APN-1236)and the local authority of Environment and Sustainable Development(S.D.S.yA.N°0155/2015),respectivelly.

    Allué CJ,Arranz A,Bava JO,Beneitez JM,Collado L,García-López JM(2010) Caracterización y cartografía fitoclimáticas del bosque nativo subantártico en la Isla Grande de Tierra del Fuego(Patagonia,Argentina).For Syst 2:189–207

    Augusto L,Ranger J,Binkley D,Rothe A(2002)Impact of several common tree species of European temperate forests on soil fertility.Ann For Sci 59:233–253

    Augusto L,Dupouey JL,Ranger J(2003)Effects of tree species on understory vegetation and environmental conditions in temperate forests.Ann For Sci 60:823–831

    Barbier S,Gosselin F,Balandier P(2008)Influence of tree species on understory vegetation diversity and mechanisms involved.A critical review for temperate and boreal forests.For Ecol Manage 254:1–15

    Bitterlich W(1984)The relascope idea.Relative measurements in forestry.CAB Londres,England

    Borchsenius F,Nielsen PK,Lawesson JE(2004)Vegetation structure and diversity of an ancient temperate deciduous forest in SW Denmark.Plant Ecol 175:121–135

    Brancaleoni L,Strelin J,Gerdol R(2003)Relationship between geomorphology and vegetation patterns in subantarctic Andean tundra of Tierra del Fuego. Polar Biol 26:404–410

    Bruce MC,Mefford MJ(1999)PC-ORD:multivariate analysis of ecological data. MJM Software Design,Oregon

    Chavez V,Macdonald SE(2012)Partitioning vascular understory diversity in mixedwood boreal forests:the importance of mixed canopies for diversity conservation.Fore Ecol Manage 271:19–26

    Correa MN(1969–1998)Flora Patagónica.Colección Científica INTA Tomo 8. Partes II,III,IVb,V,VI y VII.Buenos Aires.

    Cruz G,Ibarra M,Promis A,Cabello G(2007)Tipos de bosques,estructura y dinámica natural.In:Cruz G,Caldentey J(eds)Caracterización,silvicultura y uso de los bosques de coihue de Magallanes(Nothofagus betuloides(Mirb.) Blume)en la XII Región de Chile.Icaro Editions,Santiago,Chile

    Dale VH,Beyeler SC,Jackson B(2002)Understory vegetation indicators of anthropogenic disturbance in longleaf pine forests at Fort Benning,Georgia, USA.Ecol Indic 1:155–170

    Damascos MA(1996)Patrones morfologicos y floristicos en los bosques andinopatagónicos de Argentina.Ecol Aust 6:94–100

    Donoso CZ,Donoso P(2006)Nothofagus betuloides(Mirb)Oerst.In:Donoso Zegers C(ed)Las especies arbóreas de los bosques templados de Chile y Argentina autoescología.Marisa Cuneo Ediciones,Valdivia

    Dreiss LM,Volin JC(2013)Influence of leaf phenology and site nitrogen on invasive species establishment in temperate deciduous forest understories. For Ecol Manage 296:1–8

    Dufrêne M,Legendre P(1997)Species assemblages and indicator species:the need for a flexible asymmetrical approach.Ecol Monog 67:345–366

    Ellenberg H(1996)Vegetation Mitteleuropas mit den Alpen in okologischer Sicht. 5th edition.Ulmer,Stuttgart.

    Ferretti M,Calderisi M,Amoriello T,Bussotti F,Canullo R,Campetella G,Costantini A,Fabbio G,Mosello R(2006)Factors influencing vascular species diversity in the CONECOFOR permanent monitoring plots.Ann Inst Sper Selv 30:97–106

    Finzi AC,Canham CD,van Breemen N(1998)Canopy tree-soil interactions within temperate forests:species effects on pH an cations.Ecol Appl 8(2):447–454

    Forrester DI,Bonal D,Dawud S,Gessler A,Granier A,Pollastrini M,Grossiord C (2016)Drought responses by individual tree species are not often correlated with tree species diversity in European forests.J Appl Ecol 53:1725–1734

    Frangi JL,Barrera MD,Puigdefábregas J,Yapura PF,Arambarri AM,Richter L (2005)Ecología de los bosques de Tierra del Fuego.In:Arturi MF,Frangi JL, Goya JF(eds)Ecología y manejo de los bosques de la Argentina.Universidad Nacional de La Plata,Argentina

    Frazer GW,Fournier RA,Trofymow JA,Gall RJ(2001)A comparison of digital and film fisheye photography for analysis of forest canopy structure and gap light transmission.Agric For Meteorol 109:249–263

    Gallo EA,Lencinas MV,Martínez Pastur G(2013)Site quality influence over understory plant diversity in old-growth and harvested Nothofagus pumilio forest.For Syst 22:25–38

    Gazol A,Ibá?ez R(2010)Scale-specific determinants of a mixed beech and oak seedling–sapling bank under different environmental and biotic conditions. Plant Eco 211(1):37–48

    Gilliam FS(2007)The ecological significance of the herbaceous layer in temperate forest ecosystems.Bioscience 57:845–858

    Gracia M,Montané F,Piqué J,Retana J(2007)Overstory structure and topographic gradients determining diversity and abundance of understory shrub species in temperate forests in central Pyrenees(NE Spain).For Ecol Manage 242:391–397

    H?rdtle W,Oheimb G,von Westphal C(2003)The effects of light and soil conditions on the species richness of the ground vegetation of deciduous forests in northern Germany(Schleswig-Holstein).For Ecol Manage 182:327–338

    Hertel D,Therburg A,Villalba R(2008)Above-and below-ground response by Nothofagus pumilio to climatic conditions at the transition from the steppeforest boundary to the alpine treeline in southern Patagonia,Argentina.Plant Ecol Divers 1(1):21–33

    Huston MA(1994)Biological diversity.The coexistence of species on changing landscapes.Cambridge University Press,Cambridge,p 708

    Ivancich H,Martínez Pastur G,Roig F,Barrera M,Pulido F(2012)Changes in height growth patterns in the upper treeline forests of Tierra del Fuego in relation to climate change.Bosque 33(3):267–270

    Kreps G,Martínez Pastur G,Peri P(2012)Cambio climático en Patagonia sur. Escenarios futuros en el manejo de los recursos naturales.Ediciones INTA, Buenos Aires.

    Kudo G,Ida TY,Tani T(2008)Linkages between phenology,pollination, photosynthesis,and plant reproduction in deciduous forest understory plants.Ecology 89(2):321–331

    Kuninaga T,Hirayama K,Sakimoto M(2015)Negative canopy–understorey interaction shapes the sapling bank of Fagus crenata in a cool-temperate, conifer–hardwood mixed forest.Plant Ecol 216:1191–1202

    Lencinas MV,Martínez Pastur G,Solán R,Gallo E,Cellini JM(2008a)Forest management with variable retention impact over bryophyte communities of Nothofagus pumilio understory.Forstarchiv 79:77–82

    Lencinas MV,Martínez Pastur G,Rivero P,Busso C(2008b)Conservation value of timber quality versus associated non-timber quality stands for understory diversity in Nothofagus forests.Biodivers Conserv 17:2579–2597

    Lencinas MV,Martínez Pastur G,Gallo E,Cellini JM(2011)Alternative silvicultural practices with variable retention to improve understory plant diversity conservation in southern Patagonian forests.For Ecol Manage 262:1236–1250

    Márialigeti S,Tinya F,Bidló A,ódor P(2016)Environmental drivers of the composition and diversity of the herb layer in mixed temperate forests in Hungary.Plant Ecol 217:549–563

    Martínez Pastur G,Peri PL,Fernández MC,Stafieri G,Lencinas A(2002)Changes in understory species diversity during the Nothofagus pumilio forest management cycle.J For Res 7(3):2–10

    Martínez Pastur G,Cellini JM,Lencinas MV,Barrera M,Peri PL(2011)Environmental variables influencing regeneration of Nothofagus pumilio in a system with combined aggregated and dispersed retention.For Ecol Manage 261:178–186

    Martínez Pastur G,Jordán C,Lencinas MV,Soler R,Ivancich H,Kreps G(2012) Landscape and microenvironmental conditions influence over regeneration dynamics in old-growth Nothofagus betuloides Southern Patagonian forests. Plant Biosyst 146:201–213

    Massaccesi G,Roig F,Martínez Pastur G,Barrera M(2008)Growth patterns of Nothofagus pumilio trees along altitudinal gradients in Tierra del Fuego, Argentina.Trees 22(2):245–255

    Matteri C,Schiavone M(2002)Catálogo de los musgos(Bryophyta)de la región Fueguina en Argentina y Chile.Revista del Museo Argentino de Ciencias Naturales 4:111–138

    McCune B,Grace JB(2002)Analysis of ecological communities.MJM Software Design,Oregon

    Mittermeier RA,Mittermeier CG,Brooks TM,Pilgrim JD,Konstant WR,da Fonseca GAB,Kormos C(2003)Wilderness and biodiversity conservation.PNAS 100(18):10309–10313

    Moore D(1983)Flora of Tierra del Fuego.Anthony Nelson,England,Missouri Botanical Garden,USA

    Neirynck J,Mirtcheva S,Sioen G,Lust N(2000)Impact of Tilia platyphyllos Scop., Fraxinus excelsior L.,Acer pseudoplatanus L.,Quercus robur L.and Fagus sylvatica L.on earthworm biomass and physico-chemical properties of a loamy soil.For Ecol Manage 133:275–286

    Nordén U(1994)Influence of broad-leaved tree species on pH and organic matter content of forest topsoils in Scania,South Sweden.Scand J For Res 9:1–8

    Promis A,Cruz G,Reif A,G?rtner S(2008)Nothofagus betuloides(MIRB.)OERST 1871(Fagales:Nothofagaceae)forests in southern Patagonia and Tierra del Fuego.Ann Inst Pat 36:53–68

    Promis A,Gartner S,Reif A,Cruz G(2010)Effects of natural small-scale disturbances on below-canopy solar radiation and regeneration patterns in an old-growth Nothofagus betuloides forest in Tierra del Fuego,Chile.Allg Forst Jagdztg 181(3/4):53–64

    Roig F(1998)La vegetación de la Patagonia.In:Correa M(ed)Flora Patagónica. INTA Colección Científica,Buenos Aires

    Simonson WD,Allen HD,Coomes DA(2014)Overstorey and topographic effects on understories:evidence for linkage from cork oak(Quercus suber)forests in southern Spain.For Ecol Manage 328:35–44

    Spyreas G,Matthews JW(2006)Floristic conservation value,nested understory floras, and the development of second-growth forest.Ecol Applic 16:1351–1366

    Tejeda-Cruz C,Mehltreter K,Sosa VJ(2008)Indicadores ecológicos multitaxonómicos,p.271–278.In:Manson RH,Hernández-Ortiz V,Gallina S, Mehltreter K(eds).Agroecosistemas cafetaleros de Veracruz:Biodiversidad, manejo y conservación.Instituto de Ecología e Instituto Nacional de Ecología,México.

    Tinya F,ódor P(2016)Congruence of the spatial pattern of light and understory vegetation in an old-growth,temperate mixed forest.For Ecol Manage 381:84–92

    Tuhkanen S(1992)The climate of Tierra del Fuego from a vegetation geographical point of view and its ecoclimatic counterparts elsewhere.Acta Bot Fenn 145:64

    Van Oijen M,Rougier J,Smith R(2005)Bayesian calibration of process-based forest models:bridging the gap between models and data.Tree Physiol 25:915–927

    Veblen TT,Ashton DH,Schlegel FM,Veblen AT(1978)Influencia del estrato arbóreo sobre los estratos inferiores en un bosque,mixto,perennicaducifolio de Antillanca,Osorno,Chile.C.D.Oxf.:182.4.Bosque 2(2):88–104

    Veblen TT,Veblen AT,Schlegel FM(1979)Understorey patterns in mixed evergreen-deciduous Nothofagus forests in Chile.J Ecol 67(3):809–823

    Yu M,Sun OJ(2013)Effects of forest patch type and site on herb-layer vegetation in a temperate forest ecosystem.Fore Ecol Manage 300:14–20

    29 September 2016 Accepted:10 May 2017

    *Correspondence:lucianamestre@hotmail.com

    National Scientific and Technical Research Council(CONICET),Austral Center for Scientific Research(CADIC-CONICET),Houssay 200,9410 Ushuaia,Tierra del Fuego,Argentina

    ?The Author(s).2017 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/),which permits unrestricted use,distribution,and reproduction in any medium,provided you give appropriate credit to the original author(s)and the source,provide a link to the Creative Commons license,and indicate if changes were made.

    欧美日韩精品网址| 成人av一区二区三区在线看 | 大码成人一级视频| 国产免费一区二区三区四区乱码| 飞空精品影院首页| 国产欧美日韩精品亚洲av| 大型av网站在线播放| 99热国产这里只有精品6| 正在播放国产对白刺激| 久久人人爽人人片av| 美女午夜性视频免费| www.999成人在线观看| 精品福利永久在线观看| 91麻豆av在线| 老司机深夜福利视频在线观看 | 中文字幕色久视频| 欧美性长视频在线观看| 欧美av亚洲av综合av国产av| 99国产综合亚洲精品| 99热国产这里只有精品6| 日本av手机在线免费观看| 精品福利观看| 午夜影院在线不卡| 狠狠狠狠99中文字幕| 精品一区二区三区四区五区乱码| 99久久99久久久精品蜜桃| 久久久久久久大尺度免费视频| 人妻一区二区av| 日本a在线网址| 亚洲国产av新网站| 精品久久久精品久久久| videos熟女内射| 91成年电影在线观看| 久久久国产成人免费| 黑人欧美特级aaaaaa片| 老熟妇仑乱视频hdxx| 91精品三级在线观看| 欧美日韩成人在线一区二区| 亚洲欧美日韩高清在线视频 | 日本欧美视频一区| 日本a在线网址| 中文字幕色久视频| 欧美激情久久久久久爽电影 | 在线看a的网站| 国产国语露脸激情在线看| av国产精品久久久久影院| 在线观看www视频免费| 国产又色又爽无遮挡免| a在线观看视频网站| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 99国产精品免费福利视频| 嫁个100分男人电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 手机成人av网站| 国产成人系列免费观看| 99国产精品一区二区蜜桃av | 最黄视频免费看| 精品少妇久久久久久888优播| 久久久久久免费高清国产稀缺| 亚洲视频免费观看视频| 亚洲熟女精品中文字幕| 国产成人免费无遮挡视频| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 欧美黄色淫秽网站| 国产精品久久久久久精品电影小说| 国产在线观看jvid| 亚洲国产精品一区二区三区在线| 各种免费的搞黄视频| www.精华液| 岛国在线观看网站| 国产片内射在线| 制服人妻中文乱码| 9色porny在线观看| 亚洲 欧美一区二区三区| 国产精品一二三区在线看| 久久中文看片网| 国产色视频综合| 精品少妇一区二区三区视频日本电影| 99香蕉大伊视频| 欧美成狂野欧美在线观看| 日韩中文字幕视频在线看片| 两个人看的免费小视频| 国产精品一区二区精品视频观看| 韩国精品一区二区三区| 国产精品国产av在线观看| 19禁男女啪啪无遮挡网站| 久久国产精品男人的天堂亚洲| 欧美av亚洲av综合av国产av| 搡老岳熟女国产| 日韩欧美一区二区三区在线观看 | 久久99一区二区三区| 纵有疾风起免费观看全集完整版| 一区二区三区精品91| 国产精品秋霞免费鲁丝片| 三级毛片av免费| 大片电影免费在线观看免费| 99国产精品99久久久久| 精品国产国语对白av| 国产片内射在线| 人人澡人人妻人| 国产av国产精品国产| 一本大道久久a久久精品| 侵犯人妻中文字幕一二三四区| 婷婷成人精品国产| 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看| 国产欧美日韩精品亚洲av| 亚洲精品中文字幕一二三四区 | 国产欧美日韩一区二区精品| 亚洲av片天天在线观看| 午夜福利乱码中文字幕| 国产熟女午夜一区二区三区| 咕卡用的链子| 久久狼人影院| 成年人黄色毛片网站| 国产免费视频播放在线视频| 搡老岳熟女国产| 一级,二级,三级黄色视频| 搡老乐熟女国产| 999久久久精品免费观看国产| 捣出白浆h1v1| 日本撒尿小便嘘嘘汇集6| av电影中文网址| 国产淫语在线视频| 后天国语完整版免费观看| 精品高清国产在线一区| 日本欧美视频一区| 久久久国产成人免费| 9191精品国产免费久久| 欧美在线黄色| 一区在线观看完整版| 国产成人免费无遮挡视频| 国产野战对白在线观看| 18在线观看网站| 午夜免费鲁丝| 美女午夜性视频免费| 麻豆av在线久日| 91老司机精品| 在线观看人妻少妇| 这个男人来自地球电影免费观看| 久久 成人 亚洲| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 色视频在线一区二区三区| 18禁观看日本| 1024视频免费在线观看| av国产精品久久久久影院| videos熟女内射| 亚洲国产欧美一区二区综合| 天堂8中文在线网| 日韩制服丝袜自拍偷拍| 色精品久久人妻99蜜桃| 黄色怎么调成土黄色| 女人高潮潮喷娇喘18禁视频| 91国产中文字幕| 丝袜美足系列| 国产91精品成人一区二区三区 | 一区二区三区四区激情视频| 午夜影院在线不卡| 精品熟女少妇八av免费久了| 下体分泌物呈黄色| 性高湖久久久久久久久免费观看| 制服诱惑二区| 日本一区二区免费在线视频| 80岁老熟妇乱子伦牲交| 国产福利在线免费观看视频| 精品久久蜜臀av无| 中文字幕人妻熟女乱码| 丰满人妻熟妇乱又伦精品不卡| 在线十欧美十亚洲十日本专区| 日韩 亚洲 欧美在线| 91字幕亚洲| 啦啦啦免费观看视频1| 成年人免费黄色播放视频| 男女免费视频国产| 亚洲av男天堂| 中文字幕高清在线视频| 12—13女人毛片做爰片一| 亚洲视频免费观看视频| 天天操日日干夜夜撸| 久久久久久亚洲精品国产蜜桃av| 青青草视频在线视频观看| 曰老女人黄片| 精品人妻一区二区三区麻豆| 色94色欧美一区二区| 欧美人与性动交α欧美精品济南到| 高清欧美精品videossex| 国产精品国产av在线观看| 久久精品久久久久久噜噜老黄| 在线 av 中文字幕| 国产三级黄色录像| 91大片在线观看| 99久久人妻综合| 亚洲欧美精品自产自拍| 久久青草综合色| 天天躁狠狠躁夜夜躁狠狠躁| 别揉我奶头~嗯~啊~动态视频 | 国产男女超爽视频在线观看| 久久影院123| 韩国高清视频一区二区三区| 亚洲午夜精品一区,二区,三区| 男女边摸边吃奶| 国产成人影院久久av| 国产亚洲午夜精品一区二区久久| 久久午夜综合久久蜜桃| 激情视频va一区二区三区| 精品国产乱码久久久久久男人| 亚洲av电影在线观看一区二区三区| 亚洲精品一区蜜桃| 如日韩欧美国产精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 老司机靠b影院| 搡老岳熟女国产| 久久久久精品人妻al黑| 亚洲激情五月婷婷啪啪| 狂野欧美激情性xxxx| 夫妻午夜视频| 狂野欧美激情性bbbbbb| 久久精品aⅴ一区二区三区四区| 99国产精品一区二区蜜桃av | 精品亚洲成a人片在线观看| 亚洲av电影在线进入| 黄色片一级片一级黄色片| 欧美精品一区二区免费开放| 欧美一级毛片孕妇| 五月天丁香电影| 久久亚洲精品不卡| 色精品久久人妻99蜜桃| 欧美性长视频在线观看| 国产三级黄色录像| 国产高清videossex| av线在线观看网站| 国产精品成人在线| 亚洲欧美精品综合一区二区三区| 亚洲欧美精品自产自拍| 大片免费播放器 马上看| 国产不卡av网站在线观看| 少妇人妻久久综合中文| 美女视频免费永久观看网站| 麻豆av在线久日| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人国产一区在线观看| 亚洲国产欧美在线一区| av天堂在线播放| 中文欧美无线码| 在线十欧美十亚洲十日本专区| 国产精品九九99| 国产亚洲午夜精品一区二区久久| 天天躁狠狠躁夜夜躁狠狠躁| 菩萨蛮人人尽说江南好唐韦庄| 麻豆乱淫一区二区| 亚洲国产欧美一区二区综合| 在线观看舔阴道视频| 亚洲久久久国产精品| 国产男人的电影天堂91| 午夜福利视频在线观看免费| 精品一区二区三区四区五区乱码| 成人av一区二区三区在线看 | 国产一级毛片在线| 免费人妻精品一区二区三区视频| 手机成人av网站| 黑丝袜美女国产一区| 色视频在线一区二区三区| 一区福利在线观看| 亚洲国产av新网站| 成年av动漫网址| 丝袜人妻中文字幕| 又黄又粗又硬又大视频| 国产av精品麻豆| 我的亚洲天堂| 亚洲av美国av| 久久久久精品国产欧美久久久 | 午夜老司机福利片| 国产一级毛片在线| 日本av免费视频播放| 一边摸一边抽搐一进一出视频| 老司机影院毛片| 一级毛片女人18水好多| 人妻一区二区av| 国产亚洲午夜精品一区二区久久| 美女中出高潮动态图| 国产精品秋霞免费鲁丝片| 另类精品久久| 欧美日韩成人在线一区二区| 秋霞在线观看毛片| 国产精品 国内视频| 亚洲专区字幕在线| 精品久久久久久久毛片微露脸 | 人妻 亚洲 视频| 国产成人欧美在线观看 | 亚洲成人免费电影在线观看| 一级,二级,三级黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | 91成年电影在线观看| 天堂8中文在线网| 久久精品国产亚洲av香蕉五月 | 久9热在线精品视频| 亚洲精品一卡2卡三卡4卡5卡 | 真人做人爱边吃奶动态| 欧美激情极品国产一区二区三区| 脱女人内裤的视频| 亚洲精品美女久久av网站| 欧美精品av麻豆av| 国产不卡av网站在线观看| 2018国产大陆天天弄谢| 国产真人三级小视频在线观看| 操美女的视频在线观看| 99国产极品粉嫩在线观看| 别揉我奶头~嗯~啊~动态视频 | 99热网站在线观看| 99久久人妻综合| 香蕉国产在线看| 一区二区三区乱码不卡18| 99热国产这里只有精品6| 午夜福利影视在线免费观看| 日韩电影二区| 国产成人精品久久二区二区91| 老司机福利观看| 亚洲国产欧美日韩在线播放| 一个人免费在线观看的高清视频 | 少妇的丰满在线观看| 一级a爱视频在线免费观看| 欧美激情极品国产一区二区三区| 爱豆传媒免费全集在线观看| 国产日韩欧美视频二区| 视频区欧美日本亚洲| 免费久久久久久久精品成人欧美视频| 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频 | 巨乳人妻的诱惑在线观看| 青草久久国产| 丰满少妇做爰视频| h视频一区二区三区| 99久久国产精品久久久| 老熟妇乱子伦视频在线观看 | 亚洲av男天堂| 韩国高清视频一区二区三区| 国产有黄有色有爽视频| 伦理电影免费视频| 精品人妻一区二区三区麻豆| 国产一区二区三区av在线| 亚洲午夜精品一区,二区,三区| 99国产极品粉嫩在线观看| 99香蕉大伊视频| 国产成人欧美| 国产极品粉嫩免费观看在线| 韩国精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产精品影院久久| 国产区一区二久久| 精品亚洲成a人片在线观看| 日韩制服骚丝袜av| 日本五十路高清| bbb黄色大片| 久久久久久久国产电影| 一区二区av电影网| 精品人妻一区二区三区麻豆| 狂野欧美激情性xxxx| 丁香六月天网| 91麻豆av在线| 人妻 亚洲 视频| svipshipincom国产片| 精品视频人人做人人爽| 男女高潮啪啪啪动态图| 欧美黄色片欧美黄色片| 亚洲 欧美一区二区三区| 老汉色av国产亚洲站长工具| 日日夜夜操网爽| 亚洲精品国产精品久久久不卡| 午夜免费鲁丝| 五月开心婷婷网| 亚洲色图综合在线观看| 午夜福利在线免费观看网站| 99九九在线精品视频| 大陆偷拍与自拍| 午夜91福利影院| 精品久久久久久久毛片微露脸 | 国产在线视频一区二区| 91精品伊人久久大香线蕉| 美女高潮到喷水免费观看| 国产av国产精品国产| 欧美激情极品国产一区二区三区| 色播在线永久视频| 中文字幕另类日韩欧美亚洲嫩草| 色老头精品视频在线观看| 交换朋友夫妻互换小说| 一边摸一边做爽爽视频免费| 如日韩欧美国产精品一区二区三区| 国产一级毛片在线| 久久这里只有精品19| 乱人伦中国视频| 久久这里只有精品19| 精品人妻一区二区三区麻豆| 成人三级做爰电影| 午夜福利一区二区在线看| 一级毛片电影观看| 中文字幕人妻丝袜一区二区| 99国产精品99久久久久| 中文字幕人妻丝袜制服| 色老头精品视频在线观看| 老鸭窝网址在线观看| 亚洲五月色婷婷综合| 亚洲五月婷婷丁香| 18禁观看日本| 一本—道久久a久久精品蜜桃钙片| 在线观看人妻少妇| 国产成人啪精品午夜网站| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 国产又色又爽无遮挡免| 久久99热这里只频精品6学生| 国产一区二区三区综合在线观看| 美女福利国产在线| 天天添夜夜摸| 在线 av 中文字幕| 免费高清在线观看视频在线观看| 午夜视频精品福利| 桃花免费在线播放| 日韩 欧美 亚洲 中文字幕| 国产精品欧美亚洲77777| 国产av一区二区精品久久| 亚洲国产精品一区三区| 亚洲精品粉嫩美女一区| 天天躁日日躁夜夜躁夜夜| 男女免费视频国产| 欧美成人午夜精品| 在线观看免费午夜福利视频| 精品人妻一区二区三区麻豆| 国产男女内射视频| 国产伦理片在线播放av一区| 一级毛片精品| 黄频高清免费视频| 亚洲精品成人av观看孕妇| 欧美成狂野欧美在线观看| 18禁观看日本| av又黄又爽大尺度在线免费看| 精品人妻在线不人妻| 精品第一国产精品| 91成年电影在线观看| 欧美 亚洲 国产 日韩一| 欧美精品人与动牲交sv欧美| 婷婷色av中文字幕| 精品国产一区二区久久| 美女中出高潮动态图| 一区二区三区激情视频| 欧美日韩亚洲综合一区二区三区_| 久久久久网色| 成人手机av| 在线亚洲精品国产二区图片欧美| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| av天堂久久9| www.av在线官网国产| 一本久久精品| 麻豆av在线久日| 黄色 视频免费看| 日韩一卡2卡3卡4卡2021年| 国产xxxxx性猛交| 国产高清国产精品国产三级| 女性生殖器流出的白浆| 99国产精品一区二区三区| 中国美女看黄片| 老司机靠b影院| 两人在一起打扑克的视频| 精品一品国产午夜福利视频| 亚洲欧美日韩高清在线视频 | 日本撒尿小便嘘嘘汇集6| 亚洲欧美色中文字幕在线| 久久久精品94久久精品| 大香蕉久久网| 岛国在线观看网站| 日本av免费视频播放| 午夜免费成人在线视频| 国产精品 欧美亚洲| 亚洲国产精品999| 亚洲成人手机| 一级片免费观看大全| 亚洲精品国产区一区二| 老司机午夜十八禁免费视频| 亚洲精品久久成人aⅴ小说| 我的亚洲天堂| 免费观看人在逋| av免费在线观看网站| 国产熟女午夜一区二区三区| 国产精品久久久av美女十八| 黑人欧美特级aaaaaa片| 美女脱内裤让男人舔精品视频| 下体分泌物呈黄色| 99香蕉大伊视频| 国产成人系列免费观看| 亚洲情色 制服丝袜| 久久女婷五月综合色啪小说| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲三区欧美一区| 亚洲伊人久久精品综合| 少妇被粗大的猛进出69影院| 午夜福利视频精品| 亚洲一码二码三码区别大吗| 国产成人欧美| 91麻豆av在线| 国产在线视频一区二区| 免费日韩欧美在线观看| 欧美日韩视频精品一区| 国产成人系列免费观看| 91老司机精品| 天堂中文最新版在线下载| 天堂俺去俺来也www色官网| 性少妇av在线| 国产欧美日韩一区二区三 | 国产一级毛片在线| 亚洲一区二区三区欧美精品| 亚洲午夜精品一区,二区,三区| 黄色视频,在线免费观看| 在线观看免费高清a一片| www.999成人在线观看| 制服人妻中文乱码| 久久久欧美国产精品| 国产亚洲午夜精品一区二区久久| 亚洲精品久久成人aⅴ小说| 法律面前人人平等表现在哪些方面 | 一二三四社区在线视频社区8| 精品国产乱子伦一区二区三区 | 脱女人内裤的视频| 啦啦啦 在线观看视频| 免费一级毛片在线播放高清视频 | 女人高潮潮喷娇喘18禁视频| 夜夜夜夜夜久久久久| 波多野结衣av一区二区av| 亚洲第一欧美日韩一区二区三区 | 男人添女人高潮全过程视频| 两个人看的免费小视频| 叶爱在线成人免费视频播放| 日韩一区二区三区影片| 免费不卡黄色视频| 国产不卡av网站在线观看| 亚洲av片天天在线观看| 麻豆av在线久日| 黑人巨大精品欧美一区二区蜜桃| 午夜成年电影在线免费观看| 五月开心婷婷网| 下体分泌物呈黄色| 99热国产这里只有精品6| 久久狼人影院| 啦啦啦视频在线资源免费观看| 51午夜福利影视在线观看| 另类精品久久| 亚洲第一欧美日韩一区二区三区 | 精品国产乱码久久久久久男人| 国产视频一区二区在线看| 性少妇av在线| 欧美激情 高清一区二区三区| 真人做人爱边吃奶动态| 99精国产麻豆久久婷婷| 免费av中文字幕在线| 国产精品久久久久成人av| 少妇裸体淫交视频免费看高清 | 一区二区av电影网| 天天影视国产精品| 制服人妻中文乱码| 国产在线免费精品| 国产免费一区二区三区四区乱码| 欧美日韩亚洲综合一区二区三区_| 中文字幕精品免费在线观看视频| 国产精品 国内视频| 久久久久久免费高清国产稀缺| 免费观看a级毛片全部| 欧美日韩亚洲高清精品| 国产成人欧美在线观看 | 精品久久久精品久久久| 黄频高清免费视频| 欧美激情久久久久久爽电影 | av欧美777| 国产高清国产精品国产三级| 中文字幕人妻丝袜一区二区| 国产免费视频播放在线视频| 欧美日韩亚洲综合一区二区三区_| 免费日韩欧美在线观看| 久久天躁狠狠躁夜夜2o2o| 激情视频va一区二区三区| 国产极品粉嫩免费观看在线| 国产日韩欧美亚洲二区| 国产亚洲欧美精品永久| 宅男免费午夜| 黄色毛片三级朝国网站| 精品国产乱码久久久久久男人| a 毛片基地| 大香蕉久久成人网| 亚洲国产成人一精品久久久| 亚洲av男天堂| 老汉色av国产亚洲站长工具| av免费在线观看网站| 国产97色在线日韩免费| 在线天堂中文资源库| 一级黄色大片毛片| 爱豆传媒免费全集在线观看| 五月开心婷婷网| 亚洲中文字幕日韩| 成人三级做爰电影| 国产男女超爽视频在线观看| 亚洲精品中文字幕在线视频| 国产精品 欧美亚洲| 日韩欧美免费精品| 三上悠亚av全集在线观看| 国产精品一区二区在线不卡| 男女免费视频国产| 国产精品影院久久| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品无人区| 日本五十路高清| 久久99热这里只频精品6学生| 视频在线观看一区二区三区| 国产日韩欧美在线精品| 精品少妇内射三级|