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

    Variation of net primary productivity and its drivers in China’s forests during 2000-2018

    2020-07-16 07:19:14YuheJiGuangshengZhouTianxiangLuoYakirDanLiZhouandXiaominLv
    Forest Ecosystems 2020年2期

    Yuhe Ji,Guangsheng Zhou*,Tianxiang Luo,Yakir Dan,Li Zhou and Xiaomin Lv

    Abstract

    Keywords: Net primary production (NPP), Forest ecosystem; annual precipitation, NPP model, Fluctuation, Variability

    Background

    Vegetation net primary productivity (NPP), as a biomass increment, plays an important role in the global carbon cycle. It contributes to understanding the contribution of NPP to the carbon exchanges that take place between the biosphere and atmosphere. Vegetation NPP indicates the atmospheric carbon fixed by plants except for the carbon released by respiration. Global vegetation NPP has increased due to climate change over the last few decades (Nemani et al. 2003; Chen et al. 2019). A significant increase in NPP was observed for grasslands and shrublands after the 1980s across China (Piao et al.2001; Liu et al. 2017; Liang et al. 2015; Xu and Wang 2016; Wang et al. 2017). However, some research has found declines in global NPP over the past decade because of drought (Zhao and Running 2010; Medlyn 2011).

    There are various process models for the estimation of vegetation NPP. The Carnegie-Ames-Stanford Approach(CASA) provides a useful tool to estimate NPP by combining grid meteorological data and the satellite-derived fraction of incoming photosynthetically active radiation absorbed by vegetation (Smith et al. 2008). The Carbon Exchange between Vegetation, Soil and Atmosphere(CEVSA) model was used to explore NPP in China’s terrestrial ecosystems,and confirmed a slow increase in NPP over the period 1981-1998 (Tao et al. 2003). The Boreal Ecosystem Productivity Simulator (BEPS) was used to map the spatio-temporal variations of NPP in China for the period from 2000 to 2010 (Chen et al. 2007; Liu et al.2013). The Lund-Potsdam-Jena Dynamic Global Vegetation Model (LPJ-DGVM) was widely applied to predict NPP dynamics by combining terrestrial vegetation dynamics (Smith et al. 2008; Gao et al. 2016).Based on the eddy covariance technique, an eco-physiological process model was developed to calculate NPP directly from gross primary productivity (GPP) by the mean ratio of the NPP/GPP of different vegetation types in a local area (Wang et al.2017).All the process models can be applied in other regions. However, traditional process models have many eco-physiological parameters, and the parameters need to be modified for different research areas. It is difficult to ensure the accuracy of the parameters. Additionally, there are uncertainties in the estimation of NPP because of accumulated errors during the complex process.

    Statistics models for vegetation NPP are relatively simple but they can provide a credible output for a specific research area, although they may be difficult to apply to other regions (Adams et al. 2004). The Miami model is a classical statistics-based model that has been used widely for estimating vegetation NPP in different regions (Lieth 1973; Alexandrov and Matsunaga 2008). In China, Zhou and Zhang (1996) established a general statistical model.The model is suitable for predicting NPP for terrestrial vegetation over China, but its accuracy is not very good for estimating forest NPP. Until now, there has not been a specific statistics-based model to predict forest NPP accurately.

    Both process models and statistics models need to identify the dominant factors affecting NPP, since different dominant factors might lead to different simulation results. For example, Zhao and Running (2010) reported that global NPP had declined over the past decade, but their findings were based on outcomes from models which assumed a strong temperature dependence (Medlyn 2011).Despite recognition of the roles of temperature and precipitation in controlling NPP dynamics, there are differing opinions over their relative influence(Nemani et al.2003;Zhang et al.2017).Some observations have revealed that NPP was strongly correlated with annual precipitation at a continental scale, but interannual variability in ANPP (above-ground NPP) was not related to variability in precipitation (Knapp and Smith 2001). This conclusion is controversial, since it was reported that the interannual variation in NPP could be explained by changes in precipitation in arid and semi-arid regions,and by changes in air temperature in other regions of China (Chen et al. 2013;Liang et al. 2015; Xu and Wang 2016). In addition,vegetation structure (e.g., canopy cover and leaf area index) can change due to natural growth, climate change and human disturbance(Chen et al.2019).NPP dynamics are likely to be affected by changes in vegetation structure,but most studies ignore the effect of vegetation structure dynamics on NPP fluctuations. The effect of environment variables in NPP dynamics will most likely vary in different regions.

    This paper established a statistics-based multiple regression model to estimate forest NPP across China.The objectives were to (1) explore the long-term trends and annual fluctuations of forest NPP for different forest types at a national scale during the period 2000-2018;and (2) identify the key environmental factors controlling the variability of forest NPP. This study helps to understand the variability of forest NPP, as well as the interactions between forest productivity and environmental factors.

    Materials and methods

    Building the statistics model

    The paper assumed that the spatiotemporal patterns of forest NPP were controlled mainly by climatic, topographic and forest structure variables, although previous research has shown different relative dominance between factors affecting NPP. To establish the statistics multivariate regression model, we collected long-term observations of annual forest NPP, normalized difference vegetation index (NDVI, which was used to quantify vegetation structure), altitude, average precipitation and temperature between 1980 and 2010 at 1000 forest inventory sites (Fig.1).

    Local topography significantly affects the spatial variation of climatic variables, and each elevation increase of 100 m in China leads to an average reduction of 25 gC·m-2for annual NPP (Chen et al. 2007). The topographic effect on NPP can be described by altitude, so altitude from a digital elevation model (DEM) was selected as an important variable in establishing the statistics-based multivariate regression model. The NDVI was extracted from the annual maximum NDVI of a 16-day cycle, 2000-2010 (https://lpdaac.usgs.gov).

    Forest NPP was calculated using the traditional Volume Growth Rate Model from the data of forest inventories, 1980-2010 (Eq.1) (Luo 1996).

    where Ps, Pb, Pl and Pr are the annual net biomass of stem, branch, leaf and root in the tree layer, respectively.Pu is the annual net biomass of shrubs and herbs under the forest canopy.

    Quality control of forest NPP had been done. Firstly,forest NPP for each observation site was the average value of continuous forest inventory since 1980, and the average value reduced data uncertainty. Secondly, forest NPP for each observation site was checked by comparing with the values of the surrounding area and the values of the same forest types. Thirdly, forest NPP was checked by comparing with the results of some forest productivity models (Luo 1996; Luo et al. 2004).

    Based on the relationship between forest NPP and environmental variables (annual maximum NDVI, annual precipitation, annual temperature, and altitude) at 1000 forest inventory sites, a statistics-based multiple regression model was established (Table 1), with a multiple regression R=0.72, and a standard error of 4.34 (Eq.2),

    where NPP is forest NPP (gC·m-2·year-1), NDVI is normalized difference vegetation index (0 <NDVI ≤1), P is annual precipitation (mm), T is annual temperature (°C),and A is altitude (m).

    Simulation of annual forest NPP over China

    The statistics-based multiple regression model allowed us to simulate annual forest NPP over China. The data included climatic variables (annual temperature and annual precipitation), altitude and NDVI. Annual temperature and annual precipitation were obtained from 800 meteorological stations of the National Meteorological Information Center of the China Meteorological Administration (http://cdc.cma.gov.cn). Altitude was derived from a 90-m resolution DEM dataset(http://srtm.csi.cgiar.org/srtmdata/). NDVI with a 1-km spatial resolution was obtained from Land ProcessesDistributed Active Archive Center (https://lpdaac.usgs.gov). To reduce the residual noise caused by haze and clouds, the annual maximum NDVI was extracted from a series of 16-day compositions of NDVI in the growing season. All the input variables were interpolated or resampled into raster layers with a spatial resolution of 1 km×1 km grid cells using the Kriging method with ArcGIS10.2 (http://www.esri.com/software/arcgis).

    Table 1 The parameters of the statistics-based multiple regression model between forest net primary productivity(NPP) and environmental variables (T is annual temperature (°C), P is annual precipitation (mm),A is altitude (m),NDVI is normalized difference vegetation index (0 <NDVI ≤1))

    The distribution of annual forest NPP was obtained from overlay analysis between the Vegetation Map of the People’s Republic of China in 2008 and the simulated annual forest NPP. Prior to overlay analysis, all the forest types were merged into five dominant forest types,including evergreen broadleaf forest, evergreen needleleaf forest, broadleaf-needleleaf mixed forest, deciduous broadleaf forest, deciduous needleleaf forest.

    Accuracy assessment

    The observed forest NPP at 90 observation sites from forest inventories were used to check the accuracy of the multiple regression model. The check sites were selected following the principle of representativeness. The 90 check sites covered the typical forest types in China, and their spatial distribution could not be concentrated in one area.

    Comparison between the predicted NPP and the observed NPP provided a fitting coefficient of R2=0.536(P <0.01) (Fig.2). The results indicated that the multiple regression model had a good performance in predicting forest NPP with acceptable accuracy. Therefore, the model was capable of providing a reliable estimate of China’s forest NPP.

    Results

    Trend and fluctuation of forest NPP

    The results revealed an average forest NPP of 840 gC·m-2·year-1for all forest types in China. Evergreen forests had the highest NPP (1058 gC·m-2·year-1for the evergreen broadleaf forest, 934 gC·m-2·year-1for the evergreen needleleaf forest), because evergreen forest maintains continuous photosynthesis throughout the year. Deciduous forests had lower NPP values (759 gC·m-2·year-1for the deciduous broadleaf forest, 590 gC·m-2·year-1for the deciduous needleleaf forest) since deciduous forests have a narrower photosynthesis period.

    Variations (long-term trend and interannual fluctuation) of forest NPP in China were examined. Results revealed an overall increasing trend in forest NPP, with an average increase rate of 5.2 gC·m-2·year-1over China between 2000 and 2018. All the five major forest ecosystems showed significant increasing forest NPP, with an increase rate of 2.9-8.0 gC·m-2·year-1. The largest increase rate was found in evergreen broadleaf forest (8.0 gC·m-2·year-1), followed by evergreen needleleaf forest(6.5 gC·m-2·year-1), broadleaf-needleleaf mixed forest(6.3 gC·m-2·year-1), deciduous broadleaf forest (5.4 gC·m-2·year-1) and deciduous needleleaf forest (2.9 gC·m-2·year-1) (Fig.3).

    Interannual fluctuation in forest NPP displayed a fluctuation range of 50-269 gC·m-2·year-1among the five major forest ecosystems. NPP in the evergreen broadleaf forest had the highest fluctuation range(269 gC·m-2·year-1), followed by the evergreen needleleaf forest (202 gC·m-2·year-1). NPP in the deciduous needleleaf forest had the lowest fluctuation range(50 gC·m-2·year-1).

    The fluctuation rate (ratio of fluctuation range to average value) was calculated to reflect the degree of fluctuation. The results showed that the evergreen broadleaf forest had the largest fluctuation rate (25%), followed by the evergreen needleleaf forest (22%) and the broadleafneedleleaf mixed forest (20%). The deciduous needleleaf forest had the smallest fluctuation rate (8%) (Fig.3 and Table 2). The results indicated that, among the major forest ecosystems in China, the most unstable forest NPP occurred in the evergreen broadleaf forest, and the most stable forest NPP occurred in the deciduous needleleaf forest.

    Comparing forest NPP with its fluctuation, our results revealed an interesting phenomenon: the greater the forest NPP, the greater the fluctuation of forest NPP for the five major forest ecosystems in China.For example, evergreen broadleaf forest had the largest NPP at 1058 gC·m-2·year-1, with the largest fluctuation range, 269 gC·m-2·year-1, as well as the largest fluctuation rate (25%) among the five major forest ecosystems. Deciduous needleleaf forest had the smallest NPP at 590 gC·m-2·year-1, with the smallest fluctuation amplitude, 50 gC·m-2·year-1, as well as the smallest fluctuation rate (8%) (Fig.3 and Table 2). It is likely that this phenomenon is related to the fluctuation of dominant factors affecting NPP, but it needs further confirmation.

    Spatial distribution of forest NPP trends

    Figure 4 shows the rate of change of forest NPP during the period 2000-2018 at 1 km×1 km spatial resolution.China’s forest NPP exhibited an increasing rate of 0-161 gC·m-2·year-1, which was statistically significant over 90% of forest area during the period 2000-2018. Forest NPP over China did not exhibit a consistent increasing trend everywhere. Approximately 10% of forest area in some regions of southwest and northern China showed a slight decreasing trend in forest NPP, with a range 0-20 gC·m-2·year-1(Fig.4). Generally, there was an overall increasing trend in forest NPP, despite a slightly decreasing trend in some forest areas in China.

    Table 2 The average annual net primary productivity(NPP) and its interannual fluctuation for five main forest types over China during the period 2000-2018

    Key drivers of the variation in forest NPP

    To identify the key drivers leading to the variation in forest NPP, the fluctuation values (deviation from the mean) of NPP and the environment variables were extracted. Correlation analysis showed a significant correlation between NPP and precipitation fluctuations, with a Pearson coefficient 0.977 (P <0.01), Kendall tau-b coefficient 0.908 (P <0.01), Spearman’s rho coefficient 0.973 (P <0.01). A moderate correlation was found between NPP and temperature fluctuations, with Pearson coefficient 0.484 (P <0.01) (Table 3). The results indicated that the variation in forest NPP was much more sensitive to annual precipitation fluctuation and annual precipitation fluctuation played a pivotal role incontrolling forest NPP changes compared with other environment variables.

    Table 3 The correlations between fluctuations(deviation from the mean)of net primary productivity(NPP)and the driving factors(annual precipitation, annual temperature and normalised difference vegetation index(NDVI))

    Comparing the interannual change curves of forest NPP,annual precipitation,annual temperature,and annual maximum NDVI,more similar curves were discovered between NPP and annual precipitation. The results provide another credible proof that forest NPP fluctuation was sensitive to changes in precipitation. In contrast, temperature and NDVI made a smaller contribution to forest NPP(Fig.5).

    We also examined the spatial patterns of the changes of forest NPP, precipitation and temperature. It was discovered that the increase in forest NPP was accompanied by an increase in precipitation and temperature for most forest areas in China. The decrease in forest NPP was accompanied by a decrease in precipitation for a few forest areas in southwest China (Figs. 4, 6 and 7).

    In brief,precipitation fluctuation was a key factor leading to variations in forest NPP,and temperature was a secondary factor, judging by the correlation coefficients, the change curves and the spatial patterns of change trends.

    Discussions

    There are different estimates of China’s forest NPP, ranging from 420.1-843 gC·m-2·year-1according to previous studies (Table 4). It was difficult to judge whether one was more credible than the others because of the use of different methods. However, all (previous and our) studies conclude that forest NPP has shown an obvious increasing trend over China in recent decades(Wang et al. 2008; Mao et al. 2010; Yu et al. 2014). The increase in forest NPP is consistent with other findings that the terrestrial vegetation NPP of the northern hemisphere has increased over the past several decades(Nemani et al. 2003). Our results presented an average increase rate of 5.2 gC·m-2·year-1for China’s forest NPP in the period 2000-2018, despite different increase rates for the five major forest ecosystems in China (Fig.3).

    Previous studies have ignored the impact of forest structure on NPP, even though forest spatial structure might affect forest NPP (Fotis et al. 2018). The general statistical model of Zhou and Zhang (1996) as a typical climate productivity model also ignored the impact of forest structure, and forest NPP from the model of Zhou and Zhang was usually lower than the observed forest NPP. Our model selected vegetation structure (quantified by NDVI) as an important factor, since spatiotemporal change in forest structure was directly related to forest NPP dynamics. Our results confirmed that forest structure had significant effect on forest NPP dynamics (with Spearman’s rho coefficient 0.490, P <0.05).Therefore, vegetation structure as an input variable in our model improved the accuracy of simulation.

    There are a series of factors that may affect NPP, but they do not all have a strong influence on NPP dynamics. The dynamics of the dominant factors affecting NPP are the real causes of changes to NPP. Therefore, we extracted the fluctuation values (deviation from the mean)of NPP and environment variables prior to analysing their relationship. The relationship between these fluctuation values provided more reliable information on the factors controlling NPP dynamics.

    Our results revealed that changes in precipitation played a key role in the variation in forest NPP over China, and temperature fluctuation was a secondary factor affecting forest NPP variation. The results disagreed with previous results that air temperature was the dominant climatic factor that controlled the interannual variability in NPP throughout China, except for arid and semi-arid regions (Liang et al. 2015). The disagreement might result from different analysis methods. The results of Liang et al. (2015) came from a simple correlation analysis between annual NPP and climate conditions(i.e., mean annual air temperature and annual cumulative precipitation) from 1982 to 2010. Our results came from the correlation analysis between the fluctuation values of annual forest NPP and environmental factors(annual precipitation, annual temperature, and annual maximal NDVI) from 2000 to 2018. Furthermore, our results from the correlation analysis were confirmed by the fluctuation curves and the spatial trends of forest NPP and environmental factors.

    However, there were some uncertainties in the simulation results of our model.Our model paid more attention to natural factors, but ignored human factors. In fact,vegetation NPP in China benefited from the ecological conservation and restoration efforts, since a series of ecological protection policies had been implemented, such as Three-North Shelter Forest Program, Grain for Green Program,and Natural Forest Protection Program(Lv et al.2015;Zhu and Zheng 2019).Although our model selected vegetation structure(quantified by NDVI)as an important factor,and vegetation structure could reflect the impact of human beings on vegetation NPP to some extent, it was not enough to reflect the overall human impact.

    Table 4 Forest net primary productivity(NPP) and its trends in China in recent decades according to previous studies

    Conclusions

    This paper established a statistics-based multiple regression model to estimate annual forest NPP across China.The model had a good performance in providing a reliable estimation in forest NPP. The trends and the fluctuations of forest NPP in the period 2000-2018 were examined.The results revealed a significant increasing forest NPP for all the five major forest ecosystems, with an average increase rate of 5.2 gC·m-2·year-1over China.There was an increasing NPP trend rate of 0-161 gC·m-2·year-1over 90% of the forest area. The interannual fluctuation range of forest NPP was 50-269 gC·m-2·year-1. The evergreen broadleaf forest had the largest fluctuation(269 gC·m-2·year-1), followed by the evergreen needleleaf forest, the broadleaf-needleleaf mixed forest,the deciduous broadleaf forest and the deciduous needleleaf forest. Our results revealed that the variability in forest NPP was caused mainly by variations in precipitation, then by temperature fluctuations. Therefore, precipitation played a key role in the variation in forest NPP.

    Abbreviations

    NPP: Net primary productivity; ANPP: Above-ground net primary productivity;GPP: Gross primary productivity; CASA: Carnegie-Ames-Stanford Approach;CEVSA: Vegetation, Soil and Atmosphere model; BEPS: Boreal Ecosystem Productivity Simulator; LPJ-DGVM: Lund-Potsdam-Jena Dynamic Global Vegetation Model; DEM: Digital elevation model; NDVI: Normalized difference vegetation index

    Acknowledgments

    We gratefully acknowledge Leonie Seabrook, PhD,from Liwen Bianji, Edanz Group China(www.liwenbianji.cn/ac),for editing the English text of a draft of this manuscript.

    Authors’contributions

    Guangsheng Zhou& Yuhe Ji conceived and designed the study. Yuhe Ji wrote the paper. Tianxiang Luo, Li Zhou and Xiaomin Lv revised and perfected the design of the experiments. Yakir Dan reviewed and edited the manuscript. All authors read and approved the manuscript.

    Funding

    This study was supported by the National Natural Science Fundation of China (No.41571175, 31661143028) and the special funds for basic research and operation from the Chinese Academy of Meteorological Science(2017Y003).

    Availability of data and materials

    The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

    Ethics approval and consent to participate

    The subject has no ethic risk.

    Consent for publication

    All the data and relevant processing methods, including the automated generation of figures of evaluation results,belong to the host institution,namely, Chinese Academy of Meteorological Science.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1State Key Laboratory of Severe Weather (LASW),Chinese Academy of Meteorological Science, Beijing 100081,China.2Key Laboratory of Alpine Ecology and Biodiversity, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100101, China.3Earth& Planetary Sciences,Weizmann Institute of Science, 76100 Rehovot, Israel.

    Received: 21 September 2019 Accepted: 10 March 2020

    久久久久九九精品影院| 国产人妻一区二区三区在| av女优亚洲男人天堂| 狂野欧美白嫩少妇大欣赏| 蜜桃亚洲精品一区二区三区| av黄色大香蕉| 午夜视频国产福利| 日韩,欧美,国产一区二区三区 | 国产亚洲5aaaaa淫片| 色噜噜av男人的天堂激情| 国产成人91sexporn| 亚洲av成人av| 中国国产av一级| 91久久精品国产一区二区成人| 美女内射精品一级片tv| 亚洲高清免费不卡视频| 麻豆成人午夜福利视频| 一级av片app| 国产成人a区在线观看| 亚洲成av人片在线播放无| 久久久久久久久久久免费av| 国产爱豆传媒在线观看| 久久精品国产自在天天线| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产成人久久av| av卡一久久| 国产人妻一区二区三区在| 一卡2卡三卡四卡精品乱码亚洲| 老师上课跳d突然被开到最大视频| 亚洲三级黄色毛片| 一个人看视频在线观看www免费| 插逼视频在线观看| 日日啪夜夜撸| 亚洲伊人久久精品综合 | 久久久久网色| 国产精品一二三区在线看| 女的被弄到高潮叫床怎么办| 亚洲在线自拍视频| 亚洲,欧美,日韩| 国产av码专区亚洲av| 国产av在哪里看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | eeuss影院久久| 99久久精品热视频| 中国国产av一级| 国产亚洲av嫩草精品影院| 狂野欧美白嫩少妇大欣赏| 国产亚洲一区二区精品| 日韩一区二区三区影片| 亚洲国产精品国产精品| 97超碰精品成人国产| 精品不卡国产一区二区三区| 午夜日本视频在线| 国产成人91sexporn| 深爱激情五月婷婷| 亚洲精品国产av成人精品| 麻豆av噜噜一区二区三区| 亚洲人成网站高清观看| 国产伦理片在线播放av一区| 欧美日韩一区二区视频在线观看视频在线 | 99热6这里只有精品| 国产精品蜜桃在线观看| 免费人成在线观看视频色| 三级经典国产精品| 你懂的网址亚洲精品在线观看 | 国产精品一二三区在线看| 久久久久久久亚洲中文字幕| 欧美激情国产日韩精品一区| 日本黄大片高清| 26uuu在线亚洲综合色| 男女边吃奶边做爰视频| 免费搜索国产男女视频| 99热6这里只有精品| 国产精品麻豆人妻色哟哟久久 | 波野结衣二区三区在线| 日韩一本色道免费dvd| 草草在线视频免费看| 听说在线观看完整版免费高清| 精品久久久久久久末码| 干丝袜人妻中文字幕| 久久精品熟女亚洲av麻豆精品 | 国产精品麻豆人妻色哟哟久久 | 国产精品综合久久久久久久免费| av免费观看日本| 久久久久久国产a免费观看| 精品久久久久久久末码| 久久草成人影院| 看片在线看免费视频| 人人妻人人澡人人爽人人夜夜 | 大话2 男鬼变身卡| 国内揄拍国产精品人妻在线| 中文字幕熟女人妻在线| 亚洲国产色片| 国产av一区在线观看免费| 高清视频免费观看一区二区 | 国产 一区精品| 亚洲欧洲日产国产| 色播亚洲综合网| a级毛色黄片| 亚洲三级黄色毛片| 人体艺术视频欧美日本| 哪个播放器可以免费观看大片| 久久精品影院6| 成人无遮挡网站| 国产精品国产三级国产专区5o | 亚洲欧美日韩无卡精品| 天堂√8在线中文| 国产成人a∨麻豆精品| 久久精品国产亚洲网站| 国产大屁股一区二区在线视频| 久久久久免费精品人妻一区二区| 插逼视频在线观看| 国内揄拍国产精品人妻在线| 亚洲av一区综合| 欧美xxxx性猛交bbbb| 亚洲在线观看片| 国产精品.久久久| www日本黄色视频网| 搡女人真爽免费视频火全软件| 狂野欧美激情性xxxx在线观看| videossex国产| 国产一级毛片在线| 国产成人午夜福利电影在线观看| 最近中文字幕2019免费版| 日本色播在线视频| 久久久精品欧美日韩精品| 国产精品久久久久久精品电影小说 | av国产久精品久网站免费入址| 少妇熟女欧美另类| 成人综合一区亚洲| 小说图片视频综合网站| 婷婷色综合大香蕉| 中国国产av一级| 国产成人aa在线观看| 久久久精品大字幕| 精品不卡国产一区二区三区| 久久婷婷人人爽人人干人人爱| 偷拍熟女少妇极品色| 免费在线观看成人毛片| 天堂网av新在线| 三级男女做爰猛烈吃奶摸视频| eeuss影院久久| 69人妻影院| 99热这里只有是精品在线观看| 午夜视频国产福利| 日本免费一区二区三区高清不卡| 日韩制服骚丝袜av| 国产免费视频播放在线视频 | 男女那种视频在线观看| 国产精品1区2区在线观看.| 免费播放大片免费观看视频在线观看 | 我要看日韩黄色一级片| 能在线免费观看的黄片| 啦啦啦观看免费观看视频高清| 三级经典国产精品| eeuss影院久久| 偷拍熟女少妇极品色| 99热这里只有是精品在线观看| 看片在线看免费视频| 国产精品一区二区性色av| 色吧在线观看| 午夜免费激情av| 天堂影院成人在线观看| 精品不卡国产一区二区三区| 哪个播放器可以免费观看大片| av在线天堂中文字幕| 少妇人妻精品综合一区二区| 国产精品野战在线观看| 日韩欧美 国产精品| 我的女老师完整版在线观看| 午夜福利在线观看吧| 直男gayav资源| 看十八女毛片水多多多| 国产一区二区在线观看日韩| 岛国在线免费视频观看| 特大巨黑吊av在线直播| 亚洲在久久综合| 村上凉子中文字幕在线| 国产一区二区在线av高清观看| 2021少妇久久久久久久久久久| videos熟女内射| 久久久久国产网址| av国产免费在线观看| 欧美zozozo另类| 中文字幕亚洲精品专区| 免费看光身美女| 日本色播在线视频| 两个人视频免费观看高清| 特大巨黑吊av在线直播| 亚洲色图av天堂| 色噜噜av男人的天堂激情| 亚洲aⅴ乱码一区二区在线播放| 久热久热在线精品观看| 搞女人的毛片| 汤姆久久久久久久影院中文字幕 | 看非洲黑人一级黄片| 日日干狠狠操夜夜爽| 晚上一个人看的免费电影| 少妇熟女aⅴ在线视频| 97人妻精品一区二区三区麻豆| 嫩草影院新地址| 91精品国产九色| www.av在线官网国产| 3wmmmm亚洲av在线观看| 久久韩国三级中文字幕| 精品人妻一区二区三区麻豆| 亚洲在久久综合| 国产真实伦视频高清在线观看| 舔av片在线| 久久精品国产99精品国产亚洲性色| 欧美日韩国产亚洲二区| 美女大奶头视频| 免费观看性生交大片5| 晚上一个人看的免费电影| 国产黄片美女视频| 亚洲综合色惰| videossex国产| 国产一区有黄有色的免费视频 | 国产午夜精品论理片| 国产亚洲5aaaaa淫片| 成人漫画全彩无遮挡| 在线免费观看不下载黄p国产| 小蜜桃在线观看免费完整版高清| 永久网站在线| 国产av不卡久久| 在线免费十八禁| 成人综合一区亚洲| 中文天堂在线官网| 97在线视频观看| 亚洲欧美一区二区三区国产| 久久这里有精品视频免费| 91久久精品电影网| 精品久久久久久久久av| 老司机影院成人| 久久精品熟女亚洲av麻豆精品 | 日韩成人伦理影院| 只有这里有精品99| 久久99热这里只频精品6学生 | 99久久中文字幕三级久久日本| 69av精品久久久久久| 欧美丝袜亚洲另类| 高清视频免费观看一区二区 | 美女大奶头视频| 久久6这里有精品| 亚洲人成网站在线播| 人体艺术视频欧美日本| 亚洲av福利一区| 久久人人爽人人片av| 亚洲在久久综合| 久久婷婷人人爽人人干人人爱| 男人的好看免费观看在线视频| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| or卡值多少钱| 如何舔出高潮| 亚洲av二区三区四区| 日韩三级伦理在线观看| 亚洲人与动物交配视频| 麻豆成人av视频| 久久这里只有精品中国| 亚洲欧美清纯卡通| 久久久成人免费电影| 欧美+日韩+精品| 国产 一区精品| 欧美一区二区国产精品久久精品| 国产精品1区2区在线观看.| 精品国产三级普通话版| 最后的刺客免费高清国语| 秋霞伦理黄片| 亚洲欧洲国产日韩| 99热全是精品| .国产精品久久| 丰满乱子伦码专区| 国产一区有黄有色的免费视频 | av女优亚洲男人天堂| 国产午夜精品论理片| 免费播放大片免费观看视频在线观看 | 国产精品久久电影中文字幕| 天堂中文最新版在线下载 | 亚洲国产高清在线一区二区三| 五月玫瑰六月丁香| 亚洲美女搞黄在线观看| 中文欧美无线码| 欧美zozozo另类| 午夜视频国产福利| 亚洲精品乱久久久久久| 国产一级毛片在线| 国产精品综合久久久久久久免费| 国产精品久久视频播放| 91狼人影院| 久久精品91蜜桃| 午夜福利在线观看免费完整高清在| 亚洲四区av| 久久久久久久久久久丰满| 少妇熟女aⅴ在线视频| .国产精品久久| 亚洲精品日韩在线中文字幕| 亚洲av免费在线观看| 国产在视频线精品| 免费黄网站久久成人精品| 亚洲av不卡在线观看| 国内精品一区二区在线观看| 国产真实伦视频高清在线观看| 黑人高潮一二区| 亚洲欧美清纯卡通| 高清在线视频一区二区三区 | 亚洲国产精品合色在线| 黄片无遮挡物在线观看| 我要搜黄色片| 久久精品国产99精品国产亚洲性色| 午夜免费男女啪啪视频观看| 国产精品1区2区在线观看.| 亚洲伊人久久精品综合 | 亚洲美女搞黄在线观看| 国语自产精品视频在线第100页| 五月伊人婷婷丁香| 搡老妇女老女人老熟妇| 成人av在线播放网站| 国产日韩欧美在线精品| 欧美不卡视频在线免费观看| 国产一区二区在线观看日韩| 午夜亚洲福利在线播放| 久久久久性生活片| 中文字幕av在线有码专区| 国产免费福利视频在线观看| 亚洲国产精品sss在线观看| 乱人视频在线观看| 免费黄色在线免费观看| 观看免费一级毛片| 色网站视频免费| 国产伦精品一区二区三区视频9| 在线观看美女被高潮喷水网站| 精品无人区乱码1区二区| 亚洲在线观看片| av福利片在线观看| 只有这里有精品99| 欧美性猛交╳xxx乱大交人| 午夜激情欧美在线| 欧美一级a爱片免费观看看| 免费观看在线日韩| 日日撸夜夜添| 99久久精品国产国产毛片| 亚洲在线观看片| 免费观看在线日韩| 日韩欧美三级三区| 午夜老司机福利剧场| 日本三级黄在线观看| 欧美不卡视频在线免费观看| 国产一区二区在线观看日韩| 亚洲av一区综合| 国产免费视频播放在线视频 | 91在线精品国自产拍蜜月| 一级av片app| 欧美潮喷喷水| 国产精品一区二区三区四区免费观看| 高清午夜精品一区二区三区| 丰满人妻一区二区三区视频av| 久久午夜福利片| 欧美3d第一页| 国产精品人妻久久久久久| 蜜臀久久99精品久久宅男| 婷婷六月久久综合丁香| 国产亚洲最大av| 亚洲欧美一区二区三区国产| 国产在线男女| 伦精品一区二区三区| 国产av在哪里看| 成人漫画全彩无遮挡| 成人二区视频| 国产亚洲精品av在线| 日韩一本色道免费dvd| www.色视频.com| 我的女老师完整版在线观看| 日韩强制内射视频| 国产精华一区二区三区| 久久精品国产99精品国产亚洲性色| 天堂影院成人在线观看| 国内精品一区二区在线观看| 久久久午夜欧美精品| www.色视频.com| 91在线精品国自产拍蜜月| 亚洲av日韩在线播放| 久久亚洲国产成人精品v| 亚洲精品自拍成人| 亚洲欧美日韩东京热| 波多野结衣巨乳人妻| 1000部很黄的大片| 国产成人福利小说| 午夜爱爱视频在线播放| 黄色欧美视频在线观看| 欧美日韩综合久久久久久| av在线观看视频网站免费| 亚洲国产欧洲综合997久久,| 天天躁日日操中文字幕| 热99re8久久精品国产| 国产精品麻豆人妻色哟哟久久 | 嫩草影院入口| 国产精品蜜桃在线观看| 国产在线男女| 亚洲精品aⅴ在线观看| 搡女人真爽免费视频火全软件| 深爱激情五月婷婷| 日本一二三区视频观看| 高清毛片免费看| 边亲边吃奶的免费视频| 精品免费久久久久久久清纯| 成人美女网站在线观看视频| 久久欧美精品欧美久久欧美| 国产精品日韩av在线免费观看| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看 | 人妻系列 视频| 人人妻人人看人人澡| 精品人妻偷拍中文字幕| 亚洲av电影在线观看一区二区三区 | av国产久精品久网站免费入址| 国产中年淑女户外野战色| 欧美日韩综合久久久久久| 国产午夜精品久久久久久一区二区三区| 床上黄色一级片| 国产在视频线在精品| 国产高潮美女av| 亚洲18禁久久av| 久热久热在线精品观看| 国产精品一区二区三区四区免费观看| www.色视频.com| 禁无遮挡网站| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 亚洲高清免费不卡视频| 大香蕉97超碰在线| 久99久视频精品免费| 非洲黑人性xxxx精品又粗又长| 麻豆精品久久久久久蜜桃| 国产白丝娇喘喷水9色精品| 亚洲最大成人av| 国产免费又黄又爽又色| 精品国产三级普通话版| 女的被弄到高潮叫床怎么办| 国产女主播在线喷水免费视频网站 | 午夜激情福利司机影院| 国产高清三级在线| 99久久九九国产精品国产免费| 国产三级在线视频| 最近中文字幕高清免费大全6| 一级av片app| 两性午夜刺激爽爽歪歪视频在线观看| 男女啪啪激烈高潮av片| 亚洲精品国产成人久久av| 国产视频内射| 精品熟女少妇av免费看| 99在线人妻在线中文字幕| 精品无人区乱码1区二区| 麻豆成人av视频| 日韩中字成人| 女人十人毛片免费观看3o分钟| 99热这里只有是精品在线观看| 乱码一卡2卡4卡精品| 欧美不卡视频在线免费观看| 国产三级中文精品| 99久久人妻综合| 汤姆久久久久久久影院中文字幕 | 亚洲经典国产精华液单| 观看美女的网站| 免费人成在线观看视频色| av又黄又爽大尺度在线免费看 | 国产爱豆传媒在线观看| 日本一二三区视频观看| 插逼视频在线观看| 亚洲高清免费不卡视频| 伦理电影大哥的女人| 校园人妻丝袜中文字幕| 国内精品一区二区在线观看| 免费看美女性在线毛片视频| 男女那种视频在线观看| 91av网一区二区| 嫩草影院新地址| 美女高潮的动态| 欧美日本亚洲视频在线播放| 51国产日韩欧美| 永久免费av网站大全| 18禁在线无遮挡免费观看视频| 看十八女毛片水多多多| 69av精品久久久久久| 不卡视频在线观看欧美| 青青草视频在线视频观看| 国产成人精品一,二区| 亚洲欧美成人综合另类久久久 | 国产av在哪里看| 国产美女午夜福利| ponron亚洲| 国产真实乱freesex| 国产爱豆传媒在线观看| 精品无人区乱码1区二区| 国内少妇人妻偷人精品xxx网站| 狠狠狠狠99中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 噜噜噜噜噜久久久久久91| 看黄色毛片网站| 国产亚洲最大av| 欧美+日韩+精品| 狠狠狠狠99中文字幕| 欧美激情久久久久久爽电影| 一夜夜www| 熟妇人妻久久中文字幕3abv| 深爱激情五月婷婷| 日本一本二区三区精品| 国产高清视频在线观看网站| 免费人成在线观看视频色| or卡值多少钱| 欧美+日韩+精品| 久久精品久久精品一区二区三区| 久久久欧美国产精品| 又粗又硬又长又爽又黄的视频| 国产色爽女视频免费观看| 日本wwww免费看| 中文字幕免费在线视频6| 欧美bdsm另类| 免费看日本二区| 又粗又爽又猛毛片免费看| 国产黄色视频一区二区在线观看 | 久久久久九九精品影院| 欧美不卡视频在线免费观看| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验| 我要看日韩黄色一级片| 日韩欧美在线乱码| 中国国产av一级| 欧美日韩国产亚洲二区| 一级毛片久久久久久久久女| 中国美白少妇内射xxxbb| 色综合站精品国产| 国产免费又黄又爽又色| 搡女人真爽免费视频火全软件| 秋霞在线观看毛片| 中文天堂在线官网| 成人性生交大片免费视频hd| av在线蜜桃| 在线免费观看不下载黄p国产| 国产乱人视频| 久久99热6这里只有精品| 国产精品福利在线免费观看| 中文字幕亚洲精品专区| 久久久久精品久久久久真实原创| 国产精品久久电影中文字幕| 热99在线观看视频| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久精品电影| 亚洲国产精品久久男人天堂| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 国产精品人妻久久久影院| 国产一区有黄有色的免费视频 | 国产精品国产高清国产av| 精品久久久久久久末码| 天堂中文最新版在线下载 | 欧美激情久久久久久爽电影| 一个人观看的视频www高清免费观看| 少妇高潮的动态图| 99久久无色码亚洲精品果冻| 亚洲成人av在线免费| 国语对白做爰xxxⅹ性视频网站| 精品熟女少妇av免费看| 能在线免费观看的黄片| 亚洲精品乱码久久久久久按摩| 国产麻豆成人av免费视频| 亚洲自拍偷在线| 午夜福利成人在线免费观看| 久久这里只有精品中国| 亚洲国产精品国产精品| 日本wwww免费看| 欧美97在线视频| 综合色av麻豆| 特级一级黄色大片| 六月丁香七月| 寂寞人妻少妇视频99o| 国产乱来视频区| 亚洲成人中文字幕在线播放| 内地一区二区视频在线| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲在久久综合| 天美传媒精品一区二区| 亚洲色图av天堂| 午夜福利视频1000在线观看| 1000部很黄的大片| 三级毛片av免费| 国产亚洲5aaaaa淫片| 国产亚洲一区二区精品| 男女下面进入的视频免费午夜| 高清毛片免费看| av又黄又爽大尺度在线免费看 | 日韩精品青青久久久久久| 国产精品伦人一区二区| 成人国产麻豆网| 日本色播在线视频| 国产91av在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色综合亚洲欧美另类图片| 在线播放无遮挡| 一本一本综合久久| 热99在线观看视频| 国产伦在线观看视频一区| 亚洲国产精品成人久久小说| 热99在线观看视频| 国产伦在线观看视频一区| 国产真实乱freesex| 美女国产视频在线观看| 久久久久免费精品人妻一区二区| 国产精品嫩草影院av在线观看| 久久亚洲精品不卡| 欧美日韩国产亚洲二区| 国产精品久久久久久精品电影小说 | 亚洲aⅴ乱码一区二区在线播放| 中文天堂在线官网| 精品人妻一区二区三区麻豆|