• <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

    国产精品免费一区二区三区在线| 久久这里只有精品19| 亚洲九九香蕉| 欧美激情高清一区二区三区| 女性生殖器流出的白浆| 欧美日韩乱码在线| 少妇裸体淫交视频免费看高清 | 免费在线观看视频国产中文字幕亚洲| 欧美性长视频在线观看| 亚洲av五月六月丁香网| 妹子高潮喷水视频| 此物有八面人人有两片| 久久国产乱子伦精品免费另类| 美女免费视频网站| 色综合亚洲欧美另类图片| 国产精品日韩av在线免费观看| 欧美一区二区精品小视频在线| 久久久久九九精品影院| 精品乱码久久久久久99久播| av片东京热男人的天堂| 黄片大片在线免费观看| 亚洲国产欧洲综合997久久, | 免费av毛片视频| 国产午夜精品久久久久久| 女性生殖器流出的白浆| 久久久久久久久久黄片| 搡老岳熟女国产| 国产三级黄色录像| 久久精品91无色码中文字幕| 青草久久国产| 男人操女人黄网站| 日本一区二区免费在线视频| 久99久视频精品免费| 亚洲国产精品久久男人天堂| 中文字幕人妻熟女乱码| 精品久久久久久成人av| 99国产精品99久久久久| 欧美绝顶高潮抽搐喷水| av免费在线观看网站| 国产午夜福利久久久久久| 日韩免费av在线播放| 看免费av毛片| 久久久久久久久免费视频了| 老司机深夜福利视频在线观看| 成人国产综合亚洲| 亚洲中文字幕一区二区三区有码在线看 | 国产成人精品久久二区二区免费| 757午夜福利合集在线观看| 欧美日韩黄片免| 亚洲美女黄片视频| 黄色女人牲交| 美女高潮到喷水免费观看| 亚洲精品中文字幕一二三四区| 国产亚洲av嫩草精品影院| 久久久久久亚洲精品国产蜜桃av| 亚洲第一电影网av| 精品久久久久久久人妻蜜臀av| 岛国在线观看网站| 精品国产亚洲在线| 久久人人精品亚洲av| 午夜免费观看网址| 极品教师在线免费播放| 一边摸一边抽搐一进一小说| 热99re8久久精品国产| 国产激情久久老熟女| 50天的宝宝边吃奶边哭怎么回事| 国产免费男女视频| 妹子高潮喷水视频| 波多野结衣高清无吗| 久久99热这里只有精品18| 亚洲av日韩精品久久久久久密| 日韩av在线大香蕉| 国产熟女午夜一区二区三区| 精品日产1卡2卡| 精品第一国产精品| 男女那种视频在线观看| 欧美在线一区亚洲| 欧美日韩瑟瑟在线播放| 高清毛片免费观看视频网站| 国产成年人精品一区二区| 91成年电影在线观看| 免费在线观看完整版高清| 亚洲狠狠婷婷综合久久图片| 久久精品91无色码中文字幕| xxxwww97欧美| 国产高清videossex| 午夜日韩欧美国产| 女同久久另类99精品国产91| 欧美三级亚洲精品| 亚洲国产毛片av蜜桃av| 国产蜜桃级精品一区二区三区| 国产一区二区激情短视频| 亚洲av第一区精品v没综合| 真人一进一出gif抽搐免费| 国产亚洲欧美精品永久| 一级毛片精品| av免费在线观看网站| 午夜视频精品福利| 亚洲欧美日韩无卡精品| 中文字幕人妻丝袜一区二区| а√天堂www在线а√下载| 给我免费播放毛片高清在线观看| 狂野欧美激情性xxxx| 国产极品粉嫩免费观看在线| 亚洲一区二区三区不卡视频| 制服人妻中文乱码| 国产亚洲精品av在线| 老司机午夜十八禁免费视频| 国产又黄又爽又无遮挡在线| 黄色成人免费大全| 黄片播放在线免费| 亚洲第一电影网av| av超薄肉色丝袜交足视频| 18美女黄网站色大片免费观看| 国产欧美日韩一区二区精品| 18禁黄网站禁片午夜丰满| 哪里可以看免费的av片| 午夜a级毛片| 国产精品爽爽va在线观看网站 | 老鸭窝网址在线观看| 久99久视频精品免费| 久久草成人影院| 黑人欧美特级aaaaaa片| 国产激情欧美一区二区| 黄色视频,在线免费观看| 久久久水蜜桃国产精品网| 亚洲精品国产区一区二| 国产精品综合久久久久久久免费| 精品久久久久久,| 精品国产美女av久久久久小说| 久久精品影院6| 亚洲黑人精品在线| 桃红色精品国产亚洲av| 久久久久久亚洲精品国产蜜桃av| 日韩大尺度精品在线看网址| 午夜激情av网站| 12—13女人毛片做爰片一| 中文字幕久久专区| 亚洲自偷自拍图片 自拍| 黄网站色视频无遮挡免费观看| 波多野结衣av一区二区av| 一进一出好大好爽视频| av有码第一页| 精品福利观看| a级毛片在线看网站| 免费高清在线观看日韩| 少妇被粗大的猛进出69影院| 他把我摸到了高潮在线观看| 久久人人精品亚洲av| 国产亚洲精品第一综合不卡| 又黄又粗又硬又大视频| 侵犯人妻中文字幕一二三四区| av在线播放免费不卡| videosex国产| tocl精华| 成人18禁在线播放| 国产亚洲精品久久久久久毛片| 国产一区在线观看成人免费| 少妇裸体淫交视频免费看高清 | 热re99久久国产66热| 日本成人三级电影网站| 又黄又粗又硬又大视频| av天堂在线播放| 一a级毛片在线观看| 中亚洲国语对白在线视频| 午夜老司机福利片| 最新美女视频免费是黄的| 精品福利观看| 国产精品免费一区二区三区在线| 淫秽高清视频在线观看| 久9热在线精品视频| 正在播放国产对白刺激| 亚洲成人久久性| 亚洲精品色激情综合| 男人的好看免费观看在线视频 | 国产片内射在线| 国产亚洲精品一区二区www| 老司机靠b影院| 欧美黑人精品巨大| 1024视频免费在线观看| 亚洲精品国产一区二区精华液| 中文字幕高清在线视频| 亚洲自偷自拍图片 自拍| 一区福利在线观看| 男人的好看免费观看在线视频 | 欧美乱色亚洲激情| 无限看片的www在线观看| 亚洲av中文字字幕乱码综合 | 丁香欧美五月| 国产野战对白在线观看| 变态另类丝袜制服| 婷婷六月久久综合丁香| 久久久国产欧美日韩av| 欧美性长视频在线观看| 在线视频色国产色| 91国产中文字幕| 精品第一国产精品| 欧美国产精品va在线观看不卡| 黄色视频,在线免费观看| 国产av不卡久久| 麻豆成人午夜福利视频| 脱女人内裤的视频| 国产精品免费视频内射| 欧美不卡视频在线免费观看 | 91成人精品电影| 女性生殖器流出的白浆| 亚洲性夜色夜夜综合| 禁无遮挡网站| 精品午夜福利视频在线观看一区| xxx96com| 国产亚洲精品av在线| 久久精品91蜜桃| 国产男靠女视频免费网站| 日韩欧美三级三区| 哪里可以看免费的av片| 一个人观看的视频www高清免费观看 | 国内久久婷婷六月综合欲色啪| 日本三级黄在线观看| 国产精品一区二区精品视频观看| 亚洲第一av免费看| 亚洲av熟女| 99精品久久久久人妻精品| 一级毛片高清免费大全| 国产精品98久久久久久宅男小说| 日韩精品青青久久久久久| 99精品欧美一区二区三区四区| 99国产精品一区二区蜜桃av| 真人一进一出gif抽搐免费| 精品第一国产精品| 黑人巨大精品欧美一区二区mp4| 国产在线观看jvid| 香蕉丝袜av| 美女 人体艺术 gogo| 成人欧美大片| 国产成人欧美| 18禁黄网站禁片免费观看直播| 精品国产美女av久久久久小说| 免费人成视频x8x8入口观看| 夜夜夜夜夜久久久久| 两个人看的免费小视频| 好看av亚洲va欧美ⅴa在| 亚洲国产精品合色在线| 亚洲精品一卡2卡三卡4卡5卡| 动漫黄色视频在线观看| 在线免费观看的www视频| 女生性感内裤真人,穿戴方法视频| 国产午夜精品久久久久久| 国产真实乱freesex| 1024视频免费在线观看| 桃红色精品国产亚洲av| 亚洲人成电影免费在线| 女警被强在线播放| 亚洲人成伊人成综合网2020| 精品久久久久久久人妻蜜臀av| 好男人电影高清在线观看| 久久天堂一区二区三区四区| 免费人成视频x8x8入口观看| 99热这里只有精品一区 | 色在线成人网| 91成年电影在线观看| x7x7x7水蜜桃| 国产熟女xx| 国产不卡一卡二| 欧美日韩一级在线毛片| 欧美三级亚洲精品| 又大又爽又粗| 国产乱人伦免费视频| 嫩草影院精品99| 久久国产精品人妻蜜桃| 黄色女人牲交| 欧美中文日本在线观看视频| 亚洲av成人不卡在线观看播放网| 男人操女人黄网站| 在线观看www视频免费| 一夜夜www| 波多野结衣巨乳人妻| 男女做爰动态图高潮gif福利片| 桃色一区二区三区在线观看| 自线自在国产av| 久久久久九九精品影院| 国产精品99久久99久久久不卡| 午夜福利免费观看在线| 黄色毛片三级朝国网站| 亚洲熟妇中文字幕五十中出| 99热这里只有精品一区 | 欧美乱码精品一区二区三区| 亚洲精品美女久久av网站| 国内久久婷婷六月综合欲色啪| 无遮挡黄片免费观看| 操出白浆在线播放| 夜夜看夜夜爽夜夜摸| 最近最新中文字幕大全免费视频| 久久久久免费精品人妻一区二区 | √禁漫天堂资源中文www| 久久久久亚洲av毛片大全| 亚洲色图 男人天堂 中文字幕| 国产黄片美女视频| 成人午夜高清在线视频 | 国产精品久久久久久人妻精品电影| 成人18禁高潮啪啪吃奶动态图| 日韩精品免费视频一区二区三区| 精品日产1卡2卡| 亚洲精品一区av在线观看| 后天国语完整版免费观看| 亚洲专区国产一区二区| 亚洲精品中文字幕在线视频| 一本久久中文字幕| 麻豆成人av在线观看| 嫩草影院精品99| 99精品在免费线老司机午夜| 天堂动漫精品| xxx96com| 99热只有精品国产| 黑人欧美特级aaaaaa片| 1024香蕉在线观看| 精品久久久久久成人av| 精华霜和精华液先用哪个| 亚洲欧美激情综合另类| 日本五十路高清| 久久久国产成人免费| 亚洲熟女毛片儿| 日本 欧美在线| 天天一区二区日本电影三级| 听说在线观看完整版免费高清| 丁香六月欧美| 国产午夜精品久久久久久| 人人澡人人妻人| 两个人免费观看高清视频| 91麻豆av在线| 2021天堂中文幕一二区在线观 | 久热爱精品视频在线9| 国产又色又爽无遮挡免费看| 欧美精品啪啪一区二区三区| 亚洲第一av免费看| 精品熟女少妇八av免费久了| 婷婷亚洲欧美| 国产精品免费一区二区三区在线| 可以免费在线观看a视频的电影网站| 不卡一级毛片| 欧美性长视频在线观看| 国产黄a三级三级三级人| 色尼玛亚洲综合影院| 久久久国产欧美日韩av| 免费在线观看黄色视频的| 一级黄色大片毛片| www.精华液| 亚洲自拍偷在线| 亚洲成人久久爱视频| 一边摸一边做爽爽视频免费| 男女之事视频高清在线观看| 看片在线看免费视频| 久久精品国产亚洲av香蕉五月| 18美女黄网站色大片免费观看| 精品国内亚洲2022精品成人| 色在线成人网| 后天国语完整版免费观看| 亚洲自偷自拍图片 自拍| 国产精品久久久久久精品电影 | 国产精品1区2区在线观看.| 91成年电影在线观看| 中文在线观看免费www的网站 | 在线观看一区二区三区| 两性夫妻黄色片| 啦啦啦韩国在线观看视频| 国产三级在线视频| 日本 欧美在线| 久久久久久久久中文| 久久中文看片网| 成人国语在线视频| 久久久久久九九精品二区国产 | 欧美色欧美亚洲另类二区| 免费在线观看影片大全网站| 老鸭窝网址在线观看| 欧美乱码精品一区二区三区| 久久久水蜜桃国产精品网| 日本在线视频免费播放| 男人舔奶头视频| 十八禁人妻一区二区| 中文资源天堂在线| 亚洲精品粉嫩美女一区| 欧美成人免费av一区二区三区| 国产精品亚洲一级av第二区| a级毛片a级免费在线| 高清在线国产一区| 嫩草影视91久久| a在线观看视频网站| 超碰成人久久| 欧美性猛交黑人性爽| 人妻丰满熟妇av一区二区三区| 黄色女人牲交| 国产欧美日韩一区二区精品| 国产熟女xx| av有码第一页| 欧美最黄视频在线播放免费| 久久精品国产清高在天天线| 夜夜夜夜夜久久久久| 亚洲全国av大片| 欧美成狂野欧美在线观看| 白带黄色成豆腐渣| 免费看美女性在线毛片视频| 中文字幕久久专区| 亚洲激情在线av| 色综合亚洲欧美另类图片| 日韩有码中文字幕| 国产成人精品久久二区二区免费| 岛国视频午夜一区免费看| 久久这里只有精品19| 亚洲三区欧美一区| 嫩草影院精品99| 国产亚洲欧美98| 999久久久国产精品视频| 黄片大片在线免费观看| 一级黄色大片毛片| a级毛片在线看网站| 国产乱人伦免费视频| 91av网站免费观看| 搡老熟女国产l中国老女人| 国产精品九九99| 啦啦啦观看免费观看视频高清| 久久久久久久久久黄片| 亚洲精品国产一区二区精华液| 中文字幕人成人乱码亚洲影| bbb黄色大片| 国产免费男女视频| 欧美黑人精品巨大| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三| 99热这里只有精品一区 | 欧美性猛交╳xxx乱大交人| 亚洲在线自拍视频| 久久精品国产清高在天天线| 国产97色在线日韩免费| av免费在线观看网站| 久久人妻福利社区极品人妻图片| 精品国产乱子伦一区二区三区| 午夜免费观看网址| 高清在线国产一区| 99re在线观看精品视频| 制服诱惑二区| 久久久久国产精品人妻aⅴ院| 欧美黑人欧美精品刺激| 国产成人欧美在线观看| 久久久久久免费高清国产稀缺| 日日干狠狠操夜夜爽| 一级毛片高清免费大全| 18禁国产床啪视频网站| 精品一区二区三区四区五区乱码| 美女 人体艺术 gogo| 国产三级在线视频| 一区二区三区精品91| 99国产精品99久久久久| 欧美在线黄色| 欧美一区二区精品小视频在线| 免费高清在线观看日韩| 一本精品99久久精品77| 国产亚洲精品久久久久久毛片| 听说在线观看完整版免费高清| 久久午夜综合久久蜜桃| 精华霜和精华液先用哪个| 日韩高清综合在线| 亚洲专区国产一区二区| 久久欧美精品欧美久久欧美| 香蕉av资源在线| 成人国语在线视频| 身体一侧抽搐| 中亚洲国语对白在线视频| 国产成+人综合+亚洲专区| 亚洲av美国av| 给我免费播放毛片高清在线观看| 欧美激情久久久久久爽电影| 欧美黄色片欧美黄色片| 国产亚洲欧美98| 免费在线观看影片大全网站| 一个人观看的视频www高清免费观看 | 黄色a级毛片大全视频| 久久亚洲精品不卡| 久久久久久久精品吃奶| 国产单亲对白刺激| 久久久久亚洲av毛片大全| 日韩免费av在线播放| 欧美色视频一区免费| 亚洲三区欧美一区| 后天国语完整版免费观看| 在线观看www视频免费| 18美女黄网站色大片免费观看| 一边摸一边做爽爽视频免费| 两性夫妻黄色片| 日本免费一区二区三区高清不卡| 一进一出抽搐动态| 日韩精品青青久久久久久| 国产精华一区二区三区| 人人妻人人看人人澡| 亚洲激情在线av| 观看免费一级毛片| 亚洲av成人av| 99精品在免费线老司机午夜| 亚洲精品国产精品久久久不卡| 国产精品电影一区二区三区| 淫秽高清视频在线观看| 熟妇人妻久久中文字幕3abv| 18禁美女被吸乳视频| 欧美精品亚洲一区二区| 国产精品久久视频播放| 亚洲 欧美 日韩 在线 免费| 国产熟女午夜一区二区三区| 少妇粗大呻吟视频| 免费在线观看亚洲国产| 日本五十路高清| 精品一区二区三区四区五区乱码| 伊人久久大香线蕉亚洲五| 午夜老司机福利片| 精品一区二区三区视频在线观看免费| 搡老妇女老女人老熟妇| 18禁国产床啪视频网站| e午夜精品久久久久久久| 99久久久亚洲精品蜜臀av| 国产成人精品无人区| 亚洲中文日韩欧美视频| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品合色在线| 色精品久久人妻99蜜桃| 国产精品亚洲一级av第二区| 狠狠狠狠99中文字幕| 欧美乱码精品一区二区三区| 免费看美女性在线毛片视频| 99国产精品99久久久久| АⅤ资源中文在线天堂| 人人澡人人妻人| 99在线人妻在线中文字幕| 两个人看的免费小视频| 欧美精品啪啪一区二区三区| 久久人人精品亚洲av| 久久久国产成人免费| 国产一级毛片七仙女欲春2 | 99热这里只有精品一区 | 两个人看的免费小视频| 99国产极品粉嫩在线观看| 国产高清视频在线播放一区| 亚洲国产精品成人综合色| 国产1区2区3区精品| 久久中文看片网| 亚洲国产精品久久男人天堂| 国产爱豆传媒在线观看 | 黑丝袜美女国产一区| 大型av网站在线播放| 可以在线观看毛片的网站| 久久久久久久久久黄片| 亚洲欧洲精品一区二区精品久久久| 日本 av在线| 亚洲国产精品久久男人天堂| 亚洲男人天堂网一区| 看免费av毛片| 一级毛片高清免费大全| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 老熟妇仑乱视频hdxx| e午夜精品久久久久久久| 国产又色又爽无遮挡免费看| 欧美国产日韩亚洲一区| 久久精品影院6| 亚洲av片天天在线观看| 久久中文字幕一级| 日韩精品中文字幕看吧| 嫩草影院精品99| 国产亚洲精品一区二区www| 99国产综合亚洲精品| 真人做人爱边吃奶动态| 亚洲中文av在线| 波多野结衣巨乳人妻| 露出奶头的视频| 亚洲av片天天在线观看| 亚洲九九香蕉| 老司机在亚洲福利影院| 国产亚洲av嫩草精品影院| 黄色 视频免费看| 久久久精品国产亚洲av高清涩受| 久久久久久免费高清国产稀缺| 18禁黄网站禁片午夜丰满| 午夜福利18| 特大巨黑吊av在线直播 | 国产精品美女特级片免费视频播放器 | 成人18禁高潮啪啪吃奶动态图| 1024香蕉在线观看| 中亚洲国语对白在线视频| 一级毛片女人18水好多| 国产精品免费一区二区三区在线| 午夜久久久在线观看| 一a级毛片在线观看| 高潮久久久久久久久久久不卡| 一个人免费在线观看的高清视频| 99久久99久久久精品蜜桃| 少妇裸体淫交视频免费看高清 | 亚洲一区二区三区色噜噜| 国产精品精品国产色婷婷| 岛国在线观看网站| 一边摸一边抽搐一进一小说| 午夜两性在线视频| 一级作爱视频免费观看| 国产野战对白在线观看| 法律面前人人平等表现在哪些方面| 免费av毛片视频| 中文字幕av电影在线播放| 亚洲中文av在线| 久久久久国产精品人妻aⅴ院| 国产欧美日韩一区二区精品| 色在线成人网| 男女那种视频在线观看| 丝袜美腿诱惑在线| 手机成人av网站| 欧美日韩一级在线毛片| 久久欧美精品欧美久久欧美| 午夜成年电影在线免费观看| av有码第一页| 在线十欧美十亚洲十日本专区|