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

    Height increment of understorey Norway spruces under different tree canopies

    2014-02-24 07:33:51OlaviLaihoTimoPukkalaandErkkihde
    Forest Ecosystems 2014年1期

    Olavi Laiho,Timo Pukkalaand Erkki L?hde

    Height increment of understorey Norway spruces under different tree canopies

    Olavi Laiho1,Timo Pukkala2*and Erkki L?hde1

    Background:Stands having advance regeneration of spruce are logical places to start continuous cover forestry (CCF)in fertile and mesic boreal forests.However,the development of advance regeneration is poorly known.

    Advance regeneration;Continuous cover forestry;Growth model;Picea abies;Uneven-aged

    management

    Background

    Continuous cover forestry(CCF)and other near-natural forest management systems are gaining popularity in several countries(Brunner et al.,2006;Hasenauer,2006; Schütz et al.,2012).In Finland,public opinion is against clear-felling,and also forest landowners are often unwilling to clear-cut their forests(Valkeap?? et al.,2009).According to a recent survey,majority of forest landowners are interested in continuous cover forest management (Kumela and H?nninen,2011).

    Change in general attitude,encouraging results regarding the profitability of uneven-aged forest management (Tahvonen,2009;Tahvonen et al.,2010;Pukkala et al., 2010,2011a)and the decreasing importance of pulp and paper industries in Finland’s economy have led to increasing criticism against the prevailing silvicultural practices.As a consequence,forestry legislation is being modified(Anonym,2012),nowallowingforest landowners to use CCF without any specific justification. Also,instructions for good forestry practice(Anonym, 2006)are being updated.As a result,forest landowners now have the freedom to pursue CCF if they decide so. The sustainability of CCF depends on the amount and composition of regeneration and the survival and revival of the understorey trees.The first Finnish national forest inventory,conducted in 1921–1924,reported that on the basis of subjective assessments,only 3%of the forest area had a good-quality understorey(Ilvessalo,1956).However, a closer look at the inventory results reveals that understoreys have been utilized in forest management much more than one would deduct from their reported occurrence.For example,10-year old sapling stands covered 7%of forest land in 1921–1924,but 40 years later the percentage of 50-year old stands was 28%(Mikola,1966),instead of 7%,suggesting that a majority of 50-year-old stands originated from a released understorey.Therefore,understoreys are more common than indicated by the inventory reports but, due to their spatial heterogeneity and size variation,they are seldomly evaluated as good enough starting points foreven-aged stands.However,practical experience shows that almost all released understoreys develop into fully-stocked productive stands.Sarvas(1951)and Vuokila(1956)found that,in the 1950s,most spruce stands in Finland in fact originated from an understorey.

    Looking at the diameter distributions of the sample plots of national forest inventories instead of subjective classifications,shows that there is plenty of advance regeneration in Finnish forests.For example,L?hde et al. (1999)calculated that in mature forests growing on mineral soil sites the total number of understorey Scots pine, Norway spruce and birches(from 50 cm height to 6 cm dbh)averaged 3000 trees/ha in the 1950s(see also Sarvas,1944;Laiho et al.,2011)].Plentiful regeneration was measured also later in both Finland(L?hde,1992a,b; L?hde et al.,1999)and Sweden(Lundqvist,1993;Lundqvist and Fridman,1996;Lundqvist and Nilson,2007). Pukkala et al.(2011b)calculated that only 5%of mature Finnish stands had less than 500 understorey trees of pine,spruce and birch per hectare,and 60%of stands had at least 2000 understorey trees at the end of the 1900s.When seedlings shorter than 1.3 m were excluded,almost 50%of stands still had at least 1000 undestorey trees per hectare.

    Due to its ecological characteristics,Norway spruce does not thrive in open areas but regenerates easily and grows well under shade-intolerant pioneer species(Kalela,1949). According to Valkonen(2000),a released spruce understorey may grow almost equally well as a spruce plantation, but the released understorey is more heterogeneous than the plantation would be.In a two-storied birch-spruce stand the volume growth of overstorey birches is often higher than the growth reduction of understorey spruces (Mielik?inen and Valkonen,1995),making the two-storied stand more productive than either of the one-storied stands alone(see also Isom?ki,1979;Bergqvist,1999).

    According to official silvicultural instructions(Anonym, 2006)good-quality spruce understoreys can be used if they cover large enough continuous areas(at least 300–400 m2), are dense and even-sized,and have at least 10-cm long annual shoots.Otherwise,the understorey should be removed,after which the stand should be clear-felled,and a new stand be established by planting.Those criteria,especially the 10-cm rule,are not supported by research,which shows that even a slow-growing stagnant understorey may start growing well(Vuokila,1982;Ferlin,2002).Old and large understorey trees revive slowly(Vaartaja,1951)but may eventually grow as well as young seedlings and saplings(N?slund,1944;Sarvas,1951;Vuokila,1970).According to Kneeshaw et al.(2002),height growth does not respond immediately to increased growth space,but growth reaction occurs first in roots.Typically,it takes 4–5 years for the Norway spruce to fully adapt to the new conditions (Metslaid et al.,2005).The requirement for homogeneity is not relevant in CCF management where uneven size, growth rate and gradual revival of the understorey are benefits rather than shortcomings(L?hde et al.,2010).

    High risk of harvesting damages complicates the management of two-storied stands.According to Mets?teho (K?rn?,2006),a spruce understorey of 6000 seedlings/ha decreases the productivity of harvesting by 5%.Niemist? et al.(2012)found that retaining the understorey spruces decreased the productivity of cutting by 6–17%as compared with clear cutting where the understorey was not considered.Damage to saplings may be as high as 50% (Surakka et al.,2011).Saplings near strip roads are vulnerable to additional injuries caused by hauling and forwarding.However,the high number of spruce stands that originate from advance regeneration demonstrates that the harvesting problems are manageable.A slightly more costly harvesting,as compared to clear-felling,is compensated for by savings in stand establishment costs, immediate stocking,soil protection,aesthetic values and benefits to fauna and flora,amongst other things(Ruel et al.,2000).

    There are some Finnish studies on the revival and height increment of released Norway spruce understoreys(Koistinen and Valkonen,1993;Valkonen,2000)and on the growth of Norway spruce under birch(Mielik?inen and Valkonen,1995)and spruce(Lin et al.,2012)canopies.Although much research has been done in other countries (see Metslaid et al.,2007 for review),quantitative knowledge about the effects of various factors on the growth of understorey Norway spruce,when they are not released,is very limited.More information and better predictive tools would make it possible to evaluate spruce understoreys from the perspective of CCF management,not only as the starting point of a new even-aged stand.

    The aim of this study was to quantify the growth rate of Norway spruce(Picea abies)understorey under different canopies.In the case of abundant understorey the growth rate of the best individuals is relevant for future stand development and for CCF management.Therefore,growth variation among understorey spruces was also analyzed.A regression model was developed for describing the influence of site,canopy structure,and seedling height on the height increment of understorey spruces.The developed model was then used to illustrate the influence of various factors on the height increment of understorey.

    Methods

    Altogether 262 temporary plots were measured in different stands on mineral soils in South Finland.Stands having advance regeneration of spruce were selected following a predefined target distribution of different stand types.A matrix of seven species compositions(pine,spruce,birch,pinespruce,pine-birch,spruce-birch,pine-spruce-birch)and stand densities(<10 m2ha-1,10–30 m2ha-1,>30 m2ha-1)was created.Then,1–3 stands were measured for each class of stand composition and density,separately in two site class categories(herb-rich or better;mesic or poorer)and three different localities.The first 1–3 encountered stands per matrix cell having advance regeneration of spruce were measured in each locality and site category.A selected stand was entered,after which a random direction was chosen and 20 steps were taken to that direction,leading to plot center.

    Forest site type and the degree paludification were assessed in the field.For easier modeling,forest site type was converted into site index(dominant height at 100 years) as follows(Vuokila and V?liaho,1980):OMaT(mesotrophic herb-rich)30 m;OMT(herb-rich)27 m;MT(mesic) 24 m;VT(sub-xeric)21 m;CT(xeric)18 m.Intermediate values(e.g.MT+and MT-,i.e.,fertile and poor mesic)were given intermediate numerical values(25 m and 23 m). Afterwards,latitude of the site was obtained from a geographical map,and altitude of the stand was obtained from an elevation map.The basal area of canopy trees was measured separately for each species using relascope.The minimum,mean and maximum diameters of each species present in the canopy were measured from the relascope plot.Three trees per species were measured for diameter: the smallest,the largest,and the basal area median tree. The basal area median tree was selected visually,following the practices of Finnish compartment inventory.

    At least 20 and at most 30 conifer seedlings and saplings(from 10 cm height to 4.5 cm dbh)closest to the plot center were measured for dbh,height,and height increment of the previous growing season.This resulted in 5811 spruce seedlings and 197 pines(pine seedlings are not analyzed in this study).When plots that had no overstorey trees were also removed from the data,5601 observations remained for modeling(Table 1).About 50%of observations represented mesic growing sites, 38%represented herb-rich and better sites(OMaT and OMT),and 12%represented sub-xeric and xeric sites. Trees having dbh greater than 4.5 cm were defined as canopy trees and the remaining seedlings and saplings were defined as understorey trees.

    A mixed-effects model was fitted to the data using SPSS software.Different transformations and combinations of predictors were carefully analyzed to find a model that describes the influence of site,canopy structure and seedling height on the height increment of understorey spruces.The mixed-effect model included a random plot factor,allowing us to divide the residual variation into between-plot and within-plot components.

    Results

    Model

    The model for the height increment of understorey spruces was as follows:

    where ihijthe annual height increment of spruce seedling i in plot j(cm),h is height of the seedling(cm),SI is site index(m),P is paludification(percent coverage of peat land species among ground vegetation),G is basal area of canopy trees(m2ha-1),GSis basal area of canopy spruces(m2ha-1),D is mean diameter of canopy trees (cm),uj~N(0,σu2)is random plot factor and eij~N(0,σe2) is residual.

    The degree of explained variance(R2)and the square root of the mean of squared errors(RMSE)were calculated from back-transformed non-logarithmic growth predictions(Table 2).The Baskerville(1972)correction was used in prediction.The correction was s2(eij)/2+s2(uj)/2 for the fixed part of the mixed-effect model,and s2(eij)/2 for the full model.The residuals of the models were visually analyzed and they were found to be normally distributed with a constant variance across the ranges of predictors.

    Effect of seedling height and canopy structure

    Height increment of understorey spruces increased as a function of seedling height(Figure 1).When the basal area of canopy trees was 20 m2ha-1,a four-meter high seedling grew,on average,50 cm per year under birch and pine canopies but a one-meter-high seedling grew only 10 cm/year.

    The influence of overstorey tree species was very strong (Figure 1).In a sparse pine or birch stand(basal area 15 m2ha-1)growing on mesic site,the average height increment of a 2-m-high spruce seedlings was 19 cm/yearwhereas the predicted increment was only 11 cm/year in a spruce stand.With a stand basal area of 30 m2ha-1the average growth would be 16 cm/year under pine and birch canopies but only 6 cm/year under a spruce canopy.

    Table 1 Description of modeling data(5601 spruce seedlings)

    Mean diameter of canopy trees was also a significant predictor of height increment.With a given basal area, seedling increment improved with increasing diameter of the canopy trees(Figure 1).The result implies that seedlings grow better under a small number of large trees than under a larger number of smaller trees.

    Table 2 Model parameters.R2 and RMSE have been calculated for the back-transformed non-logarithmic height increment

    Effect of site

    Improving site index(fertility)increased the height increment of understorey spruces but the influence was not strong(Figure 2).Height increment was only 8%better on a herb-rich site(OMT)than on a sub-xeric site(VT).Paludification also affected height increment,the predicted increment being 96%higher in case of no paludification as compared to a site where the coverage of peat land species in the ground vegetation was 100%.Latitude and elevation were not significant predictors,most probably due to their small range of variation in the modeling data(330 km in y coordinate;167 m in elevation).

    Growth variation within spruce understorey

    Within-and between-stand variation in the height increment of understorey spruces was visualized by adding± s(u)or±s(e),i.e.standard deviation of random plot factor(uj)or standard deviation of residual(eij),to the non-logarithmic prediction of the mixed-effect model. Predictions with uj=s(u)are referred to as‘good stand’, uj=?s(u)as‘bad stand’,eij=s(e)as‘good seedling’,and eij=?s(e)as‘bad seedling’.About 1/6 of stands have the plot factor higher than s(u),and another 1/6 of stands have ujis less than-s(u).Similarly,within a stand,the residual is between-s(e)and s(e)for about 2/3 of seedlings.

    The visualizations(Figure 3)show that the variation unexplained by the fixed model predictors is very large, both among stands and among seedlings within a stand. A‘good seedling’grows twice as fast as a“bad seedling”, and the average growth of understorey spruces is twice as much in a“good stand”as compared to a“bad stand”when the site index,paludification,amount of hardwood understorey and canopy structure are the same.

    Discussion and conclusions

    Figure 2 Effect of site fertility and paludification on the height increment of understorey spruces when the stand basal area is 20 m2ha-1,mean diameter of canopy trees is 25 cm and seedling height is 2 m.

    Figure 3 Variation in the height increment of understorey spruces among stands and individuals on mesic site when 50% of canopy basal area is spruce.Good stand,uj=s(u);Bad stand,uj=?s(u);Good seedling,eij=s(e);Bad seedling,eij=?s(e);s(u)is standard deviation of random plot factor(uj);s(e)is standard deviation of residual(eij).

    The study analyzed the dependence of the height increment of understorey spruces on seedling height,site properties and characteristics of the main tree canopy.A regression model was fitted that describes these relationships.This model should not be used to predict future understorey development because the used dataset is insufficient for genuine growth modeling.Height increment was measured in one year only,which may bring bias to the prediction due to the weather-induced annual growth variation.In addition,selection of the plots had subjective elements and measurement of the stand characteristics was not accurate enough for growth modeling purposes.

    Results from this study indicate that 1–2 m tall spruces typically grow 15–20 cm/year under birch and pine canopies but only 5–10 cm under spruce canopies. The growth of understory trees improves when they get taller(Lieffers et al.,1996).Taking into account that spruces planted in open areas often need about ten years to reach breast height,the increment of understorey spruces is reasonably good.Given that height increment increases with increasing tree height the future increment would be slightly higher(usually 3–4%)than the past increment predicted by the model.

    The height increment of the best-growing spruces under sparse pine and birch canopies is not slower than it would be in open area(Valkonen,2000).This conclusion agrees with the result of Lieffers and Stadt(1994)who found that white spruces(Picea glauca)grow equally well under a canopy of trembling aspen(Populus tremula)as in full light if the light transmittance of the canopy is 40%or more.Comeau et al.(2003)found that the height increment of subalpine fir(Abies lasiocarpa)and white spruce was only weakly correlated with the light level under paper birch(Betula papyrifera)overstorey.

    The results of this study showed that understorey spruces grow much faster under pine and birch canopies than under spruce canopies.Increasing percentage of spruce in the main canopy decreases the growth rate of understorey.Clear beneficial effects of birch and pine admixtures on the growth of larger spruces were recently reported(Pukkala et al.2013).Earlier,positive mixture effects have been reported by several authors (Mielik?inen 1985;Pukkala et al.1994;Liang et al.2005 and Pretzsch et al.2010).The practical conclusion that can be drawn from this study is that the density of birch and pine stands is not critically important for understorey development but spruce stands must be sparse for fast understory height increment.

    It is commonly believed that birch is the best canopy for spruce understorey and birch is the best admixture in a spruce stand.However,the current study showed that this is not necessarily the case;if the effects of stand basal area and site fertility are removed,the growth is equally good or even better(Pukkala et al.,2013)under pine than under birch.The common belief of the smallercompetitive effect of birch may be partly due to lower stand basal areas and better site fertility of birch stands, as compared to pine stands.Due to fertile site and low stand basal area,understorey spruces often grow well under birch canopies although the competitive effect of pine would be fairly similar with the same basal area.

    Increasing size of canopy trees improved the height increment of understorey spruces when stand basal area was constant.The result suggests that open space below the crowns of canopy trees may improve the height increment of spruce understorey.When the vertical distance between the main canopy and understorey is large, the shading effect of canopy trees is more uniformly distributed than in stands where the crown bases are near ground(Kuuluvainen and Pukkala,1991).Therefore,the beneficial effect of large tree size in the main canopy may be related to more uniform light conditions for understorey spruces.

    A result relevant for CCF management was the improving height increment with increasing seedling height (Lieffers et al.,1996).This means that height differences among understorey spruces tend to increase when the understoreydevelops.Theunderstoreybecomes uneven-sized even when it is even-aged.Also the high within-stand variation in height increment and the positive temporal autocorrelation of annual growths of consecutive years(Miina,1993;Pasanen,1998;Metslaid et al.,2005)contribute to the differentiation of tree size.

    According to Lin et al.(2012)the height increment of the best 100 understory spruces per hectare decreases from 50 to 5 cm/year when the stand volume of a spruce-dominated stand increases from 50 to 400 m3/ha. The result agrees well with the predictions of our model. When the height of spruce seedling is 1–3 m and the basal area of overstorey birches is less than 20 m2/ha, our model gives fairly similar predictions as an earlier model(Mielik?inen and Valkonen 1995)for two-storied spruce-birch stands.Increasing overstorey basal area decreases height increment more in Mielik?inen and Valkonen(1995)than predicted by the model of this study.It is noteworthy that,similarly to this study,the model of(Mielik?inen and Valkonen 1995)also predicts that increasing tree size in overstorey improves the height increment of spruce understorey if stand basal area remains unchanged.The relationships found in this study are similar to the ones reported earlier for other spruce species(e.g.,Filipescu and Comeau 2007).Bergqvist (1999)found that diameter increment of understorey Norway spruces is reduced more than height increment by the competition caused by birch overstorey.

    In CCF stands where spruce understorey is abundant, its average growth rate is less relevant than the growth rate of the best individuals(e.g.,the best 10 or 30%,depending on the density of understorey).According to this study,the growth of the best 1/6 of seedlings is at least 50%better than the average.In sustainable CCF management,only the removed trees need to be replaced by well-growing understorey trees.If the high thinning removes 333 trees/ha,and the density of understorey is 2000 trees/ha,only 1/6 best growing understorey trees are required to take the place of the removed trees.These understorey trees grow at least 50%faster than a population-averaged model prediction.

    Between-stand variation in height increment was also high.This implies,for example,that bad experiences (slow growth)in one forest do not mean that the growth of spruce understorey is always slow.In the best stands, the growth is at least 50%faster than in an average stand and 100%faster than among the poorest stands.Therefore,overall recommendations about the utilization of spruce understorey in CCF management cannot be given only on the basis of site characteristics and overstorey stand structure.Instead,the manager needs to assess the vitality of the understorey in the field and adapt the management to the specific conditions of each stand.

    Other practical conclusions of the study,from the viewpoint of CCF management,can be summarized as follows.Species composition of the canopy is critically important for the height development of spruce understorey.An admixture of birch and pine should be maintained in the stand as long as possible.If the canopy is a pure spruce stand,it should be thinned to a low basal area for fast understorey development.Between-and within-stand variation in understorey height development is high,which should be taken into account in growth and yield prediction.Using population-averaged model prediction in simulations may lead to biased results and conclusions.

    Competing interests

    The authors declare that they have no competing interests.

    Authors’contributions

    OL and EL collected the data.TP fitted the regression model.All authors

    participated in writing.All authors read and approved the final manuscript.

    Author details

    1Joen Forest Program Consulting,Rauhankatu 41,80100 Joensuu,Finland.2Faculty of Science and Forestry,University of Eastern Finland,P.O.11180101 Joensuu,Finland.

    Received:30 May 2013 Accepted:13 August 2013

    Published:26 February 2014

    Anonym(2006)Hyv?n mets?nhoidon suositukset.Mets?talouden kehitt?miskeskus Tapio.Mets?kustannus OY,Helsinki,p 100

    Anonym(2012)Mets?nk?sittelymenetelmien monipuolistaminen–jatkoty?ryhm?n muistio.Ministry of agriculture and forestry.Memorandum 2012(7):53

    Baskerville GL(1972)Use of logarithmic regression in the estimation of plant biomass.Can J For Res 2(1):49–53

    Bergqvist G(1999)Wood volume yield and stand structure in Norway spruce understorey depending on birch shelterwood density.For Ecol Manage 122 (3):22–229

    Brunner A,Hahn K,Biber P,Skovsgaard JP(2006)Conversion of Norway spruce:A case study in Denmark based on silvicultural scenario modeling.In: Hasenauer H(ed)Sustainable Forest Management.?Springer,Berlin, Heidelberg,pp 343–372

    Comeau PG,Wang JR,Letchford T(2003)Influences of paper birch competition on growth of understory white spruce and subalpine fir following spacing. Can J For Res 33(10):1962–1973

    Ferlin F(2002)The growth potential of understorey silver fir and Norway spruce for uneven-aged forest management in Slovenia.Forestry 75(4):375–383

    Filipescu CN,Comeau PG(2007)Aspen competition affects light and white spruce growth across several sites in western Canada.Can J For Res 37:1701–1713

    Hasenauer H(2006)Concepts within tree growth modeling.In:Hasenauer H(ed) Sustainable Forest Management.?Springer,pp 3–17

    Ilvessalo Y(1956)Suomen mets?t vuosista 1921–1924 vuosiin 1951–1953. Kolmeen valtakunnan metsien inventointiin perustuva tutkimus.Enslish Summary:The forests of Finland from 1921–1924 to 1951–1953.The survey based on three national forest inventories.Comm Inst For Fenn 47(1):1–227

    Isom?ki A(1979)Kuusialikasvoksen vaikutus m?nnik?n kasvuun,tuotokseen ja tuottoon.Summary:The effect of spruce undergrowth on the increment, yield and returns of a pine stand.Folia For 392:13

    Kalela EK(1949)Ecological character of tree species and its relation to silviculture. Acta For Fenn 57(1):1–35

    K?rn? K(2006)Alikasvoksen vaikutukset ensiharvennuspuun korjuuseen. Summary:impact of undergrowth on the harvesting of first-thinning stands. Mets?tehon katsaus 18:4

    Kneeshaw DD,Williams H,Nikinmaa E,Messier C(2002)Patterns of above-and below-ground response of understory conifer release 6 years after partial cutting.Can J For Res 32:255–265

    Koistinen E,Valkonen S(1993)Models for height development of Norway spruce and Scots pine advance growth after release in Southern Finland.Silva Fenn 27(3):179–194

    Kumela H,H?nninen H(2011)Mets?nomistajien n?kemykset mets?nk?sittelymenetelmien monipuolistamisesta.Working Papers of the Finnish Forest Research Institute 203:76

    Kuuluvainen T,Pukkala T(1991)Interaction between canopy architecture and photosynthetically active radiation at different latitudes:simulation experiments and their ecological implications.Naturalia Monspeliensia Special Issue(L'arbre.Biologie et developpement):277–291

    L?hde E(1992a)Regeneration potential of all-sized spruce-dominated stands.In: Hagner M(ed)Silvicultural Alternatives.Proceedings from an Internordic Workshop,June 22–25 1992.The Swedish University of Agricultural Sciences, Department of Silviculture,Reports 35,Ume?,pp 111–116

    L?hde E(1992b)Natural regeneration of all-sized spruce-dominated stands treated by single tree selection.In:Hagner M(ed)Silvicultural Alternatives. Proceedings from an Internordic Workshop,June 22–25 1992.The Swedish University of Agricultural Sciences,Department of Silviculture,Reports 35, Ume?,pp 117–123

    L?hde E,Laiho O,Lin J(2010)Silvicultural alternatives in an uneven-sized forest dominated by Picea abies.J For Res 15:14–20

    L?hde E,Laiho O,Norokorpi Y,Saksa T(1999)Uudistuminen ja kasvatus.In: Luontaisesti syntyneiden sekametsien kehitys ja mets?nhoito.The Finnish Forest Research Institute,Research Notes 719,Helsinki,pp 32–58

    Laiho O,L?hde E,Pukkala T(2011)Uneven-vs.even-aged management in Finnish boreal forests.Forestry 84(5):547–556

    Liang J,Buongiorno J,Monserud RS(2005)Growth and yield of all-aged Douglas-fir/western hemlock forest stands:a matrix model with stand diversity effects.Can J For Res 35:2369–2382

    Lieffers VY,Stadt KJ(1994)Growth of understory Picea glauca,Calamagrostis canadensis,and Epilobium angustifolium in relation to overstory light transmission.Can J For Res 24(6):1193–1198

    Lieffers VJ,Stadt KJ,Navratil S(1996)Age structure and growth of understory white spruce under aspen.Can J For Res 26(6):1002–1007

    Lin CJ,Laiho O,L?hde E(2012)Norway spruce(Picea abies L.)regeneration and growth of understory trees under single-tree selection silviculture in Finland. Eur J For Res 131:683–691

    Lundqvist L(1993)Changes in the stand structure on permanent Picea abies plots managed with single-tree selection.Scand J For Res 8:510–517

    Lundqvist L,Fridman E(1996)Influence of local stand basal area on density and growth of regeneration in uneven-aged Picea abies stands.Scand J For Res 11:364–369

    Lundqvist L,Nilson K(2007)Regeneration dynamics in an uneven-aged virgin Norway spruce forest in northern Sweden.Scand J For Res 22:304–309

    Metslaid M,Ilisson T,Vicente M,Nikinmaa E,J?giste K(2005)Growth of advance regeneration of Norway spruce after clear-cutting.Tree Physiol 25:793–801

    Metslaid M,J?giste K,Nikinmaa E,Moser KW,Porcar-Castell A(2007)Tree variables related to growth response and acclimation of advance regeneration of Norway spruce and other coniferous species after release.For Ecol Manage 250:56–63

    Mielik?inen K(1985)Koivusekoituksen vaikutus kuusikon rakenteeseen ja kehitykseen.Summary:Effect of an admixture of birch on the structure and development of Norway spruce stands.Comm Inst For Fenn 133:1–79

    Mielik?inen K,Valkonen S(1995)Kaksijaksoisen kuusi-koivu-sekametsik?n kasvu. Mets?tieteen aikakauskirja 1995(2):81–97

    Miina J(1993)Residual variation in diameter growth in a stand of Scots pine and Norway spruce.For Ecol Manage 58:111–128

    Mikola P(1966)Alikasvosten merkitys metsien uudistamisessa.Summary: Significance of undergrowth in reforestation.Mets?taloudellinen aikakauslehti 83(1):4–7

    N?slund M(1944)Den gamla norrl?ndska granskogens reaktionsf?rm?ga efter genomhuggning.Referat:Die reaktionsf?higkeit des alten norrl?ndischen fichtenwaldes nach durchlauung.Medd Stat Skogsf?rs?ksanstalt 33:1–194

    Niemist? P,Korpunen H,Laurén A,Salom?ki M,Uusitalo J(2012)Impact and productivity of harvesting while retaining young understorey spruces in final cutting of downy birch.Silva Fenn 46(1):81–97

    Pasanen K(1998)Integrating variation in tree growth into forest planning.Silva Fenn 32(1):11–25

    Pretzsch H,Block J,Dieler J,Dong PH,Kohnle U,Nagel J,Spellmann H,Zingg A (2010)Comparison between the productivity of pure and mixed stands of Norway spruce and beach.Ann For Sci 67:712p1–712p7,DOI:10.1051/forest/ 2010037

    Pukkala T,Vettenranta J,Kolstr?m T,Miina J(1994)Productivity of a mixed Scots pine-Norway spruce stand.Scand J For Res 9:143–153

    Pukkala T,L?hde E,Laiho O(2010)Optimizing the structure and management of uneven-sized stands of Finland.Forestry 83:129–142

    Pukkala T,L?hde E,Laiho O,Salo K,Hotanen J-P(2011a)A multifunctional comparison of even-aged and uneven-aged forest management in a boreal region.Can J For Res 41:851–862

    Pukkala T,L?hde E,Laiho O(2011b)Mets?n jatkuva kasvatus.Joen Forest Program Consulting,Joensuu,229

    Pukkala T,L?hde E,Laiho O(2013)Species interactions in even-and unevenaged boreal forests.J Sust For 32:1–33

    Ruel JC,Messier C,Doucet R,Claveau Y,Comeau P(2000)Morphological indicators of growth response of coniferous advance regeneration to overstory removal in the boreal forest.For Chron 76:633–642

    Sarvas R(1944)Tukkipuun harsintojen vaikutus Etel?-Suomen yksityismetsiin. Referat:Einwirkung der S?gestammplenterungen auf die Privatw?lder Südfinnlands.Comm Inst For Fenn 33(1):1–268

    Sarvas R(1951)Tutkimuksia puolukkatyypin kuusikoista.Summary:Investigations into the spruce stands of Vaccinium type.Comm Inst For Fenn 39(1):1–82

    Schütz J-P,Pukkala T,Donoso PJ,von Gadow K(2012)Historical emergence and current application of CCF.In:Pukkala T,Gadow K(eds)Continuous Cover Forestry.Book Series Managing Forest Ecosystems Vol.24.?Springer Science +Business Media B.V,Dordrecht,pp 1–28

    Surakka H,Sirén M,Heikkinen J,Valkonen S(2011)Damage to saplings in mechanized selection cutting in uneven-aged Norway spruce stands.Scand J For Res 26(3):232–244

    Tahvonen O(2009)Optimal choice between even-and uneven-aged forestry. Natural Res Modelling 22:289–321

    Tahvonen O,Pukkala T,Laiho O,L?hde E,Niinim?ki S(2010)Optimal management of uneven-aged Norway spruce stands.For Ecol Manage 260:106–115

    Vaartaja O(1951)Alikasvosasemasta vapautettujen m?nnyntaimistojen toipumisesta ja merkityksest? mets?nhoidossa.English Summary:On the recovery of released pine advance growth and its silvicultural importance. Acta For Fenn 59(3):1–133

    Valkeap?? A,Paloniemi R,Vainio A,Vehkalahti K,Helkama K,Karppinen H, Kuuluvainen J,Ojala A,Rantala T,Rekola M(2009)Suomen mets?t jamets?politiikka–kansalaisten n?kemyksi?.University of Helsinki,Department of Forest Economics.Reports 55:36

    Valkonen S(2000)Kuusen taimikon kasvattamisen vaihtoehdot Etel?-Suomen kivenn?ismailla:Puhdas kuusen viljelytaimikko,vapautettu alikasvos ja kuusikoivusekataimikko.The Finnish Forest Research Institute,Research Papers 763:83

    Vuokila Y(1956)Etel?-Suomen hoidettujen kuusikoiden kehityksest?.Summary: On the development of managed spruce stands in southern Finland.Comm Inst For Fenn 48(1):1–138

    Vuokila Y(1970)Harsintaperiaate kasvatushakkuissa.Summary:Selection from above in intermediate cuttings.Acta For Fenn 110:1–45

    Vuokila Y(1982)Antakaa luonnon auttaa.Mets? ja Puu 99(12):6–8

    Vuokila Y,V?liaho H(1980)Viljeltyjen havumetsik?iden kasvumallit.Summary: Growth and yield models for conifer cultures in Finland.Comm Inst For Fenn 99(2):271

    Cite this article as:Laiho et al.:Height increment of understorey Norway spruces under different tree canopies.Forest Ecosystems 2014 1:4.

    10.1186/2197-5620-1-4

    *Correspondence:timo.pukkala@uef.fi

    2Faculty of Science and Forestry,University of Eastern Finland,P.O.11180101 Joensuu,Finland

    Full list of author information is available at the end of the article

    ?2014 Laiho et al.;licensee Springer.This is an Open Access article distributed under the terms of the Creative Commons

    Attribution License(http://creativecommons.org/licenses/by/2.0),which permits unrestricted use,distribution,and reproduction in any medium,provided the original work is properly cited.

    Methods:This study used regression analysis to model the height increment of spruce understorey as a function of seedling height,site characteristics and canopy structure.

    Results:An admixture of pine and birch in the main canopy improves the height increment of understorey.When the stand basal area is 20 m2ha-1 height increment is twice as fast under pine and birch canopies,as compared to spruce.Height increment of understorey spruce increases with increasing seedling height.Between-stand and within-stand residual variation in the height increment of understorey spruces is high.The increment of 1/6 fastestgrowing seedlings is at least 50%greater than the average.

    Conclusions:The results of this study help forest managers to regulate the density and species composition of the stand,so as to obtain a sufficient height development of the understorey.In pure and almost pure spruce stands, the stand basal area should be low for a good height increment of the understorey.

    观看av在线不卡| 久久久久久九九精品二区国产| 国内少妇人妻偷人精品xxx网站| 最新中文字幕久久久久| 午夜免费鲁丝| 国产精品蜜桃在线观看| 国产精品不卡视频一区二区| 国产免费一级a男人的天堂| 国产一区二区在线观看日韩| 在线观看美女被高潮喷水网站| 日韩亚洲欧美综合| 久久精品国产亚洲av天美| 亚洲在久久综合| 亚洲精品亚洲一区二区| 一级片'在线观看视频| 丰满少妇做爰视频| 乱码一卡2卡4卡精品| 日韩国内少妇激情av| 亚洲av国产av综合av卡| 深爱激情五月婷婷| 超碰97精品在线观看| 亚洲精品一区蜜桃| 国产精品一区二区三区四区免费观看| 日本午夜av视频| 中国国产av一级| 丝袜脚勾引网站| 国产色爽女视频免费观看| 亚洲av二区三区四区| 婷婷色综合大香蕉| av在线播放精品| 国产69精品久久久久777片| av视频免费观看在线观看| av在线蜜桃| 日韩中文字幕视频在线看片 | 午夜免费男女啪啪视频观看| 汤姆久久久久久久影院中文字幕| 偷拍熟女少妇极品色| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 天堂中文最新版在线下载| 亚洲国产精品专区欧美| 国产免费一区二区三区四区乱码| 成年av动漫网址| 日韩,欧美,国产一区二区三区| 汤姆久久久久久久影院中文字幕| 99九九线精品视频在线观看视频| 少妇熟女欧美另类| 亚洲性久久影院| 日本与韩国留学比较| 激情 狠狠 欧美| 国产精品一区www在线观看| 男人添女人高潮全过程视频| 亚洲精品aⅴ在线观看| 久久久久久伊人网av| 久久韩国三级中文字幕| 欧美成人午夜免费资源| 欧美区成人在线视频| 91aial.com中文字幕在线观看| 在线观看国产h片| 51国产日韩欧美| 男女无遮挡免费网站观看| 欧美区成人在线视频| 亚洲综合色惰| 中文乱码字字幕精品一区二区三区| 少妇精品久久久久久久| 欧美区成人在线视频| 色婷婷av一区二区三区视频| 插阴视频在线观看视频| 亚洲三级黄色毛片| 亚洲精品自拍成人| 黄色怎么调成土黄色| 超碰97精品在线观看| 在线精品无人区一区二区三 | 亚洲av不卡在线观看| 伊人久久精品亚洲午夜| 九九久久精品国产亚洲av麻豆| 高清不卡的av网站| 老师上课跳d突然被开到最大视频| 最近中文字幕2019免费版| 丰满人妻一区二区三区视频av| 亚洲精品国产av成人精品| 直男gayav资源| 色网站视频免费| 亚洲综合色惰| 人妻少妇偷人精品九色| 欧美人与善性xxx| 少妇人妻一区二区三区视频| 91精品国产国语对白视频| 中文字幕制服av| 国产欧美亚洲国产| 日韩三级伦理在线观看| 亚洲精品成人av观看孕妇| 国产又色又爽无遮挡免| 亚洲真实伦在线观看| 国产综合精华液| 国产精品麻豆人妻色哟哟久久| 日本免费在线观看一区| 国产亚洲午夜精品一区二区久久| 欧美3d第一页| 美女脱内裤让男人舔精品视频| 色婷婷av一区二区三区视频| 51国产日韩欧美| 久久亚洲国产成人精品v| 国产精品偷伦视频观看了| 天堂8中文在线网| 亚洲精品456在线播放app| 亚洲国产日韩一区二区| 国产 一区精品| 夫妻午夜视频| 麻豆国产97在线/欧美| 大陆偷拍与自拍| 日本av免费视频播放| 色吧在线观看| 哪个播放器可以免费观看大片| 国产视频首页在线观看| 一本—道久久a久久精品蜜桃钙片| 国产视频首页在线观看| 精品一区二区三区视频在线| 中文字幕免费在线视频6| 欧美激情国产日韩精品一区| 一级毛片黄色毛片免费观看视频| 午夜福利在线在线| a级毛片免费高清观看在线播放| 日本一二三区视频观看| 日本黄色片子视频| 最近最新中文字幕大全电影3| 日本黄色日本黄色录像| 亚洲精品自拍成人| 日本爱情动作片www.在线观看| 国产久久久一区二区三区| av.在线天堂| 久久人人爽人人爽人人片va| 秋霞在线观看毛片| 国产在线免费精品| 国产成人精品婷婷| 亚洲国产成人一精品久久久| 成年美女黄网站色视频大全免费 | 22中文网久久字幕| 人人妻人人添人人爽欧美一区卜 | 男人舔奶头视频| 天美传媒精品一区二区| 国产在线一区二区三区精| 亚洲国产高清在线一区二区三| 高清av免费在线| 少妇 在线观看| 亚洲在久久综合| www.色视频.com| 日韩免费高清中文字幕av| 精品视频人人做人人爽| 国产亚洲5aaaaa淫片| 免费观看性生交大片5| av黄色大香蕉| 精品人妻偷拍中文字幕| 校园人妻丝袜中文字幕| 黄色欧美视频在线观看| 少妇人妻一区二区三区视频| 亚洲精品日本国产第一区| 久久毛片免费看一区二区三区| 日本午夜av视频| 国产在视频线精品| 日日摸夜夜添夜夜添av毛片| 一本—道久久a久久精品蜜桃钙片| 国产美女午夜福利| 人人妻人人澡人人爽人人夜夜| 18禁裸乳无遮挡动漫免费视频| 少妇 在线观看| 日韩,欧美,国产一区二区三区| 丰满人妻一区二区三区视频av| 久久精品夜色国产| 国产成人精品福利久久| 国产精品爽爽va在线观看网站| 久久久久久久大尺度免费视频| 99视频精品全部免费 在线| 亚洲国产欧美在线一区| 国产欧美日韩一区二区三区在线 | 成人高潮视频无遮挡免费网站| 色婷婷av一区二区三区视频| 免费看av在线观看网站| 永久网站在线| 久久韩国三级中文字幕| 99久久精品一区二区三区| 亚洲精品日韩av片在线观看| 一区二区三区乱码不卡18| h视频一区二区三区| 成人黄色视频免费在线看| 最后的刺客免费高清国语| 人妻 亚洲 视频| av免费观看日本| 亚洲精品国产av成人精品| 国产精品国产av在线观看| 免费观看a级毛片全部| 亚洲va在线va天堂va国产| av专区在线播放| 国产一区二区在线观看日韩| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 亚洲精品国产成人久久av| 纯流量卡能插随身wifi吗| 亚洲国产高清在线一区二区三| 中国国产av一级| 极品少妇高潮喷水抽搐| 国产午夜精品久久久久久一区二区三区| 毛片一级片免费看久久久久| 亚洲精华国产精华液的使用体验| av又黄又爽大尺度在线免费看| 青春草国产在线视频| 老司机影院毛片| 91久久精品国产一区二区三区| 午夜精品国产一区二区电影| 中文字幕av成人在线电影| 欧美zozozo另类| 成人美女网站在线观看视频| 国产精品av视频在线免费观看| 黄色视频在线播放观看不卡| 99视频精品全部免费 在线| 舔av片在线| 国产大屁股一区二区在线视频| 久久久午夜欧美精品| 18禁在线无遮挡免费观看视频| 久久这里有精品视频免费| 99热国产这里只有精品6| 男人和女人高潮做爰伦理| 国产成人一区二区在线| 各种免费的搞黄视频| 国产视频内射| 成人毛片a级毛片在线播放| 一二三四中文在线观看免费高清| 欧美老熟妇乱子伦牲交| 有码 亚洲区| 久久久久久久国产电影| 日产精品乱码卡一卡2卡三| 精品久久久精品久久久| 国产伦在线观看视频一区| 中国三级夫妇交换| 亚洲精品国产av蜜桃| 最近最新中文字幕免费大全7| 久久久久国产精品人妻一区二区| 成人亚洲欧美一区二区av| 亚洲国产欧美人成| 丝袜喷水一区| 黄片wwwwww| 日韩 亚洲 欧美在线| h日本视频在线播放| 一级毛片电影观看| 国产成人精品一,二区| 黄片wwwwww| 狂野欧美激情性bbbbbb| 欧美xxⅹ黑人| 能在线免费看毛片的网站| 人人妻人人爽人人添夜夜欢视频 | 久久ye,这里只有精品| 新久久久久国产一级毛片| 亚洲美女黄色视频免费看| 国产老妇伦熟女老妇高清| 欧美激情极品国产一区二区三区 | 亚洲在久久综合| 少妇的逼好多水| 老熟女久久久| 精品国产乱码久久久久久小说| 国产成人aa在线观看| 国产av精品麻豆| 深爱激情五月婷婷| 18+在线观看网站| 一本—道久久a久久精品蜜桃钙片| 嘟嘟电影网在线观看| 成人毛片a级毛片在线播放| 久久午夜福利片| 麻豆成人av视频| 有码 亚洲区| 老师上课跳d突然被开到最大视频| 成人黄色视频免费在线看| 国产美女午夜福利| 交换朋友夫妻互换小说| 我要看黄色一级片免费的| 欧美激情极品国产一区二区三区 | 日本黄色日本黄色录像| 18禁在线播放成人免费| 国产精品人妻久久久久久| 亚洲自偷自拍三级| 乱码一卡2卡4卡精品| 精品久久久久久久久亚洲| 免费黄色在线免费观看| 国产女主播在线喷水免费视频网站| 亚洲成人手机| 成人高潮视频无遮挡免费网站| 91久久精品国产一区二区三区| 国产色爽女视频免费观看| 久久精品久久久久久噜噜老黄| 亚洲欧美日韩另类电影网站 | 一二三四中文在线观看免费高清| 不卡视频在线观看欧美| 一级二级三级毛片免费看| 亚洲丝袜综合中文字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 看十八女毛片水多多多| av在线播放精品| 熟女av电影| 免费久久久久久久精品成人欧美视频 | 婷婷色av中文字幕| 亚洲av欧美aⅴ国产| 丰满少妇做爰视频| 免费av不卡在线播放| 久久精品久久精品一区二区三区| 亚洲,欧美,日韩| 久久久国产一区二区| videos熟女内射| 精品国产乱码久久久久久小说| 高清欧美精品videossex| 国产又色又爽无遮挡免| 国产高潮美女av| 身体一侧抽搐| 国产国拍精品亚洲av在线观看| 国产精品嫩草影院av在线观看| 精品久久久噜噜| 婷婷色麻豆天堂久久| 亚洲国产精品成人久久小说| 边亲边吃奶的免费视频| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 日日摸夜夜添夜夜爱| 伦精品一区二区三区| 亚洲精品国产成人久久av| 日韩电影二区| 国产亚洲欧美精品永久| 少妇精品久久久久久久| 亚州av有码| 丝瓜视频免费看黄片| 精品人妻视频免费看| 51国产日韩欧美| 纵有疾风起免费观看全集完整版| 精品一区二区三区视频在线| 日韩欧美 国产精品| 亚洲综合精品二区| 插逼视频在线观看| 欧美成人一区二区免费高清观看| 黄色日韩在线| 99久久综合免费| 看免费成人av毛片| 久久6这里有精品| 一个人看视频在线观看www免费| 国产视频首页在线观看| 在线播放无遮挡| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 久久久国产一区二区| 成人美女网站在线观看视频| 国产伦精品一区二区三区视频9| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91 | 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 少妇人妻 视频| 精品一区二区三区视频在线| 国产在线男女| 欧美日韩精品成人综合77777| 春色校园在线视频观看| 国产伦理片在线播放av一区| 亚洲国产精品专区欧美| 久久精品国产亚洲网站| 国产伦在线观看视频一区| 我的女老师完整版在线观看| 97精品久久久久久久久久精品| 成人高潮视频无遮挡免费网站| 高清av免费在线| 亚洲av不卡在线观看| 欧美精品一区二区大全| 一区二区三区四区激情视频| 成人综合一区亚洲| 一区在线观看完整版| 下体分泌物呈黄色| 国产白丝娇喘喷水9色精品| av一本久久久久| av线在线观看网站| 蜜臀久久99精品久久宅男| 欧美日韩国产mv在线观看视频 | 偷拍熟女少妇极品色| 亚洲国产毛片av蜜桃av| 欧美bdsm另类| 国产乱来视频区| 看免费成人av毛片| 深爱激情五月婷婷| 精品视频人人做人人爽| 久久精品久久久久久噜噜老黄| 水蜜桃什么品种好| 成人亚洲精品一区在线观看 | 久久精品久久久久久久性| 肉色欧美久久久久久久蜜桃| 欧美zozozo另类| 少妇熟女欧美另类| 91久久精品电影网| 秋霞在线观看毛片| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| 美女主播在线视频| 欧美高清性xxxxhd video| 久久精品久久久久久噜噜老黄| 国产91av在线免费观看| 成人国产麻豆网| 国产在线一区二区三区精| 色婷婷av一区二区三区视频| 日本欧美视频一区| 一级a做视频免费观看| 国产欧美日韩精品一区二区| 中国国产av一级| 国产成人一区二区在线| 久久 成人 亚洲| 一个人看的www免费观看视频| 国产成人a∨麻豆精品| 看十八女毛片水多多多| 高清欧美精品videossex| 亚洲欧美一区二区三区黑人 | 最近最新中文字幕大全电影3| 成人国产av品久久久| 国产成人免费无遮挡视频| 性高湖久久久久久久久免费观看| 18禁裸乳无遮挡免费网站照片| 夜夜看夜夜爽夜夜摸| 中文乱码字字幕精品一区二区三区| 91久久精品电影网| 国产亚洲一区二区精品| 少妇丰满av| 精品亚洲乱码少妇综合久久| 国产男人的电影天堂91| kizo精华| 久久99热这里只频精品6学生| 免费大片黄手机在线观看| 91精品伊人久久大香线蕉| 欧美+日韩+精品| 中国三级夫妇交换| 干丝袜人妻中文字幕| 午夜福利在线观看免费完整高清在| av又黄又爽大尺度在线免费看| 欧美3d第一页| 人体艺术视频欧美日本| 色视频在线一区二区三区| 赤兔流量卡办理| 日本wwww免费看| a级毛色黄片| 国产免费福利视频在线观看| 好男人视频免费观看在线| 国产伦在线观看视频一区| 色哟哟·www| 高清在线视频一区二区三区| 欧美日韩精品成人综合77777| 人妻少妇偷人精品九色| 欧美三级亚洲精品| 99久久综合免费| 国产精品.久久久| 国产高潮美女av| 高清不卡的av网站| 久久精品熟女亚洲av麻豆精品| 三级国产精品欧美在线观看| 91在线精品国自产拍蜜月| 亚洲,欧美,日韩| 视频中文字幕在线观看| 久久久久久久久久成人| 久久人人爽av亚洲精品天堂 | 纵有疾风起免费观看全集完整版| 国产精品嫩草影院av在线观看| 亚洲第一区二区三区不卡| 青青草视频在线视频观看| 亚洲欧美成人精品一区二区| 欧美国产精品一级二级三级 | 大又大粗又爽又黄少妇毛片口| 久久人妻熟女aⅴ| 最后的刺客免费高清国语| av天堂中文字幕网| 亚洲婷婷狠狠爱综合网| 性高湖久久久久久久久免费观看| 天美传媒精品一区二区| 国产成人精品一,二区| 大码成人一级视频| videos熟女内射| 高清黄色对白视频在线免费看 | 国产精品久久久久成人av| 九草在线视频观看| 又粗又硬又长又爽又黄的视频| 欧美成人午夜免费资源| 亚洲欧美精品专区久久| 男人爽女人下面视频在线观看| 久久国产精品大桥未久av | 啦啦啦视频在线资源免费观看| 91狼人影院| 一级a做视频免费观看| 99热这里只有是精品50| 美女国产视频在线观看| 午夜精品国产一区二区电影| 91aial.com中文字幕在线观看| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在| 亚洲色图综合在线观看| 人妻一区二区av| 在线观看国产h片| h视频一区二区三区| 97在线视频观看| 男人添女人高潮全过程视频| 国产成人精品婷婷| 久久国产精品男人的天堂亚洲 | 精品一区在线观看国产| 日韩视频在线欧美| 亚洲精品色激情综合| 99热这里只有是精品50| 亚洲一级一片aⅴ在线观看| 久久精品国产a三级三级三级| 日本av免费视频播放| 99热6这里只有精品| 在线看a的网站| 春色校园在线视频观看| 国产精品一区二区三区四区免费观看| 美女xxoo啪啪120秒动态图| 久久精品熟女亚洲av麻豆精品| 亚洲av男天堂| 内地一区二区视频在线| 26uuu在线亚洲综合色| 在线亚洲精品国产二区图片欧美 | 精品视频人人做人人爽| 精品久久久久久久久av| tube8黄色片| 97超碰精品成人国产| 成人国产av品久久久| 久久精品国产自在天天线| 国产乱来视频区| 国产成人精品福利久久| 久久国产亚洲av麻豆专区| 国产欧美亚洲国产| 91久久精品国产一区二区成人| 三级经典国产精品| 97热精品久久久久久| 亚洲电影在线观看av| 日韩不卡一区二区三区视频在线| 在线亚洲精品国产二区图片欧美 | 在线观看免费视频网站a站| 久久久精品94久久精品| 久久亚洲国产成人精品v| 国产综合精华液| 黄色配什么色好看| av在线播放精品| 久久女婷五月综合色啪小说| 一区二区三区乱码不卡18| 精品熟女少妇av免费看| 国产精品嫩草影院av在线观看| 91久久精品电影网| 免费在线观看成人毛片| 汤姆久久久久久久影院中文字幕| 久久精品久久久久久久性| 中文字幕人妻熟人妻熟丝袜美| 美女主播在线视频| 国产精品蜜桃在线观看| 日本wwww免费看| 久久精品国产亚洲av天美| 成人国产麻豆网| 亚洲美女搞黄在线观看| 亚洲怡红院男人天堂| 丝瓜视频免费看黄片| 久久久亚洲精品成人影院| 亚洲精品成人av观看孕妇| 街头女战士在线观看网站| 天美传媒精品一区二区| 老熟女久久久| 午夜激情久久久久久久| 老师上课跳d突然被开到最大视频| 亚洲成人中文字幕在线播放| 欧美丝袜亚洲另类| 国内揄拍国产精品人妻在线| 丝袜喷水一区| 激情 狠狠 欧美| 网址你懂的国产日韩在线| 在线观看免费高清a一片| 人人妻人人澡人人爽人人夜夜| 18禁裸乳无遮挡免费网站照片| 岛国毛片在线播放| 国产亚洲精品久久久com| 日产精品乱码卡一卡2卡三| 一级黄片播放器| 26uuu在线亚洲综合色| 亚洲精品国产av蜜桃| 美女国产视频在线观看| 国产免费又黄又爽又色| 国产精品一区二区在线不卡| 岛国毛片在线播放| 亚洲国产最新在线播放| 三级国产精品欧美在线观看| 久久精品国产亚洲网站| 少妇人妻 视频| 男人添女人高潮全过程视频| 噜噜噜噜噜久久久久久91| 五月开心婷婷网| 一本一本综合久久| 久久精品人妻少妇| 亚洲成人一二三区av| 精品国产乱码久久久久久小说| 高清在线视频一区二区三区| 又粗又硬又长又爽又黄的视频| 国产又色又爽无遮挡免| av在线播放精品| 国产熟女欧美一区二区| 婷婷色综合www| 日韩av在线免费看完整版不卡| 男人狂女人下面高潮的视频| 精品久久久久久电影网| 老师上课跳d突然被开到最大视频| 国产精品一及| 国产大屁股一区二区在线视频| 黑人高潮一二区| 国产淫语在线视频| 噜噜噜噜噜久久久久久91| 久久久久久久久久久免费av| 男男h啪啪无遮挡| 边亲边吃奶的免费视频| 极品教师在线视频| 最近最新中文字幕免费大全7| 美女中出高潮动态图| 亚洲av福利一区| 国产成人a区在线观看|