• <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.

    一级黄色大片毛片| 欧美一区二区精品小视频在线| 国产精品一区www在线观看| 亚洲内射少妇av| 免费观看a级毛片全部| 国产精品久久久久久久电影| 成人高潮视频无遮挡免费网站| 国产色爽女视频免费观看| 亚洲欧美日韩高清专用| 日本av手机在线免费观看| 亚洲成人中文字幕在线播放| av天堂中文字幕网| 午夜视频国产福利| 成人亚洲精品av一区二区| 成人性生交大片免费视频hd| 亚洲高清免费不卡视频| 国产精品三级大全| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区激情短视频| 男插女下体视频免费在线播放| 天堂网av新在线| 美女大奶头视频| 亚洲欧美成人精品一区二区| 日本黄色片子视频| 国产精品久久视频播放| 欧美性猛交╳xxx乱大交人| 亚洲欧洲日产国产| 精品久久久久久久久亚洲| 在线观看av片永久免费下载| 精品不卡国产一区二区三区| 日本熟妇午夜| 国产日韩欧美在线精品| 欧美高清性xxxxhd video| 麻豆精品久久久久久蜜桃| 最近手机中文字幕大全| 少妇裸体淫交视频免费看高清| 亚洲国产欧美人成| 日韩欧美一区二区三区在线观看| 男人舔奶头视频| 免费观看在线日韩| 亚洲av二区三区四区| 黄色欧美视频在线观看| 欧美成人精品欧美一级黄| av女优亚洲男人天堂| 亚洲高清免费不卡视频| 91久久精品国产一区二区成人| 国产亚洲91精品色在线| 国产av在哪里看| 99视频精品全部免费 在线| 久久久久免费精品人妻一区二区| 久久久久久久久久久免费av| 免费人成在线观看视频色| 桃色一区二区三区在线观看| 乱系列少妇在线播放| 亚洲欧美中文字幕日韩二区| 99riav亚洲国产免费| av在线播放精品| 亚洲内射少妇av| 99热网站在线观看| 草草在线视频免费看| 哪个播放器可以免费观看大片| 亚洲激情五月婷婷啪啪| 不卡一级毛片| 国产黄片美女视频| 国产精品99久久久久久久久| 天堂影院成人在线观看| 桃色一区二区三区在线观看| 狂野欧美白嫩少妇大欣赏| 欧美激情久久久久久爽电影| 日日摸夜夜添夜夜添av毛片| 国产精品人妻久久久久久| 国产日韩欧美在线精品| 老师上课跳d突然被开到最大视频| 一卡2卡三卡四卡精品乱码亚洲| 国产一级毛片在线| 最近的中文字幕免费完整| 麻豆成人午夜福利视频| 麻豆国产av国片精品| 国产真实乱freesex| 男人和女人高潮做爰伦理| 最后的刺客免费高清国语| 国产黄色小视频在线观看| 综合色丁香网| 日日干狠狠操夜夜爽| 美女xxoo啪啪120秒动态图| 综合色av麻豆| 99久久中文字幕三级久久日本| 国产精品免费一区二区三区在线| avwww免费| 日韩欧美精品v在线| 国产爱豆传媒在线观看| 亚洲电影在线观看av| 国语自产精品视频在线第100页| 尾随美女入室| 亚洲欧美日韩高清在线视频| 国产大屁股一区二区在线视频| 国产精品蜜桃在线观看 | av在线天堂中文字幕| 亚洲av.av天堂| 精品国产三级普通话版| 看非洲黑人一级黄片| 婷婷色综合大香蕉| 日韩中字成人| 最新中文字幕久久久久| 国产 一区精品| 少妇人妻一区二区三区视频| 久久欧美精品欧美久久欧美| 不卡一级毛片| 精品熟女少妇av免费看| 成年女人看的毛片在线观看| av.在线天堂| 亚洲欧美日韩东京热| 内射极品少妇av片p| 搞女人的毛片| 国产成人一区二区在线| 成人午夜精彩视频在线观看| 国产色婷婷99| 99久久精品一区二区三区| 尾随美女入室| 国产精品不卡视频一区二区| 欧美变态另类bdsm刘玥| 成人一区二区视频在线观看| 亚洲国产精品国产精品| 一本精品99久久精品77| 国产 一区 欧美 日韩| 一本精品99久久精品77| 国产精品永久免费网站| 久久精品国产鲁丝片午夜精品| 成熟少妇高潮喷水视频| 成人鲁丝片一二三区免费| 啦啦啦啦在线视频资源| 啦啦啦啦在线视频资源| 中出人妻视频一区二区| 久久久国产成人精品二区| 免费看日本二区| 51国产日韩欧美| 国产成人aa在线观看| 日本撒尿小便嘘嘘汇集6| 婷婷色av中文字幕| 国产视频内射| 色视频www国产| 男人狂女人下面高潮的视频| 热99re8久久精品国产| 在线免费十八禁| 一级黄片播放器| 国产单亲对白刺激| 国产精品麻豆人妻色哟哟久久 | 久久久久免费精品人妻一区二区| 小蜜桃在线观看免费完整版高清| 亚洲人成网站高清观看| 十八禁国产超污无遮挡网站| 亚洲欧美精品综合久久99| 色尼玛亚洲综合影院| 观看免费一级毛片| 国产成人a区在线观看| 久久久欧美国产精品| 日本熟妇午夜| 精品熟女少妇av免费看| 男女啪啪激烈高潮av片| 日产精品乱码卡一卡2卡三| 久久国产乱子免费精品| 国产精品久久久久久精品电影| 亚洲一级一片aⅴ在线观看| 日本黄大片高清| 一进一出抽搐gif免费好疼| 91aial.com中文字幕在线观看| 中文亚洲av片在线观看爽| 日本黄色视频三级网站网址| 寂寞人妻少妇视频99o| 蜜臀久久99精品久久宅男| 看非洲黑人一级黄片| 国产一区亚洲一区在线观看| 久久人人爽人人爽人人片va| a级毛片a级免费在线| 欧美性猛交╳xxx乱大交人| 欧美在线一区亚洲| 久久午夜福利片| 麻豆成人午夜福利视频| a级毛色黄片| 女的被弄到高潮叫床怎么办| 婷婷色av中文字幕| 在线播放国产精品三级| 看片在线看免费视频| 国产午夜精品论理片| 黄色一级大片看看| 精品一区二区三区人妻视频| 成人毛片a级毛片在线播放| 亚洲真实伦在线观看| www.色视频.com| 亚洲精品色激情综合| 淫秽高清视频在线观看| 久久6这里有精品| 国内精品宾馆在线| 婷婷色av中文字幕| 天堂影院成人在线观看| kizo精华| 女同久久另类99精品国产91| 日产精品乱码卡一卡2卡三| 美女国产视频在线观看| 亚洲欧洲日产国产| 欧美激情在线99| 日本欧美国产在线视频| 精品人妻熟女av久视频| 亚洲av电影不卡..在线观看| 亚洲欧美精品综合久久99| 成人一区二区视频在线观看| 人妻制服诱惑在线中文字幕| 国产伦一二天堂av在线观看| 久久亚洲国产成人精品v| 夜夜看夜夜爽夜夜摸| 身体一侧抽搐| www日本黄色视频网| 国产视频首页在线观看| 天美传媒精品一区二区| 午夜福利视频1000在线观看| 亚洲精品456在线播放app| 成人综合一区亚洲| 少妇熟女aⅴ在线视频| 天堂中文最新版在线下载 | 一区福利在线观看| 久久精品综合一区二区三区| 午夜老司机福利剧场| 久久精品国产亚洲网站| 又爽又黄a免费视频| 久久这里有精品视频免费| 一级毛片aaaaaa免费看小| 内射极品少妇av片p| 人妻夜夜爽99麻豆av| 国产极品天堂在线| 国产精品av视频在线免费观看| 欧美一区二区国产精品久久精品| av在线播放精品| 精品一区二区三区视频在线| 亚洲欧美精品综合久久99| 亚洲丝袜综合中文字幕| 久久久久国产网址| 91午夜精品亚洲一区二区三区| 亚洲美女视频黄频| 色综合站精品国产| 能在线免费看毛片的网站| 国产日本99.免费观看| 老师上课跳d突然被开到最大视频| 国产午夜精品论理片| 成人毛片60女人毛片免费| 只有这里有精品99| 国产私拍福利视频在线观看| 国产精品精品国产色婷婷| 在线观看午夜福利视频| 欧美丝袜亚洲另类| 99久久成人亚洲精品观看| 精品少妇黑人巨大在线播放 | 夫妻性生交免费视频一级片| 亚洲av中文av极速乱| 国产精品av视频在线免费观看| 欧美一区二区国产精品久久精品| 免费观看人在逋| 国产片特级美女逼逼视频| 午夜福利成人在线免费观看| 国产精品电影一区二区三区| 一个人看的www免费观看视频| 精品一区二区三区人妻视频| 国产精品无大码| 狂野欧美激情性xxxx在线观看| 午夜福利视频1000在线观看| 婷婷色av中文字幕| 国产探花在线观看一区二区| 国产白丝娇喘喷水9色精品| 亚洲精品乱码久久久久久按摩| 色综合亚洲欧美另类图片| 我要搜黄色片| 国内久久婷婷六月综合欲色啪| 久久久a久久爽久久v久久| 伦精品一区二区三区| 色哟哟哟哟哟哟| 国产成人91sexporn| 国产黄色小视频在线观看| 能在线免费观看的黄片| 国产色爽女视频免费观看| 成年版毛片免费区| 久久99蜜桃精品久久| 国产一区二区在线av高清观看| 色播亚洲综合网| 啦啦啦啦在线视频资源| 在线国产一区二区在线| 亚洲欧美精品自产自拍| 日韩欧美一区二区三区在线观看| 九九热线精品视视频播放| 在线观看午夜福利视频| 午夜福利成人在线免费观看| 亚洲国产日韩欧美精品在线观看| 欧美一区二区精品小视频在线| 亚洲av.av天堂| 91久久精品国产一区二区成人| 国产日本99.免费观看| 可以在线观看毛片的网站| 天堂中文最新版在线下载 | 尾随美女入室| 此物有八面人人有两片| 欧美精品国产亚洲| 狠狠狠狠99中文字幕| 99视频精品全部免费 在线| 观看美女的网站| 欧美3d第一页| 内射极品少妇av片p| 成人欧美大片| 在线观看免费视频日本深夜| 日本色播在线视频| 99九九线精品视频在线观看视频| 国产麻豆成人av免费视频| 99热这里只有是精品50| 大香蕉久久网| 日本与韩国留学比较| or卡值多少钱| 国产成人aa在线观看| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区在线观看99 | 欧美性猛交╳xxx乱大交人| 国产老妇伦熟女老妇高清| 真实男女啪啪啪动态图| 国产熟女欧美一区二区| 色综合亚洲欧美另类图片| 中国美女看黄片| 九九爱精品视频在线观看| 国产免费男女视频| 国产av在哪里看| a级毛片免费高清观看在线播放| 国产片特级美女逼逼视频| 最近最新中文字幕大全电影3| 黄色配什么色好看| 五月伊人婷婷丁香| 日本av手机在线免费观看| 国产视频内射| 日韩成人伦理影院| 此物有八面人人有两片| 亚洲精品粉嫩美女一区| 日韩高清综合在线| 亚洲久久久久久中文字幕| 三级经典国产精品| 午夜a级毛片| 日韩av在线大香蕉| 欧美成人免费av一区二区三区| 国产精品国产高清国产av| 日韩中字成人| av视频在线观看入口| 精品日产1卡2卡| 99视频精品全部免费 在线| 婷婷六月久久综合丁香| 嫩草影院入口| 国产在视频线在精品| 夜夜看夜夜爽夜夜摸| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人人妻人人澡欧美一区二区| 成人特级av手机在线观看| 非洲黑人性xxxx精品又粗又长| 波多野结衣巨乳人妻| 日韩av不卡免费在线播放| a级毛片a级免费在线| 夜夜夜夜夜久久久久| 一个人看视频在线观看www免费| 97人妻精品一区二区三区麻豆| 搡老妇女老女人老熟妇| 蜜桃亚洲精品一区二区三区| 1000部很黄的大片| 99九九线精品视频在线观看视频| 色视频www国产| 人妻制服诱惑在线中文字幕| 午夜a级毛片| 97热精品久久久久久| 久久久精品大字幕| 寂寞人妻少妇视频99o| 综合色av麻豆| 麻豆av噜噜一区二区三区| 亚洲在久久综合| 亚洲美女视频黄频| 免费观看的影片在线观看| 久久久精品大字幕| 久久久久久九九精品二区国产| 狂野欧美激情性xxxx在线观看| 美女高潮的动态| 亚洲国产精品合色在线| 国产精品av视频在线免费观看| 性色avwww在线观看| 欧美一级a爱片免费观看看| 国产永久视频网站| 黑人欧美特级aaaaaa片| 成年av动漫网址| 99视频精品全部免费 在线| 午夜激情久久久久久久| 免费黄色在线免费观看| 看十八女毛片水多多多| 一级,二级,三级黄色视频| 一区二区三区四区激情视频| 九九爱精品视频在线观看| 最近中文字幕2019免费版| 777米奇影视久久| 三级国产精品片| 一级毛片我不卡| av线在线观看网站| 天堂中文最新版在线下载| 国产视频首页在线观看| 啦啦啦啦在线视频资源| 成人免费观看视频高清| av又黄又爽大尺度在线免费看| 免费黄色在线免费观看| 久久精品国产亚洲av涩爱| 国产一区二区在线观看日韩| 热99国产精品久久久久久7| 久久这里有精品视频免费| 男的添女的下面高潮视频| 午夜免费鲁丝| 久久精品国产a三级三级三级| 我的老师免费观看完整版| 视频中文字幕在线观看| 亚洲av在线观看美女高潮| 寂寞人妻少妇视频99o| 国产精品女同一区二区软件| 九九爱精品视频在线观看| 午夜福利在线观看免费完整高清在| 精品一区在线观看国产| 国产欧美日韩一区二区三区在线 | 一级二级三级毛片免费看| 中文字幕免费在线视频6| 国产成人精品久久久久久| 狂野欧美激情性xxxx在线观看| 九九爱精品视频在线观看| 欧美精品一区二区大全| xxxhd国产人妻xxx| 乱码一卡2卡4卡精品| 久久热精品热| 久久精品久久精品一区二区三区| 久久久国产欧美日韩av| 亚洲精品亚洲一区二区| 一二三四中文在线观看免费高清| 免费黄频网站在线观看国产| 最新的欧美精品一区二区| 大片免费播放器 马上看| 91精品一卡2卡3卡4卡| 男女无遮挡免费网站观看| 午夜福利网站1000一区二区三区| 欧美日韩视频高清一区二区三区二| videos熟女内射| 久久99蜜桃精品久久| 亚洲欧美清纯卡通| 亚洲性久久影院| 国产成人免费观看mmmm| 国产乱来视频区| 中文字幕人妻丝袜制服| 天堂8中文在线网| 一级爰片在线观看| 日本午夜av视频| 美女国产高潮福利片在线看| 一个人看视频在线观看www免费| 97超视频在线观看视频| 超碰97精品在线观看| 亚洲av日韩在线播放| 日日爽夜夜爽网站| 又黄又爽又刺激的免费视频.| 欧美日韩视频精品一区| 中文乱码字字幕精品一区二区三区| 大香蕉久久网| 伦理电影免费视频| 伊人久久国产一区二区| 午夜日本视频在线| 亚洲欧美清纯卡通| 免费大片黄手机在线观看| 午夜久久久在线观看| 老熟女久久久| 国产精品成人在线| 亚洲四区av| 伦理电影大哥的女人| 国产视频首页在线观看| 丝袜喷水一区| av播播在线观看一区| 亚洲国产精品成人久久小说| 国产免费现黄频在线看| 狂野欧美激情性xxxx在线观看| 涩涩av久久男人的天堂| 国产精品偷伦视频观看了| 欧美xxxx性猛交bbbb| 日韩熟女老妇一区二区性免费视频| 久久人人爽人人爽人人片va| 亚洲国产成人一精品久久久| 欧美精品人与动牲交sv欧美| 黄色毛片三级朝国网站| 国产不卡av网站在线观看| 久久久久视频综合| 亚洲不卡免费看| 黄片无遮挡物在线观看| 2021少妇久久久久久久久久久| 91精品一卡2卡3卡4卡| 一本色道久久久久久精品综合| 欧美成人精品欧美一级黄| 日本av免费视频播放| 一级毛片我不卡| 大陆偷拍与自拍| 国产黄片视频在线免费观看| 精品视频人人做人人爽| av福利片在线| 午夜激情福利司机影院| 欧美日韩一区二区视频在线观看视频在线| 99九九线精品视频在线观看视频| 国产高清三级在线| 欧美bdsm另类| 秋霞在线观看毛片| 亚洲国产精品一区二区三区在线| 久久 成人 亚洲| 久久av网站| 大片电影免费在线观看免费| 熟妇人妻不卡中文字幕| 视频在线观看一区二区三区| 亚洲综合精品二区| 波野结衣二区三区在线| 国产高清有码在线观看视频| 亚洲欧美成人精品一区二区| 久久久国产一区二区| 极品少妇高潮喷水抽搐| 男人操女人黄网站| 欧美精品亚洲一区二区| 全区人妻精品视频| 天美传媒精品一区二区| 国产精品人妻久久久久久| 中文精品一卡2卡3卡4更新| videos熟女内射| 欧美bdsm另类| 国产视频内射| 制服人妻中文乱码| 乱码一卡2卡4卡精品| 26uuu在线亚洲综合色| 一边摸一边做爽爽视频免费| 中文字幕人妻丝袜制服| 国产免费一区二区三区四区乱码| 不卡视频在线观看欧美| 国产男女内射视频| 亚洲av成人精品一区久久| 色婷婷av一区二区三区视频| 欧美日韩亚洲高清精品| 超色免费av| 国产精品一区www在线观看| 国产av码专区亚洲av| av在线播放精品| 欧美激情极品国产一区二区三区 | 日韩中文字幕视频在线看片| 丝瓜视频免费看黄片| 男人操女人黄网站| 亚洲精品第二区| 人成视频在线观看免费观看| 99久国产av精品国产电影| 精品国产一区二区三区久久久樱花| 一二三四中文在线观看免费高清| 欧美人与善性xxx| av天堂久久9| 日韩精品免费视频一区二区三区 | 亚洲av二区三区四区| 日本爱情动作片www.在线观看| 一级毛片aaaaaa免费看小| 大码成人一级视频| av在线老鸭窝| 一区二区av电影网| 亚洲精品乱码久久久v下载方式| 欧美精品一区二区免费开放| 最近2019中文字幕mv第一页| 久久久久久久久久久久大奶| 曰老女人黄片| 交换朋友夫妻互换小说| av女优亚洲男人天堂| 看免费成人av毛片| 日日爽夜夜爽网站| 亚洲三级黄色毛片| 国产男人的电影天堂91| 国产日韩欧美亚洲二区| 搡女人真爽免费视频火全软件| 18禁动态无遮挡网站| 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| 成人手机av| 水蜜桃什么品种好| 亚洲欧美一区二区三区黑人 | 亚洲久久久国产精品| 免费大片18禁| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 18+在线观看网站| 亚洲内射少妇av| 午夜免费观看性视频| 亚洲综合色惰| 国产精品一二三区在线看| 亚洲国产精品一区二区三区在线| 性色avwww在线观看| 人人妻人人爽人人添夜夜欢视频| 精品午夜福利在线看| 国产日韩欧美在线精品| 国产国拍精品亚洲av在线观看| 一个人免费看片子| 国产精品人妻久久久影院| 国产色婷婷99| 又粗又硬又长又爽又黄的视频| 精品久久久噜噜| 日本猛色少妇xxxxx猛交久久| 欧美变态另类bdsm刘玥| 丰满迷人的少妇在线观看| 免费看不卡的av| 欧美日韩视频高清一区二区三区二| 国产一区二区三区综合在线观看 | 亚洲精品久久午夜乱码| 国产精品.久久久| 丰满少妇做爰视频| 制服丝袜香蕉在线| 熟女人妻精品中文字幕| 亚洲av国产av综合av卡| 啦啦啦啦在线视频资源| 少妇丰满av| 99久久人妻综合| 黑人高潮一二区| 国产又色又爽无遮挡免| 人人妻人人添人人爽欧美一区卜|