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

    The Linkage between Two Types of El Ni?o Events and Summer Streamflow over the Yellow and Yangtze River Basins

    2020-04-01 02:47:16DanWANGAihuiWANGLianlianXUandXianghuiKONG
    Advances in Atmospheric Sciences 2020年2期

    Dan WANG, Aihui WANG, Lianlian XU, and Xianghui KONG

    1Nansen-Zhu International Research Center, Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China

    2University of Chinese Academy of Sciences, Beijing 100029, China

    ABSTRACT It is generally agreed that El Ni?o can be classified into East Pacific (EP) and Central Pacific (CP) types. Nevertheless,little is known about the relationship between these two types of El Ni?o and land surface climate elements. This study investigates the linkage between EP/CP El Ni?o and summer streamflow over the Yellow and Yangtze River basins and their possible mechanisms. Over the Yellow River basin, the anomalous streamflow always manifests as positive (negative)in EP (CP) years, with a correlation coefficient of 0.39 (-0.37); while over the Yangtze River basin, the anomalous streamflow shows as positive in both EP and CP years, with correlation coefficients of 0.72 and 0.48, respectively.Analyses of the surface hydrological cycle indicate that the streamflow is more influenced by local evapotranspiration (ET)than precipitation over the Yellow River basin, while it is dominantly affected by precipitation over the Yangtze River basin. The different features over these two river basins can be explained by the anomalous atmospheric circulation, which is cyclonic (anticyclonic) north (south) of 30°N over East Asia. EP years are dominated by two anticyclones, which bring strong water vapor convergence and induce more precipitation but less ET, and subsequently increase streamflow and flooding risks. In CP years, especially over the Yellow River basin, two cyclones dominate and lead to water vapor divergence and reduce moisture arriving. Meanwhile, the ET enhances mainly due to local high surface air temperature,which further evaporates water from the soil. As a result, the streamflow decreases, which will then increase the drought risk.

    Key words: summer streamflow, EP El Ni?o, CP El Ni?o, Yellow River basin, Yangtze River basin

    1. Introduction

    China has many rivers that provide rich water resources, but it is also experiencing water shortages and surpluses due to both anthropogenic and natural factors, inducing many hydroclimatological disasters. The rational utilization and management of these water resources is of great concern. In general, river basin streamflow is affected by local human activities and natural climate variabilities. For example, the construction of dams/reservoirs and irrigation systems can alter both the river flow amount and direction,while changes in climate variables will induce positive or negative feedbacks to terrestrial evaporation and runoff processes (Li and Li, 2008; Cong et al., 2009; Gardner, 2009;Huang et al., 2016; Zhou et al., 2018). Over the past several decades, natural climate variability has played a predominant role in the impact of water resources compared to local human activities. Liu et al. (2012) concluded that climate change was the dominant factor in the interannual variations of streamflow over all major river basins in East China over the past 45 years. Li et al. (2016) demonstrated that the heavy rainfall in summer over the Huaihe River valley is mainly influenced by the variation of tropical sea surface temperature.

    El Ni?o-Southern Oscillation (ENSO) is one of the prominent teleconnection patterns in the atmosphere-ocean system (Neelin et al., 1998). It is an essential factor associated with sea surface temperature anomalies (SSTAs) and has significantly influenced land climate variabilities (Qian et al.,2013; Hu et al., 2016; Tian et al., 2016). For instance, Cayan et al. (1999) found that both daily precipitation and streamflow in winter exhibit strong and systematic responses to the two phases of ENSO in the western United States. The effects of El Ni?o on streamflow have also been extensively investigated (Zhang et al., 2007; Ouyang et al., 2014). Zhang et al. (2007) indicated there is a high possibility that convective precipitation will escalate, potentially causing massive flooding in the lower reaches of the Yangtze River basin during El Ni?o years. Furthermore, the possible physical mechanisms involved in the effect of equatorial Pacific SSTAs on the terrestrial climate have been discussed as well. Wang et al. (2000) indicated that the large-scale western North Pacific (WNP) anomalous anticyclone is the key system that has bridged the equatorial central Pacific (CP) warming and weakened the East Asian winter monsoon. Lin and Lu(2009) demonstrated that the ENSO-related WNP anticyclonic anomaly in the lower troposphere may bring about increased precipitation over the Yangtze River basin, and a southward shift of the Asian jet stream in the upper troposphere may lead to decreased precipitation over the Huaihe and Yellow river basins.

    Recent studies have argued that another type of El Ni?o exists, in which the maximum SSTAs appear over the central equatorial Pacific instead of over the eastern equatorial Pacific as with traditional El Ni?o (Ashok et al., 2007; Kug et al., 2009; Ren and Jin, 2011). This type is referred to as CP El Ni?o, and the traditional El Ni?o is referred to as eastern Pacific (EP) El Ni?o. Since the 1970s, CP El Ni?o has occurred more frequently than EP El Ni?o (Ashok et al.,2007). In addition, the ratio of CP El Ni?o to EP El Ni?o events may continue to increase along with global warming,based on climate model projections (Yeh et al., 2009). To date, there have been several studies that have investigated the impacts of these two types of El Ni?o on different climate elements, such as atmospheric convection, global mean temperature, and tropical cyclones (Weng et al., 2009;Wang et al., 2013; Banholzer and Donner, 2014; Yeh et al.,2014). Previous research has also addressed the variations in streamflow under these two types of El Ni?o. In the Mississippi River basin, Liang et al. (2014) showed that EP and CP El Ni?o have opposite impacts on spring soil moisture and the soil water table. Liang et al. (2016) analyzed rivers globally and found that the streamflow in Asian rivers exhibits more complex responses to these different types of El Ni?o, and the regional precipitation is the determining factor of changes in river discharge. Yuan et al. (2016) also pointed out that CP SST shows statistically significant negative influences on the entire Yellow River basin. The above studies indicate that El Ni?o is a pivotal factor that influences the variability of streamflow.

    However, little is still known about the separate responses of streamflow in different river basins of East China to EP and CP El Ni?o, as well as the possible physical mechanisms behind those relationships. To address these issues, analyses of the relationships between the two types of El Ni?o and long-term observed streamflow over the Yellow and Yangtze river basins (the major river basins of East China) were carried out, and the possible physical mechanisms behind the relationships were investigated by analyzing the large-scale atmosphere-ocean circulation patterns during EP and CP El Ni?o years. The results of this study contribute to our understanding of land surface hydrological changes, which occur under different large-scale circulation patterns in the earth system, and therefore to improving the accuracy of land surface hydrological predictions.

    This paper is organized as follows. Section 2 describes the data and methods used in this study. Section 3 analyzes the observed streamflow and El Ni?o indices, and the relationship between summer streamflow over the Yellow and Yangtze river basins and the two types of El Ni?o events. Section 4 discusses the possible mechanisms behind the abovementioned relationships. Finally, conclusions are drawn in section 5.

    2. Data and methods

    The observed daily streamflow dataset used in this study was from the Ministry of Water Resources of China.The data were from 23 hydrological stations, which covered all China’s major river basins (e.g., the Yellow River,Yangtze River; Fig. 1), and the period of coverage was 1980-2008. This dataset was used to investigate the changes in China’s streamflow on both spatial and temporal scales (e.g., Zhang et al., 2011). Monthly observed precipitation and air temperature data were obtained from the CN05.1 dataset, with a 0.5° × 0.5° resolution (Wu and Gao,2013). It was based on daily observations at more than 2400 meteorological stations from the National Meteorological Information Center. Two daily evapotranspiration (ET) datasets were used in this study. One was from NASA’s Collection6 MODIS product, retrieved by the Penman-Monteith algorithm (Mu et al., 2007), hereafter referred to as PMMOD. The other was from the WCRP’s GEWEX Land-Flux project, which was a multi-year, global, merged benchmark synthesis ET product, produced by merging satellite data, in-situ observations, and land-surface datasets over the period 1989-2005 (Mueller et al., 2013), hereafter referred to as LandFlux-EVAL. The latter was used to support the results from the MODIS dataset. Monthly SST data were derived from NOAA’s ERSST.v4 dataset (Huang et al., 2015).In addition, to explore the possible physical mechanisms underlying the relationship between the two types of El Ni?o and streamflow, we obtained atmospheric circulation data, including monthly wind, vertical velocity, geopotential height(GPH), and surface pressure, from the NCEP-2 reanalysis products (Kanamitsu et al., 2002). Monthly specific humidity data were obtained from the NCEP-1 reanalysis products(Kalnay et al., 1996).

    Fig. 1. Locations of the Yellow and Yangtze River basins and the 23 hydrological stations (black spots), among which the black triangle symbols denote the hydrological stations of Huayuankou and Hankou. Blue lines denote the main river network.

    Several Ni?o indices (e.g., Ni?o3, and Ni?o3.4 and Ni?o4) have been defined based on sea surface temperature in different regions of the equatorial Pacific, and widely used in research (e.g., Ashok et al., 2007; Kug et al., 2009;Ren and Jin, 2011). In order to express the two different types of El Ni?o, i.e., EP and CP El Ni?o, we used the Ni?o3 index (hereafter referred to as the EP index) to represent EP El Ni?o (Trenberth, 1997) and the El Ni?o Modoki index (hereafter referred to as the CP index) to represent CP El Ni?o (Ashok et al., 2007). The EP index was obtained from NOAA’s Climate Prediction Center (http://www.cpc.ncep.noaa.gov/data/indices/ersst4.nino.mth.81-10.ascii),which was computed as the regional SSTA average over(5°S-5°N, 150°-90°W):

    The CP index was defined according to the method of Ashok et al. (2007):

    where [SSTA]C, [SSTA]Wand [SSTA]Eare the regionally averaged SSTA values over (10°S-10°N, 165°E-140°W),(15°S-5°N, 110°-70°W) and (10°S-20°N, 125°-145°E), respectively. Both the EP and CP indices were calculated from the ERSST.v4 dataset, and the relevant regions are denoted in Figs. 2c and d. To define the two types of El Ni?o,we defined typical EP (CP) El Ni?o events when the amplitude of the anomalous EP (CP) index was equal to or greater than 0.8 times its standard deviation (0.8 STD). Based on this approach, there were three EP El Ni?o years (1983,1998 and 2007) and five CP El Ni?o years (1991, 1992,1995, 2003 and 2005) from 1980 to 2008. These EP and CP El Ni?o years were in accordance with previous results(e.g., Ashok et al., 2007; Wang et al., 2013). Aside from the chosen EP and CP years, the remaining years are referred to as “normal” years.

    The Yellow and Yangtze river basins are two major river basins in East China, and two hydrological stations(e.g., Huayuankou and Hankou) are located in the lower reaches of each river basin, respectively. Therefore, these hydrological stations were selected to represent the two abovementioned river basins. The methods of linear correlation and regression analysis were used to explore the relationship between streamflow and the two types of El Ni?o years(positive index years). Composite analysis was applied to reinforce our results. All statistical significance tests for correlation and regression were performed using a two-tailed Student’s t-test. The seasonal anomalies were the departures from their respective climatic states for the study period(i.e., 1980-2008), and linear trends were extracted before analysis. For simplicity, the terms “winter” and “summer” are used to represent December-January-February (DJF) and the following June-July-August, respectively.

    Fig. 2. The winter time-series of (a) EP index and (b) EP index, and composite SSTAs in (c) EP and (d) CP El Ni?o years (units: ℃) for 1980-2008. The dash-dotted line in (a, b) denotes the corresponding 0.8 STD, and the areas of the black solid boxes in (c, d) are used to define the EP and CP El Ni?o index, respectively.

    3. Results

    3.1. Long-term summer streamflow variations

    To illustrate the evolution of streamflow, Fig. 3 shows the time series of summer streamflow anomalies for the two major river basins of East China from 1980-2008 in EP and CP El Ni?o years. The majority of the river basins’ streamflow anomalies display above-normal streamflow in EP years, but complicated variabilities in CP years, in which the streamflow varies in different river basins. Specifically,the anomalous streamflow displays as positive in EP years over the two river basins (based on the Huayuankou and Hankou hydrological stations). In contrast, in all CP years except 2005, the anomalous streamflow shows as negative at Huayuankou hydrological station, whereas it does not display such clear change at Hankou hydrological station(slightly positive). The composite anomalous streamflow in both EP and CP years is also presented in Fig. 3. It is clear that, over the Yellow River basin, the anomalous streamflow is above (below) the multi-year average, with a mean volume of 642.84 m3s-1(-276.48 m3s-1) in EP (CP) years.Meanwhile, over the Yangtze River basin, the anomalous streamflow is positive in EP years, with a mean volume of 8207.96 m3s-1, and in CP years it also shows as slightly positive, with a volume of 1111.86 m3s-1.

    Fig. 3. Annual summer streamflow anomalies in the two river basins during 1980-2008: (a) Huayuankou hydrological station; (b) Hankou hydrological station. The right-most part in the plots shows the composite volume of the streamflow in EP and CP years (units: m3 s-1). The red (blue) asterisks indicate the EP (CP) El Ni?o years. Note that the streamflow data at Huayuankou hydrological station are missing for 2008.

    Fig. 4. The mean and percentiles of summer streamflow over the Yellow and Yangtze river basins for EP, CP and normal years. The top (bottom) of the whiskers represents the 90th (10th) percentile; the top (bottom) of the boxes represents the 75th (25th) percentile; the middle line represents the 50th percentile; and the black dot in the boxes represents the mean volume of streamflow.

    Figure 4 shows the mean and percentiles of summer streamflow over the Yellow and Yangtze river basins for EP, CP and normal years. Consistent with Fig. 3, the amplitude of the streamflow volume at Huayuankou hydrological station is far smaller than that at Hankou hydrological station. At Huayuankou hydrological station (Fig. 4a), the mean volume is lowest in CP years, highest in EP years, and in-between in normal years. Moreover, the 90th percentile of the volume of streamflow (1488.74 m3s-1) in CP years is even lower than the median volume (1953.08 m3s-1) in EP years and the 75th percentile of the volume (1547.29 m3s-1)in normal years, indicating that during CP years the streamflow is smallest and there might be a high probability of hydrological drought. At Hankou hydrological station (Fig.4b), the mean streamflow volume in EP years is abundant and far greater than that in CP and normal years, and its amplitude range spans widely. The mean volume (45675.1 m3s-1) is less than the median volume (49893.0 m3s-1), implying that more than 50% of samples are higher than the average, and the 90th percentile of the volume reaches 57055.3 m3s-1. Large volumes of streamflow are likely to cause riverbanks to collapse and cause flooding over the related river basin.

    Figure 5 is a scatterplot of the normalized summer streamflow anomaly against the EP and CP indices, with the lines of best fit (determined by the least-squares method) shown for the warm and cold phases of the EP and CP indices, respectively. It can be seen from Fig. 5 that changes in streamflow are asymmetrical in (both EP and CP) La Ni?a and El Ni?o years over the Yellow River basin, but symmetrical over the Yangtze River basin. The negative Ni?o index is slightly correlated with the streamflow. The following discussions mainly focus on the warm phases. At Huayuankou hydrological station (Figs. 5a and c), the correlation coefficient of the streamflow anomaly with the EP index is 0.39,and that with the CP index is -0.37. At Hankou hydrological station (Figs. 5b and d), the correlation coefficients of the streamflow anomaly with the EP and CP indices are 0.72 and 0.48, respectively. All the correlation coefficients are statistically significant at the 90% confidence level.

    Consequently, Figs. 3-5 all show that the anomalous summer streamflow over the Yellow River basin is positive in EP years but negative in CP years, while over the Yangtze River basin it is positive in both EP and CP years,despite the composite streamflow anomaly in CP years being only slightly positive during 1980-2008.

    Fig. 5. Normalized EP index values plotted against the normalized summer streamflow anomaly at (a)Huayuankou hydrological station (over the Yellow River basin) and (b) Hankou hydrological station (over Yangtze River basin). (c, d) As in (a, b) but for the CP index. Square symbols denote the EP or CP years’streamflow anomalies. The lines of best fit for positive and negative Ni?o index (EP and CP index,respectively) values (Slope) and the correlation coefficients (R) are shown at the bottom of each plot. The correlation coefficients all pass the significance test at the 90% confidence level.

    3.2. Impacts of the two types of El Ni?o on summer precipitation and ET

    The above analyses show distinctly different changes of observational streamflow in EP and CP years over the Yellow and Yangtze river basins, respectively. In this subsection, we further analyze the connection between these two types of El Ni?o and the changes in surface hydrology elements, especially summer streamflow.

    Streamflow at a downstream station mainly depends on the rainfall over the entire upstream area in the river basin,which is triggered by evaporation, soil moisture, human activities, and so on (Gardner, 2009). Precipitation is an essential water source of streamflow and it is necessary for it to be taken into account when analyzing the variations of streamflow. Figure 6 shows the correlation between the two Ni?o indices and anomalous precipitation and the distributions of the composite anomalous precipitation in EP and CP years.Significant positive correlations between precipitation and the EP index are found primarily in the middle of mainland China between 25°N and 35°N, where the Yangtze River basin is located. The basin-average correlation coefficient is significantly positive (0.52) over the Yangtze River basin,whereas it is non-significantly negative over the Yellow River basin. The correlation coefficients between the CP index and streamflow (Fig. 6b) show that the Yellow River basin is dominated by negative correlation and the basin-average correlation coefficient is -0.31, while most of the Yangtze River basin shows non-significant positive correlation. Figures 6c and d display the distributions of composite anomalous precipitation. They are consistent with the correlation maps. During EP years, the Yangtze River basin is dominated by positive anomalous precipitation (basin average: 0.69 mm d-1), whereas slightly negative anomalous precipitation appears in the middle reaches of the Yellow River basin (0.08 mm d-1). In CP years, clear negative anomalies appear in southern China, with some of them in the region south of the Yangtze River basin. The basin-average precipitation anomaly is -0.14 mm d-1over the Yangtze River basin and 0.06 mm d-1over the Yellow River basin, which are similar to those in EP years.

    Fig. 6. (a, b) Correlation of anomalous summer precipitation with (a) EP index and (b) CP index. Dotted regions are significant at the 90% confidence levels. Black (inverted) triangles denotes (Hankou)Huayuankou hydrological station. The red curves surrounding regions are the ranges of the two river basins.(c, d) Composite maps of anomalous summer precipitation in (c) EP and (d) CP years, but with the lower anomalies (-0.3 to 0.3 mm d-1) left blank (units: mm d-1).

    Apparently, over the Yellow River basin in both EP and CP years, the area-average precipitation anomalies are very small and nearly the same, indicating that precipitation is not the only factor affecting streamflow over this basin. Actually, in arid and semiarid areas, the contributions of ET to streamflow are similarly, if not more, important than those of precipitation (Sun et al., 2017). Thus, we further discuss the changes in ET over these two river basins. Figure 7 shows the same distribution maps as Fig. 6, but based on the PM-MOD ET dataset. On the whole, the correlations between ET and the EP index are negative, but between ET and the CP index they are positive over the mainland. The correlation centers all appear in Northeast China. The Yellow River basin, which can mainly be classed as arid/semiarid, displays negative correlation with the EP index (Fig. 7a; basinaverage correlation: -0.22), indicating that in EP years the ET will be reduced and more soil water will be retained. In CP years, most parts of the Yellow River basin display nonsignificant negative correlation (Fig. 7b). Likewise, the composite distributions of ET are also in accordance with the distribution patterns of correlation. The Yellow River basin shows negative ET anomalies (basin average: -0.05 mm d-1)in EP years, and positive ET anomalies (basin average: 0.02 mm d-1) in CP years. The characteristic changes in ET over the Yangtze River basin show similar features as those over the Yellow River basin, but the area-average values are higher than those over the Yellow River basin (the anomalous ET average value is -0.18 mm d-1in EP years and 0.07 mm d-1in CP years). However, it should be noted that there are no ground observations of ET over a large area, and modelsimulated ET usually contains large uncertainties (e.g.,Wang et al., 2016). Thus, we also used LandFlux-EVAL to verify the results from the PM-MOD dataset. It was found that both ET datasets have consistent characteristics of variation, but the magnitude of LandFlux-EVAL is higher than that of PM-MOD in China. A more detailed comparison of the two ET datasets is provided in the supplemental material [Fig. S1 and Table S1 in the electronic supplementary material (ESM)].

    Moreover, we also analyzed the surface air temperature as compared with the ET (Fig. S2 in ESM). Over the humid regions, ET is more dependent on temperature than precipitation, while over the arid and semiarid regions precipitation and temperature are of equal importance for ET and other soil hydrological variables (Wang and Zeng, 2011). As a typical humid region, the Yangtze River basin usually has enough soil water for evaporation. The ET and temperature exhibit consistent changes and higher temperature would induce more ET. The Yellow River basin belongs to the arid/semiarid region and soil is usually dry and contains less moisture, which limits ET. Rising temperature would increase the atmospheric moisture demand, and then further enhance the process of ET. From our analyses, compared to normal years, the Yellow River basin displays less precipitation, lower temperatures and reduced ET in EP years, while there is less precipitation and higher temperatures and ET in CP years.

    Fig. 7. As in Fig. 6. but for ET (units: mm d-1).

    The above results indicate that the streamflow in the lower reaches is affected by precipitation, ET and temperature in the entire river basin. Over the Yellow River basin,both precipitation and ET play equal roles in the streamflow, and negative anomalous precipitation, ET and temperature will induce positive streamflow anomalies in EP years.Similar results have been obtained in a previous study (e.g.,Liu and Cui, 2011). In CP years, negative streamflow anomalies appear in the lower reaches of the Yellow River basin because of low precipitation and large ET. A positive CP index is accompanied by negative precipitation anomalies and positive ET and temperature anomalies. Over the Yangtze River basin, the streamflow in the lower reaches will increase when a positive EP index occurs with positive precipitation anomalies and negative ET anomalies, due to the river basin having a larger water input and smaller water output.Similarly, the changes of the positive CP index values usually accompany positive precipitation anomalies and negative ET anomalies, leading to positive streamflow anomalies.

    4. Possible mechanisms

    Based on the above results, it is apparent that the different characteristics of surface hydrological variables over different river basins in the two types of El Ni?o years lead to different features of streamflow. It is vital, therefore, to detect the atmospheric physical mechanisms of impact on river flows that occur during these different types of El Ni?o. The expectation is that doing so will improve our understanding of the changes in land surface hydrology with large-scale circulation patterns and assist with land surface hydrological predictions. In the following subsection, we analyze the atmospheric circulation pattern anomalies during EP and CP years, and endeavor to identify the dominant factors impacting on summer streamflow, precipitation, ET and temperature over the two river basins. More specifically, the atmospheric dynamic and thermal variations are examined through analyses of the water vapor transport,lower- to mid-tropospheric circulation patterns, and atmospheric vertical motions during the two types of El Ni?o years.

    4.1. Anomalous water vapor transport

    Moisture transport is an essential factor for precipitation and subsequently for streamflow. To examine the moisture transport during EP and CP years, we performed a correlation analysis of the vertically integrated water vapor transport flux (WVTF) anomalies from the surface to 300 hPa with the two Ni?o indices (Fig. 8). It is obvious that, in both EP and CP years (Figs. 8a-d), there are two clear moisture transport belts, located north (transport towards the east)and south (transport towards the west) of 20°N, respectively. In the mid-latitudes (50°-60°N), there is a transport belt towards the west. As a result, near 20°N (south of 30°N) is mainly occupied by an anticyclone, and near 40°N(north of 30°N) by a cyclone.

    Fig. 8. Regression of the summer vertically integrated WVTF anomalies (from the surface to 300 hPa)(vector arrows) and the convergence of WVTF (contours) with positive (a) EP and (c) CP index values, and the composite patterns in (b) EP and (d) CP years. Only those results that surpass the 80% confidence level in the regression of convergence are drawn in the maps. The letter “A” represents “anticyclone” and the letter“C” represents “cyclone”. The red curves surrounding regions indicate the ranges of the two river basins.

    In EP years, the anticyclone south of 30°N becomes strong and presents two anticyclonic centers (Figs. 8a and b): one over South China, and the other over the northern Bay of Bengal. Meanwhile, north of 30°N, a cyclonic center lies over the Japan Sea. Such a pattern exists clearly in both the composite and regression maps. The anomalous anticyclone is able to bring abundant water vapor through the Indochina Peninsula (bypassing the South China Sea) to southeastern Asia, which will then be transported to the land area of South China. Meanwhile, the anomalous cyclone over the Japan Sea shows a significant WVTF channel from the Northwest Pacific to the Mongolian Plateau, and finally arrives in the Yellow River basin. Moreover, there is significant moisture convergence over the Yangtze River basin,while over the Yellow River basin there is only a small area that shows significant moisture divergence, with other areas characterized by slight moisture convergence. Therefore,such a pattern is favorable for WVTF from the low- to midlatitude sea to the Yangtze River basin, but it is unremarkable for the Yellow River basin. In CP years, the cyclone located north of 30°N becomes robust and shows two cyclonic centers, over the west of the Japan Sea and east of the Kamchatka Peninsula, respectively (Figs. 8c and d). Also,the anticyclone appears over the West Philippines. Such a pattern would suppress the water vapor transport from low latitudes to the land area of East China. The moisture divergence is apparent over the Yellow River basin and southeast coastal area. Figure 8 also shows that there is good consistency with the distribution of both precipitation and the convergence/divergence of water vapor in the two types of El Ni?o years. For instance, in EP years, the precipitation displays a positive anomaly over the Yangtze River basin and North China, where the water vapor significantly converges(Figs. 8a and b). In CP years, the area with profound increasing precipitation is mainly located over the Huaihe River basin, where water vapor is significantly converging (Figs.8c and d). It is worth noting that, over the Yellow River basin, there is a distinct WVTF from the Mongolian Plateau and strong water vapor divergence. The above patterns reduce the precipitation over the Yellow River basin in CP years. Similarly, over South China, the water vapor anomaly shows as significantly positive and the precipitation anomaly is negative, in accordance with the obvious water vapor divergence. Furthermore, we analyzed the characteristics of vertical and horizontal motions of the atmosphere, to investigate the possible reasons of such a phenomenon.

    4.2. Anomalous atmospheric circulation changes

    Fig. 9. Similar to Fig. 8. but for the GPH (contours; units: gpm) and omega (shading; units: m s-1) 500-hPa anomalies. Dashed lines represent negative values, solid lines represent positive values. In regression maps[(a) and (c)], only regions passing the 90% confidence level are plotted (shading) and overlaid by black lines(contour).

    The spatial distributions of the regressed and composite GPH and omega anomalies are plotted for both EP and CP El Ni?o years in Fig. 9. The same patterns as those in Fig. 8 are also apparent in the GPH anomaly fields. These GPH patterns present as an anomalous high-pressure center around 20°N and an anomalous low-pressure center around 40°-50°N. The anomalous low-pressure center corresponds with ascending motion and the anomalous high-pressure center indicates sinking motion. Moreover, in the mid-latitudes the anomalous westerlies are also closely associated with ascending motion. Specifically, in EP years, the anomalous westerlies are located east of 90°E, where the Yangtze and Yellow River basins are (Figs. 9a and b); whereas, in CP years, the westerlies shift east and become narrow and weakened, resulting in the convective activity only being evident over the Huaihe River basin (Figs. 9c and d). Moreover,in the low latitudes in EP years (Figs. 9a and b), South China is dominated by sinking motion, and a high-pressure center is apparent from the anomalous GPH map, which is closely associated with the western Pacific subtropical high(WPSH). The composite WPSH (Fig. S3 in ESM) shifts westward to 115°E in EP years, as compared with CP years(130°E) and normal years (135°E). It prevails over the South China Sea and part of South China, where descending motion is dominant, and southwesterly flow prevails in the north of the WPSH. Such characteristics are consistent with the precipitation changes in South China and the water vapor transport path and its divergence process. Therefore,the changes in local precipitation are affected by both water vapor and the vertical movement.

    Combined analyses of the distributions of ET, temperature and vertical movement illustrate that they have a good and consistent relationship over land. In general, the area with rising motion corresponds to the strong convective zone, which is a favorable condition for the generation of clouds. Clouds can prevent radiation reaching the Earth’s surface, subsequently reducing the surface radiation heating.Thus, there is always negative anomalous ET and temperature where rising motion dominates; whereas, areas of sinking motion correspond to clear-sky and cloudless weather,with more solar radiation reaching the surface and the induction of positive anomalous ET and temperature. Figure 10 displays the average meridional (100°-115°E) latitude-pressure/height distribution of the rising/sinking motion over the two river basins. The area between 25°N and 35°N approximately corresponds to the location of the Yangtze River basin, and the area between 35°N and 45°N is roughly the Yellow River basin. Apparently, in EP years, the Yangtze River basin is dominated by deep rising motion, while the Yellow River basin is dominated by shallow rising motion under 500 hPa and sinking motion above 500 hPa. Accordingly, the anomalous ET and temperature over both the Yangtze and Yellow river basins are obviously negative. In CP years, the deep sinking center is over the Yellow River basin, and a shallow rising center occurs over the Yangtze River basin, corresponding to the positive ET and temperature anomalies over the Yellow River basin but slightly negative ET and temperature anomalies over the Yangtze River basin.

    Therefore, we find that EP and CP El Ni?o influence the lower reaches of both river basins’ streamflow through influencing the summer climate (e.g., precipitation, ET and temperature). Both the moisture conditions and the dynamic conditions are vitally important to the above relationships.

    5. Summary

    Fig. 10. Vertical-horizontal cross-section regression patterns averaged along 100°-115°E for summer vertical wind anomalies with (a) EP and (c) CP index values, and the composite patterns of vertical wind averaged along 100°-115°E in (b) EP and (d) CP years. Only regions passing the 90% confidence level are plotted.

    In this study, to investigate the variations in summer streamflow over the Yellow and Yangtze river basins during different types (EP and CP) of El Ni?o, we analyzed streamflow observations, precipitation, ET and large-scale atmospheric circulations based on reanalysis datasets for the period 1980-2008. We then explored the possible atmospheric physical mechanisms of the above relationships through a set of statistical analyses.

    The summer streamflow in the Yellow and Yangtze river basins shows different characteristics in the two types of El Ni?o years. Generally, the linkage of EP El Ni?o events is more evident than that of CP El Ni?o; streamflow shows a positive anomaly during EP years, but changes diversely in CP years across different river basins. For example, the streamflow anomaly appears as dramatically negative over the Yellow River basin, whereas it is slightly positive over the Yangtze River basin, in CP years. Analyses of precipitation and ET show different weights of the effects on the streamflow over the two river basins. Over the Yellow River basin, the changes in ET are more important than precipitation, while over the Yangtze river basin the change in precipitation is the main reason for streamflow changes.

    The large-scale anomalous atmospheric circulation shows obviously different characteristics in EP and CP El Ni?o years. Overall, the anomalous atmospheric circulation shows a significant cyclone north of 30°N and an anticyclone south of 30°N. Such a spatial distribution may affect the water vapor transport path and modify the distribution of sinking/ascending motion, as well as the dynamic conditions of precipitation. In EP years, the anticyclone reinforces and emerges as two anticyclonic centers: one over South China, and the other over the northern Bay of Bengal.Meanwhile, the cyclone north of 30°N lies over the Japan Sea. Such a pattern is in accordance with a special atmospheric circulation in which the water vapor transported to the Yangtze River basin is enhanced, but does not arrive at the Yellow River basin. Meanwhile, the WPSH extends to 110°E and occupies South China, leading to negative precipitation anomalies but positive ET and temperature anomalies.The westerlies at 500 hPa shift southward and the distribution of vertical movement shows ascent over the Yangtze and Yellow river basins, which is favorable for convective movement, inducing the decreasing ET and temperature.However, in CP years, the cyclone becomes robust and shows two cyclonic centers over the western Japan Sea and east of the Kamchatka Peninsula, respectively. Accordingly,the atmospheric circulation shows less water vapor entering East China, and westerlies become narrow and weakened.At the same time, the vertical movement is one of sinking motion over two river basins, especially the Yangtze River basin. The above situation will cause a decrease in precipitation but an increase in ET and temperature. Such a phenomenon is more pronounced over the Yellow River basin, which will then induce the streamflow to diminish obviously.

    In this paper, we present detailed analyses of the relationship between summer streamflow and two types of El Ni?o over the Yellow and Yangtze River basins. However, it should be pointed out that, due to the limitation of the streamflow data length, during the study period of 1980-2008, EP El Ni?o occurred only three times and CP El Ni?o only five. One might therefore be concerned that the results could be biased in terms of distinguishing the impacts from the background climatology. To make the results more robust, when we performed the correlation/regression analyses, we selected years with positive EP or CP indices, and the composite analyses were from the pre-selected cases based on the thresholds described in section 2. Therefore,we are confident that our results are sufficiently convincing.Wu et al. (2016) mentioned that hydrometeorological hazards (especially flooding) will always be a concern because of the huge economic losses and damage to property that they cause. Therefore, it is crucial to enhance our ability to predict such hazards. With this in mind, the results presented here might help us to better understand the influence of different types of El Ni?o events on land surface hydrological processes, and subsequently lead to advancements in the potential to forecast El Ni?o-related hydrometeorological disasters.

    Acknowledgements.This work was supported by the Key Project of the Ministry of Science and Technology of China (Grant No. 2016YFA0602401) and the National Natural Science Foundation of China (Grant No. 41875106). The authors would like to acknowledge the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences, for providing the China River Valley datasets, which are available for free at their website(http://www.resdc.cn). Other datasets used in this study have been acknowledged in the paper.

    Electronic supplementary material:Supplementary material is available in the online version of this article at https://doi.org/10.1007/s00376-019-9049-2.

    国产日韩欧美在线精品| 老女人水多毛片| 久久6这里有精品| 欧美日韩精品成人综合77777| 蜜臀久久99精品久久宅男| 亚洲精品乱码久久久久久按摩| 人妻夜夜爽99麻豆av| 免费看a级黄色片| 91aial.com中文字幕在线观看| 91aial.com中文字幕在线观看| 中文精品一卡2卡3卡4更新| 国产免费视频播放在线视频| 国产一级毛片在线| 国产一区有黄有色的免费视频| 久久久精品免费免费高清| 国产高清不卡午夜福利| 亚洲怡红院男人天堂| 国产国拍精品亚洲av在线观看| 五月玫瑰六月丁香| 精品午夜福利在线看| 亚洲aⅴ乱码一区二区在线播放| 麻豆成人午夜福利视频| 极品教师在线视频| 美女高潮的动态| 中文字幕亚洲精品专区| 国产久久久一区二区三区| 联通29元200g的流量卡| 亚洲av日韩在线播放| 亚洲av.av天堂| 毛片一级片免费看久久久久| 中文字幕亚洲精品专区| .国产精品久久| 久久久亚洲精品成人影院| 99精国产麻豆久久婷婷| 免费观看av网站的网址| 网址你懂的国产日韩在线| 老司机影院毛片| 91aial.com中文字幕在线观看| 欧美性猛交╳xxx乱大交人| 一本久久精品| 久久精品国产鲁丝片午夜精品| 国产欧美日韩精品一区二区| 一个人看视频在线观看www免费| 成人免费观看视频高清| 黄色一级大片看看| 成人高潮视频无遮挡免费网站| 2021天堂中文幕一二区在线观| 国产精品国产三级国产av玫瑰| 亚洲精品乱久久久久久| 中国国产av一级| 亚洲色图综合在线观看| 亚洲精品日本国产第一区| 在线观看人妻少妇| 亚洲精品aⅴ在线观看| 精品午夜福利在线看| 国产精品嫩草影院av在线观看| 真实男女啪啪啪动态图| 寂寞人妻少妇视频99o| 日韩大片免费观看网站| 久久韩国三级中文字幕| 久久国内精品自在自线图片| 日本一二三区视频观看| 一本一本综合久久| 久久99热6这里只有精品| 久久久a久久爽久久v久久| 波多野结衣巨乳人妻| 日日摸夜夜添夜夜添av毛片| av在线蜜桃| 亚洲怡红院男人天堂| 欧美另类一区| 国产精品一二三区在线看| 精品久久久久久久人妻蜜臀av| 寂寞人妻少妇视频99o| 一级片'在线观看视频| 午夜福利视频精品| 2021少妇久久久久久久久久久| 人妻系列 视频| 中文在线观看免费www的网站| 日韩国内少妇激情av| 亚洲成人av在线免费| 久久久久国产网址| 人体艺术视频欧美日本| 欧美日韩精品成人综合77777| 麻豆国产97在线/欧美| 一级毛片久久久久久久久女| 国内精品宾馆在线| 欧美成人午夜免费资源| 精品久久久久久久人妻蜜臀av| 岛国毛片在线播放| 天堂中文最新版在线下载 | www.av在线官网国产| 日产精品乱码卡一卡2卡三| 噜噜噜噜噜久久久久久91| 久久久久性生活片| 十八禁网站网址无遮挡 | 欧美最新免费一区二区三区| 免费观看在线日韩| 国产成人精品婷婷| 在线观看av片永久免费下载| 欧美精品人与动牲交sv欧美| 国产黄频视频在线观看| 亚洲综合色惰| 久久影院123| 最后的刺客免费高清国语| 免费黄频网站在线观看国产| 22中文网久久字幕| av国产精品久久久久影院| 欧美成人精品欧美一级黄| 人妻制服诱惑在线中文字幕| 欧美日韩精品成人综合77777| 国产成人a∨麻豆精品| 晚上一个人看的免费电影| 欧美一级a爱片免费观看看| 热99国产精品久久久久久7| 丰满乱子伦码专区| 69av精品久久久久久| 一级av片app| 国产高清三级在线| 99re6热这里在线精品视频| 国产永久视频网站| 我要看日韩黄色一级片| 国产一级毛片在线| 国产精品国产av在线观看| 亚洲成人av在线免费| 亚洲熟女精品中文字幕| 久久精品综合一区二区三区| 亚洲色图综合在线观看| 一区二区三区四区激情视频| av女优亚洲男人天堂| 免费看日本二区| 国产成人精品一,二区| 内射极品少妇av片p| 久久久久国产网址| 午夜福利高清视频| 看非洲黑人一级黄片| 69av精品久久久久久| 欧美极品一区二区三区四区| 中文资源天堂在线| 国产毛片a区久久久久| 亚洲精品视频女| 黑人高潮一二区| 一级片'在线观看视频| 狂野欧美激情性bbbbbb| 国产精品一区二区三区四区免费观看| av在线观看视频网站免费| 七月丁香在线播放| 国产成人福利小说| 国语对白做爰xxxⅹ性视频网站| 国产亚洲一区二区精品| 尤物成人国产欧美一区二区三区| 亚洲欧美日韩东京热| 国产精品久久久久久精品电影| 18禁裸乳无遮挡免费网站照片| 欧美高清性xxxxhd video| 波多野结衣巨乳人妻| 交换朋友夫妻互换小说| www.色视频.com| 色播亚洲综合网| 亚洲久久久久久中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 日韩一本色道免费dvd| 国产高清国产精品国产三级 | 超碰97精品在线观看| 国产久久久一区二区三区| 国产精品国产av在线观看| 日本-黄色视频高清免费观看| 男女边吃奶边做爰视频| 最近的中文字幕免费完整| 国产免费又黄又爽又色| 一个人看视频在线观看www免费| 日韩人妻高清精品专区| 九色成人免费人妻av| 一区二区三区四区激情视频| 好男人在线观看高清免费视频| 97在线人人人人妻| 欧美区成人在线视频| 精品一区二区免费观看| 日日摸夜夜添夜夜爱| 国产精品久久久久久久久免| 看黄色毛片网站| 中文天堂在线官网| 亚洲精品久久午夜乱码| 日韩一区二区视频免费看| 成人黄色视频免费在线看| 亚洲综合精品二区| 赤兔流量卡办理| 2022亚洲国产成人精品| 日韩伦理黄色片| 日韩av免费高清视频| 午夜福利视频1000在线观看| 亚洲成人一二三区av| 国产午夜精品一二区理论片| 国产午夜福利久久久久久| 少妇人妻一区二区三区视频| 大片免费播放器 马上看| 国产探花在线观看一区二区| 又大又黄又爽视频免费| 国产中年淑女户外野战色| 国产精品.久久久| 777米奇影视久久| 麻豆精品久久久久久蜜桃| 白带黄色成豆腐渣| 搡老乐熟女国产| 尤物成人国产欧美一区二区三区| 日本-黄色视频高清免费观看| 久久99蜜桃精品久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 午夜亚洲福利在线播放| 看黄色毛片网站| 亚洲国产精品专区欧美| 色婷婷久久久亚洲欧美| 欧美一级a爱片免费观看看| 天堂中文最新版在线下载 | 亚洲av成人精品一区久久| 永久免费av网站大全| 菩萨蛮人人尽说江南好唐韦庄| 日日摸夜夜添夜夜爱| 晚上一个人看的免费电影| 免费看光身美女| 我的女老师完整版在线观看| 成人二区视频| 婷婷色麻豆天堂久久| 丰满少妇做爰视频| 日本熟妇午夜| av一本久久久久| 99热这里只有精品一区| 2021天堂中文幕一二区在线观| 国产亚洲av嫩草精品影院| 另类亚洲欧美激情| 国产精品国产三级专区第一集| 我要看日韩黄色一级片| 亚洲av免费高清在线观看| 日韩视频在线欧美| 男的添女的下面高潮视频| 国产视频首页在线观看| 一级爰片在线观看| 久久热精品热| 久久久久久久午夜电影| 蜜臀久久99精品久久宅男| 99久久精品热视频| 亚洲综合精品二区| 99re6热这里在线精品视频| 久久久久久伊人网av| 大香蕉97超碰在线| 国产乱人偷精品视频| av天堂中文字幕网| 国产精品一及| 99久国产av精品国产电影| 如何舔出高潮| 老司机影院毛片| 99热国产这里只有精品6| 国产黄片视频在线免费观看| a级毛色黄片| 特级一级黄色大片| 亚洲伊人久久精品综合| 能在线免费看毛片的网站| 亚洲国产精品国产精品| 国产一区二区三区综合在线观看 | 日日摸夜夜添夜夜添av毛片| 免费看a级黄色片| 免费观看性生交大片5| 欧美极品一区二区三区四区| 熟女电影av网| 九色成人免费人妻av| 国产乱来视频区| 国内精品美女久久久久久| 国精品久久久久久国模美| 欧美变态另类bdsm刘玥| 久久久久久久久大av| 成人鲁丝片一二三区免费| .国产精品久久| 又大又黄又爽视频免费| 一区二区三区四区激情视频| 国产爽快片一区二区三区| 成人免费观看视频高清| 日韩成人av中文字幕在线观看| 天堂网av新在线| 亚洲真实伦在线观看| 欧美国产精品一级二级三级 | 国产高清有码在线观看视频| 1024香蕉在线观看| 国产亚洲欧美精品永久| 搡老岳熟女国产| 久久这里只有精品19| 久久久久久免费高清国产稀缺| 电影成人av| 2021少妇久久久久久久久久久| 90打野战视频偷拍视频| 亚洲,一卡二卡三卡| 久久久欧美国产精品| 狂野欧美激情性xxxx| 男人舔女人的私密视频| 在线观看免费午夜福利视频| 999精品在线视频| 国产黄色免费在线视频| 精品少妇一区二区三区视频日本电影 | 菩萨蛮人人尽说江南好唐韦庄| 欧美亚洲日本最大视频资源| 国产一区二区三区av在线| 国产 精品1| 肉色欧美久久久久久久蜜桃| 咕卡用的链子| av视频免费观看在线观看| 欧美人与性动交α欧美精品济南到| 亚洲欧美一区二区三区久久| 精品视频人人做人人爽| 少妇人妻久久综合中文| 黄色 视频免费看| 婷婷色麻豆天堂久久| 亚洲国产精品成人久久小说| 大话2 男鬼变身卡| av国产精品久久久久影院| 国产一区有黄有色的免费视频| 亚洲av欧美aⅴ国产| 久久热在线av| kizo精华| 中文字幕最新亚洲高清| 宅男免费午夜| 久久热在线av| 高清视频免费观看一区二区| 777米奇影视久久| 999精品在线视频| 在线精品无人区一区二区三| 国产欧美日韩综合在线一区二区| 久久精品亚洲熟妇少妇任你| 在线精品无人区一区二区三| 久久久久视频综合| 老司机深夜福利视频在线观看 | 悠悠久久av| 欧美日韩成人在线一区二区| 国产日韩欧美视频二区| 一区二区三区精品91| 久久久久精品性色| 制服丝袜香蕉在线| 亚洲国产精品成人久久小说| 母亲3免费完整高清在线观看| 老司机亚洲免费影院| 欧美最新免费一区二区三区| 精品久久久精品久久久| 美女脱内裤让男人舔精品视频| 9191精品国产免费久久| 精品国产乱码久久久久久男人| 热re99久久国产66热| 伦理电影大哥的女人| 高清不卡的av网站| 深夜精品福利| 黄色怎么调成土黄色| 国产精品国产av在线观看| e午夜精品久久久久久久| 一级a爱视频在线免费观看| av有码第一页| 国产一区二区 视频在线| 亚洲精品一二三| 男女边吃奶边做爰视频| 两个人免费观看高清视频| 日日撸夜夜添| 黄色 视频免费看| 美女大奶头黄色视频| av国产精品久久久久影院| 精品少妇黑人巨大在线播放| 久久青草综合色| 亚洲,欧美,日韩| 国产一区二区 视频在线| 国产伦理片在线播放av一区| 美女扒开内裤让男人捅视频| 熟妇人妻不卡中文字幕| 在线看a的网站| 日日摸夜夜添夜夜爱| 午夜福利一区二区在线看| 人人妻人人澡人人看| 侵犯人妻中文字幕一二三四区| 不卡视频在线观看欧美| 国产99久久九九免费精品| 黄色视频不卡| 中文字幕av电影在线播放| 国产毛片在线视频| 高清不卡的av网站| 人妻一区二区av| √禁漫天堂资源中文www| 亚洲图色成人| 中文欧美无线码| 五月天丁香电影| 久久久国产精品麻豆| 国产亚洲午夜精品一区二区久久| 亚洲av成人不卡在线观看播放网 | 久热爱精品视频在线9| 国产成人精品福利久久| 国产成人啪精品午夜网站| 婷婷成人精品国产| 亚洲七黄色美女视频| 亚洲色图综合在线观看| 人妻 亚洲 视频| 亚洲伊人久久精品综合| 亚洲精品中文字幕在线视频| 99久久精品国产亚洲精品| 热re99久久精品国产66热6| 亚洲av日韩在线播放| 久久精品亚洲熟妇少妇任你| 日韩av不卡免费在线播放| 9热在线视频观看99| 叶爱在线成人免费视频播放| 久久ye,这里只有精品| 少妇精品久久久久久久| 国产精品无大码| 国语对白做爰xxxⅹ性视频网站| kizo精华| 久久精品久久精品一区二区三区| 午夜福利,免费看| 美女扒开内裤让男人捅视频| 18禁动态无遮挡网站| 伊人亚洲综合成人网| 又黄又粗又硬又大视频| 不卡视频在线观看欧美| 欧美国产精品一级二级三级| 欧美精品高潮呻吟av久久| 少妇被粗大猛烈的视频| 亚洲三区欧美一区| 免费高清在线观看视频在线观看| 欧美日韩亚洲高清精品| 十分钟在线观看高清视频www| 免费黄频网站在线观看国产| 国产精品 欧美亚洲| 免费在线观看完整版高清| av在线播放精品| 国产男女超爽视频在线观看| 超碰成人久久| 最新的欧美精品一区二区| 性色av一级| 操出白浆在线播放| 亚洲熟女毛片儿| 男女边摸边吃奶| 亚洲精品国产av成人精品| 日韩熟女老妇一区二区性免费视频| 卡戴珊不雅视频在线播放| 欧美日韩亚洲综合一区二区三区_| 一本大道久久a久久精品| 男女免费视频国产| 一区二区三区激情视频| 午夜老司机福利片| 欧美av亚洲av综合av国产av | 久久久国产欧美日韩av| 性少妇av在线| 亚洲美女搞黄在线观看| 男人操女人黄网站| 日韩中文字幕欧美一区二区 | 看免费av毛片| 国产精品一区二区在线观看99| 午夜福利影视在线免费观看| 精品福利永久在线观看| 国产精品久久久久久精品古装| 制服诱惑二区| 亚洲av在线观看美女高潮| 亚洲一级一片aⅴ在线观看| 青春草视频在线免费观看| 狂野欧美激情性xxxx| 在线观看人妻少妇| 麻豆乱淫一区二区| 一区二区三区精品91| 18禁观看日本| 中文精品一卡2卡3卡4更新| 蜜桃在线观看..| 天美传媒精品一区二区| 精品国产乱码久久久久久男人| 97人妻天天添夜夜摸| 久久99精品国语久久久| 久久人人97超碰香蕉20202| 成年女人毛片免费观看观看9 | 高清av免费在线| 亚洲国产成人一精品久久久| 久久99精品国语久久久| 9热在线视频观看99| 在线天堂最新版资源| 高清欧美精品videossex| 久久天躁狠狠躁夜夜2o2o | 人妻一区二区av| 成人漫画全彩无遮挡| 亚洲人成网站在线观看播放| 亚洲伊人久久精品综合| 亚洲精品国产av蜜桃| 欧美日韩一级在线毛片| 在线观看www视频免费| 欧美日韩福利视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| av一本久久久久| 18禁国产床啪视频网站| 欧美另类一区| 男女无遮挡免费网站观看| e午夜精品久久久久久久| 亚洲av在线观看美女高潮| 又大又黄又爽视频免费| 亚洲国产成人一精品久久久| 只有这里有精品99| 亚洲精品第二区| a级毛片在线看网站| 国产精品嫩草影院av在线观看| 国产av精品麻豆| 国产精品嫩草影院av在线观看| 18在线观看网站| 九九爱精品视频在线观看| 一边亲一边摸免费视频| 高清视频免费观看一区二区| 午夜福利影视在线免费观看| 免费人妻精品一区二区三区视频| 97在线人人人人妻| 久久ye,这里只有精品| 成人亚洲精品一区在线观看| 国产黄色视频一区二区在线观看| 美女主播在线视频| 男人添女人高潮全过程视频| 久久精品久久精品一区二区三区| 在线看a的网站| 婷婷成人精品国产| 亚洲国产欧美在线一区| 久久久久精品人妻al黑| 五月天丁香电影| 波多野结衣一区麻豆| 在线观看一区二区三区激情| 在线观看国产h片| 久久久国产精品麻豆| 777久久人妻少妇嫩草av网站| 免费在线观看黄色视频的| 男女无遮挡免费网站观看| 久久精品久久精品一区二区三区| 亚洲av日韩在线播放| 黄色 视频免费看| 精品午夜福利在线看| 亚洲精品美女久久久久99蜜臀 | 波多野结衣av一区二区av| 中文字幕最新亚洲高清| 国产精品av久久久久免费| 午夜激情av网站| bbb黄色大片| 成人午夜精彩视频在线观看| 亚洲婷婷狠狠爱综合网| 18在线观看网站| 男女午夜视频在线观看| 成人影院久久| 在线天堂最新版资源| 啦啦啦啦在线视频资源| 成人国语在线视频| 亚洲精品美女久久久久99蜜臀 | 午夜福利乱码中文字幕| 热re99久久国产66热| 1024视频免费在线观看| av在线app专区| av一本久久久久| 不卡视频在线观看欧美| 日日摸夜夜添夜夜爱| 国产高清不卡午夜福利| 中文字幕精品免费在线观看视频| 男的添女的下面高潮视频| 亚洲av综合色区一区| 国产男人的电影天堂91| 美女高潮到喷水免费观看| a 毛片基地| 亚洲国产av新网站| 亚洲第一青青草原| 亚洲综合精品二区| 日本色播在线视频| 欧美在线黄色| 男男h啪啪无遮挡| 1024视频免费在线观看| 国产野战对白在线观看| 中国国产av一级| 男女之事视频高清在线观看 | 国产日韩欧美在线精品| 国精品久久久久久国模美| 肉色欧美久久久久久久蜜桃| 9热在线视频观看99| av在线观看视频网站免费| 热re99久久精品国产66热6| 国产国语露脸激情在线看| 极品少妇高潮喷水抽搐| 久久久国产一区二区| 成年人免费黄色播放视频| 下体分泌物呈黄色| 精品一区二区免费观看| 一本一本久久a久久精品综合妖精| 2021少妇久久久久久久久久久| 中文字幕人妻熟女乱码| 久久毛片免费看一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产精品一国产av| 日韩一区二区视频免费看| 久久精品国产综合久久久| 欧美亚洲 丝袜 人妻 在线| 搡老岳熟女国产| 美女福利国产在线| 免费人妻精品一区二区三区视频| 亚洲第一av免费看| av卡一久久| 中文欧美无线码| 久久婷婷青草| 精品国产国语对白av| 一本色道久久久久久精品综合| kizo精华| 热99国产精品久久久久久7| 国产国语露脸激情在线看| 熟女av电影| 美女国产高潮福利片在线看| 精品一区二区免费观看| 91精品三级在线观看| 国产黄色视频一区二区在线观看| 成人国产av品久久久| 欧美日韩视频高清一区二区三区二| 成年美女黄网站色视频大全免费| 久久久久久久国产电影| 久久影院123| 中文字幕另类日韩欧美亚洲嫩草| 丝袜在线中文字幕| 国产亚洲最大av| 99精国产麻豆久久婷婷| 飞空精品影院首页|