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

    Effect of disturbance on biomass, production and carbon dynamics in moist tropical forest of eastern Nepal

    2016-11-24 05:36:01TilakPrasadGautamandTejNarayanMandal
    Forest Ecosystems 2016年3期

    Tilak Prasad Gautamand Tej Narayan Mandal

    Effect of disturbance on biomass, production and carbon dynamics in moist tropical forest of eastern Nepal

    Tilak Prasad Gautam1*and Tej Narayan Mandal2

    Background:Forest biomass is helpful to assess its productivity and carbon(C)sequestration capacity.Several disturbance activities in tropical forests have reduced the biomass and net primary production(NPP)leading to climate change.Therefore,an accurate estimation of forest biomass and C cycling in context of disturbances is required for implementing REDD(Reducing Emissions from Deforestation and Forest Degradation)policy.

    Methods:Biomass and NPP of trees and shrubs were estimated by using allometric equations while herbaceous biomass was estimated by harvest method.Fine root biomass was determined from soil monolith.The C stock in vegetation was calculated by multiplying C concentration to dry weight.

    Results:Total stand biomass(Mg·ha-1)in undisturbed forest stand(US)was 960.4 while in disturbed forest stand(DS)it was 449.1.The biomass(Mg·ha-1)of trees,shrubs and herbs in US was 948.0,4.4 and 1.4, respectively,while in DS they were 438.4,6.1 and 1.2,respectively.Total NPP(Mg·ha-1·yr-1)was 26.58 (equivalent to 12.26 Mg C·ha-1·yr-1)in US and 14.91(6.88 Mg C·ha-1·yr-1)in DS.Total C input into soil through litter plus root turnover was 6.78 and 3.35 Mg·ha-1·yr-1in US and DS,respectively.

    Conclusions:Several disturbance activities resulted in the significant loss in stand biomass(53%),NPP (44%),and C sequestration capacity of tropical forest in eastern Nepal.The net uptake of carbon by the vegetation is far greater than that returned to the soil by the turnover of fine root and litter.Therefore,both stands of present forest act as carbon accumulating systems.Moreover,disturbance reflects higher C emissions which can be reduced by better management.

    Tropical forest,Disturbance,Biomass,Production,Carbon cycling,Nepal

    Background

    Tropical forests play an important role in the global C cycle.They contain about 55%of global forest C (Pan et al.2011)and account for 34%of terrestrial gross primary production(Beer et al.2010).Carbon is stored in forests mainly in biomass and soils.The C dynamics of a forest reflects environmental conditions such as climate,soil,structure,nutrient availability and disturbance(Chave et al.2001).Forest disturbances often lead to the changes in species composition,structure,stand biomass,productivity and C cycling.Therefore,a small disruption in tropical forests might result in a significant change in the global C cycle.

    Biomass and production are important parameters for understanding the functioning of a forest ecosystem.Studies on biomass help to assess the effect of disturbances on productivity,C dynamics,nutrient cycling and stability of forest stands.Primary production is usually regulated by the availability of nutrients.Decomposed fine roots and aboveground litter are the sources of nutrients in the soil.The contribution of fine roots in the C and nutrients input to soil is equivalent or even higher to that from leaf litter in tropical moist forests(Roderstein et al.2005).

    Precise quantification of forest biomass,production,and C stock demands the careful estimation of boththeabovegroundandbelowgroundaspects. Despite this requirement,most of the past studies in tropical forests are found to be limited to analyzing the aboveground systems(Chave et al.2008;Djomo et al.2011;Malhi et al.2011;Doughty et al.2013; Ngo et al.2013;Girardin et al.2014;Malhi et al. 2014)and only a few studies concern to the belowground aspects(Ibrahima et al.2010;Powers and Perez-Aviles 2013;Noguchi et al.2014).

    In Nepal,only few studies have been conducted regarding the biomass,production,and C sequestration in tropical forests(Mandal 1999;Baral et al. 2009).More information is required to understand the structure and functioning of the forests.Present studywasdesignedintropicalforestofeastern Nepal with the following specific objectives:a)to estimate the biomass and production of herb,shrub and tree species;b)to find out the effect of disturbance on biomass,production and C dynamics.

    Methods

    Study area

    The study was conducted in a Sal(Shorea robusta) dominated moist tropical forest of Sunsari District, eastern Nepal(latitude 26°41’N to 26°50’N and longitude 87°09’E to 87°21’E),within the altitude range of 220 to 370 m above msl(Fig.1).The forest lies in the catchment area of Koshi River,one of the largest rivers in Nepal.The total area occupied by the forest is 11394 ha.

    The central part(core area)of the forest is relatively undisturbed,while the peripheral part is affected by disturbance activities as removal of timber, livestock grazing,fuel-wood and litter collection,tree lopping,removal of poles for house-hold constructions and forest fires.These disturbances have caused deforestation,forest fragmentation and degradation and subsequent invasion of exotic species like Mikania micrantha,Lantana camara,and Chromolaena odorata,which adversely affect the native plant diversity in the forest.

    The climate is tropical and monsoon type with three distinct seasons:dry and warm summer(March to May), wet and warm rainy(June to October),and dry and cool winter(November to February).The mean monthly minimum and maximum air temperature during 2005–2014 ranged from 10.9°C to 25.3°C and 22.6°C to 33.2°C,respectively.The average annual rainfall for the period was 1998.6 mm(Fig.2).Relative humidity was higher in rainy season with highest value in August(92%).

    The forest is bordered by the Siwalik hill in the north and the Gangetic alluvial plains in the south.The area has been formed from soft erodible sediments of the Siwaliks and is characterized by the presence of boulder beds mixed with sand,silt,clay imparting a porous nature.The soil mainly consists of deep alfisols.

    The forest tree layer is dominated by Sal.Other main associates are Haldina cordifolia,Lagerstroemia parviflora,and Terminalia alata.Clerodendron viscosum and Murraya koenigii are some of the main shrub species while Chromolaena odorata and Achyranthes aspera are dominant herbs.

    Plant biomass estimation

    The forest was divided into 2 parts:i)the relatively undisturbed core area(treated as undisturbed stand; US),and ii)the disturbed peripheral area(treated as disturbed stand;DS).Because of the following characters the core area was considered as US:canopy area ranged between 153.8 and 226.9 m2,crown cover ranged from 70–80%,density of the tree was 466.4 individuals?ha?1and tree stumps were absent. The peripheral part of the forest is disturbed because of following characters:canopy area ranged between 28.3 and 77.5 m2,crown cover ranged from 30–40%,tree density was 234.3 individuals·ha?1and tree stump density was 70 stumps·ha?1.

    The tree biomass contained within US was estimated using 35 randomly established sample plots of 20 m×20 m size determined by species-area curve method.Same size was applied for DS as well.For the estimation of shrub and herb biomass nested quadrats of 5 m×5 m and 1 m×1 m were used,respectively.Girths of all the trees(>10 cm gbh)at breast height(1.37 cm above the soil)and of shrubs (10 cm above the ground level)present in each of the sample plots were measured.Density(individual per hectare)and basal area(girth2·(4π)?1)of the trees within plots were determined.Biomass of trees wasestimatedbyusinggirth:biomassallometric equations(Singh and Singh 1992;Mandal 1999).For estimating coarse root biomass,root:shoot ratio of 0.21 proposed for lowland tropical forests was used (Malhi et al.2009;Malhi et al.2014).

    The aboveground herbaceous biomass present in the sampling plots was harvested twice,in summer (May 2012)and at the end of rainy season(September 2012).Summer and rainy season values were averaged to obtain annual mean herbaceous biomass.

    Fine roots(<5 mm diameter)were collected from seventy randomly established sample plots,thirty five each in US and DS.Fine root biomass(FRB)was determined from soil monolith(10 cm×10 cm× 30 cm),divided into two depth ranges(upper:0–15 cm and lower:15–30 cm)at each sample plot in summer(May 2011),rainy(September 2011)andwinter(January 2012)seasons.Summer,rainy and winter season values were averaged to obtain annual mean FRB.

    Litterfall and litter mass

    For the estimation of litterfall one litter trap(1 m× 1 m)was fixed on the forest floor at each of the seventysamplingplots.Litterfallwascollectedat monthly intervals from April 2011 to March 2012 and categorized into leaf and non-leaf components. Litter mass accumulated at each sampling plot was collected once every season from one 1 m×1 m plot. The turnover rate of litter was calculated according to Jenny et al.(1949).

    Net production

    Using the allometric equations,the aboveground biomass (AGB)of different components of marked trees/shrubs in permanent sampling plots was computed for 2011(B1) and 2012(B2)from respective girth measurements.The net changes in biomass(ΔB=B2–B1)of components yielded annual biomass increments which were summed to get the net AGB accretion in the trees/shrubs.Aboveground herbaceous net production was estimated as thedifferences between maximum and minimum biomass values through the year.

    The annual leaf fall was added to the foliage biomass accumulation to represent leaf production,and wood and miscellaneous litterfall values were added to the biomass accumulation in twigs(Singh et al. 1994).Coarse root production was estimated as a fraction of stem production(Malhi et al.2009;Malhi et al.2014).Fine root production(FRP)was estimated as the differences between maximum(rainy season values)and minimum(summer season values) biomass values.The fine root turnover was calculated as a ratio of its production and annual mean biomass(Srivastava et al.1986).

    Carbon estimation in vegetation and litter

    Samples of different tree,shrub and herb(aboveground)components of all species were collected from each sampling plot.The litter(leaf and nonleaf)and fine root(<2 and 2–5 mm diameter)samples were also collected.All the samples were oven dried at 80°C to constant weight,powdered and used separately for C analyses.

    Carbon present in plant materials was estimated by ash content method.Carbon concentrations were assumed to be approximately 50%of ash free weight (McBrayer and Cromack 1980).In this method oven dried plant components(stem,branch,root,leaf, and litter)were burnt separately in electric furnace at 400°C.Ash content(inorganic elements in the form of oxides)left after burning was weighed and carbon concentration was calculated by using the following equation:

    The C stock in vegetation was calculated by multiplying C concentration(a conversion factor)to dry weight. A conversion factor of 0.470 was used to convert aboveground and coarse root biomass to C.To obtain C stock in fine roots a conversion factor of 0.430 and 0.455 was used for<2 and 2–5 mm size class,respectively.The conversion factors of 0.453 and 0.468 were used for leaf and non-leaf litter,respectively.

    Soil sampling and carbon estimation

    Soil samples were collected from 70 randomly selected pots(35 each in US and DS).Soil was collected from three pits(10 cm×10 cm×30 cm each) from each plot.The soils of three pits were mixed and pooled as one replicate(Singh et al.2001).Bulk density(BD)was determined by using a metallic tube while soil organic carbon(SOC)by dichromate oxidation method(Kalembasa and Jenkinson 1973).Carbon content in soil(Mg C·ha?1soil)was calculated using the formula:

    where SOC is soil organic carbon in%and BD is bulk density(g·cm?3).

    Statistical analysis

    Statistical tests were carried out in MS excel 2007 and SPSS(IBM Statistics,ver.20)packages.All the data were normally distributed except that of basal area of trees which were log transformed before analysis.Pearson’s correlation coefficients were obtained between normal and linear data of biomass,net production,stand density and basal area of trees.Students T-test was performed to compare the mean biomass and production of trees between US and DS.

    Results

    Plant biomass

    The biomass and C stocks of trees,shrubs,herbs and fine roots are summarized in Tables 1 and 2.The estimated total biomass for the US was 960.4 Mg·ha–1(equivalent to 452.06 Mg C·ha–1),while for DS it was 449.1 Mg·ha–1(equivalent to 211.33 Mg C·ha–1).

    The total biomass(Mg·ha–1)of the tree layer in US and DS was found to be 948.0 and 438.4,respectively. The tree biomass in the US and DS was significantly different(P<0.001 at 95%confidence level).Of the total tree biomass,83%was aboveground and 17%belowground(excluding fine root)in both the stands.The maximum contribution to tree AGB was made by bole in both stands(64%each)and minimum by foliage(2% each in both stands).

    The shrub biomass(Mg·ha–1)showed increasing trend from 4.4 at US to 6.1 at DS.Contribution of shrub biomass to the total stand biomass was higher in DS (1.4%)than US(0.5%).The aboveground herbaceous biomass contributed 0.1%in US and 0.3%in DS.The fine roots shared 0.7–0.8%to total stand biomass.The annual FRB was lower in DS by 49.5%.

    Table 1 Oven-dried stand biomass(Mg·ha?1)and carbon stock(Mg C·ha?1)in moist tropical forest of eastern Nepal

    Litterfall

    The total annual litterfall(Mg·ha?1·yr?1)was 11.8 in US and 5.4 in DS.Leaves accounted for 69%(US) to 76%(DS)of total litterfall while non-leaf litter formed the rest.The total litter mass(Mg·ha?1) comprised 6.7 in US and 3.6 in DS.The turnover rate(per yr)of the total litter ranged from 0.79 in DS to 0.83 in US.

    Net production

    The total NPP(Mg·ha–1·yr–1)of the forest was 26.6 (equivalent to an annual C sequestration of 12.26 Mg C·ha–1·yr–1)in US and 14.9(i.e.6.88 Mg C·ha–1·yr–1)in DS(Table 3).Among the different life forms:tree,shrub, and herb comprised 72%,2%,and 6%of NPP in US and 67%,5%,and 9%in DS,respectively;while rest 20%NPP in US and 19%in DS were contributed by stand fine root.

    The contribution in NPP by different components of trees was in the order leaf>bole>twig>coarse root>branch,in both stands(Table 4).The most prominent C sink was found to be leaf that accounted for 42–43%of the total NPP.For shrub,maximum contribution to NPP was made by stem in both stands and minimum by foliage.The NPP of herbs(aboveground)ranged between 1.3 and 1.7 Mg·ha–1·yr–1in DS and US,respectively(Table 3).

    Total aboveground NPP(Mg·ha–1·yr–1)combining tree,shrub and herb were 19.93 in US and 10.97 in DS. In belowground parts,percentage allocation of NPP was up to 25%of total NPP(26.58 Mg·ha–1·yr–1)in US and 26%(3.94 Mg·ha–1·yr–1)in DS(Table 5).Of this,both tree and shrub’s coarse root contributed about 21 and 26%to total belowground NPP in US and DS,respectively;whereas fine root comprised nearly 79%of belowground NPP in US and 74%in DS.

    Table 2 Carbon stock(Mg C·ha?1)estimates for undisturbed and disturbed stands of moist tropical forest of eastern Nepal

    Consideringbothdepths,averageturnoverrate (per yr)of smaller fine roots(0–2 mm diameter) was 0.75 in US and 0.86 in DS,whereas it was 0.74 and 0.76 for larger size class(2–5 mm)in US and DS,respectively.Thebiomassaccumulationratio (biomass/net production)decreased to 30 due to the consequences of forest disturbance in DS,while this value was higher(36)in US.Net production of trees was significantly higher in US than DS(P<0.001). Positive correlations were found between stand density,basal area,biomass and NPP of trees in both stands(Table 5).

    Table 3 Net primary production(NPP;Mg·ha?1·yr?1)and carbon allocation(Mg C·ha?1·yr?1)in different components of vegetation in moist tropical forest of eastern Nepal

    Table 4 Net production(Mg·ha?1·yr?1)of trees and shrubs and their distribution(%)in different components in moist tropical forest of eastern Nepal

    Table 5 Pearson's correlations between stand density,basal area,biomass and net production of trees in undisturbed and disturbed stands of moist tropical forest in eastern Nepal

    Carbon budget and flux

    The dry matter values for standing crops,net production,litterfall etc.of both undisturbed and disturbed stands were converted to C(Figs.3 and 4).Aboveground C storage in both stands was 82%of the total C in vegetation and 64–69%of that stored in stand(vegetation plus soil).The forest received a C input through NPP of 12.26 Mg·ha–1·yr–1in US and 6.88 Mg·ha–1·yr–1in DS. Of this,69%was associated with aboveground and 31% with root NPP in US while in DS,contributions from aboveground and root parts to total NPP was 67 and 33%,respectively.

    Wood(bole plus branch)and roots accounted for 47%(in US)to 53%(in DS)of the total C input,the remaining 53%in US and 47%in DS being utilized forthe extension of the photosynthetic parts(tree and shrub foliage plus herbaceous shoots).The input from the foliage compartment to the litter compartment was 3.90 Mg·ha–1·yr–1(equivalent to 100%of herbaceous aboveground production plus 86%of foliage production) in US while it was 2.17 Mg·ha–1·yr–1for DS.About 14% (0.51 Mg·ha–1·yr–1in US and 0.26 Mg·ha–1·yr–1in DS)of foliage production did not find its way into leaf litterfall.

    Carbon addition to the soil through the turnover of roots and aboveground litter amounted to 1.98 and 4.80 Mg·ha–1·yr–1,respectively in US(fine root and litter turnover were 86 and 83%per yr,respectively). Similarly,C input due to the turnover of roots and aboveground litter was 1.05 and 2.30 Mg·ha–1·yr–1, respectively in DS(fine root and litter turnover were 82 and 79%per yr,respectively).Evidently,contribution of belowground plant parts to soil C is substantial,representing 29 and 32%of the total input in US and DS,respectively.

    Discussion

    Biomass and production

    The average biomass of the vegetation in both US and DS was 704.8 Mg?ha–1(331.7 Mg C?ha–1).More than double biomass in US as compared to DS may be related to the presence of more trees of larger girth classes.The lower biomass in DS may be the result of lower density of trees (234 individuals?ha–1)compared to US(466 individuals?ha–1).It may also be associated with disturbance activities like tree felling and removing for timber,firewood collection, lopping,grazing,and selective logging.

    Biomass allocated to shrub was higher in DS than US.It may be the result of their ability to utilize the space and resource created by disturbances.The higher biomass of herbs in US as compared to DS may be due to the nutrient rich soil at this stand(Gautam and Mandal 2013).The different causes are put forward by different workers regarding the accumulation of biomass in the forest.Mean annual precipitation explains 55%of the variation in biomass in seasonally dry tropical forests(Becknell et al. 2012).The available nutrients,soil,land use history,species composition and stand age are also responsible for the remaining variations(Powers et al.2009).

    Most of the C stock in the forest was associated with AGB(64–69%).Belowground biomass and soil up to the depth of 30 cm contained 31 to 36%C stock.Past study also found 33%C in top soil of 1 m in primary tropical forest of Singapore(Ngo et al.2013).Globally, tropical forests store about 50%C in AGB and next 50%within 1 m soil(Dixon et al.1994).The contribution of these pools to the total C stocks varies among the sites.For example,an African moist tropical forest had more than three times as much C in AGB as in soil to 1 m depth(Djomo et al.2011),while a Peruvian montane forest had twice as much C in soil as in AGB(Gibbon et al.2010).Two Asian forests,a tropical seasonal forest in China(Lü et al.2010)and a lowland dipterocarp forest in Malaysia(Saner et al.2012)store twice as much C in biomass as in soil.The contribution of leaf fall to total litterfall in the present study(69%at US and 76%at DS)is in the range of past studies in tropical forests(72–74%)(Girardin et al.2014;Malhi et al.2014).

    The data obtained for NPP are comparable with the different tropical forests of world(Table 6).Although having higher biomass,the NPP of present forest is relatively low indicating its maturity.Generally,in old growth tropical forests the NPP remains either stationary or decreases with age.The higher NPP of some old tropical forests may be linked to forest dynamics.The more dynamic forests favor faster growing trees and species with lower biomass and maintenance costs(Malhi et al.2009).

    Aboveground NPP of US(8.54 Mg C·ha?1·yr–1,excluding herbs)is comparable with the tropical forests of world while quite low than the eastern Amazonian forest(Doughty et al.2013).Moreover,aboveground NPP of DS(4.54 Mg C·ha?1·yr–1,excluding herb)is near to paleotropical forests,and mean of the world tropical forests.As far as the belowground NPP is concerned(1.77–2.98 Mg C·ha?1·yr–1),it is comparable to the Amazonian tropical forest(Table 6).

    Disturbance affects the NPP either by increasing resource availability or through the changes in functional properties of community.The disturbance results in canopy opening,which increases availability of light, favor new leaf production,and finally above-ground allocation increases(Arag?o et al.2009).But in the present study,percentage allocation to aboveground parts was almost same in both US and DS.It may be the result of severe lopping activities in DS,which is also indicated by less percentage of NPP allocation to twig and leaf,and higher percentage allocation to bole and branch in DS. Relatively higher percentage of NPP in coarse root of DS indicated higher allocation to root to collect more water and nutrients from nutrient deficient soil.

    Among different components of trees,leaf accounted highest percentage(43%in US and 42%in DS)of tree NPP.It was almost similar to the earlier findings for tropical forests(Hertel et al.2009;Doughty et al.2013). A global dataset of NPP allocation in tropical forests reported 34%in canopy,39%in wood,and 27%in fine roots(Malhi et al.2011).Present data of US,excluding herb NPP,also showed almost same pattern of allocation (33%canopy,46%woody tissue including coarse root, and 21%fine roots).

    Higher allocation of NPP in shrubs of DS as compared to that of US may be due to open canopy created as a result of logging and lopping of trees.The higher NPP of herbs in US might be associated with higher moisture and nutrients due to accumulation of more litter on the soil.Fine root,although sharing a very small percentage to total plant biomass,is one of the major belowground components of the NPP.In present study,FRB shared less than 1%to total biomass but contributed 19–20% to total NPP.The estimated FRP in this study was within the ranges of 1.1–6.2 Mg·ha?1·yr–1for global tropical forests(Vogt et al.1986).

    Carbon budget and flux

    Aboveground C storage in both stands was 82%of the total stored in vegetation and 64–69%of that stored in the stand(vegetation plus soil).These values compare with the aboveground storage of 83%vegetation C and 51%stand C in a dry tropical forest in India(Singh and Singh 1991).As in other tropical forests,relatively higher contribution of aboveground parts to the stand C in present case may be due to the presence of large sized trees with wide canopy.

    The contribution of aboveground(74–75%)and root(25–26%)parts to stand NPP in both forest stands is comparable to 72%NPP for aboveground parts(Singh and Singh 1991).The higher contribution of wood parts to the total C input as compared to foliage in DS may be due to the lopping of branches and twigs.On the other hand,closed canopy of trees in US accounted for higher input of foliage C compared to wood.

    The transfer of C from foliage compartment to the litter compartment involved 100%ofherbaceous aboveground production plus 86%of foliage production while the remaining amount,i.e.,about 14%of foliage production did not find its way into leaf litterfall.It happens due to the losses during senescence (Singh and Singh 1991).

    Table 6 Net primary production(NPP)of different forests of the World

    Litterfall and FRB are important vectors of nutrient recycling in forest ecosystems.Their turnover is usually determined by species,age groups,canopy cover, weather conditions and biotic factors.Total C input into soil through litter plus root turnover was 6.78 and 3.35 Mg·ha–1·yr–1in US and DS,respectively, suggesting substantial retention of C in the vegetation over in annual cycle(45%in US and 51%in DS).It indicates that the present forest is C-accumulating system,acting as a significant global C sink,as other wet tropical forests(Pan et al.2011).

    Conclusions

    The present study indicated that various types of anthropogenic disturbances have altered the structure and functioning of the economically important Saldominated forest.Although in having higher biomass,the NPP of present forest is relatively low indicating its maturity.Several disturbance activities like lopping,fodder collection,litter removal,grazing etc.result in the significant loss in stand biomass (53%)and net production(44%).Due to forest disturbance C stock and C sequestration capacity are reduced which reflects the higher C emissions.From the management point of view,better policy to reduce the C emission through vegetation should be formulated as per the objective of REDD+.In spite of this,both stands of this forest appear to act as C-accumulating system,an important global C sink.

    Abbreviations

    AGB:aboveground biomass;DS:disturbed forest stand;NPP:net primary production;US:undisturbed forest stand.

    Competing interests

    The authors declare that they have no competing interests.

    Authors’contributions

    TP carried out field and laboratory works,analyzed data,designed and drafted the manuscript.TN designed and supervised the research.Both authors read and approved the final manuscript.

    Authors’information

    1)TP participated in this research as a doctoral student.At present,he is working as Associate Professor in the Department of Botany,Mahendra Morang Adarsha Multiple Campus,Tribhuvan University,Biratnagar,Nepal. 2)TN is a professor of Botany in Post Graduate Campus,Tribhuvan University,Biratnagar,Nepal.

    Acknowledgements

    The authors are grateful to the Head,Department of Botany,Post Graduate Campus,Tribhuvan University,Biratnagar,Nepal for providing laboratory and library facilities.First author is thankful to the University Grants Commission, Nepal for Scholarship.Our heartfelt thanks go to Mr.KP Bhattarai who helped in data collection and laboratory analysis.

    Author details

    1Department of Botany,Mahendra Morang Adarsha Multiple Campus, Tribhuvan University,Biratnagar,Nepal.2Department of Botany,Post Graduate Campus,Tribhuvan University,Biratnagar,Nepal.

    References

    Arag?o LEOC,Malhi Y,Metcalfe DB,Silva-Espejo JE,Jiménez E,Navarrete D, Almeida S,Costa ACL,Salinas N,Phillips OL,Anderson LO,Alvarez E,Baker TR,Goncalvez PH,Huamán-Ovalle J,Mamani-Solórzano M,Meir P, Monteagudo A,Pati?o S,Pe?uela MC,Prieto A,Quesada CA,Rozas-Dávila A, Rudas A,Silva JA Jr,Vásquez R(2009)Above-and below-ground net primary productivity across ten Amazonian forests on contrasting soils. Biogeosciences 6:2759-2778

    Baral SK,Malla R,Ranabhat S(2009)Above-ground carbon stock assessment in different forest types of Nepal.Banko Jankari 19:10-14

    Becknell JM,Kissing Kucek L,Powers JS(2012)Aboveground biomass in mature and secondary seasonally dry tropical forests:A literature review and global synthesis.For Ecol Manag 276:88-95

    Beer C,Reichstein M,Tomelleri E,Ciais P,Jung M,Carvalhais N,Rodenbeck C, Altaf Arain M,Baldocchi D,Bonan GB,Bondeau A,Cescatti A,Lasslop G, Lindroth A,Lomas M,Luyssaert S,Margolis H,Oleson KW,Roupsard O, Veenendaal E,Viovy N,Williams C,Woodward FI,Papale D(2010)Terrestrial gross carbon dioxide uptake:global distribution and covariation with climate. Science 329:834-838

    Bombelli A,Henry M,Castaldi S,Adu-Bredu S,Arneth A,de Grandcourt A,Grieco E,Kutsch WL,Lehsten V,Rasile A,Reichstein M,Tansey K,Weber U,Valentini R(2009)The sub-Saharan Africa carbon balance,an overview.Biogeosci. Discuss 6:2085-2123.

    Chave J,Olivier J,Bongers F,Chatelet P,Forget P-M,van der Meer P,Norden N, Riera B,Charles-Dominique P(2008)Above-ground biomass and productivity in a rain forest of eastern South America.J Trop Ecol 24:355-366

    Chave J,Riera B,Dubois MA(2001)Estimation of biomass in a Neotropical forest of French Guiana:Spatial and temporal variability.J Trop Ecol 17:79-96

    Dixon RK,Solomon AM,Brown S,Houghton RA,Trexler MC,Wisniewski J(1994) Carbon pools and flux of global forest ecosystems.Science 263:185-190

    Djomo AN,Knohl A,Gravenhorst G(2011)Estimations of total ecosystem carbon pools distribution and carbon biomass current annual increment of a moist tropical forest.For Ecol Manag 261:1448-1459

    Doughty CE,Metcalfe DB,da Costa MC,de Oliveira AAR,Neto GFC,Silva JA, Arag?o LEOC,Almeida SS,Quesada CA,Girardin CAJ,Halladay K,da Costa AC,Malhi Y(2013)The production,allocation and cycling of carbon in a forest on fertile terra preta soil in eastern Amazonia compared with a forest on adjacent infertile soil.Plant Ecol Diversity 7:41-53

    Gautam TP,Mandal TN(2013)Soil characteristics in moist tropical forest of Sunsari district,Nepal.Nepal J Sci Tech 14:35-40

    Gibbon A,Silman MR,Malhi Y,Fisher JB,Meir P,Zimmermann M,Dargie GC, Farfan WR,Garcia KC(2010)Ecosystem carbon storage across the grasslandforest transition in the high Andes of Manu National Park,Peru.Ecosystems 13:1097-1111

    Girardin CAJ,Espejob JES,Doughty CE,Huasco WH,Metcalfe DB,Durand-Baca L, Marthews TR,Aragao LEOC,Farfán-Rios W,García-Cabrera K,Halladay K, Fisher JB,Galiano-Cabrera DF,Huaraca-Quispe LP,Alzamora-Taype I,Eguiluz-Mora L,Salinas-Revilla N,Silman MR,Meir P,Malhi Y(2014)Productivity and carbon allocation in a tropical montane cloud forest in the Peruvian Andes. Plant Ecol Diversity 7:107-123

    Hertel D,Moser G,Culmsee H,Erasmi S,Horna V,Schuldt B,Leuschner C(2009) Below-and above-ground biomass and net primary production in a paleotropical natural forest(Sulawesi,Indonesia)as compared to neotropical forests.For Ecol Manag 258:1904-1912

    Ibrahima A,Mvondo ZEA,Ntonga JC(2010)Fine root production and distribution in the tropical rainforests of south-western Cameroon:effects of soil type and selective logging.iForest:Biogeosc For 3:130-136

    Jenny H,Gessel SP,Bingham FT(1949)Comparative study of decomposition rates of organic matter in temperate and tropical region.Soil Sci 68:419-432

    Kalembasa SJ,Jenkinson DS(1973)A comparative study of titremetric and gravimetric methods for the determination of organic carbon in soil.J Sci Food Agri 24:1085-1090

    Lü XT,Yin JX,Jepsen MR,Tang JW(2010)Ecosystem carbon storage and partitioning in a tropical seasonal forest in Southwestern China.For Ecol Manag 260:1798-1803

    Malhi Y,Aragao LEO,Metcalfe DLB,Paiva R,Quesada CA,Almeida S,Anderson L, Brando P,Chambers JQ,Costa ACL,Hutyra L,Oliveira P,Patino S,Pyle EH, Robertson AL,Teixeira LM(2009)Comprehensive assessment of carbonproductivity,allocation and storage in three Amazonian forests.Glob Change Biol 15:1255-1274

    Malhi Y,Doughty C,Galbraith D(2011)The allocation of ecosystem net primary productivity in tropical forests.Philos Trans R Soc B:Biol Sci 366:3225-3245

    Malhi Y,Farfán Amézquita F,Doughty CE,Silva-Espejo JE,Girardin CAJ,Metcalfe DB,Arag?o LEOC,Huaraca-Quispe LP,Alzamora-Taype I,Eguiluz-Mora L, Marthews TR,Halladay K,Quesada CA,Robertson AL,Fisher JB,Zaragoza-Castells J,Rojas-Villagra CM,Pelaez-Tapia Y,Salinas N,Meir P,Phillips OL (2014)The productivity,metabolism and carbon cycle of two lowland tropical forest plots in south-western Amazonia,Peru.Plant Ecol Diversity 7: 85-105

    Mandal T(1999)Ecological analysis of recovery of landslide damaged Sal forest ecosystem in Nepal Himalaya.Dissertation,Banaras Hindu University,Varanasi

    McBrayer JF,Cromack KJ(1980)Effect of snowpack on oak-litter breakdown and nutrient release in a Minnesota forest.Pedobiologia 20:47-54

    Ngo KM,Turner BL,Muller-Landau HC,Davies SJ,Larjavaara M,Nik Hassan NF, Lum S(2013)Carbon stocks in primary and secondary tropical forests in Singapore.For Ecol Manag 296:81-89

    Noguchi H,Suwa R,de Souza CAS,da Silva RP,dos Santos J,Higuchi N,Kajimoto T,Ishizuka M(2014)Examination of vertical distribution of fine root biomass in a tropical moist forest of the Central Amazon,Brazil.Japan Agri Res Quarterly 48:231-235

    Pan Y,Birdsey RA,Fang J,Houghton R,Kauppi PE,Kurz WA,Phillips OL, Shvidenko A,Lewis SL,Canadell JG,Ciais P,Jackson RB,Pacala SW,McGuire AD,Piao S,Rautiainen A,Sitch S,Hayes D(2011)A large and persistent carbon sink in the world's forests.Science 333:988-993

    Powers JS,Becknell JM,Irving J,Pèrez-Aviles D(2009)Diversity and structure of regenerating tropical dry forests in Costa Rica:Geographic patterns and environmental drivers.For Ecol Manag 258:959-970

    Powers JS,Perez-Aviles D(2013)Edaphic factors are a more important control on surface fine roots than stand age in secondary tropical dry forests.Biotropica 45:1-9

    Roderstein M,Hertel D,Leuschner C(2005)Above and below-ground litter production in three tropical montane forests in southern Ecuador.J Trop Ecol 21:483-492

    Saner P,Loh YY,Ong RC,Hector A(2012)Carbon stocks and fluxes in tropical lowland dipterocarp rainforests in Sabah,Malaysian Borneo.PloS One 7.doi: 10.1371/journal.pone.0029642

    Singh JS,Singh SP(1992)Forests of Himalaya.Gyanodaya Prakashan,Nainital

    Singh KP,Mandal TN,Tripathi SK(2001)Patterns of restoration of soil physciochemical properties and microbial biomass in different landslide sites in the sal forest ecosystem of Nepal Himalaya.Ecol Eng 17:385-401

    Singh L,Singh JS(1991)Storage and flux of nutrients in a dry tropical forest in India.Ann Botany 68:275-284

    Singh SP,Adhikari BS,Zobel DB(1994)Biomass,productivity,leaf longevity,and forest structure in the Central Himalaya.Ecol Monogr 64:401-421

    Srivastava SK,Singh KP,Upadhyay RS(1986)Fine root growth dynamics in teak (Tectona grandis Linn.F.).Can J For Res 16:1360-1364

    Vogt KA,Grier CC,Vogt DJ(1986)Production,turnover and nutrient dynamics of above and belowground detritus of world forests.Adv Ecol Res 15:303-377

    *Correspondence:tilakg024@gmail.com

    1Department of Botany,Mahendra Morang Adarsha Multiple Campus, Tribhuvan University,Biratnagar,Nepal

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

    ?2016 Gautam and Mandal.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.

    10 February 2016 Accepted:18 April 2016

    黄片大片在线免费观看| 777久久人妻少妇嫩草av网站| 村上凉子中文字幕在线| 国产精品一及| 国内揄拍国产精品人妻在线| 欧美一区二区国产精品久久精品 | 午夜成年电影在线免费观看| 啦啦啦韩国在线观看视频| 熟妇人妻久久中文字幕3abv| 波多野结衣高清无吗| 国产成人精品无人区| 非洲黑人性xxxx精品又粗又长| 级片在线观看| 国产精品乱码一区二三区的特点| 99久久久亚洲精品蜜臀av| 久久99热这里只有精品18| 一个人免费在线观看电影 | 午夜精品一区二区三区免费看| 少妇熟女aⅴ在线视频| 一区二区三区高清视频在线| 波多野结衣高清作品| 母亲3免费完整高清在线观看| 国产日本99.免费观看| 国内精品久久久久精免费| 亚洲专区国产一区二区| 男女做爰动态图高潮gif福利片| 日本免费a在线| 日本一区二区免费在线视频| 99国产精品99久久久久| 97碰自拍视频| 99国产综合亚洲精品| 久9热在线精品视频| 在线播放国产精品三级| 国产黄片美女视频| 男女做爰动态图高潮gif福利片| 精品免费久久久久久久清纯| 欧美黄色淫秽网站| 一本综合久久免费| 国产亚洲精品第一综合不卡| 我要搜黄色片| 性欧美人与动物交配| 亚洲午夜理论影院| 人人妻,人人澡人人爽秒播| 黄色视频不卡| 国产成人一区二区三区免费视频网站| 丁香六月欧美| 亚洲最大成人中文| 亚洲片人在线观看| 19禁男女啪啪无遮挡网站| 日本精品一区二区三区蜜桃| 亚洲成av人片免费观看| 精品国产乱码久久久久久男人| 久久中文字幕人妻熟女| 亚洲片人在线观看| 亚洲av成人不卡在线观看播放网| 三级国产精品欧美在线观看 | av视频在线观看入口| 久久久久久免费高清国产稀缺| 在线观看舔阴道视频| 怎么达到女性高潮| 黄色a级毛片大全视频| 日韩欧美在线二视频| 97超级碰碰碰精品色视频在线观看| 欧美一级毛片孕妇| 国产av又大| 日本 欧美在线| 欧美成人一区二区免费高清观看 | 午夜激情福利司机影院| 亚洲一码二码三码区别大吗| 精品免费久久久久久久清纯| 十八禁网站免费在线| 亚洲av美国av| 九九热线精品视视频播放| 国产一区二区三区在线臀色熟女| 欧美zozozo另类| 欧美乱妇无乱码| 热99re8久久精品国产| 欧美日韩亚洲国产一区二区在线观看| 欧美3d第一页| 精品无人区乱码1区二区| 欧美中文综合在线视频| 国产激情欧美一区二区| 亚洲人成网站在线播放欧美日韩| 午夜免费成人在线视频| АⅤ资源中文在线天堂| 免费看十八禁软件| 国产av又大| 欧美3d第一页| 久久热在线av| 91在线观看av| 日韩精品中文字幕看吧| 熟女电影av网| 欧美日韩黄片免| 美女大奶头视频| 人人妻,人人澡人人爽秒播| 日韩精品青青久久久久久| 亚洲成人免费电影在线观看| 亚洲专区中文字幕在线| 欧美又色又爽又黄视频| 日韩欧美国产一区二区入口| 日韩 欧美 亚洲 中文字幕| 女生性感内裤真人,穿戴方法视频| 亚洲av熟女| 三级毛片av免费| 可以免费在线观看a视频的电影网站| 亚洲性夜色夜夜综合| 午夜免费成人在线视频| 国产精品,欧美在线| 亚洲精品中文字幕一二三四区| 亚洲免费av在线视频| 夜夜躁狠狠躁天天躁| 欧美不卡视频在线免费观看 | 国产午夜精品久久久久久| 这个男人来自地球电影免费观看| 后天国语完整版免费观看| 九色国产91popny在线| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品999在线| 欧美日韩亚洲综合一区二区三区_| av福利片在线| 成人18禁高潮啪啪吃奶动态图| 少妇熟女aⅴ在线视频| 三级国产精品欧美在线观看 | 亚洲专区中文字幕在线| 亚洲国产精品999在线| 天堂√8在线中文| 亚洲精品国产精品久久久不卡| 91麻豆精品激情在线观看国产| 欧美日韩乱码在线| 日本撒尿小便嘘嘘汇集6| 国产片内射在线| 亚洲av五月六月丁香网| 欧美午夜高清在线| 免费观看人在逋| 欧美日韩福利视频一区二区| 久久午夜综合久久蜜桃| 一级黄色大片毛片| 日韩欧美一区二区三区在线观看| 狠狠狠狠99中文字幕| av福利片在线| 可以在线观看毛片的网站| 每晚都被弄得嗷嗷叫到高潮| 免费电影在线观看免费观看| 哪里可以看免费的av片| 一二三四社区在线视频社区8| 黄频高清免费视频| 免费观看人在逋| 1024视频免费在线观看| 日韩av在线大香蕉| 色综合站精品国产| 精品久久久久久久久久久久久| 午夜成年电影在线免费观看| 日日摸夜夜添夜夜添小说| 婷婷精品国产亚洲av| 97人妻精品一区二区三区麻豆| 亚洲精品美女久久av网站| 又大又爽又粗| www国产在线视频色| 欧美高清成人免费视频www| 成人av一区二区三区在线看| 88av欧美| 在线观看免费日韩欧美大片| 高清在线国产一区| 手机成人av网站| 国产三级在线视频| 久久久久久国产a免费观看| 国产高清有码在线观看视频 | 狠狠狠狠99中文字幕| 午夜福利在线观看吧| 久久人妻av系列| 动漫黄色视频在线观看| 亚洲18禁久久av| 亚洲精华国产精华精| 精品欧美国产一区二区三| 久久久久久九九精品二区国产 | 国产精品影院久久| av在线播放免费不卡| 亚洲av第一区精品v没综合| 国产亚洲欧美98| 亚洲在线自拍视频| 欧美成人午夜精品| 少妇被粗大的猛进出69影院| 亚洲一码二码三码区别大吗| 欧美3d第一页| 一本综合久久免费| 欧美又色又爽又黄视频| 亚洲无线在线观看| 欧美日韩一级在线毛片| 中文字幕久久专区| 亚洲国产高清在线一区二区三| 香蕉av资源在线| 久久精品国产亚洲av高清一级| 最近最新中文字幕大全电影3| 国产精品综合久久久久久久免费| 母亲3免费完整高清在线观看| 丁香六月欧美| 18禁国产床啪视频网站| 精品午夜福利视频在线观看一区| 老司机靠b影院| 日本免费a在线| 99热这里只有是精品50| 搡老岳熟女国产| 后天国语完整版免费观看| 特级一级黄色大片| 欧美日韩中文字幕国产精品一区二区三区| 国内久久婷婷六月综合欲色啪| 午夜福利高清视频| 1024视频免费在线观看| 国产伦一二天堂av在线观看| 免费电影在线观看免费观看| 两人在一起打扑克的视频| 久热爱精品视频在线9| 淫秽高清视频在线观看| 夜夜躁狠狠躁天天躁| 一级a爱片免费观看的视频| 午夜精品在线福利| 黄色成人免费大全| 性色av乱码一区二区三区2| 亚洲人成网站高清观看| 男女下面进入的视频免费午夜| 在线观看日韩欧美| 中文亚洲av片在线观看爽| 中文字幕熟女人妻在线| 午夜精品一区二区三区免费看| 亚洲国产精品sss在线观看| 日本一二三区视频观看| 免费在线观看视频国产中文字幕亚洲| 亚洲av成人av| 欧美国产日韩亚洲一区| 99国产极品粉嫩在线观看| 一级片免费观看大全| 国产1区2区3区精品| 一区福利在线观看| 国产视频内射| 又粗又爽又猛毛片免费看| 久久精品91无色码中文字幕| 757午夜福利合集在线观看| 亚洲精品美女久久久久99蜜臀| 老司机福利观看| 91字幕亚洲| 色精品久久人妻99蜜桃| 少妇熟女aⅴ在线视频| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3| 日韩欧美精品v在线| 最近在线观看免费完整版| 老司机在亚洲福利影院| 我要搜黄色片| 亚洲成av人片在线播放无| 亚洲精品色激情综合| 99久久综合精品五月天人人| 999久久久国产精品视频| 日日干狠狠操夜夜爽| 可以在线观看毛片的网站| 国产av一区二区精品久久| 美女 人体艺术 gogo| 黑人巨大精品欧美一区二区mp4| 国内毛片毛片毛片毛片毛片| 看片在线看免费视频| 丰满的人妻完整版| 亚洲av电影不卡..在线观看| 亚洲乱码一区二区免费版| 男女视频在线观看网站免费 | 免费人成视频x8x8入口观看| 亚洲自偷自拍图片 自拍| 哪里可以看免费的av片| 黄色a级毛片大全视频| 亚洲国产看品久久| 很黄的视频免费| 精品国产超薄肉色丝袜足j| 久久人妻福利社区极品人妻图片| 一级毛片精品| 波多野结衣高清作品| 午夜成年电影在线免费观看| av在线播放免费不卡| 亚洲人成伊人成综合网2020| 国产高清视频在线播放一区| 好男人在线观看高清免费视频| av免费在线观看网站| 中文字幕熟女人妻在线| 中文字幕av在线有码专区| 美女黄网站色视频| 欧美大码av| 99re在线观看精品视频| 老司机午夜福利在线观看视频| 999久久久国产精品视频| 男插女下体视频免费在线播放| 色老头精品视频在线观看| 天天添夜夜摸| 99久久99久久久精品蜜桃| 亚洲国产欧洲综合997久久,| 国产高清videossex| 国产在线观看jvid| 成年免费大片在线观看| 最好的美女福利视频网| 亚洲人成伊人成综合网2020| 男人舔奶头视频| 久久精品91无色码中文字幕| 国产激情久久老熟女| 午夜福利18| 岛国在线免费视频观看| www日本在线高清视频| 亚洲性夜色夜夜综合| 窝窝影院91人妻| 美女高潮喷水抽搐中文字幕| 久久久久久久久免费视频了| 草草在线视频免费看| 国产亚洲欧美98| 国产精品综合久久久久久久免费| 国产成人精品久久二区二区91| 男女那种视频在线观看| 午夜免费观看网址| 亚洲免费av在线视频| 国产黄片美女视频| 亚洲五月天丁香| 超碰成人久久| a级毛片a级免费在线| 国产激情欧美一区二区| 国产黄a三级三级三级人| 在线播放国产精品三级| 亚洲激情在线av| svipshipincom国产片| 99re在线观看精品视频| 亚洲精品国产一区二区精华液| 在线视频色国产色| 欧美一级毛片孕妇| 亚洲国产看品久久| 国产激情久久老熟女| 啦啦啦免费观看视频1| 亚洲免费av在线视频| 亚洲精品在线观看二区| 亚洲欧美激情综合另类| 欧美日本亚洲视频在线播放| 国产区一区二久久| 亚洲成人精品中文字幕电影| 狂野欧美激情性xxxx| 99精品在免费线老司机午夜| 亚洲,欧美精品.| 法律面前人人平等表现在哪些方面| 美女免费视频网站| 成人亚洲精品av一区二区| 国产在线观看jvid| 欧美日韩亚洲国产一区二区在线观看| 国产一区在线观看成人免费| 欧美成人免费av一区二区三区| 亚洲成人精品中文字幕电影| 国产v大片淫在线免费观看| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 亚洲精品在线观看二区| 久久久久久亚洲精品国产蜜桃av| 国产野战对白在线观看| 欧美中文日本在线观看视频| 亚洲国产精品999在线| 欧美丝袜亚洲另类 | 久久久久久国产a免费观看| 国产激情欧美一区二区| 久久久国产精品麻豆| 波多野结衣巨乳人妻| 亚洲熟妇熟女久久| 国产成人av激情在线播放| 亚洲精品色激情综合| 亚洲人成网站高清观看| 99精品欧美一区二区三区四区| 91大片在线观看| 少妇被粗大的猛进出69影院| 久久国产精品人妻蜜桃| 天堂√8在线中文| 看免费av毛片| 久久久久久大精品| 亚洲国产日韩欧美精品在线观看 | 天堂√8在线中文| 怎么达到女性高潮| 久久久久久久午夜电影| 国产熟女xx| 99国产综合亚洲精品| 欧美日韩国产亚洲二区| 午夜亚洲福利在线播放| 亚洲一区中文字幕在线| 老鸭窝网址在线观看| 又黄又爽又免费观看的视频| 午夜福利在线观看吧| 一级黄色大片毛片| 香蕉国产在线看| 黄色片一级片一级黄色片| 欧美日韩国产亚洲二区| 两个人视频免费观看高清| 亚洲一区中文字幕在线| 又大又爽又粗| 我的老师免费观看完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 大型黄色视频在线免费观看| aaaaa片日本免费| 神马国产精品三级电影在线观看 | 国内揄拍国产精品人妻在线| 午夜久久久久精精品| 宅男免费午夜| 看免费av毛片| 嫩草影院精品99| 搡老熟女国产l中国老女人| 免费看美女性在线毛片视频| 婷婷亚洲欧美| 丝袜人妻中文字幕| 每晚都被弄得嗷嗷叫到高潮| 久久久久性生活片| 美女大奶头视频| 熟女少妇亚洲综合色aaa.| 两个人免费观看高清视频| 欧美大码av| 国产高清视频在线观看网站| 看黄色毛片网站| 欧美日韩亚洲综合一区二区三区_| 久久久久免费精品人妻一区二区| 欧美成人性av电影在线观看| 精品少妇一区二区三区视频日本电影| 国产亚洲av嫩草精品影院| 亚洲成人久久性| 午夜激情福利司机影院| 久久久国产欧美日韩av| 成熟少妇高潮喷水视频| 丰满的人妻完整版| 成人一区二区视频在线观看| 俄罗斯特黄特色一大片| 午夜激情av网站| 99久久无色码亚洲精品果冻| 久久香蕉精品热| 欧美色视频一区免费| 国产高清视频在线播放一区| 久久久久久亚洲精品国产蜜桃av| 日本 av在线| 中文在线观看免费www的网站 | 好看av亚洲va欧美ⅴa在| 性色av乱码一区二区三区2| 欧美日韩精品网址| 午夜精品在线福利| 婷婷精品国产亚洲av在线| 久久久久国产一级毛片高清牌| avwww免费| 国产aⅴ精品一区二区三区波| 欧美黑人欧美精品刺激| 中文字幕高清在线视频| 欧美成人性av电影在线观看| 国产精品久久久久久久电影 | 亚洲片人在线观看| 国产区一区二久久| 亚洲 欧美 日韩 在线 免费| 色播亚洲综合网| 亚洲欧美日韩东京热| or卡值多少钱| 亚洲欧美日韩高清在线视频| 成人欧美大片| 特级一级黄色大片| www.999成人在线观看| 国产又黄又爽又无遮挡在线| 亚洲 欧美一区二区三区| 最好的美女福利视频网| 1024手机看黄色片| 悠悠久久av| 毛片女人毛片| 亚洲电影在线观看av| 欧美日本亚洲视频在线播放| 一本久久中文字幕| 午夜激情av网站| 国产精品久久久久久人妻精品电影| 欧美日本视频| 九色成人免费人妻av| 欧洲精品卡2卡3卡4卡5卡区| 五月伊人婷婷丁香| 成年免费大片在线观看| 日韩欧美 国产精品| 欧美乱妇无乱码| 人妻夜夜爽99麻豆av| 亚洲色图 男人天堂 中文字幕| 日韩精品中文字幕看吧| 日本 欧美在线| 在线观看66精品国产| 成人午夜高清在线视频| 国产精品野战在线观看| 亚洲精品色激情综合| 成人18禁高潮啪啪吃奶动态图| 国产在线观看jvid| 真人一进一出gif抽搐免费| 国产一级毛片七仙女欲春2| 高清毛片免费观看视频网站| 久久久久久免费高清国产稀缺| 俺也久久电影网| 身体一侧抽搐| 欧美日韩福利视频一区二区| 此物有八面人人有两片| 国产成人一区二区三区免费视频网站| 亚洲av成人不卡在线观看播放网| 久久这里只有精品19| 日本五十路高清| 国产日本99.免费观看| 国产三级中文精品| bbb黄色大片| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久末码| 欧美绝顶高潮抽搐喷水| 我要搜黄色片| 日本熟妇午夜| 欧美日韩乱码在线| 成人永久免费在线观看视频| 国产av不卡久久| 亚洲人成网站高清观看| 亚洲精品美女久久久久99蜜臀| 久久午夜综合久久蜜桃| 久久久久久久久久黄片| 最近最新中文字幕大全电影3| 久久99热这里只有精品18| 91麻豆av在线| 丁香欧美五月| 成年女人毛片免费观看观看9| 亚洲真实伦在线观看| 18禁黄网站禁片午夜丰满| 国产高清视频在线播放一区| 欧美色视频一区免费| 亚洲精品中文字幕在线视频| 18美女黄网站色大片免费观看| 精品福利观看| 男女之事视频高清在线观看| 国产午夜精品论理片| 最近视频中文字幕2019在线8| 窝窝影院91人妻| 中文字幕人妻丝袜一区二区| 亚洲欧洲精品一区二区精品久久久| 他把我摸到了高潮在线观看| 高清在线国产一区| 国产蜜桃级精品一区二区三区| 久久久久久久久免费视频了| 性色av乱码一区二区三区2| 亚洲国产高清在线一区二区三| 在线看三级毛片| 久久精品aⅴ一区二区三区四区| 国产精品电影一区二区三区| 免费搜索国产男女视频| 又粗又爽又猛毛片免费看| 国产精品 欧美亚洲| 亚洲av日韩精品久久久久久密| 国产黄片美女视频| 中文字幕精品亚洲无线码一区| 日韩欧美一区二区三区在线观看| 亚洲国产欧美网| 又大又爽又粗| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 搡老妇女老女人老熟妇| 悠悠久久av| 亚洲男人的天堂狠狠| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美精品综合一区二区三区| 午夜老司机福利片| 他把我摸到了高潮在线观看| 搡老熟女国产l中国老女人| 桃红色精品国产亚洲av| 日韩欧美精品v在线| 怎么达到女性高潮| 亚洲国产精品成人综合色| 亚洲 国产 在线| 精品久久久久久久久久久久久| 欧美性猛交╳xxx乱大交人| 欧美三级亚洲精品| 亚洲第一欧美日韩一区二区三区| 日韩三级视频一区二区三区| 亚洲午夜精品一区,二区,三区| 欧美日韩福利视频一区二区| 少妇被粗大的猛进出69影院| 国产午夜精品论理片| 国产av又大| av欧美777| 国产亚洲精品综合一区在线观看 | 亚洲午夜精品一区,二区,三区| 18禁黄网站禁片午夜丰满| 国产精品久久久久久精品电影| 国产1区2区3区精品| 免费看十八禁软件| 日韩欧美精品v在线| 国产99久久九九免费精品| 亚洲国产欧美一区二区综合| 女人高潮潮喷娇喘18禁视频| 亚洲欧美日韩高清在线视频| 亚洲aⅴ乱码一区二区在线播放 | 熟女少妇亚洲综合色aaa.| 九九热线精品视视频播放| 亚洲熟女毛片儿| 国产精品亚洲av一区麻豆| 九色成人免费人妻av| 变态另类成人亚洲欧美熟女| 国产亚洲欧美98| 国产欧美日韩精品亚洲av| 欧美精品亚洲一区二区| 免费在线观看成人毛片| 天堂动漫精品| 国产精品98久久久久久宅男小说| 夜夜夜夜夜久久久久| 老汉色av国产亚洲站长工具| 曰老女人黄片| x7x7x7水蜜桃| 在线免费观看的www视频| 欧美乱色亚洲激情| 久久香蕉激情| 久久99热这里只有精品18| 欧美乱色亚洲激情| 女警被强在线播放| 日韩欧美三级三区| 成人国产综合亚洲| 免费在线观看成人毛片| 亚洲精品美女久久久久99蜜臀| 国语自产精品视频在线第100页| 精品无人区乱码1区二区| 国产高清视频在线观看网站| 婷婷精品国产亚洲av| 久久久国产成人精品二区| 琪琪午夜伦伦电影理论片6080| 日韩中文字幕欧美一区二区|