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

    Effects of reservoirs on seasonal discharge of Irtysh River measured by Lepage test

    2014-03-15 07:15:46FengHUANGZiqiangXIALidanGUOFuchengYANG
    Water Science and Engineering 2014年4期

    Feng HUANG*, Zi-qiang XIA, Li-dan GUO, Fu-cheng YANG

    1. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, P. R. China

    2. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, P. R. China

    3. Sichuan Electric Power Design and Consulting Co., Ltd., Chengdu 610016, P. R. China

    Effects of reservoirs on seasonal discharge of Irtysh River measured by Lepage test

    Feng HUANG*1,2, Zi-qiang XIA1,2, Li-dan GUO1,2, Fu-cheng YANG3

    1. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, P. R. China

    2. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, P. R. China

    3. Sichuan Electric Power Design and Consulting Co., Ltd., Chengdu 610016, P. R. China

    The Irtysh River is an international river partially joining the territories of China, Kazakhstan, and Russia. Cascade reservoirs have been constructed in the upper reaches of the river and their effects on the seasonal discharge of the middle and lower reaches were analyzed considering the mean and dispersion of the seasonal discharge. The Lepage test, which is a nonparametric, two-sample test for detecting location and dispersion, was used to measure the significance of difference between the pre-dam and post-dam seasonal discharge. The results show that the reservoirs’ effects on the seasonal discharge varied with the season. In the middle reaches of the river, the summer and autumn discharge decreased significantly and their inter-annual variabilities also decreased significantly. The summer and autumn precipitation over the Irtysh River Basin changed little before and after the operation of the reservoir, which indicates that the discharge changes mainly due to water storage of the reservoirs. The reservoirs store water in summer and autumn and store more water in a wet year, which leads to the reduction of the mean and dispersion of the summer and autumn discharge. The winter discharge increased significantly because the reservoirs released water for power generation. The spring discharge changed slightly. In the lower reaches, only the winter discharge increased significantly, and the other seasonal discharge changed slightly. The reservoirs’ effects on the seasonal discharge are more significant in the middle reaches than in the lower reaches.

    seasonal discharge; effects of reservoirs; Lepage test; Irtysh River

    1 Introduction

    One of the main hydrologic impacts of dams is the decrease in streamflow variability due to peak flow damping and minimum flow enhancement, which lead to flow homogenization downstream of dams. The homogenization also occurs with respect to the inter-annual variability of streamflow (Ye et al. 2003; Yang et al. 2004a; Lu and Siew 2006; Lajoie et al. 2007; Poff et al. 2007; Assani et al. 2011; Zhang et al. 2012), where the trend becomes almostuniform downstream of dams. From an ecological standpoint, more uniform streamflow and variability cause homogeneity of aquatic fauna, leading to significant loss of biodiversity and the loss of native species (Moyle and Mount 2007; Guo et al. 2011; Yu et al. 2014).

    Three reservoirs were constructed in the upper reaches of the Irtysh River in the 1950s and 1980s. Considering the influences of the reservoirs, the seasonal discharge series were divided into pre-dam and post-dam sub-series, and the significance of the difference between the two sub-series was measured by the Lepage test (Lepage 1971). The Lepage test has been demonstrated to be more statistically powerful than other similar tests, such as the Student’s t test, the chi-square test, and the Wilcoxon-Mann-Whitney test (Hirakawa 1974). This method has few underlying assumptions and has been successfully applied in many hydrological and meteorological studies (Yonetani and McCabe 1994; Chen et al. 2009; Yang et al. 2009; Zhang et al. 2009; Liu et al. 2011).

    The hydrological regime of the Irtysh River has changed dramatically over recent decades, probably due to the influence of human activities and climate changes (McClelland et al. 2004). Yang et al. (2004b) found that the streamflow decreases in summer and increases in winter in the upstream reaches of the Irtysh River Basin. Huang et al. (2012) found abrupt decreases in annual runoff and abrupt decreases in variability and concentration degree of intra-annual runoff in the middle reaches of the Irtysh River. Generally, much attention has been paid to the change trend and the intra-annual variability in previous studies. In this study, the seasonal discharge of the Irtysh River was analyzed, with a focus on its inter-annual variability. The results can add to our understanding of the impacts of multi-year regulating reservoirs on the Irtysh River and provide a vital basis for river management.

    2 Study area and data

    The Irtysh River (Fig. 1), the largest tributary of the Ob River, is 4 248 km long and drains an area of approximately 160 × 104km2. It originates in the Altai Mountains in the Sinkiang Uyghur Autonomous Region, in China, flows west through Zaysan Lake and northwest across Eastern Kazakhstan, and finally joins the Ob River in Khanty-Mansiysk City, in Russia. The range of the mainstream basin in Kazakhstan and China is 59.6 × 104km2, of which 5.6 × 104km2is the headwater area (Hrkal et al. 2006). The Irtysh River has two main tributaries: the Ishim River and the Tobol River, which are both on the left side. The Ishim River is 2 450 km long and drains an area of 14.4 × 104km2. The Tobol River is 1 591 km long and drains an area of 42.6 × 104km2.

    The Irtysh River Catchment has various natural conditions. The basin can be divided into two completely different parts in terms of morphology: the headwater area and a flat country area. The headwater area lies in the Altai Mountains, where local rocks are crystalline Paleozoic complex metamorphic rocks with some magmatic intrusions. The flat country area of the river catchment is part of the southern tip of the Siberian platform. Different types ofclimate can be observed in the headwater area. This region is an important water source for the whole catchment, and long-term annual total precipitation varies between 1 500 and 2 000 mm. The region has more evenly distributed precipitation than the Siberian platform. The flat country area, which is in the middle reaches of the Irtysh River Catchment, shows characteristics of a continental climate with extreme temperature values in the winter and summer. Generally, spring is from April to May, which is the snowmelt period, summer is from June to August, autumn is from September to October, and winter is from November to March of the next year.

    Large cascade reservoirs (Table 1) have been constructed in the upper reaches of the Irtysh River (Fig. 1). The Ust-Kamenogorsk Reservoir is the re-regulating reservoir of the Bukhtarma Reservoir, whose effective volume is only 1.7 × 108m3, and, thus, its influence on the discharge regime is relatively smaller compared with the Bukhtarma Reservoir. Therefore, the study period was segmented by the year 1960, when the Bukhtarma Reservoir was constructed. The Shul’binsk Reservoir was constructed in 1989. Thus, the post-dam discharge regime has mainly been influenced by the Bukhtarma Reservoir before 1989 and by cascade reservoirs after 1989.

    Table 1 Reservoirs in upper reaches of Irtysh River

    Fig. 1 Irtysh River Basin

    Two hydrometric stations with long-term observed data, one in a middle reach and the other in a lower reach of the Irtysh River, were selected (Fig. 1). The Omsk Hydrometric Station is located at 55.0°N and 73.4°E and drains an area of about 76.9 × 104km2. The multi-year mean runoff of the Omsk Hydrometric Station is about 277 × 108m3. The Tobolsk Hydrometric Station is located at 58.2°N and 68.2°E and drains an area of about 150 × 104km2. The multi-year mean runoff of the Tobolsk Hydrometric Station is about 676 × 108m3. Monthly discharge data of the two hydrometric stations from 1936 to 1999 were collected. To investigate the reservoirs’ effects on the seasonal discharge, the discharge series were divided into two sub-series: the pre-dam series (1936 to 1960) and the post-dam series (1961 to 1999).

    The Irtysh River has different water sources in different seasons. It is mainly supplemented by precipitation in summer and autumn. Twelve meteorological stations (Fig. 1) with long-term observed precipitation data were selected (Table 2). The monthly precipitation data from summer and autumn were collected and analyzed with the Lepage test.

    Table 2 Meteorological stations of Irtysh River Basin

    3 Methods

    The Lepage test is a nonparametric, two-sample test for detecting location and dispersion, and for measuring the significant difference between two samples, even if the distributions of parent populations are unknown. The Lepage test is a sum of the squares of the standardized Wilcoxon’s and Ansari-Bradley’s statistics:

    where W and A are Wilcoxon’s and Ansari-Bradley’s statistics, respectively, E is the mean value, and V is variance value. HK is the index of the Lepage test. It follows the chi-square distribution with two degrees of freedom. If HK exceeds 5.99 (4.21, 9.21), then the difference between the two samples corresponds to a significant confidence level of 95% (90%, 99%). The Lepage test is calculated as follows: Let x= (x1, x2,...xn1) and y= (y1, y2,...yn2) be two independent sam ples with sizes n1and n2, respectively. XY= (x1,x2,...xn1,y1,y2,...yn2) is the combined sample of x and y. Assume that ui=1 if the ith smallest data in XY belongs to x and ui=0 if it belongs to y.

    The statistics in Eq. (1) can be derived based on the following equations:

    If n1+n2is even,E(A) and V (A) are estimated as

    Otherwise,

    Statistical features of the segments divided by the change point are represented by the mean and the coefficient of variation (Cv). Furthermore, the empirical frequency curves (EFCs) of the segments are drawn to display the changes of the mean and the dispersion.

    The EFC is a simple and effective method of summarizing the distribution of discharge for a given catchment. The shape of the EFC is determined by the precipitation pattern, size, and physiographic characteristics of a catchment. The shape of the EFC can also be influenced by water resources development and land use. The EFC is widely used to measure the discharge regime as it provides an easy way of displaying the complete range of discharge. It can also be used to assess the changes in the discharge regime affected by human activities and climate changes, by considering discharge changes in percentiles. Each value of discharge (Q) has a corresponding exceedance probabilityp, and an EFC is simply a plot of Qp, which is the pth quantile or percentile of Q versus exceedance probability p, with p defined as

    The quantile Qpis a function of observed discharge, and since this function depends upon observations, it is often termed the empirical quantile function.

    4 Results and discussion

    4.1 Lepage test results of seasonal discharge change

    The seasonal discharge series were divided into two sub-series by the year 1960, when the Bukhtarma Reservoir began to operate. The Lepage test was used to measure whether the post-dam seasonal discharge changed significantly. The results are shown in Table 3.

    For the Omsk Hydrometric Station, the discharge changed more significantly in summerthan in other seasons. Both the mean and the standard deviation of the summer discharge series decreased after the operation of the reservoir. The mean of the summer discharge decreased from 1 730 m3/s to 1 041 m3/s. The autumn discharge also changed significantly and the pattern was similar with that of the summer discharge. However, the mean of the winter discharge increased by 31%, from 368 m3/s to 483 m3/s, after the operation of the reservoir. The spring discharge changed slightly after the operation of the reservoir.

    The Omsk Hydrometric Station is located in a middle reach of the Irtysh River, where the annual runoff is nearly equal to the effective volume of the Bukhtarma Reservoir, so its seasonal discharge was dramatically changed by the operation of the reservoir. The Tobolsk Hydrometric Station is located in a lower reach, where the water sources are more complex, so the seasonal discharge was not significantly influenced by the operation of the reservoir. The Lepage test shows that only the winter discharge significantly increased, by 25%, from 717 m3/s to 893 m3/s. The Irtysh River is mainly supplemented by the groundwater in the winter, so the water released from the reservoirs increases the winter discharge directly and effectively.

    Table 3 Statistical features of seasonal discharge and Lepage test results

    4.2 EFC analysis

    The EFCs of the pre-dam and post-dam seasonal discharge are shown in Fig. 2 and Fig. 3, which describe the changes of the mean and the dispersion of the seasonal discharge. For the pre-dam and post-dam sub-series, the EFCs of the seasonal discharge that changed slightly are nearly coincident. For the Omsk Hydrometric Station, the post-dam EFC of the summer discharge is remarkably lower than the pre-dam EFC, indicating a dramatic decrease in the mean of the summer discharge. Furthermore, the dispersion of the summer discharge also decreased. The summer discharge ranged from 698 m3/s to 2 923 m3/s during the period from 1936 to 1960 and from 586 m3/s to 1 676 m3/s during the period from 1961 to 1999. The dispersion of the autumn discharge also obviously decreased. The autumn discharge ranged from 375 m3/s to 1 760 m3/s during the period from 1936 to 1960 and from 472 m3/s to 910 m3/sduring the period from 1961 to 1999. For the winter discharge at the Omsk and Tobolsk hydrometric stations, the post-dam EFCs are high above the pre-dam EFCs, indicating a significant increase in the mean of the winter discharge.

    Fig. 2 EFC changes of seasonal discharge at Omsk Hydrometric Station

    Fig. 3 EFC changes of seasonal discharge at Tobolsk Hydrometric Station

    For the purposes of flood control, irrigation, power generation, and some other reasons, the reservoirs store water in summer and autumn and release water in winter, which reduces the intra-annual variability of discharge. Furthermore, the reservoirs store more water in wet years to maximize the power generation and reduce the inter-annual variabilities of the summer and autumn discharge.

    Feng HUANG et al. Water Science and Engineering, Oct. 2014, Vol. 7, No. 4, 363-372 369

    4.3 Lepage test results of summer and autumn precipitation changes

    Precipitation is the main water sources of the Irtysh River in summer and autumn. If the summer and autumn precipitation changed little before and after the operation of the reservoir, this can indicate that the discharge changed mainly due to the operation of the reservoir. The summer and autumn precipitation was analyzed with the Lepage test. The results are shown in Table 4. At the meteorological station 36665, the autumn precipitation increased significantly. In general, the summer and autumn precipitation over the Irtysh River Basin changed only slightly and would not have led to the seasonal discharge reduction. Therefore, the reduction in the summer and autumn discharge may be mainly attributed to the operation of the reservoir. It should be noted that some other reasons, e.g., increased agricultural and industrial water use, may also have contributed to the discharge reduction. However, their influences are much smaller compared with those of large reservoirs.

    Table 4 Statistical features of summer and autumn precipitation and Lepage test results

    5 Conclusions

    The Lepage test was used to evaluate the effects of reservoirs on the seasonal discharge of the Irtysh River. The summer and autumn precipitation was analyzed to separate the reservoirs’ effects from the natural influences. The main conclusions are as follows:

    (1) For the Omsk Hydrometric Station, located in a middle reach, the summer and autumn discharge decreased significantly and their inter-annual variabilities also decreased significantly. The winter discharge significantly increased and the spring discharge changed slightly.

    (2) For the Tobolsk Hydrometric Station located in a lower reach, only the winter discharge increased significantly due to the operation of the reservoir, and the discharge in other seasons changed slightly.

    (3) The summer and autumn precipitation changed little after the beginning of operation of the reservoir. This phenomenon indicates that the discharge reduction is mainly attributable to the operation of the reservoir.

    (4) The effects of the reservoirs on the seasonal discharge varied with the season. The reservoirs’ effects are greater on the seasonal discharge in the middle reaches than in the lower reaches. The water storage of the reservoirs significantly reduced the mean and dispersion of the summer and autumn discharge in the middle reaches, but only slightly affected in the lower reaches. The reservoirs’ effects on the winter and spring discharge of the middle reaches are coincident with those of the lower reaches. The reservoirs increase the winter discharge but only slightly change the spring discharge.

    Assani, A. A., Landry, R., Daigle, J., and Chalifour, A. 2011. Reservoirs effects on the interannual variability of winter and spring streamflow in the St-Maurice River Watershed (Quebec, Canada). Water Resources Management, 25(14), 3661-3675. [doi:10.1007/s11269-011-9875-1]

    Chen, Y. D., Zhang, Q., Xu, C. Y., Yang, T., Chen, X. H., and Jiang, T. 2009. Change-point alterations of extreme wate levels and underlying causes in the Pearl River delta, China. River Research and Applications, 25(9), 1153-1168. [doi:10.1002/rra.1212]

    Guo, W. X., Wang, H. X., Xu, J. X., and Xia, Z. Q. 2011. Ecological operation for Three Gorges Reservoir. Water Science and Engineering, 4(2), 143-156. [doi:10.3882/j.issn.1674-2370.2011.02.003]

    Hirakawa, K. 1974. The comparison of powers of distribution-free two-sample tests. TRU Mathematics, 10(1974), 65-82.

    Hrkal, Z., Gadalia, A., and Rigaudiere, P. 2006. Will the river Irtysh survive the year 2030? Impact of long-term unsuitable land use and water management of the upper stretch of the river catchment (North Kazakhstan). Environmental Geology, 50(5), 717-723. [doi:10.1007/s00254-006-0244-y]

    Huang, F., Xia, Z. Q., Li, F., Guo, L. D., and Yang, F. C. 2012. Hydrological Changes of the Irtysh River and the Possible Causes. Water Resources Management, 26(11), 3195-3208. [doi:10.1007/s11269-012-0067-4]

    Lajoie, F., Assani, A. A., Roy, A. G., and Mesfioui, M. 2007. Impacts of dams on monthly flow characteristics. The influence of watershed size and seasons. Journal of Hydrology, 334(3-4), 423-439. [doi:10.1016/j.jhydrol.2006.10.019]

    Lepage, Y. 1971. A combination of Wilcoxon’s and Ansari-Bradley’s statistics. Biometrika, 58(1), 213-217. [doi:10.1093/biomet/58.1.213]

    Liu, Y., Huang, G., and Huang, R. H. 2011. Inter-decadal variability of summer rainfall in Eastern Chinadetected by the Lepage test. Theoretical and Applied Climatology, 106(3-4), 481-488. [doi:10.1007/s00704-011-0442-8]

    Lu, X. X., and Siew, R. Y. 2006. Water discharge and sediment flux changes over the past decades in the Lower Mekong River: Possible impacts of the Chinese dams. Hydrology and Earth System Sciences, 10(2), 181-195. [doi:10.5194/hess-10-181-2006]

    McClelland, J. W., Holmes, R. M., Peterson, B. J., and Stieglitz, M. 2004. Increasing river discharge in the Eurasian Arctic: Consideration of dams, permafrost thaw, and fires as potential agents of change. Journal of Geophysical Research, 109(D18102), 1-12. [doi:10.1029/2004JD004583]

    Moyle, P. B., and Mount, J. F. 2007. Homogenous rivers, homogenous faunas. Proceedings of the National Academy of Sciences of the United States of America, 104(14), 5711-5712. [doi:10.1073/pnas.070 1457104]

    Poff, N. L., Olden, J. D., Merritt, D. M., and Pepin, D. M. 2007. Homogenization of regional river dynamics by dams and global biodiversity implications. Proceedings of the National Academy of Sciences of the United States of America, 104(14), 5732-5737. [doi:10.1073/pnas.0609812104]

    Yang, D. Q., Ye, B. S., and Kane, D. L. 2004a. Streamflow changes over Siberian Yenisei River Basin. Journal of Hydrology, 296(1-4), 59-80. [doi:10.1016/j.jhydrol.2004.03.017]

    Yang, D. Q., Ye, B. S., and Shiklomanov, A. 2004b. Discharge characteristics and changes over the Ob River watershed in Siberia. Journal of Hydrometeorology, 5(4), 595-610. [doi:10.1175/1525-7541(2004) 005〈0595:dcacot〉2.0.co;2]

    Yang, T., Chen, X., Xu, C. Y., and Zhang, Z. C. 2009. Spatio-temporal changes of hydrological processes and underlying driving forces in Guizhou region, Southwest China. Stochastic Environmental Research and Risk Assessment, 23(8), 1071-1087. [doi:10.1007/s00477-008-0278-7]

    Ye, B. S., Yang, D. Q., and L.Kane, D. 2003. Changes in Lena River streamflow hydrology: Human impacts versus natural variations. Water Resources Research, 39(7), SWC 8-1-14. [doi:10.1029/2003wr001991]

    Yonetani, T., and McCabe, G. J. 1994. Abrupt changes in regional temperature in the conterminous United States, 1895-1989. Climate Research, 4(1), 13-23.

    Yu, Z. B., Yang, T., and Schwartz, F. W. 2014. Water issues and prospects for hydrological science in China. Water Science and Engineering, 7(1), 1-4. [doi:10.3882/j.issn.1674-2370.2014.01.001]

    Zhang, Q., Xu, C. Y., and Yang, T. 2009. Variability of water resource in the Yellow River Basin of past 50 years, China. Water Resources Management, 23(6), 1157-1170. [doi:10.1007/s11269-008-9320-2]

    Zhang, W., Mu, S. S., Zhang, Y. J., and Chen, K. M. 2012. Seasonal and interannual variations of flow discharge from Pearl River into sea. Water Science and Engineering, 5(4), 399-409. [doi:10.3882/j.issn.1674-2370.2012.04.004]

    (Edited by Fang WANG)

    This work was supported by the Fundamental Research Funds for the Central Universities (Grant No. 2013/B13020312), the Ministry of Water Resources’ Special Funds for Scientific Research on Public Causes, for the People’s Republic of China (Grant No. 201001052), and the Innovative Project of Scientific Research for Postgraduates in Ordinary Universities of Jiangsu Province (Grant No. CXZZ11_0433).

    *Corresponding author (e-mail: hfeng0216@163.com)

    Received Sep. 17, 2013; accepted Mar. 6, 2014

    9191精品国产免费久久| 亚洲中文av在线| 国产成人精品在线电影| 久久精品人人爽人人爽视色| 国产亚洲精品久久久久5区| 国产精华一区二区三区| 国产国语露脸激情在线看| 老司机午夜十八禁免费视频| 性欧美人与动物交配| 脱女人内裤的视频| 老司机靠b影院| 91麻豆精品激情在线观看国产| 变态另类丝袜制服| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜久久久在线观看| 看黄色毛片网站| 免费在线观看影片大全网站| 在线天堂中文资源库| 天天添夜夜摸| av网站免费在线观看视频| 亚洲人成电影观看| 欧美久久黑人一区二区| 亚洲国产高清在线一区二区三 | 国产乱人伦免费视频| 女人高潮潮喷娇喘18禁视频| 黑丝袜美女国产一区| 青草久久国产| 久久热在线av| 国产99白浆流出| 91字幕亚洲| 亚洲久久久国产精品| 最近最新免费中文字幕在线| 制服人妻中文乱码| 欧美日韩黄片免| 1024香蕉在线观看| 久久久国产成人免费| 最新美女视频免费是黄的| 午夜激情av网站| 99re在线观看精品视频| 久久久久久久久久久久大奶| 午夜久久久久精精品| 国产乱人伦免费视频| 在线观看免费视频网站a站| 狠狠狠狠99中文字幕| 国产单亲对白刺激| 一边摸一边抽搐一进一出视频| 一进一出抽搐gif免费好疼| 欧美老熟妇乱子伦牲交| 久久性视频一级片| 91老司机精品| 好男人电影高清在线观看| 国产精品二区激情视频| 波多野结衣一区麻豆| 久久午夜综合久久蜜桃| 成人特级黄色片久久久久久久| 日韩精品免费视频一区二区三区| 色播亚洲综合网| 国产日韩一区二区三区精品不卡| 亚洲国产精品合色在线| 国产免费av片在线观看野外av| 欧美激情高清一区二区三区| 亚洲av电影不卡..在线观看| 美女国产高潮福利片在线看| 国产欧美日韩一区二区三| 国产精品久久久久久人妻精品电影| 午夜老司机福利片| 亚洲专区字幕在线| 久久精品aⅴ一区二区三区四区| 狠狠狠狠99中文字幕| 美女 人体艺术 gogo| 真人做人爱边吃奶动态| 无遮挡黄片免费观看| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 两个人看的免费小视频| 老汉色∧v一级毛片| 99国产精品免费福利视频| 久久九九热精品免费| 人人澡人人妻人| 欧美日韩黄片免| 在线观看一区二区三区| 色综合亚洲欧美另类图片| 涩涩av久久男人的天堂| 一a级毛片在线观看| 成在线人永久免费视频| 中文字幕色久视频| 久热爱精品视频在线9| 黄色女人牲交| 18禁裸乳无遮挡免费网站照片 | 伦理电影免费视频| 国产精品国产高清国产av| 男人舔女人下体高潮全视频| 十八禁人妻一区二区| 国产aⅴ精品一区二区三区波| 亚洲午夜精品一区,二区,三区| 精品久久久久久,| 黄频高清免费视频| 91麻豆av在线| 精品福利观看| 美女高潮喷水抽搐中文字幕| 99香蕉大伊视频| 人妻丰满熟妇av一区二区三区| 在线观看免费视频日本深夜| 美女高潮喷水抽搐中文字幕| 精品一区二区三区四区五区乱码| 好看av亚洲va欧美ⅴa在| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人一区二区免费高清观看 | 午夜免费观看网址| 波多野结衣一区麻豆| 亚洲欧美日韩高清在线视频| 亚洲国产欧美日韩在线播放| 久久精品人人爽人人爽视色| 人妻久久中文字幕网| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 一边摸一边抽搐一进一小说| 麻豆一二三区av精品| 神马国产精品三级电影在线观看 | 日韩欧美一区二区三区在线观看| 这个男人来自地球电影免费观看| 精品国内亚洲2022精品成人| 亚洲欧美日韩无卡精品| 午夜免费激情av| 中文字幕久久专区| 国产男靠女视频免费网站| 国产伦人伦偷精品视频| 久久久久久久久免费视频了| 人人妻人人澡欧美一区二区 | 国产乱人伦免费视频| 亚洲,欧美精品.| 成人手机av| 日日摸夜夜添夜夜添小说| 国产精品影院久久| 午夜福利在线观看吧| 国产在线精品亚洲第一网站| 欧美+亚洲+日韩+国产| 一夜夜www| www.熟女人妻精品国产| 韩国精品一区二区三区| 中亚洲国语对白在线视频| 中文字幕最新亚洲高清| 神马国产精品三级电影在线观看 | 久久精品影院6| 日韩有码中文字幕| 老司机深夜福利视频在线观看| 久久人人精品亚洲av| a级毛片在线看网站| 久久影院123| 纯流量卡能插随身wifi吗| 天天添夜夜摸| 黑人欧美特级aaaaaa片| 久久久国产欧美日韩av| 精品欧美国产一区二区三| 在线国产一区二区在线| 丰满的人妻完整版| 日本撒尿小便嘘嘘汇集6| 欧美精品啪啪一区二区三区| 19禁男女啪啪无遮挡网站| 两人在一起打扑克的视频| 超碰成人久久| 国产精品二区激情视频| 一进一出抽搐动态| 国产精品野战在线观看| 日韩 欧美 亚洲 中文字幕| 黑丝袜美女国产一区| 很黄的视频免费| 欧美激情久久久久久爽电影 | 亚洲精品av麻豆狂野| 嫩草影视91久久| 狠狠狠狠99中文字幕| 精品第一国产精品| 国产又爽黄色视频| 在线观看日韩欧美| 一二三四在线观看免费中文在| 国产蜜桃级精品一区二区三区| 在线观看午夜福利视频| 国产亚洲精品av在线| 国产一区在线观看成人免费| 欧美乱妇无乱码| 91在线观看av| 亚洲成av片中文字幕在线观看| 国产欧美日韩一区二区精品| www日本在线高清视频| 1024香蕉在线观看| 欧美黄色淫秽网站| 少妇粗大呻吟视频| 一级片免费观看大全| 一级毛片女人18水好多| 在线观看免费视频网站a站| 久9热在线精品视频| 亚洲男人天堂网一区| 亚洲一卡2卡3卡4卡5卡精品中文| 我的亚洲天堂| 国产亚洲精品第一综合不卡| 桃红色精品国产亚洲av| 不卡av一区二区三区| 欧美老熟妇乱子伦牲交| 一卡2卡三卡四卡精品乱码亚洲| 国产精品1区2区在线观看.| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美激情综合另类| 亚洲电影在线观看av| 久久久精品欧美日韩精品| 亚洲国产中文字幕在线视频| 日本五十路高清| 母亲3免费完整高清在线观看| 校园春色视频在线观看| 欧美+亚洲+日韩+国产| 久久久久久久精品吃奶| 亚洲欧美日韩高清在线视频| 国产精品香港三级国产av潘金莲| 一级,二级,三级黄色视频| 国产精品电影一区二区三区| 男人舔女人下体高潮全视频| 桃色一区二区三区在线观看| 成人特级黄色片久久久久久久| 久久香蕉国产精品| 成人三级做爰电影| 欧美精品啪啪一区二区三区| 天天躁夜夜躁狠狠躁躁| 91成人精品电影| 嫩草影院精品99| 欧美中文综合在线视频| 91字幕亚洲| 国产精品 国内视频| 国产1区2区3区精品| 一区福利在线观看| 久久久久久亚洲精品国产蜜桃av| 国产精品国产高清国产av| 一二三四社区在线视频社区8| 国产麻豆成人av免费视频| 757午夜福利合集在线观看| 亚洲,欧美精品.| 久久久国产成人精品二区| 久久香蕉精品热| 国产亚洲精品久久久久5区| 欧美一级a爱片免费观看看 | 黄网站色视频无遮挡免费观看| 亚洲精品粉嫩美女一区| 两人在一起打扑克的视频| 一边摸一边做爽爽视频免费| 国产高清有码在线观看视频 | 亚洲人成电影观看| 岛国在线观看网站| 丝袜美腿诱惑在线| 在线国产一区二区在线| 日韩三级视频一区二区三区| 两个人看的免费小视频| 久热这里只有精品99| 日韩一卡2卡3卡4卡2021年| 黄网站色视频无遮挡免费观看| 久久久久久大精品| 91老司机精品| 久久精品国产清高在天天线| 校园春色视频在线观看| 亚洲国产精品999在线| 在线观看日韩欧美| 999久久久精品免费观看国产| 午夜福利免费观看在线| 50天的宝宝边吃奶边哭怎么回事| 美女午夜性视频免费| 免费在线观看黄色视频的| 免费无遮挡裸体视频| 午夜福利免费观看在线| 丁香六月欧美| www.999成人在线观看| 亚洲avbb在线观看| 一级a爱片免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 亚洲专区字幕在线| 淫秽高清视频在线观看| 午夜免费鲁丝| 精品一品国产午夜福利视频| 女人精品久久久久毛片| 亚洲精品在线观看二区| 在线观看免费日韩欧美大片| 变态另类成人亚洲欧美熟女 | 欧美性长视频在线观看| 国产一区二区激情短视频| 国产精品 欧美亚洲| 欧美日本亚洲视频在线播放| 日本vs欧美在线观看视频| 色哟哟哟哟哟哟| 一级,二级,三级黄色视频| 中文字幕精品免费在线观看视频| 女性生殖器流出的白浆| 50天的宝宝边吃奶边哭怎么回事| av视频免费观看在线观看| 很黄的视频免费| 国产成人欧美| 国产色视频综合| 成年版毛片免费区| 满18在线观看网站| 久久精品亚洲熟妇少妇任你| 性色av乱码一区二区三区2| 一级毛片精品| 长腿黑丝高跟| 欧美一级a爱片免费观看看 | 97人妻精品一区二区三区麻豆 | 欧美色欧美亚洲另类二区 | 国产欧美日韩综合在线一区二区| 侵犯人妻中文字幕一二三四区| 免费高清在线观看日韩| 日日摸夜夜添夜夜添小说| 黄片小视频在线播放| 人妻久久中文字幕网| 神马国产精品三级电影在线观看 | 国产亚洲精品久久久久5区| 免费观看精品视频网站| 在线永久观看黄色视频| 久久伊人香网站| 国产精品久久久av美女十八| 女人被狂操c到高潮| 男人舔女人的私密视频| 99久久国产精品久久久| ponron亚洲| 久久九九热精品免费| 国语自产精品视频在线第100页| 国产亚洲欧美在线一区二区| 精品国内亚洲2022精品成人| 亚洲精品国产精品久久久不卡| 欧洲精品卡2卡3卡4卡5卡区| 婷婷六月久久综合丁香| 精品久久久久久成人av| 国产真人三级小视频在线观看| 久久香蕉激情| 亚洲人成网站在线播放欧美日韩| 搡老妇女老女人老熟妇| 亚洲国产日韩欧美精品在线观看 | 亚洲成av人片免费观看| 日韩中文字幕欧美一区二区| 99riav亚洲国产免费| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 美女扒开内裤让男人捅视频| 欧美黑人精品巨大| 色尼玛亚洲综合影院| 一a级毛片在线观看| 国产区一区二久久| 人人妻人人澡欧美一区二区 | 最好的美女福利视频网| 精品一品国产午夜福利视频| 国产亚洲欧美98| 亚洲 国产 在线| 99香蕉大伊视频| 久久婷婷成人综合色麻豆| cao死你这个sao货| 色婷婷久久久亚洲欧美| 日日干狠狠操夜夜爽| av电影中文网址| 免费高清在线观看日韩| 中文字幕av电影在线播放| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 国产精品影院久久| 热99re8久久精品国产| 悠悠久久av| tocl精华| 免费在线观看影片大全网站| 国产极品粉嫩免费观看在线| 自线自在国产av| 黄色片一级片一级黄色片| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆成人av免费视频| 老熟妇仑乱视频hdxx| 国产av又大| 动漫黄色视频在线观看| 成人手机av| 久久香蕉精品热| 亚洲国产欧美一区二区综合| 午夜两性在线视频| 好男人在线观看高清免费视频 | 日韩欧美一区二区三区在线观看| 亚洲男人的天堂狠狠| 一区二区三区高清视频在线| 国产av精品麻豆| 在线播放国产精品三级| 日韩中文字幕欧美一区二区| or卡值多少钱| 日日夜夜操网爽| 亚洲av成人av| 久久国产亚洲av麻豆专区| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美在线二视频| 美女午夜性视频免费| 校园春色视频在线观看| 欧美午夜高清在线| 50天的宝宝边吃奶边哭怎么回事| 天堂动漫精品| 中文字幕高清在线视频| 久久中文字幕人妻熟女| 免费在线观看视频国产中文字幕亚洲| 国产精品精品国产色婷婷| 美国免费a级毛片| 少妇裸体淫交视频免费看高清 | 国产精品亚洲av一区麻豆| 午夜福利18| 老司机午夜福利在线观看视频| 免费看a级黄色片| 非洲黑人性xxxx精品又粗又长| 午夜精品国产一区二区电影| 日本vs欧美在线观看视频| 村上凉子中文字幕在线| www.精华液| 婷婷六月久久综合丁香| 精品久久久久久久久久免费视频| 88av欧美| 99国产精品99久久久久| 99在线人妻在线中文字幕| 日本在线视频免费播放| 精品国产美女av久久久久小说| 精品久久蜜臀av无| 欧美+亚洲+日韩+国产| 男人的好看免费观看在线视频 | 手机成人av网站| 日本 欧美在线| 中文亚洲av片在线观看爽| 亚洲五月婷婷丁香| 亚洲全国av大片| 桃色一区二区三区在线观看| av在线播放免费不卡| 国产av一区在线观看免费| 亚洲男人天堂网一区| 国产精品电影一区二区三区| 丝袜美腿诱惑在线| 国产精品一区二区精品视频观看| 又大又爽又粗| 纯流量卡能插随身wifi吗| 很黄的视频免费| 亚洲国产精品999在线| 国产成人啪精品午夜网站| 国产精品久久久人人做人人爽| 久久久久精品国产欧美久久久| 少妇被粗大的猛进出69影院| 精品一区二区三区av网在线观看| 日本精品一区二区三区蜜桃| 亚洲国产欧美网| 香蕉丝袜av| 欧美精品啪啪一区二区三区| 亚洲人成电影观看| 久久人人97超碰香蕉20202| 成人国语在线视频| av中文乱码字幕在线| 亚洲人成电影免费在线| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久精品国产欧美久久久| 1024香蕉在线观看| 麻豆一二三区av精品| 欧美日韩乱码在线| 精品不卡国产一区二区三区| 国产精品美女特级片免费视频播放器 | 电影成人av| 精品国产亚洲在线| 亚洲欧美激情综合另类| 宅男免费午夜| 亚洲情色 制服丝袜| 国产黄a三级三级三级人| 亚洲国产日韩欧美精品在线观看 | 久久久国产成人精品二区| 一区二区日韩欧美中文字幕| 精品国产乱子伦一区二区三区| 亚洲精品粉嫩美女一区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美av亚洲av综合av国产av| 伦理电影免费视频| av有码第一页| 岛国在线观看网站| 午夜精品在线福利| 99国产极品粉嫩在线观看| 国产亚洲欧美精品永久| 精品乱码久久久久久99久播| 99国产精品99久久久久| 亚洲aⅴ乱码一区二区在线播放 | 曰老女人黄片| 大香蕉久久成人网| 国产男靠女视频免费网站| 欧美中文综合在线视频| 变态另类丝袜制服| 亚洲最大成人中文| 给我免费播放毛片高清在线观看| 久久人人97超碰香蕉20202| 少妇的丰满在线观看| 国产欧美日韩一区二区三| 日日摸夜夜添夜夜添小说| 男女之事视频高清在线观看| 久久精品91蜜桃| 极品人妻少妇av视频| 亚洲片人在线观看| 可以在线观看毛片的网站| 搡老妇女老女人老熟妇| 动漫黄色视频在线观看| 最新在线观看一区二区三区| 1024香蕉在线观看| 久久精品国产综合久久久| cao死你这个sao货| 女人被躁到高潮嗷嗷叫费观| 大码成人一级视频| 国产精品一区二区三区四区久久 | 两性午夜刺激爽爽歪歪视频在线观看 | 午夜a级毛片| 人妻丰满熟妇av一区二区三区| 极品人妻少妇av视频| 中文字幕人成人乱码亚洲影| 午夜免费鲁丝| 一夜夜www| 在线观看www视频免费| 非洲黑人性xxxx精品又粗又长| 搡老岳熟女国产| 9色porny在线观看| 国内毛片毛片毛片毛片毛片| 亚洲熟女毛片儿| 三级毛片av免费| 亚洲精品久久成人aⅴ小说| 高潮久久久久久久久久久不卡| 一个人免费在线观看的高清视频| 岛国在线观看网站| 久久香蕉激情| 国产亚洲精品第一综合不卡| 制服人妻中文乱码| 极品教师在线免费播放| 国语自产精品视频在线第100页| 午夜精品在线福利| 国产精品久久久av美女十八| 亚洲在线自拍视频| 精品免费久久久久久久清纯| 亚洲一区二区三区色噜噜| 狂野欧美激情性xxxx| 免费在线观看完整版高清| 精品国产国语对白av| 亚洲国产欧美一区二区综合| 久久人人97超碰香蕉20202| 国产成人啪精品午夜网站| av有码第一页| 亚洲欧美精品综合一区二区三区| 99精品久久久久人妻精品| 一a级毛片在线观看| 99精品在免费线老司机午夜| 69av精品久久久久久| 一区二区三区高清视频在线| 国产精品久久电影中文字幕| 久久久久久久午夜电影| 国产精品1区2区在线观看.| 三级毛片av免费| 在线永久观看黄色视频| 欧美日本视频| 禁无遮挡网站| av天堂在线播放| 精品欧美一区二区三区在线| 在线观看免费视频网站a站| 日韩欧美在线二视频| 无遮挡黄片免费观看| 久久久久国产一级毛片高清牌| 国产精品久久久久久亚洲av鲁大| 欧美日韩一级在线毛片| 波多野结衣高清无吗| 十分钟在线观看高清视频www| 成人国产一区最新在线观看| 十分钟在线观看高清视频www| 国产亚洲精品第一综合不卡| 老司机福利观看| 午夜免费观看网址| 精品国产乱子伦一区二区三区| 亚洲久久久国产精品| 99精品在免费线老司机午夜| 久久久久久大精品| 欧美 亚洲 国产 日韩一| 亚洲一码二码三码区别大吗| 久久国产精品影院| 制服诱惑二区| 亚洲av成人不卡在线观看播放网| 美女 人体艺术 gogo| 亚洲国产毛片av蜜桃av| 国产成人啪精品午夜网站| 日韩欧美三级三区| 黄色毛片三级朝国网站| 久久中文字幕人妻熟女| 成年版毛片免费区| 国产精华一区二区三区| 好男人电影高清在线观看| 国产高清videossex| 国产精品精品国产色婷婷| 国产99白浆流出| 成熟少妇高潮喷水视频| 亚洲无线在线观看| 亚洲国产高清在线一区二区三 | 多毛熟女@视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产真人三级小视频在线观看| 亚洲欧美日韩高清在线视频| 免费搜索国产男女视频| 免费看a级黄色片| 国内精品久久久久精免费| 亚洲aⅴ乱码一区二区在线播放 | 夜夜看夜夜爽夜夜摸| 亚洲成国产人片在线观看| 国产黄a三级三级三级人| 麻豆成人av在线观看| 少妇的丰满在线观看| 美女高潮到喷水免费观看| 国产亚洲精品久久久久久毛片| 国产亚洲欧美98| 亚洲精品美女久久av网站| 亚洲欧美精品综合一区二区三区| x7x7x7水蜜桃| www.999成人在线观看| 欧美精品亚洲一区二区| 国产区一区二久久| 亚洲精品国产精品久久久不卡| 国语自产精品视频在线第100页| 午夜老司机福利片| 99精品欧美一区二区三区四区| 欧美国产精品va在线观看不卡| 国产熟女午夜一区二区三区| 午夜福利一区二区在线看|