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

    Biomass and net primary productivity of mangrove communities along the Oligohaline zone of Sundarbans,Bangladesh

    2018-01-03 09:20:12MdKamruzzamanShamimAhmedandAkiraOsawa
    Forest Ecosystems 2017年4期

    Md.Kamruzzaman,Shamim Ahmedand Akira Osawa

    Biomass and net primary productivity of mangrove communities along the Oligohaline zone of Sundarbans,Bangladesh

    Md.Kamruzzaman1,2*,Shamim Ahmed2and Akira Osawa1

    Background:The article presents the first estimates of biomass and productivity for mangrove forests along the Oligohaline zone of the Sundarbans Reserve Forest(SRF),Bangladesh.This study was conducted overone year from March 2016 to April 2017.Stand structure,above and below-ground biomass changes,and litterfall production were measured within a 2100 m2sample plot.Methods:All trees in the study plots were numbered and height(H)and diameter at breast height(DBH)were measured.Tree height(H)and DBH for each tree were measured in March 2016 and 2017.We apply the above and belowground biomass equation for estimating the biomass of the mangrove tree species(Chave et al.Oecologia 145:87?99,2005;Komiyama et al.J Trop Ecol 21:471–477,2005).Litterfall was collected using 1-mm mesh litter traps with collection area of 0.42 m2.Net Primary Production(NPP)was estimated by the summation method of Ogawa Primary productivity of Japanese forests:productivity of terrestrial communities,JIBP synthesis(1977)and Matsuura and Kajimoto Carbon dynamics of terrestrial ecosystem:Systems approach to global environment(2013).Results:Heritiera fomes has maintained its dominance of the stand and also suffered the highest tree mortality(2.4%)in the suppressed crown class.The total above-ground biomass(AGB)and below-ground biomass(BGB)of the studied stand was 154.8 and 84.2 Mg·ha?1,respectively.Among the total biomass of the trees,64.8%was allocated to AGB and 35.2%to BGB.In case of species-wise contribution of biomass allocation,Avicennia officinalis showed the highest score and Aglaia cucullata the lowest.Mean annual total litterfall was 10.1 Mg·ha?1·yr?1,with the maximum litterfall in winter or dry season and late summer or rainy season.The mean AGB increment and above-ground net primary productivity(AGNPP)were 7.1 and 17.2 Mg·ha?1·yr?1,respectively.Total net primary productivity(NPP)was estimated to be 21.0 Mg·ha?1·yr?1over the observed period.The results in the Sundarbans mangrove forests exhibited that mangrove communities with similar height and diameter produced different biomass production with the different basal area.The present analysis revealed that the root biomass was large enough and the mean ratio of above?/below-ground biomass was estimated to be 1.84.Conclusions:Mangrove communities growing at the oligohaline zone of the Sundarbans,Bangladesh showed high biomass and net primary production indicating their ecological and conservation significance that may be considered in future decision making process for the area as well as in understanding the role of Sundarbans mangrove forest on mitigating the effect of global warming.

    Primary production,Litterfall,Monthly pattern,Biomass carbon,Species wise contribution,Sundarbans

    Background

    Mangrove forestsplay an importantrole in the functioning of tropical coastal ecosystems(Day et al.1987),as well as unique wetland ecosystems in intertidal coastal regions of the tropics and subtropics(Lugo and Snedaker 1975;Nagarajan et al.2008).Mangroves are also important contributors of nutrients to estuarine and inshore productivity through litterfall.Litterfall is a valuable indicator of their productivity and it is the most easily measured component of total net primary productivity(Day et al.1987;Mackey and Smail 1995).Biomass and net primary productivity of mangrove forest have been studied previously in different mangrove forests across the world(e.g.,Putz and Chan 1986;Day et al.1987,1996;Saintilan 1997;Komiyama et al.2000;Kamruzzaman et al.2017)and the purposes of their biomass estimation and productivity are mainly related to consideration of ecosystem management and evaluation of carbon stock(Tamai et al.1986,Komiyama et al.1987,2000).

    The Sundarbans,the world’s largest continuous patch of mangrove forest,is located in the estuary of the river Ganges-Brahmaputra.The forest is distributed over two neighboring countries, Bangladesh and India. In Bangladesh,the forest covers 6017 km2(21o30′–22o30′N,89o00′–89o55′E).It is the first mangrove forest in the world which has been brought under scientific management since 1982,thereafter no apparent forest harvest and catastrophic disturbance have occurred.As recently as 200 years ago the Sundarbans extended across most of southwestern Bangladesh(Islam 2001).From the available literature it is known that during the period 1873–1933,the forest was reduced from 7599 to 6000 km2(Curtis 1933;Blasco 1977;Iftekhar and Saenger 2008)and in the recent years the boundary is almost same as 6017 km2.Chaffey et al.(1985)divided the Sundarbans mangrove forest into three ecological zones viz.fresh water(Oligohaline),moderately saline water(Mesohaline)and salt water(Polyhaline)zones.Based on salinity distribution,three salinity zones:Oligohaline(salinity <2 dsm?1),Mesohaline(salinity 2–4 dsm?1),Polyhaline(>4 dsm?1)zones could be recognized(Siddiqi 2001).Changes in salinity might be responsible for the spatial distribution of plant communities(Ahmed et al.2011).

    Our study was conducted in a mangrove community within the Oligohaline zone of the Sundarbans,Bangladesh.A few studies have been conducted,such as allometry,biomass and carbon storage of the Sundarbans Reserve Forest(SRF),Bangladesh with respect to its biodiversity and role in carbon sink among the other forest ecosystems in the tropics(Iftekhar and Saenger 2008;Ahmed et al.2011;Rahman et al.2015;Hossain et al.2015).However,no studies were found in literature that evaluated net primary productivity of the mangrove communities in this area.

    Thus,the objective of this study was to determine net primary productivity of mangrove forests along the Oligohaline zone of SRF.Its relationship to different parameters were also examined and presented.For the majority of mangrove studies,dead tree biomass was excluded and litterfall was measured just for a short period(Christensen 1978;Day et al.1987;Hernádez et al.2011).However,in this study,we examined 2 years of data for each of the components and took into account tree mortality during calculation of net primary productivity.

    Methods

    Study site

    The study was conducted in Dhangmari,Karamjol,and Ghagramari areas of SRF,Bangladesh,over 1 year from March 2016 to April 2017.This area receives regular tidal inundation through the River Passur.The mean annual rainfall in the region varies from about 2000 mm in the east and 1600 mm in the west.The rainy season is from June to September.October to February is considered winter and a drier season.Higher temperatures(26 °C–34 °C)occur during March and June and lower temperatures(12 °C–25 °C)during December and February.The annual relative humidity varies from 70 to 80%(Rahman and Asaduzzaman 2010).We established 21 plots(10 m×10 m),which are located in the Oligohaline zone of the Sundarbans,covering 2100 m2where all were considered around 200 m away from the shore for avoiding destruction of the plot due to river erosion and damages due to storms.All trees in the study plots were numbered and height(H)and diameter at breast height(DBH)were measured.Tree height(H)and DBH for each tree were measured in March 2016 and 2017.Some error in measuring tree growth was likely because only 1 year of DBH change was recorded.

    Biomass estimation

    Although allometric relationships have been developed between DBH and weight for a number of mangrove species around the world,biomass production depends on the interaction between edaphic,climatic and topographic factors of the specified area and the species that occurred on the specific area.Unfortunately there is only one developed allometric equation for biomass estimation of the studied mangrove species encountered in this study.The above-ground biomass of each individual of E.agallocha was estimated using allometric relationships between DBH and biomass,established by Hossain et al.(2015).So we tentatively apply the following aboveground biomass equation for the other mangrove tree species(Chave et al.2005).

    where AGB=aboveground biomass,ρ=wood density,D=DBH.The wood density data were obtained from Global Wood Density Database(Chave et al.2009).

    Below-ground biomass of mangrove tree species was estimated using common allometric relationships between DBH and biomass,established by Komiyama et al.(2005).

    where BGB=belowground biomass,ρ=wood density,D=DBH.

    Individual tree’s biomass calculation for each species and its increment for each tree were estimated using the same allometric equation.The common allometric relationship for root weight was derived from the relationship between below-ground and above-ground weight of trees,so wood density is the same for both above-ground and belowground weight estimation.

    Litterfall collection

    Litterfall was collected using 1-mm mesh litter traps with collection area of 0.42 m2.Two litter traps were placed in each 10 m × 10 m plot at>1 m high from the ground to avoid tidal water,i.e.,a total of 42 litter traps were placed within the 21 plots.The litter traps were emptied monthly;the collected litterfall was kept in a cotton bag and carried to the laboratory where it was separated into leaves,stipules,branches,flower buds,flowers,and propagules.Individual litterfall components were dried at 80°C for 48 h,desiccated at room temperature,and then weighed using a digital balance(EK-600H,A&D Co.,Ltd.,Tokyo,Japan).

    Net primary production

    Net Primary Production(NPP), ΔPn(Mg·ha?1·yr?1),was estimated by the summation method of Ogawa(1977)and Matsuura and Kajimoto(2013)as follows:

    where Δy(Mg·ha?1·yr?1)is biomass increment disregarding tree mortality,ΔL(Mg·ha?1·yr?1)is total dead matter,which is sum of total litter production and amount of dry matter due to tree mortality,and ΔG(Mg·ha?1·yr?1)is amount of grazing.We do not account for direct herbivory in this study,therefore,ΔG is negligibly small in the present study area,and ΔG is ignored in the current study(Kamruzzaman et al.2017).

    Results

    Community structure and composition

    Mangrove species and their structural characteristics in the study area along the oligohaline zone of SRF were presented in Table 1.A total of six mangrove speciesbelonging to five families were observed in the tree layer:Heritiera fomes,Excoecaria agallocha, Xylocarpus mekongensis,Bruguiera sexangula,Avicennia officinalis,and Aglaia cucullata.The mean density,height,DBH,and wood density of H.fomes,E.agallocha,X.mekongensis,B.sexangula,A.officinalis,and A.cucullata were indicated in Table 2.Tree density decreased with increasing age of the stand.The mortality rate of H.fomes was the highest than other species.The highest mortality was suffered by small size H.fomes(2.4%)trees in the suppressed crown class.No mortality was occurred in the rest of the species during the study period.Figure 1 showed the relationship between mean DBH and mean height of different plots of the studied stands.Tree height of the mangrove stands revealed positive correlation with tree DBH.This pattern means that tree height is increased with increasing DBH.

    Table 1 Stand structure data of mangrove communities along the Oligohaline zone of Sundarbans,Bangladesh

    Aboveground and belowground biomass

    Table 3 showed the biomass production and proportions of biomass apportioned to the different components of the trees.The mean values of above-ground and belowground biomass of the stand were 154.8±3.5 and 84.2 ± 1.9 Mg·ha?1·yr?1,respectively,which accounted for 64.8%and 35.2%,respectively,of the total tree biomass.Figure 2 showed the species-wise contribution in biomass allocation among the different components.Of these,A.officinalis accounted for 40.3%and 36.6%,which was thehighest contributor both in above-ground biomass and below-ground biomassproduction,respectively.The above-and below-ground biomass was plotted against structural parameters of the stand such as mean height,mean DBH,and basal area(Fig.3).Correlation analysis between biomass of the studied species with their respective mean height,mean DBH,and basal area exhibited the arrangement of the communities according to their degree of preponderance.The communities with high above-and below-ground biomass had greater values of basal area.The total biomass of the studied plots ranged from 22.8 Mg·ha?1(Plot 20)to 1135.0 Mg·ha?1(Plot 16).The total basal area of the plots ranged from 3.9 m2·ha?1(Plot 20)to 85.5 m2·ha?1(Plot 16).The extreme variation in basal area and biomass is partly the result of edge effects associated with the very small plot size.This pattern also indicates that basal area may be a predictor of total biomass of trees in the mangrove stands.Figure 4 showed a relationship between AGB/BGB and tree density.There is no correlation between them.

    Table 2 Specific density,species wise DBH,H and wood density in different plot location in mangrove communities along the Oligohaline zone of Sundarbans

    Fig.1 Relationship between mean DBH and height of trees in the study area(R2=0.79);●:2016;□:2017

    Litterfall

    Litterfall were observed throughout the year,and we were so far unable to separate the litterfall components at the species level.The study area consists of six mangrove species and out of these X.mekongensis and E.agallocha are deciduous species and the rest are evergreen species.Figure 5 showed the seasonality of the litterfall of all the mangrove species.Leaf-fall was the dominant component of litterfall throughout the year.There were definite peaks in each litterfall component. However,litterfall amounts of leaf,branch,and reproductive organs did not show any clear monthly patterns ofthe mixed mangrove species.Because leaf shedding and leaf emergence times were mixed among the evergreen and deciduous species.For example,the leaf shedding time of X.mekongensis started in December and continued until March,i.e.,winter and spring seasons,whereas leaf shedding of E.agallocha started from February and continued until May,i.e.,spring season.The highest peak of leaf litterfall was in February,i.e.,winter season and the next highest peak was observed in July,i.e.,in the late summer or rainy season.Branch litterfall was high in June-July(rainy season),and this may be due to strong rain and storms at that time.Reproductive organ litterfall was highest in March,i.e.,in the spring season when almost all studied mangrove species started their flowering.Young,tender,green,and yellow leaves contributed the most to total litterfall.Mean values of litterfall in leaf,branch,and reproductive organ(±SE)were 7.1 ± 0.05,1.3 ± 0.02 and 1.7 ± 0.02 Mg·ha?1·yr?1,respectively,which represented 70.2%,12.9%and 16.9%of the total litterfall,respectively.

    Net primary productivity

    The mean biomass increment and mean litterfall was 7.1 and 10.1 Mg·ha?1·yr?1,respectively for the mangrove community(Table 4).One-year measurement is usually not sufficient to obtain a reliable estimate of long term average mortality rate.Therefore,tree mortality is negligibly small in the present study and we simply summed the biomass increment of all the surviving trees and then added litterfallto get the above-ground net primary productivity for each smallplot.Above-ground net primary production(AGNPP)was 17.2 Mg·ha?1·yr?1.Mean below ground increment was estimated to be 3.8 Mg·ha?1·yr?1.Similarly,belowground production includes increment of roots of living trees plus turn of fine roots.Most of this turnover is the mortality of fine roots of living trees,same as the litterfall of above-ground components.But in this study,we were unable to measure the fine root production,so our assessment of below-ground production did not include fine root turnover.NPP for both aboveground and belowground combined was 21.0 Mg·ha?1·yr?1during the study period.AGNPP contributed 81.9% of the total NPP and belowground contributed 18.1%of the total NPP.

    Table 3 Biomass and its contribution in the mangrove forest along the Oligohaline zone of Sundarbans,Bangladesh

    Discussion

    Fig.2 Species-wise contribution in total biomass allocation of mangrove forest along the oligohaline zone of the Sundarbans,Bangladesh;solid dark column:AGB;open column:BGB

    Stand density(2305 ha?1)of the present study area was comparable to mature riverine mangroves or mixed mangrove stands(3310–3917 ha?1)at French Guiana(Fromard et al.1998).They reported higher density at pioneer stage of mangrove species than its mature stage.In the present study we observed that communities with small sized trees had very high density;while communities with large sized trees had lower density.This is in agreement with Fromard et al.(1998)suggesting that density was the most discriminating factor for early development stage of mangroves-a young stand matures by decreasing the number of individuals.The present study showed that mangrove communities with similar mean height and mean diameter produced different biomass production with different basal area(Fig.2).For example,mean height(6.5 m)and mean DBH(6.4 cm)of a particular stand showed a large difference in biomass production(127–279 Mg·ha?1)and(84–279 Mg·ha?1),respectively.This variation is due to the differences in basal area of the stand.These results are generally consistent with the findings of Fromard et al.(1998),who observed that large biomass differences(180–315 t·ha?1)for Avicennia germinans populations with similar height and diameter but with different basal area.

    Aboveground and belowground biomass

    This study showed that the mean above-ground biomass of 154.8 Mg·ha?1was within the documented range of 40.70–279.03 Mg·ha?1in the mangrove forestsofEastSumatra,Indonesia(Kusmanaetal.1992),or of 123.5–383.5 Mg·ha?1in the mangrove forests of Dominican Republic(Sherman et al.2003).The present value of AGB was almost identical to that(159 Mg·ha?1) ofRhizophora apiculata Bl.dominated mangrove forestin southern Thailand(Christensen 1978)and(163 Mg·ha?1)of B.gymnorrhiza dominated forestalong theOkukubiRiver,Okinawa Island,Japan(Kamruzzaman et al.2017).This value was higher than those reported for subtropical mixed mangrove forest at Ishigaki Island,southern Japan (78.6 Mg·ha?1;Suzuki and Tagawa 1983).The present value of accumulation of mean aboveground biomass was within the converted range ofAGB from aboveground biomasscarbon for Sundarbans mangrove forests(55.3–153.4 Mg·ha?1C;Chandaetal.2016).Compared to theprevious studies,the presentstudy indicatesthatbiomass productivity of the mangrove species in the present study was relatively high among the mangrove forests in the tropical and subtropical areas,showing that mean annual wood production of 7.1 Mg·ha?1·yr?1.A similar observation was made by Putz and Chan(1986),who reported that the average stem production over 31 years for a Rhizophora-Bruguiera forest in Malaysia was 6.7 Mg·ha?1·yr?1.The mean value of AGB increment was higher than the 5.9 Mg·ha?1·yr?1recorded on R.stylosa dominated mangrove forest at Missionary Bay,Hinchinbrook Channel,Australia(Clough 1998).It is important to mention that AGB of the studied mangrove communities is accounted for mainly due to stem biomass and is also a permanent indicator of biomass increment.Small branches and leaves mainly contributed in the litterfall production.

    Fig.3 a-f Relationships between biomass and main structural parameters of mangrove communities;●:2016;□:2017

    Fig.4 Relationship between AGB/BGB and density;●:2016;□:2017

    The mean above/below-ground biomass ratio of the present stand was 1.84.This value is higher than the records for Ceriops tagal(Perr.)C.B.Rob.at southern Thailand(1.05;Komiyama et al.2000)and for B.gymnorrhiza at Okinawa Island,Japan(1.3;Kamruzzaman et al.2017).This mean value is similar with records for other mangrove species at New South Wales(1.3–2.0;Saintilan 1997).The present value of above?/belowground ratio is lower than trees of the terrestrial inland forest of tropical areas(5.1–10.7;Ogawa et al.1965;Hozumi et al.1969).It indicates that a major portion of biomass is allocated to the underground parts in the present studied mangrove species than in terrestrial plants so the above?/below-ground ratio in Sundarbans mangroves is generally lower than terrestrial plants.Mangroves are usually coping a saline environment with the stress of high water tables but physiologically dry condition for the plants and deficiency of oxygen(Ball 1988;Havanond and Maxwell 1996),so a large portion of biomass is allocated to the underground parts of the mangrove species to adapt to the harsh environment.A low above/below-ground ratio of biomass was also reported in a larch ecosystem that grows on permaforest.Thislarch speciesmustalso copewith aharsh environment of limited soil nutrients and low soil temperature(Kajimoto et al.2010),showing similarity to the mangrove communities of the present study as well as to the mangrove forest at Okinawa Island,Japan(Kamruzzaman et al.2017).

    Fig.5 Seasonal changes of litterfall during the study period;solid dark column:leaf;open column:branch;vertical striping column:reproductive organs.The bar indicates standard error

    Table 4 Mean annual net primary production and its components of mangrove forest(Mg·ha?1·yr?1)along the Oligohaline zone of the Sundarbans,Bangladesh

    Litterfall

    This is the first report of litterfall production of the Sundarbans mangrove forest,Bangladesh.The mean total litterfall production(10.1 Mg·ha?1·yr?1)of the present mangrove forest was within the documented range for mangrove forests at Laguna de Términos,Mexico(8.4–12.5 Mg·ha?1·yr?1;Day et al.1987),in Unguja Island,Zanzibar,Tanzania(16 Mg·ha?1·yr?1;Shunula and Whittick 1999),and at Central Java,Indonesia(7.1–10.4;Sukardjo 1996).The mean total litterfall of the present study was higher than those recorded on southeastern Mexican mangrove forest(3.1–7.9 Mg·ha?1·yr?1;Day et al.1996).Mangrove leavesareshed continuouslythroughouttheyear especially in case of a community of mixed mangrove species.Day etal.(1987)reported that higher productivity in mangroves is related to higher input of nutrients.Thus,regular tidal inundation,lower salinity level,and regular nutrient input may be the factorsrelated to highermean productivity in the present study.This study showed that maximum leaf litterfall of the investigated species occurred in winter or dry season and late summer or rainy season.In many mangrove areas,similar findings were observed that peak litterfall occurs during the rainy season(Leach and Burgin 1985;Day et al.1987)and Flores-Verdugo et al.(1987)reported high litterfall at the end ofthe dry season.In the tropics,where temperatures are always favorable for tree growth,seasonal development is often not correlated with climate.In tropical rain forest many trees flush and flower at the same specie-specific time each year,but others do so at irregular intervals(Borchert 2017).The tropical climate could cause multimodal peaks of leaf litterfall in mangroves in this region.

    Net primary productivity

    This is the first report of net primary productivity of the Sundarbans Mangrove forest,Bangladesh.The present stand had above-ground biomass of 154.8 Mg·ha?1·yr?1,and above-ground net primary production during the study period was estimated to be 17.2 Mg·ha?1·yr?1.Those values were lower to those of AGB(159 Mg·ha?1·yr?1)and AGNPP(27 Mg·ha?1·yr?1)for R.apiculata dominated mangrove forest at Southern Thailand(Christensen 1978).Our estimates of AGNPP(17.2 Mg·ha?1·yr?1)were near the high end of the range of values reported in mangrove forests at Florida,USA(8.1–26.1 Mg·ha?1·yr?1;Ross et al.2001).The present value of NPP was higher than those reported for R.apiculata dominated forest at Matang mangrove forest,Malaysia(17.7 Mg·ha?1·yr?1;Putz and Chan 1986).The study area is located in the oligohaline zone of the Sundarbans mangrove forest and the area is flushed frequently by tides that may cause higher NPP of this mangroves.Similar findings were observed by Day et al.(1989),who reported that mangrove forest flushed frequently by tides and exposed to high nutrient concentrations have higher net primary production.There are three ecological zones in the Sundarbans such as oligohaline,mesohaline,and polyhaline zone based on the degree of salinity and floristic composition.The floristics composition of Sundarbans is defined by the distribution of three species:H.fomes,E.agallocha,and C.decandra.All three species occur throughout the Sundarbans but in different proportions depending on salinity.H.fomes is the characteristics of the oligohaline zone,E.agallocha of the mesohaline zone,and C.decandra of the polyhaline zone(Rahman and Islam 2015).Depending on the site of mangroves along the salinity gradient of an estuary and with distance inland from the shore,the properties of mangrove communities vary within an environmental setting(Chen and Twilley 1999).So in the present study,mangrove communities in the oligohaline zone may have higher productivity and turnover than other ecological zones of Sundarbans,Bangladesh.

    When comparing our estimated rates of litterfall and biomass increment both in AGB and BGB,it becomes evident that litterfall production only amount up to 48.1% of the total mangrove NPP.There are very few reports that describe the contribution of litterfall to total NPP,and reported that litterfall represents up to 30%of the overall NPP(Alongi et al.2005).NPP of the mangrove forest along the oligohaline zone of Sundarbans,Bangladesh,is conservative and the estimation was lower because we do not have data regarding coarse root and fine root production and their contribution is not included.As well as we do not have any data regarding vegetative damage including leaf and reproductive organs by direct consumption of the herbivory animals.The ratio of AGNPP and litterfall of the present study was 1.7:1.Our results agree with the summarization of Teas(1979),who reported that the assumption of total AGNPP for mangroves is three times as large as the amount of total litterfall.

    Conclusion

    It is not possible to clearly determine how much difference there are in the production of mangrove forests among different studies,due to differences in methodology,inter-annual variation,sitecondition,and stage of mangrove development.Moreover,there is no standard technique to compare the production capability of the mangrove stands all over the world.It is concluded that,the results of the first study at Sundarbans,Bangladesh indicatethatthebiomass accumulation (239 Mg·ha?1) and productivity(21.0 Mg·ha?1·yr?1)ofthe Sundarbans mangrove forest,Bangladesh contribute significantly to global carbon budget because higher biomass accumulation,higher productivity and large area of the remaining mangrove communities spanning two countries(Bangladesh and India)are related to higher carbon turnover of the stand in absolute terms.The unique characteristics of the Sundarbans mangrove forest also provide important information to decision makers on the strategy for sustainable managementofthe mangrove forests.

    This work was supported by a Grant-in-Aid for Scientific Research(JSPS KAKENHI 15F15389)from the Japan Society for the promotion of Science to Professor Akira Osawa and Md.Kamruzzaman.The authors are grateful to the Forest Department,Govt.of the People’s Republic of Bangladesh for establishing the research plots in the SRF and for assisting the research team in collecting data.The authors are grateful to Mrs.Minhaj-Uz-Siraj,Kalan Basak,and Sumanto Paul for their invaluable assistance during field data collection.

    Author’s contributions

    Md.K developed the concept of the article.Md.K also established the plots,forest inventory and acquired the data.SA helped in data collection.Md.K worked on the manuscript preparation and AO edited the draft manuscript and all authors read and approved the final manuscript.

    Competing interests

    The authors declare that they have no competing interests.

    Ahmed A,Aziz A,Khan AZMNA,Islam MN,Iqbal KF,Nazma IMS(2011)Tree diversity as affected by salinity in the Sundarban mangrove forests,Bangladesh.Bangladesh J Bot 40:197–202

    Alongi DM,Clough BF,Robertson AI(2005)Nutrient-use efficiency in arid-zone forests of the mangroves Rhizophora stylosa and Avicennia marina.Aquat Bot 82:121–131

    Ball MC(1988)Salinity tolerance in mangroves Aegiceras corniculatum and Avicennia marina.I.Water use in relation to growth,carbon partitioning,and salt balance.Aust J Plant Physiol 15:447–464

    Blasco F(1977)Outlines of ecology,botany and forestry of the mangals of the Indian subcontinent.In:Chapman VJ(ed)Wet coastal ecosystems,ecosystems of the world,vol 1.Elsevier,Amsterdam,pp 241–260

    Borchert R(2017)The phenology of tropical trees.http://www.borchert.faculty.ku.edu.Accessed 18 May 2017

    Chaffey DR,Miller FR,Sandom JH(1985)A forest inventory of the Sundarban,Bangladesh.Project Report 140,Land Resources Development Centre,Surrey,England.196

    Chanda A,Mukhopadhyay A,Ghosh T,Akhand A,Mondal P,Ghosh S,Mukherjee S,Wolf J,Lázár AN,Rahman MM,Salehin M,Chowdhury SM,Hazra S(2016)Blue carbon stock of the Bangladesh Sundraban mangroves:what could be the scenario after a century?Wetlands 36:1033–1045

    Chave J,Andalo C,Brown S,Cairns MA,Chambers JQ,Eamus D,F?lster H,Fromard F,Higuchi N,Kira T,Lescure JP,Nelson BW,Ogawa H,Puig H,Riéra B,Yamakura T(2005)Tree algometry and improved estimation of carbon density and balance in tropical forests.Oecologia 145:87–99

    Chave J,Coomes DA,Jansen S,Lewis SL,Swenson NG,Zanne AE(2009)Towards a worldwide wood economics spectrum.Ecol Lett 12:351–366

    Chen R,Twilley RR(1999)Patterns of mangrove forest structure and oil nutrient dynamics along the shark river estuary,Florida.Estuaries 22:955–970

    Christensen B(1978)Biomass and primary production of Rhizophora apiculata Bl.In a mangrove in southern Thailand.Aquat Bot 4:43–52

    Clough B(1998)Mangrove forest productivity and biomass accumulation in Hinchinbrook Channel,Australia.Mangr Salt Marsh 2:191–198

    Curtis SJ(1933)Working plan for the forests of the Sundarbans Division for the period from 1stApril 1931 to 31stMarch 1951,Voulme III Part of Appendix III Description of the Compartments and their Histories.Government of Bengal

    Day JW,Conner WH,Ley-Lou F,Day RH,Navarro AM(1987)The productivity and composition of mangrove forests,Laguna De Términos,Mexico.Aquat Bot 27:267–284

    Day JW,Coronado-Molina C,Vera-Herrera FR,Twilley R,Rivera-Monroy VH,Alvarez-Guillen H,Day R,Conner W(1996)A 7 year record of aboveground net primary production in a southern Mexican mangrove forest.Aquat Bot 55:39–60

    Day JW,Hall C,Kemp M,Yáňez-Arancibia A(1989)Estuarine ecology.Wiley Interscience,New York,p 558

    Flores-Verdugo FJ,Day JW,Briseňo-Dueňas R(1987)Structure,litterfall,decomposition,and detritus dynamics of mangroves in a Mexican coastal lagoon with ephemeral inlet.Mar Ecol Prog Ser 35:51–56

    Fromard F,Puig H,Mougin E,Marty G,Betoulle JL,Cadamuro L(1998)Structure,above-ground biomass and dynamics of mangrove ecosystems:new data from French Guiana.Oecologia 115:39–53

    Havanond S,Maxwell GS(1996)Strategies for mangrove restoration.FORTROP 10:21–36

    Hernádez CMA,Zaragoza CG,Iriarte-Vivar S,Flores-Verdugo FJ,Casasola PM(2011)Forest structure,productivity and species phenology of mangroves in the la Mancha lagoon in the Atlantic coast of Mexico.Wetl Ecol Manag 19:273–293

    Hossain M,Siddique MRH,Saha S,Abdullah SMR(2015)Allometric models for biomass,nutrients and carbon stock in Excoecaria agallocha of the Sundarbans,Bangladesh.Wetl Ecol Manag 23:765–777

    Hozumi K,Yoda K,Kokawa S,Kira T(1969)Production ecology of tropical rain forest in southern Cambodia I.Plant biomass.Nat Life Se Asia 6:1–51

    Iftekhar MS,Saenger P(2008)Vegetation dynamics in the Bangladesh Sundarbans mangroves:a review of forest inventories.Wetl Ecol Manag 16:291–312

    Islam MS(2001)Sea-level changes in Bangladesh:the last ten thousand years.Asiatic Society of Bangladesh,Dhaka

    Kajimoto T,Osawa A,Usoltsev VA,Abaimov AP(2010)Biomass and productivity of Siberian Larch forest ecosystems.In:Osawa A,Zyryanova OA,Matsuura Y,Kajimoto T,Wein RW(eds.)Permafrost Ecosystems:Siberian Larch forests.Ecological studies,Springer 209:99–122

    Kamruzzaman M,Osawa A,Deshar R,Sharma S,Mouctar K(2017)Species composition,biomass,and net primary productivity of mangrove forest in Okukubi River,Okinawa Island,Japan.Regional Studies in Marine Science 12:19–27

    Komiyama A,Havanond S,Srisawatt W,Mochida Y,Fujimoto K,Ohnishi T,Ishihara S,Miyagi T(2000)Top/root biomass of a secondary mangrove(Ceriops tagal(Perr.)C.B.Rob.)forest.Forest Ecol Manag 139:127–134

    Komiyama A,Ogino K,Aksomkoae S,Sabhasri S(1987)Root biomass of a mangrove forest in southern Thailand I.Estimation by the trench method and the zonal structure of root biomass.J Trop Ecol 3:97–108

    Komiyama A,Poungparn S,Kato S(2005)Common allometric equations for estimating the tree weight of mangroves.J Trop Ecol 21:471–477

    Kusmana C,Sabiham S,Abe K,Watanabe H(1992)An estimation of above ground tree biomass of a mangrove forest in East Sumatra,Indonesia.Tropics 1:243–257

    Leach GJ,Burgin S(1985)Litter production and seasonality of mangroves in Papua New Guinea.Aquat Bot 23:215–224

    Lugo AE,Snedaker SC(1975)Properties of a mangrove forest in southern Florida.In:Walsh G,Snedaker S,Teas H(eds)Proceedings of the international symposium on biology and management of mangroves.University of Florida,Gainesville,pp 170–212

    Mackey AP,Smail G(1995)Spatial and temporal variation in litterfall of Avicennia marina(Forssk.)Vierh.In the Brisbane River,Queensland,Australia.Aquat Bot 52:133–142

    Matsuura Y,Kajimoto T(2013)Measurement and analysis of Carbon pool and flow by the ecological approach.In:Oikawa T,Yamamoto S(eds)Carbon dynamics of terrestrial ecosystem:Systems approach to global environment.Kyoto University press,p 414

    Nagarajan B,Pandiarajan C,Krishnamoorthy M,Sophia P(2008)Reproductive fitness and success in mangroves:implication on conservation.In:Sengupta M,Dalwani R(eds)Proceedings of Tall 2007:the 12thWorld Lake Conference,pp 29–33

    Ogawa H(1977)Principles and methods of estimating primary production in forests.In:Shidei T,Kira T(eds)Primary productivity of Japanese forests:productivity of terrestrial communities,JIBP synthesis.University of Tokyo Press,pp 29–37

    Ogawa H,Yoda K,Ogino K,Kira T(1965)Comparative ecological studies on three main type of forest vegetation in Thailand.Nat Life Se Asia 4:49–80

    Putz FE,Chan HT(1986)Tree growth,dynamics,and productivity in a mature mangrove forest in Malaysia.Forest Ecol Manag 17:211–230

    Rahman MR,Asaduzzaman M(2010)Ecology of Sundarban,Bangladesh J Sci Foundation 8:35–47

    Rahman MM,Islam SA(2015)Phenophase of five mangrove species of the Sundarbans of Bangladesh.IJBSC Res 4:77–82

    Rahman MM,Khan MNI,Hoque AKF,Ahmed I(2015)Carbon stock in the Sundarbans mangrove forest:spatial variations in vegetation types and salinity zones.Wetl Ecol Manag 23:269–283

    Ross MS,Ruiz PL,Telesnicki GJ,Meeder JF(2001)Estimating above-ground biomass and production in mangrove communities of Biscayne National Park,Florida(USA).Wetl Ecol Manag 9:27–37

    Saintilan N(1997)Above-and below-ground biomasses of two species of mangrove on the Hawkesbury River estuary,new South Wales.Mar Freshw Res 48:147–152

    Sherman RE,Fahey TJ,Martinez P(2003)Spatial patterns of biomass and aboveground net primary productivity in a mangrove ecosystem in the Dominican Republic.Ecosystems 6:384–398

    Shunula JP,Whittick A(1999)Aspects of litter production in mangroves from Unguja Island,Zanzibar,Tanzania.Estuar Coast Shelf Sci 49:51–54

    Siddiqi NA(2001)Mangrove forestry in Bangladesh.Institute of Forestry and Environmental Sciences.University of Chittagong,Chittagong

    Sukardjo S(1996)The relationship of litterfall to basal area and climatic variables in the Rhizophora Mucronata Lamarck plantation at Tritih,central java,Indonesia.J Southeast Asian Stud 34:424–432

    Suzuki E,Tagawa H(1983)Biomass of a mangrove forest and a sedge marsh on Ishigaki Island,South Japan.Jap J Ecol 33:231–234

    Tamai S,Nakasuga T,Tabuchi R,Ogino K(1986)Standing biomass of mangrove forests in southern Thailand.J Jpn Fo Soc 68:384–388

    Teas HJ(1979)Silviculture with saline water.In:Hollander A,Aller JC,Epstein E,Pietro AS,Zaborsky OR(eds)The baseline concept.Plenum Press,New York,pp 117–161

    *Correspondence:zaman_dulu@yahoo.com

    1Laboratory of Ecosystem Production and Dynamics,Graduate School of Global Environmental Studies,Kyoto University,Kyoto 606-8501,Japan

    2Forestry and Wood Technology Discipline,Khulna University,Khulna 9208,Bangladesh

    ?The Author(s).2017Open AccessThis 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.

    6 June 2017 Accepted:29 August 2017

    亚洲精品美女久久久久99蜜臀 | 99国产综合亚洲精品| 亚洲成国产人片在线观看| 午夜视频国产福利| 国产精品嫩草影院av在线观看| 亚洲综合色惰| 熟女电影av网| 久久久精品免费免费高清| 亚洲欧洲日产国产| 亚洲欧美成人精品一区二区| 国产成人精品一,二区| www日本在线高清视频| 国产日韩一区二区三区精品不卡| 超色免费av| 亚洲精品美女久久av网站| av在线app专区| 中文字幕精品免费在线观看视频 | 久久久国产一区二区| 免费人妻精品一区二区三区视频| 日本爱情动作片www.在线观看| www.av在线官网国产| 国产一区二区三区综合在线观看 | 日韩一本色道免费dvd| 高清视频免费观看一区二区| 一级黄片播放器| 久久午夜福利片| 中文字幕最新亚洲高清| 青春草亚洲视频在线观看| 免费av中文字幕在线| 在线观看三级黄色| 一级毛片 在线播放| 国产男女内射视频| 中文精品一卡2卡3卡4更新| 日本午夜av视频| 中国三级夫妇交换| 久久午夜综合久久蜜桃| 久久青草综合色| 美女福利国产在线| 熟女人妻精品中文字幕| 欧美日韩av久久| 午夜老司机福利剧场| 亚洲 欧美一区二区三区| 亚洲av电影在线观看一区二区三区| 视频区图区小说| 日本vs欧美在线观看视频| 国产 一区精品| 国产精品偷伦视频观看了| 99国产综合亚洲精品| 国产成人精品福利久久| 午夜影院在线不卡| 国产精品一二三区在线看| 免费黄色在线免费观看| 看免费av毛片| 日本爱情动作片www.在线观看| 午夜久久久在线观看| 如日韩欧美国产精品一区二区三区| 亚洲激情五月婷婷啪啪| 视频在线观看一区二区三区| 成年人免费黄色播放视频| 国产成人午夜福利电影在线观看| 成人综合一区亚洲| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线| 大话2 男鬼变身卡| 午夜福利网站1000一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产精品久久久久久久电影| 满18在线观看网站| 99国产综合亚洲精品| 波野结衣二区三区在线| av在线老鸭窝| 婷婷色av中文字幕| 亚洲欧美中文字幕日韩二区| 婷婷色综合www| 97人妻天天添夜夜摸| 欧美亚洲日本最大视频资源| 国产深夜福利视频在线观看| 女的被弄到高潮叫床怎么办| 亚洲欧美清纯卡通| 美女内射精品一级片tv| 欧美精品国产亚洲| 日本色播在线视频| 女人久久www免费人成看片| 国产一区二区激情短视频 | 亚洲第一区二区三区不卡| 国产毛片在线视频| 亚洲丝袜综合中文字幕| 亚洲国产精品999| 80岁老熟妇乱子伦牲交| 国产精品不卡视频一区二区| 国产欧美日韩一区二区三区在线| 国产乱来视频区| 一边摸一边做爽爽视频免费| 亚洲欧美日韩另类电影网站| 美国免费a级毛片| 久久 成人 亚洲| 9191精品国产免费久久| 丰满少妇做爰视频| 大片电影免费在线观看免费| 最近最新中文字幕免费大全7| 国产一区二区三区av在线| 香蕉国产在线看| 少妇人妻 视频| 亚洲精品,欧美精品| 欧美精品一区二区大全| 男女边摸边吃奶| 免费大片黄手机在线观看| 天天影视国产精品| 涩涩av久久男人的天堂| 成年动漫av网址| 99国产精品免费福利视频| 亚洲综合色惰| 成人毛片a级毛片在线播放| 久久这里只有精品19| 久久久久久久精品精品| 麻豆乱淫一区二区| 午夜视频国产福利| 热99久久久久精品小说推荐| 欧美xxⅹ黑人| 搡女人真爽免费视频火全软件| 一级毛片电影观看| 久久国内精品自在自线图片| 999精品在线视频| 考比视频在线观看| 国产av精品麻豆| av福利片在线| 中国国产av一级| 最后的刺客免费高清国语| 中文字幕制服av| 欧美激情 高清一区二区三区| 国产免费一级a男人的天堂| 极品人妻少妇av视频| 久久精品久久精品一区二区三区| 黄色视频在线播放观看不卡| av在线app专区| 国产成人精品福利久久| 国产黄频视频在线观看| 高清av免费在线| 一本久久精品| 午夜精品国产一区二区电影| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 熟妇人妻不卡中文字幕| 日韩成人伦理影院| 欧美xxⅹ黑人| 国产精品秋霞免费鲁丝片| 久久精品久久精品一区二区三区| a级毛色黄片| 亚洲人成网站在线观看播放| 精品福利永久在线观看| 69精品国产乱码久久久| www.熟女人妻精品国产 | 久久久久久久久久人人人人人人| 欧美亚洲 丝袜 人妻 在线| 国产亚洲午夜精品一区二区久久| 最后的刺客免费高清国语| 高清毛片免费看| 人人妻人人爽人人添夜夜欢视频| freevideosex欧美| 夫妻午夜视频| 亚洲av免费高清在线观看| 伊人久久国产一区二区| 看免费成人av毛片| 黄片无遮挡物在线观看| 国产精品欧美亚洲77777| 亚洲精品日本国产第一区| 丁香六月天网| 老司机影院毛片| 又大又黄又爽视频免费| 国产成人a∨麻豆精品| 成年av动漫网址| 亚洲精品美女久久久久99蜜臀 | 香蕉国产在线看| 亚洲av.av天堂| 国产片内射在线| 午夜免费鲁丝| 最近2019中文字幕mv第一页| 成人综合一区亚洲| 久久人妻熟女aⅴ| 爱豆传媒免费全集在线观看| 久久97久久精品| 久热久热在线精品观看| a级片在线免费高清观看视频| av网站免费在线观看视频| 国产欧美另类精品又又久久亚洲欧美| 99久久综合免费| 久久久久久久大尺度免费视频| 十八禁网站网址无遮挡| 午夜av观看不卡| 丰满迷人的少妇在线观看| 秋霞在线观看毛片| 国产成人一区二区在线| 亚洲国产欧美在线一区| 狠狠精品人妻久久久久久综合| 制服丝袜香蕉在线| 免费观看性生交大片5| 一级爰片在线观看| 看非洲黑人一级黄片| 在现免费观看毛片| 精品少妇黑人巨大在线播放| 大片电影免费在线观看免费| 国产乱人偷精品视频| a级毛片在线看网站| 国产成人a∨麻豆精品| 精品人妻在线不人妻| 国产一区二区三区av在线| 又黄又粗又硬又大视频| 五月开心婷婷网| 男女边吃奶边做爰视频| 国产色爽女视频免费观看| 91国产中文字幕| 性色avwww在线观看| a级片在线免费高清观看视频| 亚洲av男天堂| 国产欧美日韩一区二区三区在线| 汤姆久久久久久久影院中文字幕| 日韩欧美精品免费久久| 国产一区二区三区综合在线观看 | 亚洲欧美成人精品一区二区| 女人精品久久久久毛片| 晚上一个人看的免费电影| 国产福利在线免费观看视频| 最后的刺客免费高清国语| 欧美激情 高清一区二区三区| 女人精品久久久久毛片| 成年av动漫网址| 1024视频免费在线观看| 五月天丁香电影| 久久久久精品人妻al黑| www.熟女人妻精品国产 | 91久久精品国产一区二区三区| 成年av动漫网址| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 黄色配什么色好看| 国产精品一区二区在线不卡| 日本91视频免费播放| 亚洲熟女精品中文字幕| 欧美成人午夜精品| 中文天堂在线官网| av不卡在线播放| 国产午夜精品一二区理论片| 日韩中文字幕视频在线看片| 亚洲精品一区蜜桃| 日日啪夜夜爽| 视频区图区小说| 少妇熟女欧美另类| 国产精品 国内视频| 成人无遮挡网站| 久久av网站| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产av蜜桃| 亚洲一区二区三区欧美精品| 高清欧美精品videossex| 欧美日本中文国产一区发布| 黑人欧美特级aaaaaa片| 18+在线观看网站| 我的女老师完整版在线观看| 免费黄网站久久成人精品| 视频在线观看一区二区三区| 国产 精品1| 国产高清不卡午夜福利| 久久 成人 亚洲| 亚洲精品一二三| xxx大片免费视频| 久久久久久久国产电影| 哪个播放器可以免费观看大片| 亚洲国产色片| 天堂中文最新版在线下载| 国产免费视频播放在线视频| 国产一区二区三区综合在线观看 | 女人久久www免费人成看片| 亚洲欧洲精品一区二区精品久久久 | 国产乱人偷精品视频| 91aial.com中文字幕在线观看| 亚洲欧美中文字幕日韩二区| 全区人妻精品视频| 日韩欧美精品免费久久| 亚洲av电影在线进入| 观看美女的网站| 两个人免费观看高清视频| 一本色道久久久久久精品综合| 如何舔出高潮| 中文欧美无线码| 大香蕉97超碰在线| 狠狠婷婷综合久久久久久88av| 啦啦啦中文免费视频观看日本| 久久久欧美国产精品| 欧美精品国产亚洲| 日本欧美视频一区| 99久久中文字幕三级久久日本| 三上悠亚av全集在线观看| 亚洲精品,欧美精品| 国产 精品1| 日本-黄色视频高清免费观看| 丝袜美足系列| 五月开心婷婷网| 天堂8中文在线网| 秋霞伦理黄片| √禁漫天堂资源中文www| 涩涩av久久男人的天堂| 七月丁香在线播放| 久久精品国产鲁丝片午夜精品| 哪个播放器可以免费观看大片| 秋霞在线观看毛片| 这个男人来自地球电影免费观看 | 在线观看免费高清a一片| 亚洲,欧美,日韩| 国产亚洲最大av| 国产色婷婷99| 一本色道久久久久久精品综合| 久久精品国产亚洲av天美| 婷婷色综合大香蕉| 国产一级毛片在线| 欧美日韩视频精品一区| 亚洲精品456在线播放app| 最新的欧美精品一区二区| a级毛片黄视频| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 人体艺术视频欧美日本| 亚洲精品国产av成人精品| 免费av不卡在线播放| 日本av手机在线免费观看| 一区二区日韩欧美中文字幕 | 久久精品国产a三级三级三级| 观看av在线不卡| 免费看av在线观看网站| 91aial.com中文字幕在线观看| 日韩一区二区三区影片| 久久99精品国语久久久| 人妻系列 视频| a级毛片在线看网站| 国产片特级美女逼逼视频| 九九在线视频观看精品| 一本大道久久a久久精品| 夫妻性生交免费视频一级片| 性色avwww在线观看| 夫妻性生交免费视频一级片| 精品久久久精品久久久| 狠狠婷婷综合久久久久久88av| 丰满少妇做爰视频| 精品熟女少妇av免费看| 亚洲欧美清纯卡通| 寂寞人妻少妇视频99o| 中文字幕精品免费在线观看视频 | 免费黄色在线免费观看| 久久99精品国语久久久| 欧美亚洲日本最大视频资源| 26uuu在线亚洲综合色| 欧美成人精品欧美一级黄| 国产精品无大码| 99热这里只有是精品在线观看| 精品人妻偷拍中文字幕| 国产日韩欧美视频二区| 如日韩欧美国产精品一区二区三区| 天堂俺去俺来也www色官网| 久久久亚洲精品成人影院| av女优亚洲男人天堂| 成人毛片60女人毛片免费| 国产亚洲精品久久久com| av.在线天堂| 婷婷色综合www| 国产精品久久久久久久电影| 99久久精品国产国产毛片| 中文字幕人妻丝袜制服| 在现免费观看毛片| a级毛片在线看网站| 国产黄色视频一区二区在线观看| 最新的欧美精品一区二区| 国产日韩欧美视频二区| 男人操女人黄网站| 内地一区二区视频在线| 国产成人免费观看mmmm| 又黄又爽又刺激的免费视频.| 各种免费的搞黄视频| 国产精品久久久av美女十八| 欧美精品国产亚洲| 少妇的逼水好多| 亚洲欧洲精品一区二区精品久久久 | 久久影院123| 尾随美女入室| 国产av国产精品国产| 国产一区二区激情短视频 | 高清在线视频一区二区三区| 色吧在线观看| 一级爰片在线观看| 免费av中文字幕在线| 少妇被粗大的猛进出69影院 | 青春草视频在线免费观看| 一边亲一边摸免费视频| 精品人妻在线不人妻| 国产亚洲最大av| 欧美激情 高清一区二区三区| 久久综合国产亚洲精品| 亚洲中文av在线| 嫩草影院入口| 国产精品久久久久久久电影| 丝袜脚勾引网站| 寂寞人妻少妇视频99o| 国产片特级美女逼逼视频| 22中文网久久字幕| 久久久精品区二区三区| 国产极品粉嫩免费观看在线| 精品亚洲成国产av| 有码 亚洲区| 免费看不卡的av| 久久综合国产亚洲精品| 女性生殖器流出的白浆| 又黄又粗又硬又大视频| 日韩一本色道免费dvd| 韩国av在线不卡| 91aial.com中文字幕在线观看| 久热久热在线精品观看| 日本vs欧美在线观看视频| 亚洲国产毛片av蜜桃av| 黑人猛操日本美女一级片| 久久久欧美国产精品| 人人澡人人妻人| 大香蕉久久网| 99久久精品国产国产毛片| 男女下面插进去视频免费观看 | 久久久久精品久久久久真实原创| 老女人水多毛片| 亚洲av男天堂| 国产xxxxx性猛交| 国国产精品蜜臀av免费| 99re6热这里在线精品视频| 高清不卡的av网站| 人体艺术视频欧美日本| 成人无遮挡网站| videosex国产| 午夜免费鲁丝| 久久亚洲国产成人精品v| 精品一区二区三区四区五区乱码 | 国内精品宾馆在线| 人人妻人人澡人人爽人人夜夜| 人人妻人人添人人爽欧美一区卜| 成人午夜精彩视频在线观看| 一本—道久久a久久精品蜜桃钙片| 三上悠亚av全集在线观看| 国产精品偷伦视频观看了| 最近的中文字幕免费完整| 亚洲精品aⅴ在线观看| 亚洲精品国产av成人精品| 亚洲经典国产精华液单| 校园人妻丝袜中文字幕| 国产成人午夜福利电影在线观看| 少妇的逼好多水| 少妇精品久久久久久久| 国产麻豆69| 在线亚洲精品国产二区图片欧美| 一区二区三区精品91| 人妻少妇偷人精品九色| 日本-黄色视频高清免费观看| 精品久久久精品久久久| 欧美最新免费一区二区三区| av国产精品久久久久影院| 亚洲精品一二三| av一本久久久久| 两个人免费观看高清视频| 国产男女内射视频| 成年动漫av网址| 黑人欧美特级aaaaaa片| 老司机影院毛片| 亚洲国产色片| 90打野战视频偷拍视频| 日日摸夜夜添夜夜爱| 少妇高潮的动态图| 久久久久网色| 亚洲精品aⅴ在线观看| 久久久精品区二区三区| 日本vs欧美在线观看视频| 日韩精品免费视频一区二区三区 | 久热这里只有精品99| 色视频在线一区二区三区| 国产高清不卡午夜福利| 国产精品99久久99久久久不卡 | 2022亚洲国产成人精品| 欧美成人午夜免费资源| 有码 亚洲区| 精品少妇黑人巨大在线播放| 春色校园在线视频观看| 亚洲精品视频女| 国产精品蜜桃在线观看| 女人精品久久久久毛片| 18+在线观看网站| 91精品国产国语对白视频| 大话2 男鬼变身卡| 另类精品久久| 亚洲av.av天堂| 欧美变态另类bdsm刘玥| 少妇猛男粗大的猛烈进出视频| 最近最新中文字幕大全免费视频 | 中国三级夫妇交换| 18禁国产床啪视频网站| 精品少妇内射三级| a级毛片在线看网站| 中文字幕亚洲精品专区| 大片电影免费在线观看免费| 免费看光身美女| 欧美 亚洲 国产 日韩一| 欧美亚洲 丝袜 人妻 在线| 亚洲图色成人| www.熟女人妻精品国产 | 美女国产高潮福利片在线看| 亚洲天堂av无毛| 久久久久久久亚洲中文字幕| 久久影院123| 黑丝袜美女国产一区| 热99国产精品久久久久久7| 久久99蜜桃精品久久| av在线app专区| 亚洲欧洲精品一区二区精品久久久 | 热99国产精品久久久久久7| 日本与韩国留学比较| 日日啪夜夜爽| 久久午夜综合久久蜜桃| 99国产综合亚洲精品| 欧美成人午夜精品| 日日撸夜夜添| 蜜桃在线观看..| 精品国产国语对白av| 亚洲一区二区三区欧美精品| 国产乱人偷精品视频| 亚洲精品国产色婷婷电影| 国产精品一区二区在线观看99| 亚洲欧美一区二区三区国产| 久久这里只有精品19| 亚洲人成网站在线观看播放| 18禁在线无遮挡免费观看视频| 一级,二级,三级黄色视频| 亚洲精品国产av蜜桃| 有码 亚洲区| 五月开心婷婷网| 97人妻天天添夜夜摸| 日韩不卡一区二区三区视频在线| 成人午夜精彩视频在线观看| 美女内射精品一级片tv| 国产极品天堂在线| 日本91视频免费播放| 黄片播放在线免费| 狠狠婷婷综合久久久久久88av| 欧美3d第一页| 国产成人精品婷婷| 亚洲,欧美,日韩| 日本欧美国产在线视频| 五月伊人婷婷丁香| 午夜福利乱码中文字幕| 精品少妇内射三级| 国产极品天堂在线| 国产色爽女视频免费观看| 人人澡人人妻人| 亚洲国产精品一区二区三区在线| 欧美日韩亚洲高清精品| 成人国产av品久久久| 精品99又大又爽又粗少妇毛片| 精品久久国产蜜桃| 女性生殖器流出的白浆| 少妇 在线观看| 精品一区二区免费观看| 人妻一区二区av| 热99国产精品久久久久久7| 日韩中字成人| 欧美变态另类bdsm刘玥| 日韩欧美精品免费久久| 如日韩欧美国产精品一区二区三区| 免费少妇av软件| 免费观看性生交大片5| 男的添女的下面高潮视频| 高清毛片免费看| 久久午夜综合久久蜜桃| 午夜福利在线观看免费完整高清在| 免费看av在线观看网站| 高清av免费在线| 亚洲经典国产精华液单| 99热这里只有是精品在线观看| 久久久久视频综合| 精品人妻熟女毛片av久久网站| 9热在线视频观看99| 日韩制服骚丝袜av| 人人妻人人爽人人添夜夜欢视频| 亚洲精品日本国产第一区| 少妇人妻 视频| a级片在线免费高清观看视频| 大香蕉97超碰在线| 久久精品久久久久久噜噜老黄| 午夜免费男女啪啪视频观看| 最后的刺客免费高清国语| 午夜av观看不卡| 久久这里有精品视频免费| 日韩欧美精品免费久久| 哪个播放器可以免费观看大片| 久久这里有精品视频免费| 精品久久久精品久久久| 亚洲熟女精品中文字幕| 久久久精品免费免费高清| 一本久久精品| 男人操女人黄网站| 成人午夜精彩视频在线观看| 午夜福利影视在线免费观看| 美女主播在线视频| 久久久欧美国产精品| 国语对白做爰xxxⅹ性视频网站| 日韩一本色道免费dvd| 久久久久久久大尺度免费视频| 日韩三级伦理在线观看| 18禁国产床啪视频网站| 99久久综合免费| 国产精品 国内视频| 纯流量卡能插随身wifi吗| 9热在线视频观看99| 欧美亚洲 丝袜 人妻 在线| 久久久久久久大尺度免费视频| 中国三级夫妇交换| 一区二区日韩欧美中文字幕 |