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

    Influence of cascade reservoirs on spatiotemporal variations of hydrogeochemistry in Jinsha River

    2021-09-09 07:55:30QiushengYuanPeifangWangJuanChenChaoWangShengLiuXunWang
    Water Science and Engineering 2021年2期

    Qiu-sheng Yuan,Pei-fang Wang,Juan Chen,Chao Wang,Sheng Liu,Xun Wang

    Key Laboratory of Integrated Regulation and Resources Development of Shallow Lakes,Ministry of Education,Hohai University,Nanjing 210098,China

    Abstract River hydrogeochemistry offers necessary guidance for effective water environmental management.However,the influence of cascade reservoirs on river hydrogeochemistry remains unknown.In this study,the Jinsha River,the headwaters of the Yangtze River of China,was selected to investigate the spatiotemporal variations of hydrogeochemistry after the construction of six cascade reservoirs.Major ions,total dissolved solids,electrical conductivity,and pH values of sampled water in the upper natural reaches and lower reservoir-regulated reaches were analyzed in both flood and dry seasons.The results of Piper diagram and Gibbs plots showed that the hydrogeochemistry of the Jinsha River was naturally controlled by both evaporation-crystallization and carbonate weathering processes,but it was also artificially affected by reservoirs.The impoundment of cascade reservoirs affected the hydrodynamic condition of the river.The river flow in the flood season was reduced by approximately 24.5%,altering the proportions of water sources and leading to notable hydrogeochemical alterations in reservoir-regulated reaches.Conversely,river hydrogeochemistry generally remained unchanged in the dry season,owing to the insignificant effect of cascade reservoirs on river flow.In contrast to what has been observed in previous studies of individual reservoirs,the cumulative influence of cascade reservoirs on the Jinsha River flow regime did not cause abrupt hydrogeochemical changes between the upstream and downstream areas of each reservoir.Moreover,the water quality assessments revealed that the impoundment of cascade reservoirs improved downstream irrigational water quality,with lower Na+ratio values in the flood season.This study provides the earliest evaluation of cascade reservoir influence on the hydrogeochemistry of the Jinsha River.The findings of this study can be used as a reference for scientific guidelines for future environmental management of cascade reservoirs in large rivers.

    Keywords:River hydrogeochemistry;Spatiotemporal variations;Controlling factors;Cascade reservoirs;Jinsha River

    1.Introduction

    Large rivers play a significant and irreplaceable role in the global water cycle and geochemical cycle involving continents and oceans(Immerzeel et al.,2013).Water quality of large rivers is a crucial and sensitive issue worldwide because it is always affected by a variety of natural and artificial factors(Jin et al.,2015,2020;Pant et al.,2018).In general,river hydrogeochemistry is primarily controlled by the evaporationcrystallization,atmospheric precipitation,and rock weathering processes(Gibbs,1970;Martin and Meybeck,1979),and it is secondarily influenced by surface runoff,groundwater discharge,and anthropogenic disturbance(Qu et al.,2019;Wang et al.,2018;Wu et al.,2008).Study of river hydrogeochemistry provides important guidance for water quality management,because it not only enables comprehensive understanding of the composition and characteristics of riverine chemicals,but also offers necessary indicators for aquatic environment protection(Pant et al.,2018).Due to the importance of river hydrogeochemical attributes in water quality management,the chemical signatures and relevant controlling factors of many large rivers worldwide,including the Mekong River(Huang et al.,2009),Yangtze River(Chen et al.,2002),and Amazon River(Stallard and Edmond,1987),have been thoroughly studied.These studies have provided important information on the sources,patterns,and controlling mechanisms of riverine chemicals in global water cycling(Meybeck and Ragu,2012;Viers et al.,2009).

    Over the past half century,many cascade reservoirs have been constructed in large rivers for various purposes,including flood control,power generation,navigation,drinking water,and irrigation(Winemiller et al.,2016).However,cascade reservoirs have simultaneously transformed lotic river systems into serial,discontinuous,and slow-flow reservoir systems(Shi et al.,2017),resulting in considerable modification of geomorphology,temperature,and flow regimes,as well as spontaneous cycles of riverine chemicals(Ellis and Jones,2013;Chao et al.,2016).It has been reported that individual reservoirs constructed in large rivers have the potential to significantly alter river hydrogeochemistry,owing to the enrichment of anthropogenic contamination within reservoirs or the erosion of released water in downstream areas of reservoirs(Kelly,2010;Gao et al.,2013;Wang et al.,2018).For example,Wang et al.(2018)found that after the impoundment of the Three Gorges Reservoir in the Yangtze River of China,the concentrations of anions(e.g.,SO24-,Cl-,and NO-3)and cations(e.g.,Ca2+,Mg2+,and Na+)sharply increased because of the modification of hydrological and geochemical processes as well as the mixing patterns of different water sources.In addition,Gao et al.(2013)proposed that the impoundment of the Wujiangdu Reservoir could significantly change riverine chemical weathering,which in turn could modify the river hydrogeochemical composition downstream.Recent studies have only focused on the influence of individual reservoirs.However,the impact of cascade reservoirs is still unknown.The repetitive and accumulative effects of cascade reservoirs on river flow,sediment transport,biogenic substance,and plankton may significantly alter river hydrogeochemistry through the dam blockage,intermittent discharge,and reservoir backwater processes,which differ from the conditions in individual reservoirs(Chen et al.,2020;Huang et al.,2018;Winemiller et al.,2016).It has been demonstrated that cascade reservoirs have cumulative impacts on river hydrodynamic process and flow regime,and their effect is notably greater than those related to individual reservoirs(Huang et al.,2018).For example,Huang et al.(2010)concluded that the cumulative influence of multiple cascade reservoirs on the upper Yangtze River flow regime is stronger than the total impacts of all reservoirs under independent operation.Therefore,it can be hypothesized that cascade reservoirs may produce cumulative effects on river hydrogeochemistry as they cumulatively impact river flow regime.

    The Jinsha River is a large river in China that originates on the Qinghai-Tibet Plateau and flows through the arid and semiarid alpine regions.This river is regulated by typical cascade reservoirs(Li et al.,2013).As the headwaters of the Yangtze River,the Jinsha River serves as a vital resource for the downstream drinking water and agricultural irrigation.Owing to the social and economic importance of this river,there is immense concern about its water quality(Yuan et al.,2019).It is therefore essential to investigate the hydrogeochemistry of the Jinsha River,because of its important role in sustaining human livelihood and maintaining ecological balance.Previous studies have investigated the specific signatures of natural and anthropogenic disturbances to the hydrogeochemistry of the Jinsha River(Qu et al.,2015;Wu et al.,2008).However,no published works have focused on the spatial variations of river hydrogeochemistry as influenced by six cascade reservoirs successively put into operation between 2010 and 2014(Li et al.,2018).Recently,the seasonal regime of Jinsha River flow has been significantly affected by the cascade reservoirs,exhibiting a reduced flood-season river flow and an increased dry-season river flow(Huang et al.,2018).However,the way in which the six cascade reservoirs affect river hydrogeochemistry in different seasons remains poorly understood.

    This study hypothesized that(1)the operation of cascade reservoirs affects the spatial and seasonal variations of river hydrogeochemistry,and(2)water quality is significantly influenced by cascade reservoirs.To test these hypotheses,the hydrogeochemical composition of the Jinsha River was investigated using water samples from the upper natural reaches and lower reservoir-regulated reaches in both flood and dry seasons.Furthermore,the quality of drinking water was assessed using the World Health Organization guidelines(WHO,2011),and the quality of irrigation water was assessed through the percent sodium(P)and sodium adsorption ratio(R)(Thomas et al.,2014).This study intended to provide an integrated understanding of the influence of cascade reservoirs on the spatiotemporal variations of hydrogeochemistry in the Jinsha River.The findings of this study are expected to contribute to future environmental management of cascade reservoirs in large rivers.

    2.Materials and methods

    2.1.Study area and field sampling

    Fig.1.Sampling sites in natural reaches(N1 through N12)and reservoir-regulated reaches(LY-U/D,AH-U/D,JAQ-U/D,LKK-U/D,LDL-U/D,and GYY-U/D)and lithologic map of Jinsha River Basin extracted from Wu et al.(2008).

    The Jinsha River is 2 290 km long,with an average slope of 1.44 m/km and a drainage area of 362 000 km2(Fig.1).The river originates from the Tongtian River,which is located in an arid and cold area with intense evaporation and an annual precipitation below 500 mm.However,the lower Jinsha River region falls within the India Monsoon climate zone and is wetter than the upstream region.The annual precipitation varies significantly with altitude.The upper Jinsha River from the Benzilan to Shigu hydrological stations has an annual precipitation of approximately 600 mm.By contrast,the annual precipitation in the middle and lower reaches from the Huiyi to Yibin stations is 900 to 1 300 mm(Wu et al.,2008).Precipitation and runoff are mainly concentrated from May to October in the flood season,a period that accounts for over 90%of the annual total amount(data source from the Yangtze River Water Resources Commission,http://www.cjw.gov.cn).The lithology of the Jinsha River Basin is complicated(Fig.1).The headwaters arise from the Tanggula and Bayan Har fold belts,which are made up of clastic rocks(ClR)and carbonate rocks(CaR),as well as quaternary fluvial deposits(QFD)entrained with volcanic rocks and evaporite rocks.In the upper reaches,CaR,ClR,ophiolitic m′elanges(OpM),and low-grade metamorphic rocks(LMR)are widely distributed.Moreover,CaR,ClR,QFD,continental basalts,high-grade metamorphic rocks(HMR),and granitoids(Gra)are commonly found in the middle and lower river basin(Wu et al.,2008).Over the past decade,six mega cascade reservoirs were constructed and put into operation in the middle Jinsha River reaches,including the Liyuan(LY),Ahai(AH),Jinanqiao(JAQ),Longkaikou(LKK),Ludila(LDL),and Guanyinyan(GYY)reservoirs(Fig.1).Detailed information on these reservoirs is presented in the supplementary data.In this study,the upstream reaches of the Jinsha River were defined as natural reaches because they maintained their original appearance with negligible influence from human activities.Given that the downstream river was greatly influenced by the cascade reservoirs,it was defined as a regulated reach.The average discharges at the Shigu(upstream of all six reservoirs)and Panzhihua(downstream of all six reservoirs)hydrological stations before(in 2007 and 2008)and after(in 2015 and 2016)the operation of cascade reservoirs are shown in Fig.2.

    Considering the potential influence of climate change,water samples were collected in the flood season(May)and the dry season(December)of 2017.All samples were collected in a week without rainfall effects to ensure relatively consistent aquatic conditions.Moreover,most of the sampling areas were in the high mountains,with steep gorges,where various uncertain geological disasters usually occur in the middle and late flood season.Thus,the sampling time was selected to occur in the early flood season(May).The selection of sampling locations was based on the hydro-climatic variations,in terms of altitude,distance,topography,and reservoir sites without the influence of water confluence due to large tributaries.Of the total 24 total sampling sites,12 were in the natural reaches,and 12 were situated across six reservoirs(Fig.1).Detailed information on the locations of these sampling sites is available in the supplementary data.All water samples were collected in the middle of the river channel with 2-L wide mouth polyethylene bottles.The sampling equipment was rinsed with river water twice before sampling.Samples were taken against the river flow direction at a depth of approximately 50 cm.

    2.2.Water physicochemical analysis

    Water pH and electrical conductivity(Ec)were determined in situ with a calibrated multiparameter water quality analyzer(HQ40d,Hach Co.,USA).Water samples were subsequently filtrated through 0.2-μm millipore nitrocellulose filters(Shanghai Xin Ya Purification Equipment Co.,Ltd.,China).The samples for cation analysis were acidified to pH<2.0 using high-purity HNO3.These filtrated water samples were stored at 4°C for chemical analysis.

    Concentrations of major cations,such as Na+,K+,Ca2+,and Mg2+,and dissolved silicon(Si),were measured with an inductively coupled plasma-optical emission spectrometer(ICP-OES,Prodigy,USA).The detection limits of ICP-OES were 0.01μg/L for Si and 0.05μg/L for other elements.Concentrations of major anions,such as SO24-,Cl-,HCO-3,and NO-3,were detected by an ion chromatography system(Dionex ICS-900,Thermo Scientific,USA)with a detection limit of 0.1μg/L.The analysis followed standard operating procedures,with blank reagent as well as duplicate samples,and the analytical error of all major ions was no higher than 3%.For quality control,standard samples of known concentrations were periodically analyzed.The estimates of total dissolved solid concentration(ρ(TDS))were obtained using the following formula(Wu et al.,2008):

    2.3.Water quality assessment

    Riverine water quality for drinking purposes was assessed by comparing the values of major water quality index values with the standard values from the World Health Organization(WHO,2011).In addition,P and R were used to evaluate the irrigation water quality(Thomas et al.,2014).The estimates of P and R were mathematically expressed as follows:

    Fig.2.Average discharges at Shigu and Panzhihua hydrological stations and discharge differences before and after operation of cascade reservoirs(data from the Bureau of Hydrology,Ministry of Water Resources of the People's Republic of China).

    Detailed information on the standards of irrigation water quality is presented in the supplementary data.

    2.4.Data analysis

    Descriptive statistical analysis was performed to evaluate the spatiotemporal variation of each variable.The variations in different seasons were compared through paired t-tests for the upper natural reaches,the lower reservoir-regulated reaches,and the entire river.Furthermore,Piper(1944)and Gibbs(1970)plots as well as ionic relationship analysis(Qu et al.,2019)were used to clarify the controlling mechanisms of river hydrogeochemistry.

    3.Results and discussion

    3.1.General hydrogeochemistry of entire Jinsha River

    Table 1 presents the geochemical composition of the Jinsha River in flood and dry seasons.Theρ(TDS)value for the entire Jinsha River water samples(501 mg/L)was five times higher than the average of global rivers,owing to the influence of the arid/semiarid alpine environment in southwestern China,and was similar to the conditions in the upper Mekong River and Yarlung Tsangpo River,as shown in Table 2(Meybeck and Ragu,2012;Qu et al.,2015,2017).

    The average concentrations of major cations throughout the Jinsha River could be ranked in the following descending order:ρ(Na+)>ρ(Ca2+)>ρ(Mg2+)>ρ(K+).Na+accounted for 48.79%of the total cation budget of the Jinsha River,followed by Ca2+and Mg2+,accounting for 35.60%and 14.20%,respectively(Table 1).The content of Na+for the entire Jinsha River water was over ten times higher than the global average(Meybeck and Ragu,2012).This is attributed to the predominant influence of intense evaporation at the headwaters(Wu et al.,2008;Qu et al.,2015).Furthermore,the concentrations of Ca2+and Mg2+were three and five times higher than the global averages,respectively.This implies that carbonate weathering was another factor that controlled the variation of hydrogeochemistry along the Jinsha River(Meybeck,2003).

    The average concentrations of major anions followed the pattern ofρ(HCO-3)>ρ(Cl-)>ρ(SO2-4)>ρ(NO-3)(Table 1).These results are consistent with previous findings,and it was determined that intense rock weathering and evaporite dissolution were responsible for the relatively higher concentrations of HCO-3,Cl-,and SO2-4than those of global averages(Meybeck,2003).The silicate weathering was relatively less intense in the Jinsha River,and the average dissolved silica concentration(3.41 mg/L)was lower than the average of global rivers(Meybeck,2003).Compared to the aforementioned geochemical processes,silicate weathering was less pronounced(Pant et al.,2018).These findings are similar to those of previous studies on the Jinsha River and surroundingrivers in southwestern China(Wu et al.,2008;Huang et al.,2009;Qu et al.,2015,2017).

    Table 1 Geochemical composition of Jinsha River in flood and dry seasons.

    Table 2 Geochemical composition of Jinsha River and surrounding large rivers in China.

    3.2.Influence of natural processes on river hydrogeochemistry

    3.2.1.pH,total dissolved solids,and electrical conductivity

    The Jinsha River Basin is carbonate-dominated,leading to a mildly alkaline river water condition with a pH value ranging from 7.85 to 9.43(Fig.3 and supplementary data).Due to the effect of surface runoff,the pH value increases along the Jinsha River(Wu et al.,2008;Qu et al.,2015).The paired ttest between the pH values in different seasons revealed that the pH value in the dry season was 0.28 higher than that in the flood season(p<0.05)(see supplementary data).This was probably due to the strong riverine microbial respiration(high air temperature)in summer and the lesser influence of acid rain(low precipitation)in winter(see supplementary data).

    The upper natural reaches of the Jinsha River were characterized by notably higherρ(TDS)and Ecvalues and greater spatial variation than in the lower reservoir-regulated reaches.The highestρ(TDS)and Ecvalues were observed at the headwaters of the Jinsha River(sample N1),and the lowest values were recorded at the AH reservoir(sample AH-D)(Fig.3 and supplementary data).The higherρ(TDS)and Ecvalues in the upper Jinsha River were due to intense evaporation that significantly exceeds precipitation.They were also caused by excessive hydraulic erosion and chemical weathering in the arid/semiarid alpine environment.Moreover,the substantially increased river flow and precipitation along the Jinsha River Basin played a pivotal role in diluting the ionic concentration(Fig.2 and supplementary data).Simultaneously,the degree of hydraulic erosion and chemical weathering was reduced,owing to the increased river flow and widened river channel.The relatively low concentration of suspended solids in downstream areas indirectly reduced the CaR dissolution(Chen et al.,2002).Thus,theρ(TDS)value tended to decrease along the Jinsha River(Fig.3).However,theρ(TDS)and Ecvalues of sample N2 in the dry season were abnormally lower than those of surrounding sampling sites.Glaciers were observed in most areas of the upstream region,located in a valley with an elevation higher than 3 000 m(see supplementary data).As reported in a previous study on hydrogeochemistry of the Gandaki River,the dilution of glacier meltwater may potentially decrease ionic concentrations(Pant et al.,2018).

    The seasonal variations ofρ(TDS)and Ecbetween the upper natural reaches and lower regulated reaches of the Jinsha River were significantly different(p<0.05)(see supplementary data).In the flood season,higherρ(TDS)and Ecvalues in the upper natural reaches were closely correlated with the strong evaporation-crystallization processes(high air temperature in May)(see supplementary data).The significantly increased river flow and precipitation in the flood season may decreaseρ(TDS)and Ecthrough the dilution effect(see supplementary data),but they played a less important role than the effect of evaporation processes.However,these natural phenomena suddenly disappeared in the lower reaches because of the effect of cascade reservoirs.This issue is discussed in section 3.3.

    3.2.2.Concentrations of major ions

    Concentrations of major ions exhibited a decreasing trend along the Jinsha River except those of Si and NO-3.This is consistent with previous findings on the spatial variations of ρ(TDS)and Ec.The Si and NO-3concentrations were relatively lower than those of other chemicals,and they were fairly consistent along the entire river(Fig.3 and supplementary data).Several elevated NO-3concentrations were recorded in the lower reservoir-regulated reaches,such as the downstream area of the LKK Reservoir and upstream area of the GYY Reservoir.This was due to the anthropogenic nitrogenous fertilizer contamination from riparian agriculture drainage in the dry-hot valley of the Jinsha River(Pant et al.,2018;Zhang et al.,2006).Most ionic concentrations in the lower Jinsha River were less than half of those in the upper reaches,whereas the river flow was almost doubled downstream(Figs.2 and 3).Due to the rain shadow climate and intense evaporate dissolution caused by high solar radiation in the upper arid/semiarid alpine regions,the Na+and Clconcentrations in the upstream region were much higher than those in the downstream region(Qu et al.,2015,2019).In addition,water recharge from the surrounding saline lakes and hot springs also contributed to the elevated Na+and Clconcentrations(Qu et al.,2019).

    In the flood season,the water ionic concentration changed from the patterns ofρ(Na+)>ρ(Ca2+)>ρ(Mg2+)>ρ(K+)andρ(Cl-)>ρ(HCO-3)>ρ(SO2-4)>ρ(NO-3)in the upstream region to the patterns ofρ(Ca2+)>ρ(Na+)>ρ(Mg2+)>ρ(K+)andρ(HCO-3)>ρ(SO2-4)>ρ(Cl-)>ρ(NO-3)in the downstream region.However,the ionic abundance of river water in the dry season followed the unchanged patterns of ρ(Na+)>ρ(Ca2+)>ρ(Mg2+)>ρ(K+)andρ(HCO-3)>along the entire Jinsha River(Fig.3 and supplementary data).As mentioned earlier,the Jinsha River Basin is CaR-dominated,and therefore,the concentration of HCO-3was the highest among all anions(Wu et al.,2008).Although river water with elevated Cl-concentrations discharged from the headwaters,surface runoff from the CaR-dominated lithology increased the downstream HCO-3concentrations.The coupled effect of high precipitation and intense surface runoff resulted in the high river flow in the flood season.This was responsible for the lower concentrations of Na+,K+,Cl-,,andand the higher concentrations of Ca2+and HCO-3along the Jinsha River in the flood season than in the dry season.

    The spatial and seasonal variations of major ions demonstrated that the hydrogeochemistry of the Jinsha River was mainly controlled by natural processes according to different climate,geology,and hydrology.However,the variation patterns of these hydrogeochemical processes in upper natural reaches were not similar to those in lower reservoir-regulated reaches,where river flow regime was influenced after the operation of cascade reservoirs(Huang et al.,2018).

    3.3.Impact of cascade reservoirs on river hydrogeochemistry

    3.3.1.Comparison of natural and reservoir-regulated conditions

    Fig.3.Spatiotemporal variations of hydrogeochemistry in upper natural reaches(NR)and lower reservoir-regulated reaches(RR)of Jinsha River.

    Fig.4 shows the hydrogeochemistry at the Shigu(upstream area of all six reservoirs)and Panzhihua(downstream area of all six reservoirs)hydrological stations in flood and dry seasons for the time periods before(in 2005)and after the operation of cascade reservoirs(in 2017).In the figure,F-2005 and D-2005 represent the flood and dry seasons in 2005,respectively,and F-2017 and D-2017 represent the flood and dry seasons in 2017,respectively.Data for 2005 were from Wu et al.(2008)and data for 2017 were from this study.The change rate of ionic concentrations between Shigu and Panzhihua hydrological stations was used to investigate the influence of cascade reservoirs on river hydrogeochemistry.

    After the construction of cascade reservoirs,the concentration change rates of Na+,K+,and Cl-were largely altered in the flood season,from-41.4%to-18.9%in 2005 to-85.8%to-72.7%in 2017,whereas the change rates for these chemicals in the dry season increased from-63.6%to-23.2%in 2005 to 16.5%to 26.7%in 2017(Fig.4).This phenomenon shows a remarkable difference from the conditions before the construction of cascade reservoirs as reported in Wu et al.(2008).Furthermore,the change rates of HCO-3concentration increased in the flood season and decreased in the dry season under the cascade reservoir regulations(Fig.4).In addition to the effects of reservoirs,human activities also induced the changes of ionic concentrations.Abnormally high concentrations ofin the Panzhihua region were largely ascribed to the effect of coal consumption in Sichuan Province(Chen et al.,2002;Wu et al.,2008),where heavy metal pollution is also severe,as reported in previous studies(Yuan et al.,2019,2021).Meanwhile,the concentrations of NO-3were notably increased in both the Shigu and Panzhihua areas after the construction of reservoirs(from 0.20 to 0.90 mg/L in 2005 to 1.47 to 2.28 mg/L in 2017)because of the growing impact of agricultural nitrogenous fertilizers(Zhang et al.,2006).The population density in the lower reservoirregulated reaches is over 100 per km2higher than in the upper natural reaches,and the Sichuan Basin is a large ricefarming center(Wu et al.,2008).Therefore,the change rates of SO42-and NO-3cannot accurately reflect the impacts of cascade reservoirs,and they were much more significantly influenced by other anthropogenic activities.

    Fig.4.Hydrogeochemistry of ionic concentrations at Shigu and Panzhihua hydrological stations and change rates of ionic concentrations between these two stations in flood and dry seasons before(in 2005)and after(in 2017)operation of cascade reservoirs.

    Cascade reservoirs limit the natural transfer of river water,decrease the longitudinal velocity of water flow,and thus facilitate seasonal variations of river flow(Huang et al.,2018;Li et al.,2018;Chao et al.,2016;Wang et al.,2018).Due to the construction of cascade reservoirs along the Jinsha River,the monthly river flow regimes also significantly changed,causing reduced river flow in the flood season and increased river flow in the dry season(Huang et al.,2018).Before the operation of cascade reservoirs(2007-2008),the variation of the Jinsha River flow ranged from(800.48±60.05)m3/s at the Shigu hydrological station to(1 223.14±139)m3/s at the Panzhihua hydrological station in May.However,after the operation of cascade reservoirs(2015-2016),the variation of the river flux changed from(785.33±23.14)m3/s at the Shigu hydrological station to(849.36±182.12)m3/s at the Panzhihua hydrological station.The river flow at the Shigu hydrological station showed no significant alteration.By contrast,the river flow at the Panzhihua hydrological station was notably reduced by approximately 24.5%in the flood season,as calculated based on the average discharge data in May before and after the operation of reservoirs(Fig.2).This further altered river hydrogeochemistry.In the flood season,river water supplied from the upper reaches and additionally added water from tributaries converged and was partially trapped within reservoirs.Therefore,the concentrations of Na+and Cl-were abruptly diluted to one third of their original contents. However, the HCO-3concentration increased,which may have been due to the effect of surface runoff from the CaR-dominated lithology(Wang et al.,2018).Conversely,in the dry season,river hydrogeochemistry was more significantly influenced by the upstream water,owing to the slight river flow increase of approximately 1.1%,regulated by the cascade reservoirs(Fig.2).Therefore,major ionic compositions were relatively consistent from the upper natural reaches to the lower reservoir-regulated reaches.In addition,the arid and cold climate drastically reduced the amount of water supplied from tributaries along the Jinsha River Basin in the dry season(Wu et al.,2008).

    3.3.2.Comparison between upstream and downstream regions of each reservoir

    The major ionic concentrations showed no abrupt changes between the upstream and downstream regions of each reservoir along the Jinsha River,except K+in the LDL Reservoir and NO-3in the GYY Reservoir(Fig.3 and supplementary data).This is an interesting phenomenon,different from the hydrogeochemical variations associated with individual reservoirs.Generally,if the river water has been notably polluted by anthropogenic activities(e.g.,industrial and agricultural inputs),the impoundment of individual reservoirs may increase ionic concentrations,for example,the hydrogeochemistry within the Three Gorges Reservoir on the Yangtze River(Wang et al.,2018).On the other hand,due to the release of hungry/clean water from individual reservoirs,the intense weathering of carbonate or silicate may increase the downstream ionic concentrations,for example,the downstream hydrogeochemistry of the Wujiangdu Reservoir on the Wujiang River(Gao et al.,2013).

    However,the hydrogeochemistry along cascade reservoirs reported in this study contrasts with what has been reported for individual reservoirs.No significant change of major ionic concentrations between the upstream and downstream areas of each reservoir occurred(Fig.3),and only K+and NO-3had elevated concentrations in some reservoirs(see supplementary data),likely due to the occasional agricultural pollution of nitrogen and potassium fertilizers(Zhang et al.,2006).Therefore,it can be surmised that the influence of cascade reservoirs on river hydrogeochemistry has a different mechanism in comparison with that of individual reservoirs.The cascade reservoirs transformed a lotic ecosystem into a serial,discontinuous,and slow-flow ecosystem(Li et al.,2018;Huang et al.,2018).The water storage regions of these reservoirs were connected end to end,which limited the downstream hydraulic erosion of each reservoir.Compared to individual reservoirs,the post-dam rock weathering along the cascade reservoirs may be relatively mild,with no clearly elevated concentration of chemical materials.This unexpected cumulative effect of cascade reservoirs on the hydrogeochemistry of the Jinsha River warrants further investigations with long-term field observations.

    3.4.Hydrogeochemical facies and their controlling factors

    3.4.1.Hydrogeochemical facies in Jinsha River

    The influence of cascade reservoirs on the hydrogeochemical facies of the Jinsha River was evaluated through the ionic milli-equivalent percentages in the Piper diagram(Fig.5(a)).There are six distinct types:(1)Ca-HCO3,(2)Na-Cl,(3)mixed Ca-Mg-Cl,(4)mixed Ca-Na-HCO3,(5)Ca-Cl,and(6)Na-HCO3(Khadka and Ramanathan,2013).In the flood season,the sample clusters of the upper natural reaches were close to the Na++K+apex and Clapex.By contrast,water samples in the lower reservoirregulated reaches abruptly shifted to the Ca2+apex and HCO-3apex.Water in the upper Jinsha River was classified as the Na-Cl and mixed Ca-Mg-Cl types,whereas all water from the lower reservoir-regulated reaches was classified as the Ca-HCO3type.In contrast,all the samples were clustered together with consistent ionic equivalent percentage patterns in the dry season.Except for the headwaters(the Na-Cl type),the hydrogeochemical facies of the entire river were regulated by the mixed Ca-Mg-Cl type(Fig.5(b)and supplementary data).

    The spatial heterogeneity of ionic equivalent percentages along the Jinsha River demonstrated that the hydrogeochemical facies were primarily shaped by natural processes,as discussed in section 3.2,including upstream evaporite dissolution and carbonate weathering across the entire river basin.Furthermore,as the hydrogeochemical facies between natural and reservoir-regulated reaches tended to diverge in the flood season but converge in the dry season,the seasonal differences were also visible(Fig.5).This was the consequence of cascade reservoir regulations.

    3.4.2.Factors controlling river hydrogeochemistry

    Fig.5.Piper diagram characterizing river hydrogeochemistry in upper natural reaches and lower reservoir-regulated reaches of Jinsha River in flood and dry seasons.

    The mechanisms that control river hydrogeochemistry include natural processes (e.g., groundwater leaching,evaporation-crystallization,atmospheric precipitation,and rock weathering)and anthropogenic activities(Gibbs,1970;Martin and Meybeck,1979;Viers et al.,2009).The natural mechanisms can be inferred from the Gibbs diagram and identified by three end-members(Fig.6):(1)precipitation,with aρ(TDS)less than 10 mg/L and aρ(Na+)/[(ρ(Na+)+ρ(Ca2+)]weight ratio of 0.5-1.0;(2)rock weathering,with aρ(TDS)of 70-300 mg/L and aρ(Na+)/[(ρ(Na+)+ρ(Ca2+)]weight ratio of less than 0.5;and(3)evaporation,with aρ(TDS)higher than 300 mg/L and aρ(Na+)/[(ρ(Na+)+ρ(Ca2+)]weight ratio of 0.5-1.0(Gibbs,1970;Qu et al.,2015).In Fig.6,all sampling points were bounded within dotted lines to reveal the mechanisms that control river hydrogeochemistry.

    The Gibbs diagram illustrates that evaporation was the major controlling factor of hydrogeochemistry in the upper Jinsha River,and rock weathering became the dominant ionic source in the lower reservoir-regulated reaches in the flood season.By contrast,the evaporation process played a primary role in affecting the hydrogeochemistry of the Jinsha River in the dry season.Generally,the natural ionic sources are mainly derived from bedrocks rather than atmospheric precipitation.The results of the Gibbs plot in this study agreed with previous observations of the Jinsha River and surrounding rivers(Wu et al.,2008;Qu et al.,2015).However,the operation of cascade reservoirs affected river hydrogeochemistry in the lower reaches,especially in the flood season.This is a new finding that has not been reported in previous studies.

    3.4.3.Chemical weathering

    The ionic equivalent ratios of various hydrogeochemical attributes were used to further determine the way chemical weathering of bedrock affected the hydrogeochemistry of the Jinsha River(supplementary data).The average molar ratio of ρ(Na+)/ρ(Cl-)at the headwaters of the Jinsha River(approximately 1.14)was notably higher than that in marine water(supplementary data;Meybeck,2003).This demonstrates that the ionic contributions of precipitation are negligible in rivers within arid and semiarid alpine regions(Singh et al.,2014).The dissolved load in the headwaters mainly came from the dissolution of ClR,CaR,and QFD,which were intercalated with volcanic rocks and evaporites(e.g.,halite,soft sulphate minerals,and sulfide oxidation).

    In the flood season,the notably elevated average molar ratios of[ρ(Ca2+)+ρ(Mg2+)]/[ρ(Na+)+ρ(K+)]and ρ(HCO-3)/[ρ(Na+)+ρ(K+)]were approximately 4.04 and 3.13,respectively.This finding confirmed that the lower reservoir-regulated reaches were dominated by the weathering of ClR,CaR,and HMR minerals(supplementary data;Qu et al.,2019).In addition,the proton-producing reaction is essential for carbonate weathering,including sulfide oxidation and carbonation.The relative importance of these two actions can be explained by theρ(HCO-3)/[ρ(SO2-4)+ρ(HCO-3)]ratio(C ratio)(Pant et al.,2018).The mean C ratio of the lower Jinsha River(approximately 0.77)was significantly higher than 0.5(supplementary data),indicating the vital role of CO2dissolution and carbonate in proton-producing processes,as reported in previous studies(Qu et al.,2019;Fearnside and Pueyo,2012).Furthermore,the increased ratio ofρ(Ca2+)/along the Jinsha River confirmed the predominant function of carbonate weathering as discussed elsewhere(Singh et al.,2014).

    Fig.6.Spatiotemporal variations ofρ(Na+)/[(ρ(Na+)+ρ(Ca2+)]as function ofρ(TDS)in Gibbs graph(Gibbs,1970)from natural and regulated reaches of Jinsha River in flood and dry seasons(data in 2005 and 2012 are from Wu et al.(2008)and Qu et al.(2015),respectively).

    All these results demonstrated that the impoundment of cascade reservoirs notably altered the hydrogeochemical composition of the water in the lower Jinsha River,especially in the flood season.The river flow supplied from the upper reaches and supplemented by surrounding tributaries converged and was partially trapped within the reservoirs.This effect altered the proportions of water sources and caused the lower Jinsha River water to be more significantly affected by the process of carbonate weathering than evaporation.

    3.5.Environmental implication:Water quality for drinking and irrigation

    The Jinsha River is an important water source for the Yangtze River and plays a key role in sustaining downstream human livelihoods and maintaining ecological equilibrium(Li et al.,2013).Furthermore,agriculture,a significant income source in southwestern China,is directly affected by the water quality of regional precipitation and rivers(Thomas et al.,2014).In the Jinsha River Basin,there are many economic crops,such as flue-cured tobacco,as well as various fruits,such as mango,banana,pomegranate,and loquat(Zhang et al.,2006).Thus,a comprehensive evaluation of drinking and irrigation water quality is necessary.

    The drinking water quality of the Jinsha River basically fell within the safety thresholds of the WHO guidelines(WHO,2011).However,a few samples exceeded the maximum permissible limits,including sample N1 with regard toρ(TDS)and samples N1,N2,and N3 with regard toρ(Na+)andρ(Cl-)(Fig.3).Due to the dissolution of evaporites and groundwater transport from surrounding springs and saline lakes at the headwaters,all the water samples exceeding the WHO guidelines were collected from the upper arid and semiarid alpine reaches in the flood season.These results are consistent with the observations of other rivers on the Qinghai-Tibet Plateau(Qu et al.,2019).Water quality from the lower reservoir-regulated reaches was well within the WHO guidelines,although the construction of reservoirs had a significant influence on the ionic compositions.

    Water quality for agricultural irrigation depends on the type and concentration of dissolved salts,particularly sodium(Na+)(Thomas et al.,2014).Elevated concentrations of Na+in irrigation water may lead to the displacement of Ca2+and Mg2+,thereby reducing soil permeability.Ultimately,crop yield is impacted,leading to calcium deficiency,deflocculation,and tilth impairment.Therefore,two important indices,P and R,were used to assess the irrigation water quality(Thomas et al.,2014).The calculated values and categories of P and R are presented in Fig.7 and supplementary data.Except for a few samples in the upper arid and semiarid alpine reaches of the Jinsha River,the overall water quality was safe for agricultural irrigation.The P and R values were higher in the upstream region due to the impacts of evaporite dissolution,especially in the flood season.Moreover,after the construction of cascade reservoirs,the P and R values in the lower reaches tended to decrease in the flood season and increase in the dry season(Fig.7).Obviously,the operation of cascade reservoirs reduced the water flow in the Jinsha River in the flood season.The reservoir-trapped water from upper reaches and surrounding tributaries decreased the concentrations of Na+to one third of its original contents,which potentially improved the water quality for agricultural irrigation.

    Fig.7.Irrigation water quality assessments through P and R in flood and dry seasons along Jinsha River.

    In addition to the influence of cascade reservoirs on drinking and irrigation water,the alteration of hydrogeochemistry in the Jinsha River may also affect zooplankton,phytoplankton,and other biological communities(Chen et al.,2020).The present study provided fundamental knowledge for further investigations into the ecological and environmental impacts of cascade reservoirs.

    4.Conclusions

    This study provided an integrated understanding on the influence of cascade reservoirs on the hydrogeochemistry of the Jinsha River.The impoundment of cascade reservoirs affected the hydrodynamic conditions.The river flow in the flood season was reduced by approximately 24.5%.It altered the proportions of water sources and caused notable hydrogeochemical alterations in reservoir-regulated reaches.The concentrations of Na+and Cl-were diluted to one third of their original contents,whereas the increased concentration of HCO-3was attributed to the effects of surface runoff from the CaR-dominated lithology.In contrast to that of individual reservoirs,the cumulative influence of cascade reservoirs on river flow regime did not cause abrupt hydrogeochemical changes between the upstream and downstream areas of each reservoir.The reservoir-trapped water from upper reaches and surrounding tributaries decreased the concentrations of Na+,potentially improving the downstream irrigation water quality in the flood season.All these findings provide important implications for further water environment management of cascade reservoirs in large rivers.

    Declaration of competing interest

    The authors declare no conflicts of interest.

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.wse.2021.06.008.

    永久免费av网站大全| 国产伦人伦偷精品视频| 男女边摸边吃奶| 操美女的视频在线观看| 国产av国产精品国产| 国产精品99久久99久久久不卡 | 国产欧美日韩综合在线一区二区| 日日摸夜夜添夜夜爱| 黄网站色视频无遮挡免费观看| 人妻 亚洲 视频| 国产高清不卡午夜福利| 99热网站在线观看| 无限看片的www在线观看| 免费黄网站久久成人精品| 亚洲美女搞黄在线观看| 欧美日韩视频高清一区二区三区二| 桃花免费在线播放| 国产国语露脸激情在线看| 在线观看一区二区三区激情| 亚洲国产毛片av蜜桃av| 久久精品国产a三级三级三级| 女人高潮潮喷娇喘18禁视频| 亚洲美女搞黄在线观看| 亚洲国产成人一精品久久久| 国产精品久久久久久精品古装| 婷婷成人精品国产| 国产精品久久久久久精品古装| 国产精品99久久99久久久不卡 | 纵有疾风起免费观看全集完整版| 国产精品三级大全| 日韩中文字幕欧美一区二区 | 视频在线观看一区二区三区| 国产精品欧美亚洲77777| 亚洲一区中文字幕在线| 久久久欧美国产精品| 美女大奶头黄色视频| 色综合欧美亚洲国产小说| 9191精品国产免费久久| 国产日韩欧美在线精品| 日韩人妻精品一区2区三区| 爱豆传媒免费全集在线观看| 国产老妇伦熟女老妇高清| 亚洲国产毛片av蜜桃av| av网站免费在线观看视频| 精品一区二区免费观看| 少妇人妻久久综合中文| 国产av一区二区精品久久| 午夜免费观看性视频| 一级爰片在线观看| av在线老鸭窝| 日韩视频在线欧美| 18禁国产床啪视频网站| 国产精品一区二区在线观看99| 精品国产乱码久久久久久小说| av国产精品久久久久影院| 如何舔出高潮| 国产一区二区激情短视频 | 大陆偷拍与自拍| 亚洲专区中文字幕在线 | 大片免费播放器 马上看| 狠狠精品人妻久久久久久综合| 国产一区二区三区综合在线观看| 在线观看免费午夜福利视频| 欧美乱码精品一区二区三区| 国产xxxxx性猛交| 亚洲七黄色美女视频| 狂野欧美激情性xxxx| 亚洲自偷自拍图片 自拍| 天天躁狠狠躁夜夜躁狠狠躁| 成人毛片60女人毛片免费| 久久久久国产一级毛片高清牌| 无限看片的www在线观看| 嫩草影视91久久| 人人妻人人澡人人爽人人夜夜| 精品国产一区二区久久| av不卡在线播放| 国产成人系列免费观看| 久久99一区二区三区| 啦啦啦在线观看免费高清www| 99久国产av精品国产电影| 在线观看一区二区三区激情| 超色免费av| 人人澡人人妻人| 999精品在线视频| 国产在线免费精品| av又黄又爽大尺度在线免费看| 亚洲精品日本国产第一区| 精品一区二区三卡| 超色免费av| 午夜免费鲁丝| 涩涩av久久男人的天堂| 国产女主播在线喷水免费视频网站| 日日撸夜夜添| 七月丁香在线播放| 男女床上黄色一级片免费看| 国产成人午夜福利电影在线观看| 91国产中文字幕| 久久人人爽av亚洲精品天堂| www.自偷自拍.com| 永久免费av网站大全| 夫妻午夜视频| 国产精品蜜桃在线观看| 国产成人精品在线电影| 久久久久视频综合| 啦啦啦啦在线视频资源| 老司机影院毛片| 国产精品亚洲av一区麻豆 | 国产一区二区三区综合在线观看| 少妇人妻久久综合中文| 欧美国产精品va在线观看不卡| 美女中出高潮动态图| 久久久亚洲精品成人影院| 搡老乐熟女国产| 热re99久久精品国产66热6| 日本黄色日本黄色录像| 爱豆传媒免费全集在线观看| 9热在线视频观看99| 国产成人午夜福利电影在线观看| 日韩伦理黄色片| av国产久精品久网站免费入址| 国产精品麻豆人妻色哟哟久久| 满18在线观看网站| 午夜免费鲁丝| 搡老岳熟女国产| 日韩一本色道免费dvd| 青春草亚洲视频在线观看| 亚洲精品中文字幕在线视频| 狠狠婷婷综合久久久久久88av| 大片电影免费在线观看免费| 婷婷成人精品国产| 亚洲欧美一区二区三区黑人| 青青草视频在线视频观看| 免费黄网站久久成人精品| 肉色欧美久久久久久久蜜桃| 久久久久久久久免费视频了| av在线老鸭窝| 免费人妻精品一区二区三区视频| 丰满少妇做爰视频| 宅男免费午夜| 免费在线观看黄色视频的| 天天躁日日躁夜夜躁夜夜| 男人添女人高潮全过程视频| 欧美人与性动交α欧美精品济南到| 午夜日韩欧美国产| 国产成人午夜福利电影在线观看| 日韩一卡2卡3卡4卡2021年| 欧美 亚洲 国产 日韩一| 午夜免费观看性视频| 国产精品一国产av| 天天躁日日躁夜夜躁夜夜| 一区福利在线观看| 亚洲国产av影院在线观看| 欧美国产精品va在线观看不卡| 亚洲精华国产精华液的使用体验| 操美女的视频在线观看| 黄色视频不卡| 精品少妇一区二区三区视频日本电影 | 视频区图区小说| 日日撸夜夜添| 欧美乱码精品一区二区三区| av在线app专区| 国产免费福利视频在线观看| 国产麻豆69| 国产日韩一区二区三区精品不卡| 丝袜在线中文字幕| 国产成人精品无人区| 久久人人爽av亚洲精品天堂| 国产日韩欧美亚洲二区| 国产黄频视频在线观看| 老司机深夜福利视频在线观看 | 无遮挡黄片免费观看| 亚洲欧美精品自产自拍| 色吧在线观看| 在线天堂中文资源库| 交换朋友夫妻互换小说| 国产国语露脸激情在线看| 欧美另类一区| av卡一久久| 免费黄色在线免费观看| 在线观看人妻少妇| 成人亚洲欧美一区二区av| 午夜av观看不卡| 国产免费一区二区三区四区乱码| 精品人妻在线不人妻| 青草久久国产| 咕卡用的链子| 亚洲欧美成人精品一区二区| 中文字幕色久视频| 熟妇人妻不卡中文字幕| 大香蕉久久网| 色婷婷av一区二区三区视频| 性色av一级| 国产熟女欧美一区二区| 久久久久久久精品精品| 成年动漫av网址| 午夜av观看不卡| 一本大道久久a久久精品| 久久毛片免费看一区二区三区| 少妇被粗大猛烈的视频| 最近2019中文字幕mv第一页| 香蕉丝袜av| 国产视频首页在线观看| 免费黄频网站在线观看国产| 中文字幕人妻熟女乱码| 少妇人妻精品综合一区二区| 香蕉丝袜av| 久久热在线av| avwww免费| 精品国产乱码久久久久久小说| 亚洲综合精品二区| 国产一区二区三区综合在线观看| 男女国产视频网站| 国产成人系列免费观看| 日本一区二区免费在线视频| 水蜜桃什么品种好| 欧美日韩一区二区视频在线观看视频在线| 观看美女的网站| 亚洲欧美色中文字幕在线| 久久精品亚洲av国产电影网| 男人爽女人下面视频在线观看| 日本黄色日本黄色录像| av一本久久久久| 女的被弄到高潮叫床怎么办| 亚洲精品日韩在线中文字幕| av国产久精品久网站免费入址| 国产一卡二卡三卡精品 | av不卡在线播放| 捣出白浆h1v1| 亚洲精品乱久久久久久| 热99国产精品久久久久久7| 亚洲综合色网址| 人人妻人人澡人人看| 少妇猛男粗大的猛烈进出视频| 欧美黑人精品巨大| 十分钟在线观看高清视频www| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 精品一区二区三区av网在线观看 | 欧美亚洲 丝袜 人妻 在线| 国语对白做爰xxxⅹ性视频网站| 一本一本久久a久久精品综合妖精| 大码成人一级视频| 亚洲av在线观看美女高潮| 丝袜人妻中文字幕| 午夜激情av网站| 午夜福利在线免费观看网站| 日本爱情动作片www.在线观看| 在线观看免费视频网站a站| 日本wwww免费看| 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 亚洲精品av麻豆狂野| 国产av一区二区精品久久| 桃花免费在线播放| 亚洲av成人不卡在线观看播放网 | 男女下面插进去视频免费观看| 中文字幕av电影在线播放| 欧美黑人精品巨大| 成人国产av品久久久| 亚洲精品日韩在线中文字幕| 丝袜脚勾引网站| 90打野战视频偷拍视频| 久久精品久久精品一区二区三区| 免费在线观看完整版高清| 精品久久蜜臀av无| 亚洲成人免费av在线播放| 制服诱惑二区| 免费黄色在线免费观看| 狠狠婷婷综合久久久久久88av| 国产av国产精品国产| 久久毛片免费看一区二区三区| 亚洲成国产人片在线观看| 丝袜在线中文字幕| 这个男人来自地球电影免费观看 | 久久性视频一级片| 最近的中文字幕免费完整| 99国产综合亚洲精品| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看| 一级毛片 在线播放| 飞空精品影院首页| 男人操女人黄网站| 亚洲三区欧美一区| 汤姆久久久久久久影院中文字幕| 欧美精品亚洲一区二区| 免费观看人在逋| 美女福利国产在线| 在线天堂中文资源库| 久久久久人妻精品一区果冻| 久热这里只有精品99| 少妇的丰满在线观看| 天天躁夜夜躁狠狠躁躁| 99精品久久久久人妻精品| 亚洲熟女毛片儿| 高清黄色对白视频在线免费看| 精品国产国语对白av| 国产一区二区三区综合在线观看| 亚洲欧洲精品一区二区精品久久久 | 伦理电影大哥的女人| 9色porny在线观看| 亚洲精品一二三| 菩萨蛮人人尽说江南好唐韦庄| √禁漫天堂资源中文www| 国产精品秋霞免费鲁丝片| 中文字幕av电影在线播放| 亚洲一区中文字幕在线| 九草在线视频观看| 男女国产视频网站| 在线观看人妻少妇| 别揉我奶头~嗯~啊~动态视频 | 欧美日韩视频精品一区| 国产精品久久久久成人av| 99re6热这里在线精品视频| 精品久久蜜臀av无| 九九爱精品视频在线观看| 久久久久久久久久久久大奶| 中文字幕高清在线视频| 考比视频在线观看| 免费不卡黄色视频| 一级黄片播放器| 精品酒店卫生间| 国产精品免费大片| 九九爱精品视频在线观看| 国产极品天堂在线| av电影中文网址| 亚洲第一青青草原| 99久久精品国产亚洲精品| 秋霞在线观看毛片| 亚洲中文av在线| 亚洲欧洲精品一区二区精品久久久 | 在线观看一区二区三区激情| av福利片在线| 一级黄片播放器| 制服丝袜香蕉在线| 亚洲精品久久午夜乱码| 精品一品国产午夜福利视频| 亚洲天堂av无毛| 欧美精品亚洲一区二区| 美女脱内裤让男人舔精品视频| 91老司机精品| 天天影视国产精品| 日本av免费视频播放| 老司机亚洲免费影院| 亚洲精品美女久久av网站| 夫妻午夜视频| 乱人伦中国视频| 日韩大码丰满熟妇| 又黄又粗又硬又大视频| 亚洲,欧美,日韩| 熟女少妇亚洲综合色aaa.| 亚洲精品国产av蜜桃| 男女边摸边吃奶| 咕卡用的链子| 亚洲在久久综合| 国产精品欧美亚洲77777| 日韩制服丝袜自拍偷拍| 男女边摸边吃奶| 一本大道久久a久久精品| 亚洲精品一二三| 国产成人91sexporn| 亚洲婷婷狠狠爱综合网| 丝袜喷水一区| 街头女战士在线观看网站| 无限看片的www在线观看| 丝袜喷水一区| 久久精品亚洲av国产电影网| 免费观看性生交大片5| 国产深夜福利视频在线观看| 成人黄色视频免费在线看| 啦啦啦视频在线资源免费观看| 无遮挡黄片免费观看| 中文字幕人妻熟女乱码| 一二三四在线观看免费中文在| 国产精品二区激情视频| 久久久久国产一级毛片高清牌| 久久久久久人人人人人| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| 欧美黄色片欧美黄色片| 日韩大片免费观看网站| 建设人人有责人人尽责人人享有的| 成人亚洲欧美一区二区av| 黑人欧美特级aaaaaa片| 欧美在线黄色| 日本一区二区免费在线视频| 亚洲精品在线美女| 国产女主播在线喷水免费视频网站| 精品少妇内射三级| 中文字幕人妻丝袜一区二区 | 国产欧美日韩一区二区三区在线| 男女边摸边吃奶| 色婷婷久久久亚洲欧美| 黑丝袜美女国产一区| 国产精品欧美亚洲77777| 十八禁高潮呻吟视频| 午夜日本视频在线| 午夜日韩欧美国产| 欧美日本中文国产一区发布| 熟女少妇亚洲综合色aaa.| 乱人伦中国视频| 一区在线观看完整版| 麻豆精品久久久久久蜜桃| 男的添女的下面高潮视频| 国产精品女同一区二区软件| 亚洲伊人色综图| 欧美日韩成人在线一区二区| 好男人视频免费观看在线| www.自偷自拍.com| 免费高清在线观看视频在线观看| 亚洲精品国产av成人精品| kizo精华| 欧美在线一区亚洲| 中文字幕高清在线视频| 久久婷婷青草| 国产熟女午夜一区二区三区| 日日爽夜夜爽网站| avwww免费| 少妇人妻 视频| 亚洲精品在线美女| 亚洲中文av在线| 国产男女内射视频| 亚洲欧美精品综合一区二区三区| 国产免费现黄频在线看| av一本久久久久| 少妇精品久久久久久久| 在线观看免费午夜福利视频| 国产精品免费视频内射| 国产精品.久久久| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇一区二区三区视频日本电影 | 亚洲国产精品成人久久小说| 久久久欧美国产精品| 亚洲天堂av无毛| 99久国产av精品国产电影| 午夜福利在线免费观看网站| 99久国产av精品国产电影| 亚洲国产欧美一区二区综合| 亚洲成人一二三区av| 看非洲黑人一级黄片| 日日啪夜夜爽| 国产激情久久老熟女| 久久天堂一区二区三区四区| 亚洲精品av麻豆狂野| 99国产综合亚洲精品| av女优亚洲男人天堂| 超色免费av| 亚洲国产成人一精品久久久| 欧美黑人精品巨大| 欧美日韩亚洲国产一区二区在线观看 | 国产精品一二三区在线看| 欧美激情高清一区二区三区 | 欧美日韩一区二区视频在线观看视频在线| 毛片一级片免费看久久久久| 美女高潮到喷水免费观看| 天美传媒精品一区二区| 亚洲成人国产一区在线观看 | 久久97久久精品| 熟女av电影| 午夜激情久久久久久久| 少妇精品久久久久久久| e午夜精品久久久久久久| 在线亚洲精品国产二区图片欧美| 婷婷色麻豆天堂久久| 免费日韩欧美在线观看| 成人免费观看视频高清| 亚洲精品久久午夜乱码| 免费在线观看完整版高清| √禁漫天堂资源中文www| 午夜影院在线不卡| 女人精品久久久久毛片| 乱人伦中国视频| 人人妻人人澡人人看| 97在线人人人人妻| 国产精品人妻久久久影院| 久久久欧美国产精品| 国产av国产精品国产| av.在线天堂| 女人爽到高潮嗷嗷叫在线视频| 亚洲成色77777| 在线观看免费高清a一片| 宅男免费午夜| 日韩中文字幕视频在线看片| 五月天丁香电影| 男女下面插进去视频免费观看| 亚洲欧美日韩另类电影网站| 日韩不卡一区二区三区视频在线| 国产精品香港三级国产av潘金莲 | 水蜜桃什么品种好| 美女主播在线视频| 国产xxxxx性猛交| 国产精品女同一区二区软件| 人体艺术视频欧美日本| 国产亚洲av高清不卡| 日本一区二区免费在线视频| 男人添女人高潮全过程视频| 日本爱情动作片www.在线观看| 国产有黄有色有爽视频| 日韩一区二区三区影片| 国产精品秋霞免费鲁丝片| 日韩精品免费视频一区二区三区| 精品少妇久久久久久888优播| 亚洲成人av在线免费| 一级a爱视频在线免费观看| 哪个播放器可以免费观看大片| 伊人久久大香线蕉亚洲五| 婷婷色综合大香蕉| 天天躁夜夜躁狠狠躁躁| 亚洲美女视频黄频| 熟妇人妻不卡中文字幕| 国产极品粉嫩免费观看在线| 成年av动漫网址| 人人妻人人澡人人看| 九九爱精品视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产区一区二| 日韩大码丰满熟妇| 久久精品久久久久久久性| 久久精品亚洲av国产电影网| 超碰97精品在线观看| 91aial.com中文字幕在线观看| 国产一级毛片在线| 亚洲天堂av无毛| 久久久久久人人人人人| 国产精品一国产av| 亚洲av综合色区一区| 男人操女人黄网站| 晚上一个人看的免费电影| 最黄视频免费看| 欧美日韩亚洲高清精品| 少妇被粗大猛烈的视频| 伊人久久大香线蕉亚洲五| 99久久精品国产亚洲精品| 国产成人午夜福利电影在线观看| 久久精品国产综合久久久| 亚洲av中文av极速乱| 国产深夜福利视频在线观看| 精品一区在线观看国产| 国产男女超爽视频在线观看| 日韩一区二区三区影片| 最新在线观看一区二区三区 | 亚洲国产日韩一区二区| 男人爽女人下面视频在线观看| 中文字幕人妻熟女乱码| 秋霞伦理黄片| 国产免费一区二区三区四区乱码| 国产精品国产三级专区第一集| 97在线人人人人妻| 亚洲情色 制服丝袜| 51午夜福利影视在线观看| 丝袜在线中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久成人aⅴ小说| 国产成人免费无遮挡视频| 日韩一本色道免费dvd| 自线自在国产av| 国产亚洲最大av| 国产av精品麻豆| 亚洲精品在线美女| 啦啦啦在线观看免费高清www| 在线天堂中文资源库| 水蜜桃什么品种好| 精品一区二区免费观看| 尾随美女入室| 欧美国产精品va在线观看不卡| 美女主播在线视频| 久热这里只有精品99| 久久精品国产综合久久久| 久久精品亚洲av国产电影网| 亚洲第一av免费看| 国产成人免费无遮挡视频| 午夜日本视频在线| 无遮挡黄片免费观看| 久久99热这里只频精品6学生| 纵有疾风起免费观看全集完整版| 亚洲精品国产av蜜桃| 亚洲精品国产色婷婷电影| 国产有黄有色有爽视频| 老司机靠b影院| √禁漫天堂资源中文www| 国产高清不卡午夜福利| 久久久久久免费高清国产稀缺| 自拍欧美九色日韩亚洲蝌蚪91| 女性生殖器流出的白浆| 少妇 在线观看| 女人精品久久久久毛片| 国产一区二区三区综合在线观看| 97人妻天天添夜夜摸| 极品人妻少妇av视频| 熟女av电影| 久久久久网色| 免费在线观看黄色视频的| 免费看不卡的av| 不卡av一区二区三区| 波多野结衣一区麻豆| 久久婷婷青草| 女人爽到高潮嗷嗷叫在线视频| 免费观看a级毛片全部| 久久精品aⅴ一区二区三区四区| 天天躁日日躁夜夜躁夜夜| 日本色播在线视频| 在线观看免费视频网站a站| 欧美日韩国产mv在线观看视频| 亚洲熟女精品中文字幕| 成年动漫av网址| 天天躁日日躁夜夜躁夜夜| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲人成网站在线观看播放| 老司机深夜福利视频在线观看 | 两性夫妻黄色片| h视频一区二区三区| av电影中文网址| 欧美日韩福利视频一区二区| 18禁国产床啪视频网站| 亚洲一区中文字幕在线| av片东京热男人的天堂|