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

    Hydrological response to climate change and human activities:A case studyof Taihu Basin,China

    2020-08-18 09:05:44JunWuZhiyongWuHejunLinHipingJiMinLiu
    Water Science and Engineering 2020年2期

    Jun Wu ,Zhi-yong Wu *,He-jun Lin ,Hi-ping Ji ,Min Liu

    a Bureau of Hydrology,Taihu Basin Authority of Ministry of Water Resources,Shanghai 200434,China

    bCollege of Hydrology and Water Resources,Hohai University,Nanjing 210098,China

    Abstract Climate change and human activities have changed a number of characteristics of river flow in the Taihu Basin.Based on long-term time series of hydrological data from 1986 to 2015,we analyzed variability in precipitation,water stage,water diversion from the Yangtze River,and net inflow into Taihu Lake with the Mann-Kendall test.The non-stationary relationship between precipitation and water stage was first analyzed for the Taihu Basin and the Wuchengxiyu(WCXY)sub-region.The optimized regional and urban regulation schemes were explored to tackle high water stage problems through the hydrodynamic model.The results showed the following:(1)The highest,lowest,and average Taihu Lake water stages of all months had increasing trends.The total net inflow into Taihu Lake from the Huxi(HX)sub-region and the Wangting Sluice increased significantly.(2)The Taihu Lake water stage decreased much more slowly after 2002;it was steadier and higher after 2002.After the construction of Wuxi urban flood control projects,the average water stage of the inner city was 0.16-0.40 m lower than that of suburbs in the flood season,leading to the transfer of flooding in inner cities to suburbs and increasing inflow from HX into Taihu Lake.(3)The regional optimized schemes were more satisfactory in not increasing the inner city flood control burden,thereby decreasing the average water stage by 0.04-0.13 m,and the highest water stage by 0.04-0.09 m for Taihu Lake and the sub-region in the flood season.Future flood control research should set the basin as the basic unit.Decreasing diversion and drainage lines along the Yangtze River can take an active role in flood control.? 2020 Hohai University.Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    Keywords:Hydrological response;Climate change;Human activities;Flood control;Mann-Kendall test;Taihu Basin

    1.Introduction

    The hydrological cycle of a catchment is a complex process that is influenced by climate,physical characteristics of the catchment,and human activities(Liu and Xu,2015;Wang et al.,2011b;Zhu et al.,2015b;Xia et al.,2018;Hasan and Wyseure,2018).Climate determines the main inputs to the hydrological cycle,such as rainfall and radiation.Human activities have persistent influences on runoff distribution.In the past sixty years,both global climate change and intensive human activities have caused remarkable influences on hydrological process all over the world(Liu et al.,2016).A number of studies have shown that elements of the hydrological cycle,such as precipitation and water stage,have changed both in time and space,due to hydraulic engineering construction(Gu et al.,2017).Urbanization on river systems has been widely recognized as the most significant influence of all human activities(Shi et al.,2010),leading to hydrological regime variation and uncertainties in flood control.With impervious areas associated with buildings and transportation infrastructure increasing(Gregory,2006;Zhang et al.,2017),both surface runoff and flow velocity into rivers have increased significantly,leading to frequent regional precipitous flood occurrence (Yuan et al., 2006).

    With the background of global climate change (Todorov et al., 2018), the distribution of water resources and the frequency of floods and droughts for the Taihu Basin in recent years has undergone a lot of change,elevating the discrepancy between water supply and demand (Zhou et al., 2013). Most rivers in the Taihu Basin have been interrupted by sluices,pump stations,and floodgates(Zhao et al.,2011).There are 68 sluices along the Yangtze River,186 sluices surrounding Taihu Lake, and nine sluices along Hangzhou Bay to defend against flooding and meet water resources demand (Yin et al., 2009).

    According to Li et al. (2013), urbanization in China has been in a state of acceleration since 1978. Following a basinwide catastrophic flood in 1991, the first round of comprehensive regulation projects were begun. The three main basin flood control projects(the Wangyu River,the Taipu River,and Lake Levee) were completed by the end of 1999, while regional flood control projects were completed by the end of 2002. The Taihu Basin Authority (TBA) began to carry out water diversion from the Yangtze River to the Taihu Basin through the Wangyu River and to supply water downstream through the Taipu River in 2002. To ensure flood prevention security in the cities, including Suzhou, Huzhou, Jiaxing,Wuxi,and Changzhou,urban flood control projects have been built around Taihu Lake since 2003 with a flood protection standard of 100 years or 200 years.The flood control capacity has been further enhanced since 2002, when the complex system of flood control infrastructure was formed, including dikes, large sluice gates, and pumping stations. This study addressed the effects of human activities on the basin, and on regional and urban flood control projects across the basin.Therefore, 2002 can be considered an important demarcation point in this study or any study of hydrological regime variation in the Taihu Basin.

    Research in the Taihu Basin has focused on hydrological simulation(Liu et al.,2013;Wang and Yang,2007),flood risk management (Gong and Lin, 2009), water quality (Ye et al.,2017), and ecological problems (Liu et al., 2018). Research on hydrological response to urbanization has been a topic of wide concern in recent years (Zhang et al., 2018). Zhu et al.(2015a) found that river connectivity in Shanghai degraded significantly and land use changed with more building areas and fewer agricultural areas. Yang et al. (2014) analyzed precipitation differences between cities and suburbs in the Taihu Basin and revealed that the maximum daily precipitation and the number of rainstorm days increased significantly.Deng et al.(2015)investigated temporal and spatial change of river systems in the Taihu Basin, and demonstrated that changes in river density, water surface ratio, and main river area-to-length ratio in the rapid urbanization period were much greater than those in the slow urbanization period.Wang et al. (2016) assessed the contributions of precipitation and human activity to water stage increase in the plain river network region of the Taihu Basin, demonstrating that human activities have played more and more important roles in extreme water stage changes since the late 1980s.

    By the end of 2015, the population of Taihu Basin reached 59.97 million and the GDP was 6 688 billion Yuan, representing about 4.4% and 9.9% of the nation's totals, respectively. Due to its rapid social-economic development, the natural hydrological cycle in the Taihu Basin has been converted to an artificially dominated natural-social water cycle,and the Taihu Basin suffers high vulnerability to natural disasters (Zhao and Wen, 2012). Therefore, it is important to evaluate the impacts of climate change and human activities,summarize new characteristics of potential threats to storage and drainage capacity, and propose countermeasures to deal with the challenging situation of hydrological regime variation in the basin.

    In this study, both climate change (precipitation) and changes in human activities (water diversion and drainage along the Yangtze River, inflow and outflow of Taihu Lake,urban flood control projects, and urbanization) were considered to identify hydrological response at different spatiotemporal scales. First, the Mann-Kendall test was used to detect the trends of sub-region precipitation of the Taihu Basin, the water stage of Taihu Lake, water diversion and drainage along the Yangtze River, and inflow and outflow of Taihu Lake. Second, possible reasons for Taihu Lake water stage change and increase of net inflow into Taihu Lake were examined after considering the operation of Wuxi urban flood control projects. Third, with the aim of increasing basin and regional drainage ability and tackling high water stage problems in the Jiangnan Canal, optimized regulation schemes of regional and urban projects were determined through hydrodynamic models. These can provide scientific suggestions for further research, such as impacts of urbanization on the hydrological cycle, water resources availability, and sustainable development of water safety in the Taihu Basin.

    2. Study area and data

    2.1. Study area

    The Taihu Basin (30°28′N to 32°15′N, and 119°11′E to 121°53′E) is located in the Yangtze River Delta, with the Yangtze River at the north, the Qiantang River at the south,and the East Sea at the east. The western parts of the Taihu Basin are mountainous areas(Peng et al.,2016).The total area of the Taihu Basin is 36 895 km2, with portions in Jiangsu Province covering 19 399 km2, portions in Zhejiang Province covering 12 095 km2, portions in Shanghai City covering 5 176 km2,and portions in Anhui Province covering 225 km2.The hydrographic network of the Taihu Basin is complicated(Deng et al., 2016), with rivers and lakes accounting for 17%of the total area of the basin. The elevation of low land and polder around the Taihu Lake is no more than 5 m.About 80%of the Taihu Basin is plains, while the remaining 20% is occupied by low hills and mountains. As the third largest freshwater lake in China, with a water area of 2 338 km2,Taihu Lake is a flood storage and regulation center, whose average water depth is 1.89 m and water volume is 4.428 billion m3. The average lake bottom elevation is 1 m above mean sea level. The longest river in the Taihu Basin is the Jiangnan Canal,with a total length of 203 km,flowing through Changzhou, Wuxi, and Suzhou.

    The average annual precipitation of the Taihu Basin is 1 218.1 mm, mainly concentrated in the flood season from May to September (Wang et al., 2011a). The normal Meiyu period with intensive precipitation from mid-June to early July often lasts for 25 d, accounting for about 19.8% of the annual precipitation, but large differences of precipitation amount exist in various years. For example, the Meiyu did not appear in 1958 and 1978, whereas the precipitation amounts of the Meiyu period in 1991 and 1999 were much greater than the average precipitation.

    2.2. Data

    Daily precipitation and evaporation of hydrometric stations,water stage stations of Taihu Lake and Wuxi, and the water diversion and drainage amounts along the Yangtze River were excerpted from theAnnual Hydrological Reports. The inflow and outflow data of Taihu Lake were obtained from the Bureau of Hydrology, Taihu Basin Authority of Ministry of Water Resources.

    (1)Precipitation data were from 1986 to 2015.As shown in Fig. 1, based on the topography and water features, the whole basin can be divided into seven sub-regions with 106 hydrometric stations to calculate average precipitation,including the Zhexi (ZX) sub-region with 23 stations, the Huxi (HX) subregion with 19 stations, the Taihu (TH) sub-region with eight stations, the Wuchengxiyu (WCXY) sub-region with 12 stations, the Yangchengdianmao (YCDM) sub-region with 13 stations,the Hangjiahu(HJH)sub-region with 17 stations,and the Pudongpuxi (PDPX) sub-region with 14 stations. The arithmetic mean method was applied in calculating the average precipitation of each sub-region, while the areaweighted mean method was applied in calculating the average precipitation of the Taihu Basin based on sub-region area. In Fig. 1, the green shadow represents polder areas greater than 33.3 km2, and the blue line represents the Jiangnan Canal, flowing through Changzhou, Wuxi, and Suzhou.

    Fig. 1. Taihu Basin and its sub-regions.

    (2) Water stage data were from 1986 to 2015. The Taihu Lake water stage is the average value of five stations surrounding Taihu Lake: Dapukou, Wangtingtai, Xishan, Xiaomeikou, and Jiapu. The Wuxi suburb water stage station is located in the Jiangnan Canal close to urban flood control projects, while the Wuxi inner station is located in the urban area, which is considered to be protected by the urban flood control projects.Suffering from severe land subsidence caused by extensive groundwater exploitation since the 1980s, the maximum cumulative subsidence was over 1 m for some stations in the Taihu Basin (Wu et al., 2009). However, the water stage in theAnnual Hydrological Reportswas observed based on a stationary datum. The occurrence of land subsidence leads to height datum distortion. Thus, the water level ought to be modified. All the water stage data were based on the Zhenjiang Wusong elevation system and modified before trend analysis. Daily observed flood events with different magnitudes of precipitation and peak feature(single-or multipeak) of Taihu Lake and Wuxi water stage stations were selected to evaluate urbanization effects.

    (3) The water diversion data, including inflow and outflow of Taihu Lake, were from 1986 to 2015. The water diversion project from the Yangtze River to the Taihu Basin through the Wangyu River and the supply downstream through the Taipu River were launched in 2002. Water diversion and drainage amount were from 14 sluices along the Yangtze River, including four sluices in HX, four sluices in WCXY, the Wangyu Sluice, and five sluices in YCDM.Inflow and outflow of Taihu Lake were from 23 tour gauging discharge stations, including the Taipu Sluice, the Wangting Sluice, three stations (Changxing, Yangjiabu, and Hangchangqiao) in ZX, six stations (Zhihugang, Wujingang,Yapugang,Wuxiqiao,Chengdonggangqiao,and Dagangqiao)in HX, five stations (Daxuanhe pump station, Dushan Sluice,Meilianghu pump station, Wulihu Sluice, and Yanhu Sluice)in WCXY, four stations (Tongkeng Sluice, Xujiangdaqiao,Guajingkou, and Lianhuqiao) in YCDM, and three stations(Chengbei Sluice, Huanlou Sluice, and Tuanjieqiao) in HJH.All the stations surrounding Taihu Lake are shown in Fig. 2.

    (4) Land use data were acquired by digitizing the 1985,2000, and 2010 land cover maps (1:10 000), which were monitored by the interpreted thematic mapper(TM).The land cover maps were classified into four categories: cropland,forest, urban area, and water area. All scenes adopted from Landsat TM by manual interpretation included radiometric and geometric corrections. The land cover map of 2010 was adopted in the hydrodynamic model.

    (5)Evaporation and discharge data from 1999,2000,2003,2008, 2009, and 2010 were used in the model. Daily evaporation and discharge data from the Gangkou hydrometric station, Deqing hydrometric station, Fushi Reservoir, and Qingshan Reservoir were used in the Xin'anjiang model.Daily evaporation data from the Hangchangqiao Station and the Taipu Sluice were used in the river network hydrodynamic model.

    3. Methods

    3.1.Nonparametric Mann-Kendall test for monotonic trend

    The Mann-Kendall test is a simple but useful method to identify changes in long-term hydrological series (Mann,1945), especially for complicated basins where it is difficult for mathematical modeling to provide satisfactory results(Guo et al.,2018a).In the Mann-Kendall test,the null hypothesisH0states that there is no trend in the population data set(x1,x2,…,xn)(Guo et al.,2018b).The alternative hypothesisH1of a twosided test is that a monotonic trend exists in the data set. The statisticSis computed as follows:

    For sample sizes larger than 10, the statisticsSare approximately normally distributed, with meanE(S) and varianceV(S) as follows:

    The standardized statistics (Z) are formulated as follows:

    Positive and negativeZvalues indicate upward and downward trends, respectively. A value ofZ1-α/2is determined using a standard normal distribution table.The existing trend is considered to be statistically significant if |Z|>Z1-α/2.In this study,the significance level α was set to be 0.05 and 0.10 withZ1-α/2values equal to 1.96 and 1.65,respectively.

    3.2.Runoff model and hydrodynamic model in Taihu Basin

    In order to study the effects of water diversion projects and urban flood control projects on hydraulic structure regulation,a water quantity numerical model for the Taihu Basin was adopted.It was composed of rainfall-runoff simulation and the river network hydrodynamic model. Developed by Hohai University (Guo et al., 2016), the rainfall-runoff simulation provides the lateral inflow and discharge boundary for the river network hydrodynamic model. The water quantity numerical model is composed of six parts:(1)rainfall-runoff simulation,(2) unsteady flow simulation of river networks, (3) generalization of river networks and lakes, (4) generalization and simulation of hydraulic structures, (5) simulation of boundary conditions, and (6) calibration and verification.

    Fig. 2. Locations of water stage and discharge stations surrounding Taihu Lake.

    (1) Rainfall-runoff simulation: The runoff model includes runoff yield and concentration. Both hilly regions and plain river networks exist in the Taihu Basin, and the traditional Xin'anjiang model(Zhao,1992)was applied to hilly regions of the Zhexi sub-region, whose discharge was offered to plain areas as inputs. Each sub-region in plain areas was classified into four kinds of underlying surfaces: water surface, paddy fields, arid land, and urban construction land. Therefore,different runoff models and confluence processing patterns were adopted based on different characteristics of confluence inside and outside of polder areas. For the water surface,runoff could be obtained using rainfall minus evaporation.For paddy fields, irrigation programs of different growth phases were considered, as well as water demand coefficients, infiltration rate, and irrigation-drainage mode. For arid land, the Xin'anjiang model of three layers of evaporation was used.For urban construction land, runoff could be obtained using net rainfall multiplied by the runoff coefficient.Drainage modulus was considered in the polder area confluence, and a plain region confluence unit hydrograph was used outside of polder areas.

    (2) Unsteady flow simulation of river networks: Basic equations of one-dimensional unsteady flow in the aqueduct included the continuous equation and momentum equation(Chen et al.,2016).At each river network node,conservation of mass and momentum was available. The river network was dispersed with the four-point implicit difference method,and basic equations were resolved with the iteration method.

    (3)Generalization of river networks and lakes:As shown in Fig.3,1 482 rivers,4 275 cross-sections,and 1 164 nodes were generalized in the hydrodynamic model, including 138 storage nodes and 215 control buildings.The whole generalized area for the Taihu Basin was 36847.8 km2(which is very close to the real basin area), including the plain regions of 28534.4 km2,river networks and lakes of 1 083.3 km2, and hilly regions of 7 230.1 km2.

    Fig. 3. Generalized rivers and nodes for hydrodynamic model in Taihu Basin.

    (4) Generalization and simulation of hydraulic structures:The hydraulic structures include water gates, ship locks, culverts, and pump stations. The nodes were set at upstream and downstream of hydraulic structures, and discharge was calculated with the weir formula.

    (5) Simulation of boundary conditions: The boundary conditions included precipitation, evaporation, and tidal level. For each sub-region, the initial water stage of each cross-section was assigned the same value, while the initial discharge was 0. The boundary condition along the Yangtze River and Hangzhou Bay was the tidal level. For tidal stations with observation data, the unit hydrograph was adopted to calculate integral tidal level based on daily semidiurnal tide. For tidal stations lacking observation data, the tidal level interpolation formulas were adopted through relative distance with observed stations.

    (6) Calibration and verification: Due to limitations of the evaporation data, only data from the years 1999, 2000, 2003,2008,2009,and 2010 were collected.In order to keep consistent with typical scenarios of existing basin and regional planning,the data of years 1999,2000,2008,and 2009 were chosen in the calibration period to represent extraordinarily abundant,abundant, normal, and scarce precipitation, respectively. The 2003 and 2010 data used in the verification period are two typical scenarios representing scarce precipitation and normal-topartially abundant precipitation. The Xin'anjiang model parameters that needed to be calibrated for hilly regions of the Zhexi sub-region included four modules: evapotranspiration,runoff generation, runoff separation, and runoff routing. The parameters of hydrodynamic models for plain regions included roughness, the free flow coefficient, and the submerged flow coefficient.The model parameters are listed in Table 1.

    Nash-Sutcliffe efficiency (NSE) was used to evaluate model performance. WhenNSEapproaches 1.0, the model simulates the measured data perfectly. WhenNSEis negative,the model is a worse predictor than the measured mean value.The equation forNSEis as follows:

    whereZoiandZsiare the observed and simulated water levels at timei,is the mean observed data over the simulation period, andnis the total number of observations. For the calibration period at the daily level in this study,NSEwas 0.98. For the validation period,NSEwas 0.97.

    Table 1 Parameters of Xin'anjiang model and river network hydrodynamic model.

    Table 2 Land use changes in Taihu Basin from 1985 to 2010.

    4. Results

    From 1985 to 2010, the most noticeable change in land use was conversion of cropland to urban area. In 1985,cropland occupied the largest proportion of land (64%),followed by forest (14%), water area (13%), and urban area(10%). In 2000, the proportion of cropland decreased to 53%, while the proportion of urban area increased to 19%.In 2010, the proportion of cropland decreased to 48%,whereas the proportion of urban area increased to 24%. The land use change is shown in Table 2. In terms of the change rate, the urban area has undergone the most significant changes, by 151% (from 3 570 km2in 1985 to 8 943 km2in 2010), followed by cropland, with a decreasing rate of 25%(from 23 482 km2in 1985 to 17 675 km2in 2010). The change was more significant from 1985 to 2000. As shown in Fig. 4, small towns expanded around the former urban center, while the transport infrastructure was well developed, causing the urban boundary to blur and city growth to become scattered.

    4.1. Trends of precipitation, water stage, diversion,inflow, and outflow

    The increasing trends of annual precipitation and precipitation in the flood season (from May to September) over the Taihu Basin were not statistically significant. As shown in Fig.5,monthly basin precipitation decreased significantly at a significance level of 0.05 in April,May,and September,withZvalues of -2.54, -2.33, and -2.82, respectively. The subregion precipitation in the non-flood season (January to March,October to December)demonstrated increasing trends,but only precipitation in January increased significantly at a significance level of 0.05.Monthly sub-region precipitation in April, May, and September decreased significantly at a significance level of 0.1, except in HX.

    The increasing trend of the highest annual Taihu Lake water stage(HWS)was not significant,with aZvalue of 1.45.However, theZvalue of the lowest annual Taihu Lake water stage (LWS) reached 5.08, and theZvalue of the average annual Taihu Lake water stage (AWS) reached 2.57, which indicated that LWS and AWS increased significantly at a significance level of 0.05. As shown in Fig. 6(a), the HWS,LWS,and AWS of all months had increasing trends,especially in January, February, March, April, August, and December.Apart from this, the LWS in May, July, September, and November also increased significantly, while the AWS in September and November increased significantly (Fig. 6(b)).Considering the 11 key projects completed and water diversion begun in 2002,the HWS,LWS,and AWS variations were compared before 2001 and after 2002. As shown in Fig. 6(b),from January to April and in December,the HWS increased by more than 0.18 m, the LWS increased by more than 0.19 m,and the AWS increased by more than 0.17 m. The largest increasing extent of the HWS, LWS, and AWS reached more than 0.29 m in March, whereas the least increasing extent reached more than 0.04 m in July.

    Fig. 4. Land use changes in Taihu Basin from 1985 to 2010.

    Fig. 5. Trend analysis of precipitation over Taihu Basin from 1954 to 2015.

    Monthly net inflow into the Taihu Lake were defined as the daily accumulative amount of water diversion minus drainage. As shown in Fig. 7(a) and (b), net inflow from HX into Taihu Lake had significant increasing trends for all months,and inflow from the Wangting Sluice into Taihu Lake had insignificant increasing trends in non-flood seasons except for April. Compared with net inflow before 2001, net inflow from HX after 2002 increased by 120 × 106m3to 402 × 106m3, while the net inflow from the Wangyu Sluice increased by 5 × 106m3to 185 × 106m3. However, net inflow from ZX into Taihu Lake decreased by 36×106m3to 202 × 106m3in the flood season, and decreased by 1×106m3to 76×106m3in the non-flood season except in October.As shown in Fig.7(c)and(d),the net outflow of the Taipu Sluice had decreasing trends in the non-flood season except in March. In April, May, June, July, September,October, and November, the net outflow of YCDM had decreasing trends,whereas that of HJH had increasing trends.Compared with net outflow before 2001, net outflow of HJH after 2002 increased by 32×106m3to 155×106m3,and the Taipu Sluice decreased by 49×106m3to 72×106m3in the non-flood season except in March.The net outflow of YCDM decreased by 11 × 106m3to 121 × 106m3from April to November, and that of WCXY increased by 11 × 106m3to 34 × 106m3in the non-flood season.

    The flood movement of the Taihu Basin changed a lot.During the period from 2002 to 2015, total water diversion from the Yangtze River through the Wangyu River reached 272.4 × 108m3, with net inflow into Taihu Lake being 124.2 × 108m3, and water moving downstream through the Taipu River reaching 182.8 × 108m3. Compared with the annual net flow (differences between inflow and outflow)before 2001, the annual net inflow from HX increased by 24.9 × 108m3, the value from the Wangting Sluice increased by 7.1 × 108m3, and the value from ZX decreased by 7.1 × 108m3. However, net outflow of HJH increased by 11.4×108m3,the value of WCXY increased by 1.9×108m3,the Taipu Sluice decreased by 3.0 × 108m3, and the value of YCDM decreased by 3.7 × 108m3.

    4.2.Possible reasons for Taihu Lake water stage increase and flood risks in Taihu Basin

    As mentioned above,the year of 2002 can be considered an important demarcation point to study hydrological regime variation in the Taihu Basin.In order to eliminate the influence of water diversion on the Taihu Lake water stage, we selected 218 flood events from 1986 to 2015 to study the relationships between the Taihu Lake water stage increase and precipitation under low net water diversion conditions. The Taihu Lake water stage increase has a strong relationship with precipitation,and the correlation coefficient was 0.86 before 2001,but 0.83 after 2002. Similarly, in order to eliminate the influence of precipitation on the Taihu Lake water stage increase, we selected 233 flood events to study the relationships between the Taihu Lake water stage increase and net water diversion from the Yangtze River under low-precipitation conditions.The correlation coefficient was 0.30 before 2001, but 0.39 after 2002.

    Fig. 6. Trends and increase of HWS, LWS, and AWS.

    Fig. 7. Trends and variation of net inflow and net outflow of different sub-regions.

    For each period of water diversion from the Yangtze River into Taihu Lake after 2002, the Taihu Lake water stage demonstrated a very slow decreasing trend, taking five to six days to complete a 0.01-m decrease. However, for the same period before 2001,it took only two to three days to complete a 0.01-m decrease. Therefore, the Taihu Lake water stage decreased much more slowly after 2002, becoming steadier and generally higher after 2002.Therefore,the increase of the Taihu Lake water stage is still dominated by precipitation,while the decrease is affected by water diversion from the Yangtze River after 2002.

    For the period before 2001, we selected 116 heavy flood events with an average basin precipitation of 57.8 mm, while the Taihu Lake water stage increased from 3.10 m to 3.22 m on average. For the period after 2002, we selected 39 flood events with an average precipitation of 55.9 mm, while the water stage increased from 3.27 m to 3.38 m on average.Compared with water stage increase for the same precipitation before 2001, shown in Table 3, water stage increased more after 2002 when precipitation exceeded 70 mm, while the water stage increased more significantly when precipitation exceeded 120 mm. For the same 0.01 m of the Taihu Lake water stage increase after 2002,5 mm of average precipitation was needed under the 50-mm precipitation scenario, 4 mm of average precipitation was needed under the 100-mm scenario,and only 3.8 mm of precipitation was needed for the 150-mm scenario.

    The Wuxi urban flood control projects were put in use in 2007,and can be used as an example to analyze the effects of human activities on water stage in the WCXY sub-region.The Wuxi urban flood control projects are located in the WCXY sub-region, downstream of HX and upstream of YCDM.Before construction of urban flood control projects, the water stage of Wuxi inner cities was the same as that of suburbs.As shown in Fig. 8, the annual average water stage of the inner city was 0.16-0.40 m lower than that of suburbs in flood seasons after construction. We selected 36 intensive flood events to explore non-stationary relationships between area precipitation and regional average water stage differences in the WCXY sub-region. Before intensive precipitation occurred, the average water stage difference between inner cities and suburbs was 0.24 m. However, after precipitation,the average water stage difference increased to 0.46 m on average. With more intensive precipitation, the water stage of suburbs increased more quickly,and more differences between inner cities and suburbs arose. As shown in Fig. 9(a), with more precipitation, the maximum water stage differences ΔZmaxbetween inner cities and suburbs grew larger.A 10-mm precipitation increase may lead to a 0.054-m ΔZmaxincrease on average. Since an initial water stage difference ΔZintmayhave already existed before intensive precipitation occurred,we also built relationships between precipitation and the real increase of water stage differences. The real increase of water stage differences ΔZrealis defined as ΔZmaxminus ΔZint. As shown in Fig.9(b),a 10-mm precipitation increase may lead to 0.046 m of ΔZreal. Therefore, flooding in inner cities was transferred to suburbs through urban flood control projects,leading to higher and higher water stages in suburbs,which are especially severe in the Jiangnan Canal.Instead of flowing into WCXY, flooding upstream was forced from HX into Taihu Lake, increasing the inflow from HX into Taihu Lake.

    Table 3 Taihu Lake water stage increase for different precipitation scenarios before 2001 and after 2002.

    Fig. 8. Average water stage differences of Wuxi inner cities and suburbs and precipitation distribution.

    The peak lag time of inflow from HX into Taihu Lake was defined as the time duration between precipitation peak and inflow peak. Based on the flood fluctuation process under low net water diversion along the Yangtze River,the peak lag time of HX inflow into Taihu Lake was one to two days before 2001,while the peak lag time decreased to only one day after 2002.For the same inflow into Taihu Lake of 16×106m3,16 mm of precipitation was needed before 2001, while only 11 mm precipitation was needed after 2002. With a shorter peak lag time, more inflow was concentrated in shorter time intervals.Thus,the efficiency of flow concentration increased after 2002.

    Possible reasons for Taihu Lake water stage increase can therefore be divided into three categories:

    (1) The water diversion increase along the Yangtze River:The HWS, LWS, and AWS of April and May increased,although monthly precipitation of the same period decreased,which may have been due to the water diversion of the Wangyu Sluice from the Yangtze River. The HWS, LWS, and AWS in August and September demonstrated significantly increasing trends, although the precipitation of the same period decreased significantly, which may have been caused by the water diversion of the Wangyu Sluice and HX from the Yangtze River storing water resources in the non-flood season.The HWS, LWS, and AWS in November and December increased significantly, which may have been due to the precipitation increase as well as water diversion of the Wangyu Sluice and WCXY from the Yangtze River in order to satisfy water requirements and improve the regional water environment.

    (2) Increase of net inflow into Taihu Lake: The increase of multi-year average inflow was more significant than that of outflow, and it was particularly significant for HX and the Wangting Sluice. Compared with the multi-year average annual net inflow before 2001,the annual net inflow after 2002 increased by 18.1 × 108m3. The increase of net inflow from HX into Taihu Lake may be due to increasing water diverted from the Yangtze River and frequent use of Wuxi urban flood control projects. The Taihu Lake water stage decreased much more slowly after 2002,causing the Taihu Lake water stage to be steadier and higher after 2002. Flooding in Wuxi inner cities was transferred to suburbs through urban flood control projects,leading to a higher water stage in the Jiangnan Canal.Instead of flowing into WCXY, flooding upstream was forced from HX into Taihu Lake,increasing inflow of HX into Taihu Lake.

    Fig. 9. Relationships between precipitation and maximum water stage differences and between precipitation and real increase of water stage differences.

    (3)Land cover change with more impervious areas:With a significant decrease in cropland and an increase in urban area,both surface runoff and flow velocity from HX into Taihu Lake increased significantly.More inflow from HX into Taihu Lake was concentrated in a shorter time interval and the efficiency of flow concentration increased after 2002, leading to a quicker and higher Taihu Lake water stage increase,as well as a heavy flood control burden for both the basin and regions.

    5.Discussion on optimized flood control strategy in Taihu Basin

    With the drainage power of urban flood control projects along the Jiangnan Canal(especially in Changzhou,Wuxi,and Suzhou)increasing,more water in inner cities has been moved to suburbs, leading to higher water stages and severe flood risks in the Jiangnan Canal. Therefore, a certain incompatibility has already emerged in projects’ regulation mode. In order to increase basin and regional drainage and tackle high water stage problems in the Jiangnan Canal, the optimized regulation schemes of regional and urban projects were explored through water quantity numerical model computation. The control of annual average water stage increase(ΔZ)was the basis for the optimized regulation scheme.

    One basic scheme and six optimized schemes with 1991 design precipitation were investigated in this study. The heavy precipitation in 1991 was caused by the Meiyu, and mainly concentrated in the northern part of the Taihu Basin.The homogenous frequency enlargement method was adopted to calculate the 1991 design precipitation of with a 50-year return period. For basic schemes, the regulation modes of regional and urban flood control projects are currently in use. As shown in Table 4, for urban optimized scheme 1(CO1), the water stage was increased by 50%ΔZon the drainage line of inner cities on the Jiangnan Canal.For urban optimized scheme 2 (CO2), the water stage was increased by 100%ΔZon drainage line of inner cities on the Jiangnan Canal. For regional optimized scheme 1 (RO1) along the Yangtze River,the water stage was reduced by 50%ΔZfor the regional drainage line. For regional optimized scheme 2(RO2), the water stage was reduced by 50%ΔZfor the regional diversion line,while the water stage was reduced by 100%ΔZfor the regional drainage line. The urban and regional optimized scheme 1 (CRO1) was a combination scheme of CO1 and RO1, while the urban and regional optimized scheme 2 (CRO2) was a combination scheme of CO2 and RO2.

    As shown in Table 5, variation between basic and optimized schemes for average and highest water stage in the flood season was compared, with HX, WCXY, and YCDM representing regional water stage. The CO1 and CO2 had no clear effects on reducing the average Taihu Lake and regional water stages in the flood season. RO1 decreased the stage by 0.04-0.06 m, RO2 decreased the stage by 0.08-0.13 m, the effect of CRO1 was basically identical to that of RO1,and the effect of CRO2 was basically identical to that of RO2. The water stage decrease in WCXY was more significant,followed by the decreases in HX and YCDM.For reducing the average suburb water stage, RO1 was better than CO1, and RO2 was better than CO2, while CRO1 included comprehensive effects of CO1 and RO1, and CRO2 included the comprehensive effects of CO2 and RO2.The decrease of the suburb water stage was 0.05-0.06 m for RO1, while the decrease was 0.11-0.13 m for RO2, and more significant in Changzhou.The decrease of the suburb water stage was 0.06-0.08 m for CRO1, while the decrease was 0.13-0.18 m for CRO2, and more significant in Suzhou. For CRO1 and CRO2, the suburb water stage decrease was more significant than the city water stage increase.

    For reducing the highest Taihu Lake and regional water stage, the combined urban and regional optimized scheme CRO1 (CRO2) was more effective than corresponding regional optimized scheme RO1 (RO2), followed by corresponding urban optimized scheme CO1 (CO2). Changes of water stages in WCXY were more significant than those in HX and Taihu Lake for CRO1 and CRO2. The urban optimized schemes decreased the suburb water stage by 0.01-0.04 m,the regional optimized schemes decreased by 0.03-0.14 m,and urban and regional optimized schemes decreased by 0.06-0.15 m. However, for RO1, RO2, CRO1, and CRO2 in Suzhou, the suburb water stage decrease was more significant than the city water stage increase.

    Therefore, research on sustainable development of water safety should set the basin as the basic unit, and orchestrate flood control relationships among basin, region, and city. In order to achieve more harmonious regulation,the principles of a part being subordinate to the whole should be followed.From the aspect of the basin, regulation schemes should bemodified by adoption of more flexible methods based on regional and urban water regime variations. From the aspects of regional and city needs, the optimized regulation scheme should take overall consideration of highlighting flood control priorities and sharing risks, coordinating reasonable drainage lines of inner cities,and diversion and drainage lines along the Yangtze River.

    Table 4 Basic and optimized computation schemes.

    Table 5 Variations of average and highest water stages of optimized schemes.

    6. Conclusions

    In this study,precipitation,water stage,net water diversion along the Yangtze river,and inflow and outflow of Taihu Lake using long time series data were analyzed. The optimized regulation schemes of regional and urban projects were explored through a hydrodynamic model. The following conclusions were made:

    (1)The HWS,LWS,and AWS of all months had increasing trends.The HWS increased by more than 0.18 m,the LWS by more than 0.19 m, and the AWS by more than 0.17 m in the non-flood season after 2002. Total net inflow into Taihu Lake increased by 18.1×108m3,with HX and the Wangting Sluice increasing significantly after 2002. Monthly HWS, LWS, and AWS increased, although precipitation in the same period decreased, which may be due to the water diversion of the Wangyu Sluice,the WCXY sub-region,and the HX sub-region from the Yangtze River.

    (2) The Taihu Lake water stage decreased much more slowly after 2002, becoming steadier and higher after 2002.More inflow was concentrated in a shorter time interval,leading to higher inflow concentration efficiency into Taihu Lake.After construction of Wuxi urban flood control projects,the average water stage of the inner city was 0.16-0.40 m lower than that of suburbs in flood seasons. Flooding in inner cities was transferred to suburbs, leading to a higher water stage in suburbs of the WCXY sub-region. Instead of flowing into the WCXY sub-region, flooding upstream was forced from the HX sub-region into Taihu Lake, increasing HX inflow into Taihu Lake.

    (3) Urban and regional optimized schemes were more significant than single regional or urban optimized schemes.All six optimized schemes demonstrated advantages in water stage decrease in WCXY, Taihu Lake, and suburbs along the Jiangnan Canal (Changzhou, Wuxi, and Suzhou), among which regional optimized schemes acquired more satisfying effects of without increasing inner city flood control burden.Regional schemes decreased the average Taihu Lake and regional water stages by 0.04-0.13 m, and decreased the highest water stage by 0.04-0.09 m in the flood season.Future research on sustainable development of water safety should set the basin as the basic unit, and orchestrate flood control relationships between the basin, region, and city. Decreasing diversion and drainage lines along the Yangtze River could take an active role in flood control.

    Declaration of competing interest

    The authors declare no conflicts of interest.

    一二三四社区在线视频社区8| 啦啦啦韩国在线观看视频| 好男人在线观看高清免费视频| 久久99热这里只有精品18| 亚洲精品在线观看二区| 色老头精品视频在线观看| 国产精品,欧美在线| 日本五十路高清| 亚洲自拍偷在线| 免费av毛片视频| www.999成人在线观看| 久久精品国产清高在天天线| 男女做爰动态图高潮gif福利片| 美女高潮喷水抽搐中文字幕| 麻豆国产97在线/欧美| 18美女黄网站色大片免费观看| 最近在线观看免费完整版| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产亚洲欧美在线一区二区| 一本一本综合久久| 欧美在线一区亚洲| 国产精品永久免费网站| 一个人看的www免费观看视频| ponron亚洲| 宅男免费午夜| 男插女下体视频免费在线播放| 亚洲人与动物交配视频| 18+在线观看网站| 国产一级毛片七仙女欲春2| 99久久精品热视频| 一区二区三区免费毛片| 最新美女视频免费是黄的| 亚洲午夜理论影院| 欧美又色又爽又黄视频| 精品一区二区三区视频在线观看免费| 精品一区二区三区人妻视频| 国产av在哪里看| 女生性感内裤真人,穿戴方法视频| 波多野结衣高清作品| a级毛片a级免费在线| 性色av乱码一区二区三区2| 舔av片在线| 久久久久九九精品影院| 国产成人av教育| 亚洲欧美精品综合久久99| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 深爱激情五月婷婷| 天天躁日日操中文字幕| 欧美日韩一级在线毛片| 亚洲久久久久久中文字幕| 国产高清三级在线| 搞女人的毛片| 精品不卡国产一区二区三区| www.www免费av| 国产一区二区三区视频了| 嫩草影视91久久| 久久精品国产自在天天线| 久久婷婷人人爽人人干人人爱| 中文字幕人成人乱码亚洲影| 亚洲成人中文字幕在线播放| 欧美bdsm另类| 国产高清视频在线观看网站| 俄罗斯特黄特色一大片| 青草久久国产| 日本撒尿小便嘘嘘汇集6| 日韩欧美国产在线观看| 无限看片的www在线观看| 在线天堂最新版资源| 最好的美女福利视频网| 午夜a级毛片| 女生性感内裤真人,穿戴方法视频| 久久6这里有精品| 别揉我奶头~嗯~啊~动态视频| 色av中文字幕| 在线国产一区二区在线| 每晚都被弄得嗷嗷叫到高潮| 啦啦啦观看免费观看视频高清| 高清日韩中文字幕在线| 成人高潮视频无遮挡免费网站| 老司机在亚洲福利影院| 欧美三级亚洲精品| ponron亚洲| 亚洲av成人不卡在线观看播放网| 欧美日韩国产亚洲二区| 国产精品亚洲av一区麻豆| 国产欧美日韩一区二区精品| 欧美zozozo另类| 嫁个100分男人电影在线观看| 麻豆国产av国片精品| 亚洲最大成人手机在线| 天堂av国产一区二区熟女人妻| 亚洲天堂国产精品一区在线| 日本一本二区三区精品| 国产熟女xx| 午夜视频国产福利| 动漫黄色视频在线观看| 精品电影一区二区在线| 男人舔奶头视频| 丰满的人妻完整版| 最近视频中文字幕2019在线8| 久久精品国产亚洲av香蕉五月| 国产精品日韩av在线免费观看| 青草久久国产| 欧美色视频一区免费| 久久精品人妻少妇| av国产免费在线观看| 可以在线观看毛片的网站| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av中文字字幕乱码综合| 成人国产综合亚洲| 神马国产精品三级电影在线观看| 午夜福利18| 五月玫瑰六月丁香| 18禁在线播放成人免费| 欧美日韩乱码在线| 国产蜜桃级精品一区二区三区| eeuss影院久久| 色吧在线观看| 男人舔女人下体高潮全视频| 热99在线观看视频| 一个人免费在线观看的高清视频| 男女视频在线观看网站免费| 成人av在线播放网站| svipshipincom国产片| 国产69精品久久久久777片| 亚洲人成网站高清观看| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 国产成人福利小说| 欧美在线一区亚洲| 99热6这里只有精品| 最近最新中文字幕大全电影3| 色老头精品视频在线观看| 亚洲成人精品中文字幕电影| 听说在线观看完整版免费高清| 国产欧美日韩一区二区精品| 日韩欧美国产在线观看| 欧美日韩精品网址| eeuss影院久久| 欧美性感艳星| 一边摸一边抽搐一进一小说| 欧美成人免费av一区二区三区| 啦啦啦观看免费观看视频高清| 国产精品99久久99久久久不卡| 变态另类成人亚洲欧美熟女| 色噜噜av男人的天堂激情| a级毛片a级免费在线| 国产成人欧美在线观看| 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| 观看美女的网站| 欧美黄色淫秽网站| 亚洲在线观看片| 欧美日韩瑟瑟在线播放| 国产精品久久久久久亚洲av鲁大| 成年女人永久免费观看视频| 久久久久九九精品影院| 天堂影院成人在线观看| 变态另类成人亚洲欧美熟女| 黄色女人牲交| 精品免费久久久久久久清纯| 欧美色视频一区免费| 蜜桃久久精品国产亚洲av| 天天一区二区日本电影三级| 午夜福利高清视频| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 国产成人av激情在线播放| 熟女电影av网| 久久久久久久久大av| 日日干狠狠操夜夜爽| 国产成人欧美在线观看| 少妇的逼水好多| 国产高清视频在线播放一区| 免费高清视频大片| 免费观看人在逋| 亚洲avbb在线观看| 国产亚洲精品av在线| 久久精品国产亚洲av涩爱 | 日韩av在线大香蕉| 免费看日本二区| 人人妻人人看人人澡| 亚洲午夜理论影院| 十八禁网站免费在线| 亚洲精品粉嫩美女一区| 真人一进一出gif抽搐免费| 法律面前人人平等表现在哪些方面| 又黄又粗又硬又大视频| 国产乱人视频| 黄色女人牲交| av在线天堂中文字幕| 丰满乱子伦码专区| av视频在线观看入口| 国产淫片久久久久久久久 | 亚洲激情在线av| 99在线人妻在线中文字幕| 日韩欧美一区二区三区在线观看| av福利片在线观看| 欧美av亚洲av综合av国产av| 十八禁网站免费在线| 精品免费久久久久久久清纯| 99热6这里只有精品| 99热这里只有是精品50| 亚洲 欧美 日韩 在线 免费| 欧美bdsm另类| www.www免费av| 欧美另类亚洲清纯唯美| 91在线观看av| 久久久久亚洲av毛片大全| 九九热线精品视视频播放| 午夜福利在线观看吧| 精品久久久久久成人av| 成人高潮视频无遮挡免费网站| 国产精品影院久久| 欧美性猛交╳xxx乱大交人| 女人被狂操c到高潮| 亚洲无线观看免费| 日韩精品青青久久久久久| 欧美黄色片欧美黄色片| 亚洲欧美日韩卡通动漫| 在线看三级毛片| 免费一级毛片在线播放高清视频| 久久国产精品影院| www日本黄色视频网| www.www免费av| 天堂影院成人在线观看| 日韩精品中文字幕看吧| 中文字幕人成人乱码亚洲影| 亚洲七黄色美女视频| 国产成人欧美在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲最大成人中文| 欧美+亚洲+日韩+国产| 亚洲熟妇熟女久久| 国产真实乱freesex| 日韩欧美国产在线观看| 国产精品国产高清国产av| 日本五十路高清| 男女做爰动态图高潮gif福利片| 偷拍熟女少妇极品色| 国产97色在线日韩免费| 小蜜桃在线观看免费完整版高清| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 亚洲aⅴ乱码一区二区在线播放| 国产日本99.免费观看| 真实男女啪啪啪动态图| 成人国产综合亚洲| 国产精品,欧美在线| 久久精品91无色码中文字幕| 色综合欧美亚洲国产小说| 亚洲va日本ⅴa欧美va伊人久久| 久久久久亚洲av毛片大全| 琪琪午夜伦伦电影理论片6080| 中亚洲国语对白在线视频| 99久国产av精品| 少妇高潮的动态图| 午夜a级毛片| 麻豆成人午夜福利视频| xxx96com| 好男人电影高清在线观看| 成人午夜高清在线视频| 国内精品久久久久久久电影| 免费在线观看日本一区| 伊人久久精品亚洲午夜| 91av网一区二区| 国产黄色小视频在线观看| 亚洲成人久久爱视频| 啦啦啦观看免费观看视频高清| av天堂在线播放| 中国美女看黄片| 最近最新中文字幕大全电影3| 嫩草影视91久久| 日本a在线网址| 日韩欧美国产在线观看| 亚洲精品在线观看二区| 90打野战视频偷拍视频| 亚洲成人久久爱视频| 日韩国内少妇激情av| 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| netflix在线观看网站| av国产免费在线观看| 黄片大片在线免费观看| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 深夜精品福利| 可以在线观看的亚洲视频| 亚洲中文日韩欧美视频| 日本五十路高清| 免费高清视频大片| 欧美色视频一区免费| 9191精品国产免费久久| 久久人人精品亚洲av| 69人妻影院| 日本在线视频免费播放| 日日夜夜操网爽| 久久久久亚洲av毛片大全| 日韩高清综合在线| 亚洲国产日韩欧美精品在线观看 | 露出奶头的视频| 最新在线观看一区二区三区| 内射极品少妇av片p| 国产成人系列免费观看| 91在线精品国自产拍蜜月 | 黄色女人牲交| 成人性生交大片免费视频hd| 国产午夜精品论理片| 国产精品,欧美在线| 一级黄片播放器| 婷婷精品国产亚洲av| 国产高清视频在线观看网站| 亚洲五月婷婷丁香| 久久久精品大字幕| 国产一区二区三区视频了| 国产三级黄色录像| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 国产国拍精品亚洲av在线观看 | 亚洲国产日韩欧美精品在线观看 | 日本一本二区三区精品| 欧美日韩综合久久久久久 | 高清在线国产一区| 亚洲美女黄片视频| 欧美在线一区亚洲| 亚洲在线自拍视频| 亚洲中文字幕一区二区三区有码在线看| 俄罗斯特黄特色一大片| 国产精品一区二区免费欧美| 午夜福利在线观看免费完整高清在 | 老鸭窝网址在线观看| 狂野欧美白嫩少妇大欣赏| 国产成+人综合+亚洲专区| 在线视频色国产色| www日本在线高清视频| 波多野结衣高清无吗| 丝袜美腿在线中文| 高潮久久久久久久久久久不卡| 成人三级黄色视频| 性色avwww在线观看| 老司机深夜福利视频在线观看| 亚洲av成人av| 久9热在线精品视频| 人人妻人人看人人澡| eeuss影院久久| 亚洲精品久久国产高清桃花| 精品人妻一区二区三区麻豆 | 免费在线观看影片大全网站| 亚洲精品国产精品久久久不卡| 国产av一区在线观看免费| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 国产综合懂色| 亚洲精品在线观看二区| 美女cb高潮喷水在线观看| 欧美乱色亚洲激情| 久久精品国产亚洲av香蕉五月| 毛片女人毛片| 午夜视频国产福利| 日本 av在线| 97超级碰碰碰精品色视频在线观看| 白带黄色成豆腐渣| 99在线人妻在线中文字幕| 国产97色在线日韩免费| 精品福利观看| 国产一区二区激情短视频| 熟女少妇亚洲综合色aaa.| 小说图片视频综合网站| 91麻豆精品激情在线观看国产| av女优亚洲男人天堂| 国产精品国产高清国产av| 男女床上黄色一级片免费看| 国产亚洲欧美98| 亚洲精品在线美女| 欧美激情在线99| 脱女人内裤的视频| 女人高潮潮喷娇喘18禁视频| 国产亚洲精品久久久com| 在线国产一区二区在线| 久久婷婷人人爽人人干人人爱| 日本撒尿小便嘘嘘汇集6| 国产午夜福利久久久久久| 最近最新中文字幕大全免费视频| 免费观看精品视频网站| 国产免费av片在线观看野外av| 在线看三级毛片| 麻豆久久精品国产亚洲av| 免费人成视频x8x8入口观看| 啦啦啦韩国在线观看视频| 岛国视频午夜一区免费看| 成人18禁在线播放| 村上凉子中文字幕在线| 中亚洲国语对白在线视频| 美女被艹到高潮喷水动态| 久久精品国产自在天天线| 97超级碰碰碰精品色视频在线观看| 国产极品精品免费视频能看的| 嫩草影院精品99| 国产又黄又爽又无遮挡在线| 欧美又色又爽又黄视频| 精品久久久久久久末码| 国产成+人综合+亚洲专区| 亚洲狠狠婷婷综合久久图片| 色av中文字幕| 伊人久久大香线蕉亚洲五| 极品教师在线免费播放| 精品久久久久久久久久久久久| 91av网一区二区| 非洲黑人性xxxx精品又粗又长| 全区人妻精品视频| 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 午夜福利高清视频| 男女视频在线观看网站免费| 香蕉久久夜色| 亚洲熟妇中文字幕五十中出| 内地一区二区视频在线| 久久精品91蜜桃| 午夜福利免费观看在线| 国产精品永久免费网站| 国产三级中文精品| 一个人观看的视频www高清免费观看| 亚洲在线自拍视频| 日本 av在线| 蜜桃亚洲精品一区二区三区| 国产一区二区三区视频了| 叶爱在线成人免费视频播放| 精品久久久久久久久久免费视频| 亚洲精品日韩av片在线观看 | 国产日本99.免费观看| 国产老妇女一区| av视频在线观看入口| 成人精品一区二区免费| 国产欧美日韩一区二区三| 国产亚洲欧美98| 又粗又爽又猛毛片免费看| 亚洲最大成人手机在线| 黄色成人免费大全| 一个人免费在线观看电影| 窝窝影院91人妻| 亚洲精品色激情综合| 亚洲男人的天堂狠狠| 日韩成人在线观看一区二区三区| 国产美女午夜福利| 亚洲中文字幕一区二区三区有码在线看| 性色av乱码一区二区三区2| 国产亚洲精品综合一区在线观看| 国产av不卡久久| 男女那种视频在线观看| 91麻豆精品激情在线观看国产| 精品久久久久久成人av| 欧美日本亚洲视频在线播放| 黄色片一级片一级黄色片| 人妻夜夜爽99麻豆av| 禁无遮挡网站| 国产精品久久久久久久电影 | 在线视频色国产色| 久久久久九九精品影院| 午夜免费男女啪啪视频观看 | 久久6这里有精品| 在线a可以看的网站| av女优亚洲男人天堂| 免费看a级黄色片| 99久久综合精品五月天人人| 99热只有精品国产| 日韩大尺度精品在线看网址| 免费看美女性在线毛片视频| 麻豆成人av在线观看| 成人av在线播放网站| 村上凉子中文字幕在线| 精品人妻一区二区三区麻豆 | or卡值多少钱| 色老头精品视频在线观看| 男人的好看免费观看在线视频| 大型黄色视频在线免费观看| 色av中文字幕| 亚洲精品美女久久久久99蜜臀| 在线播放无遮挡| 夜夜爽天天搞| 老熟妇仑乱视频hdxx| 桃色一区二区三区在线观看| 久久这里只有精品中国| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 国产伦一二天堂av在线观看| 十八禁人妻一区二区| 日韩精品青青久久久久久| 亚洲av电影不卡..在线观看| 88av欧美| 欧美+日韩+精品| 久久久久国内视频| 国产乱人视频| 精品久久久久久成人av| 亚洲内射少妇av| 日本熟妇午夜| 真人一进一出gif抽搐免费| 午夜福利18| 欧美大码av| 国产一区二区三区在线臀色熟女| 国语自产精品视频在线第100页| 99在线视频只有这里精品首页| 国产午夜精品久久久久久一区二区三区 | 亚洲成av人片在线播放无| 两个人视频免费观看高清| 九九在线视频观看精品| 中文字幕精品亚洲无线码一区| 少妇人妻一区二区三区视频| 欧美大码av| 天天添夜夜摸| 亚洲专区国产一区二区| 国产美女午夜福利| 久久这里只有精品中国| 一级毛片高清免费大全| 91久久精品国产一区二区成人 | 色播亚洲综合网| 99久久久亚洲精品蜜臀av| 两人在一起打扑克的视频| avwww免费| 两个人看的免费小视频| 亚洲性夜色夜夜综合| 女同久久另类99精品国产91| 少妇人妻精品综合一区二区 | 国产av不卡久久| 天天一区二区日本电影三级| 欧美色视频一区免费| 一级a爱片免费观看的视频| 两个人看的免费小视频| 一个人免费在线观看电影| 亚洲精品在线美女| www国产在线视频色| 18禁国产床啪视频网站| 亚洲一区二区三区不卡视频| 午夜老司机福利剧场| 国产av在哪里看| 中文在线观看免费www的网站| 欧美性猛交黑人性爽| 亚洲精品一卡2卡三卡4卡5卡| 热99re8久久精品国产| 久久99热这里只有精品18| 人人妻人人看人人澡| 免费看十八禁软件| 2021天堂中文幕一二区在线观| 最近最新中文字幕大全电影3| 国内毛片毛片毛片毛片毛片| 一区二区三区激情视频| 老司机午夜福利在线观看视频| 88av欧美| 在线十欧美十亚洲十日本专区| 麻豆国产av国片精品| 一a级毛片在线观看| 国产伦精品一区二区三区四那| 一个人观看的视频www高清免费观看| 中文字幕av成人在线电影| 色尼玛亚洲综合影院| 中文字幕熟女人妻在线| 亚洲国产高清在线一区二区三| 国产成人欧美在线观看| 天堂√8在线中文| 老熟妇乱子伦视频在线观看| 三级毛片av免费| 亚洲激情在线av| 欧美+日韩+精品| 麻豆国产av国片精品| av欧美777| 精品久久久久久久人妻蜜臀av| 一区二区三区免费毛片| 2021天堂中文幕一二区在线观| 亚洲人成网站高清观看| 久久6这里有精品| 好男人电影高清在线观看| 精品人妻一区二区三区麻豆 | 亚洲国产欧美网| 桃红色精品国产亚洲av| 97人妻精品一区二区三区麻豆| 国产精品电影一区二区三区| 亚洲精华国产精华精| 免费看a级黄色片| 欧美一级a爱片免费观看看| 亚洲aⅴ乱码一区二区在线播放| 亚洲专区中文字幕在线| 精品人妻一区二区三区麻豆 | 亚洲成av人片在线播放无| 国产精华一区二区三区| 身体一侧抽搐| 午夜福利成人在线免费观看| 久久久久久大精品| 亚洲电影在线观看av| 国产精品久久久久久久久免 | 免费搜索国产男女视频| 精品人妻1区二区| 午夜精品一区二区三区免费看| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美国产一区二区入口| 国产精品一区二区免费欧美| 丝袜美腿在线中文| 女人高潮潮喷娇喘18禁视频| 国产色婷婷99| 免费搜索国产男女视频| 最好的美女福利视频网| 99热精品在线国产| 欧美成人一区二区免费高清观看| 精品久久久久久,| 99久久精品一区二区三区| 午夜福利在线观看吧| 亚洲国产精品sss在线观看| 国产激情偷乱视频一区二区| 高清毛片免费观看视频网站| 国产一区在线观看成人免费| 99久久99久久久精品蜜桃| 亚洲欧美日韩高清专用| 精品国产三级普通话版| 好看av亚洲va欧美ⅴa在| 99精品久久久久人妻精品| 高清日韩中文字幕在线| 国内精品一区二区在线观看|