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

    Upper- and Lower-tropospheric Circulation Anomalies Associated with Interannual Variation of Pakistan Rainfall during Summer

    2020-10-15 10:09:26RiyuLUSaadiaHINAandXiaoweiHONGLASGInstituteofAtmosphericPhysicsChineseAcademyofSciencesBeijing100029China
    Advances in Atmospheric Sciences 2020年11期

    Riyu LU, Saadia HINA, and Xiaowei HONGLASG, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

    2Key Laboratory of Meteorological Disaster (KLME), Ministry of Education &Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD), Nanjing University of Information Science and Technology, Nanjing 210044, China

    3Department of Environmental Sciences and Engineering, Government College University, Faisalabad 38000, Pakistan

    4Climate Change Research Center (CCRC), Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China

    5College of Earth and Planetary Sciences, University of the Chinese Academy of Sciences, Beijing 100049, China

    ABSTRACT This study investigated the large-scale circulation anomalies, in both the upper and lower troposphere, associated with the interannual variation of rainfall in Pakistan during summer, using the station observation data in this country and circulation data of the NCEP?NCAR reanalysis from 1981 to 2017. Results showed that the upper- and lower-tropospheric circulation anomalies associated with monthly rainfall variability exhibit similar features from June to August, so analyses were performed on June?August circulation and Pakistan rainfall data. The analyzed results indicated that summer rainfall in Pakistan is enhanced when there is an anticyclonic anomaly to the northwest of Pakistan in the upper troposphere and easterly anomalies along the southern foothills of the Himalayas in the lower troposphere, and vice versa. These upper- and lower-tropospheric circulation anomalies were found to be related, but show unique features. The upper-tropospheric anticyclonic anomaly is closely related to the Silk Road Pattern along the Asian westerly jet, while the lower-tropospheric easterly anomalies are related to the cyclonic anomaly to the south of Pakistan, i.e., intensified South Asian monsoon trough. The results presented here suggest that the interannual variability of summer rainfall in Pakistan is a combined result of upper- and lower-tropospheric circulation anomalies, and of extratropical and tropical circulation anomalies.

    Key words: rainfall, Pakistan, interannual variability, Silk Road Pattern, monsoon trough

    1. Introduction

    Rainfall, especially summer rainfall, is crucial for agriculture, water resources, food security and industry in Pakistan(Webster et al., 2011; Latif et al., 2017). Pakistan receives less rainfall in comparison with other regions of South Asia,so it is particularly sensitive to variability in rainfall. Therefore, various studies have been performed to understand the spatiotemporal features and variability of rainfall in Pakistan (e.g., Hussain and Lee, 2013; Adnan et al., 2016;Hina and Saleem, 2019). The period from June to September is considered as the monsoon season in Pakistan, and this season yields nearly 60% of the total annual rainfall(e.g., Hussain and Lee, 2016). During summer, rainfall is much greater in the northern part of Pakistan than the south-ern part (e.g., Hanif et al., 2013).

    As a synoptic disturbance, Pakistan rainfall is greatly affected by large-scale circulations in both the lower and upper troposphere. Syed et al. (2010) found that an upper-tropospheric anticyclonic anomaly appears to the northwest of Pakistan a few days before the active phase of rainfall, and is centered at the northernmost region of Pakistan on the onset day of the active phase. In addition, Latif and Syed(2016) indicated that the upper-tropospheric anticyclonic anomaly over West Asia shows its strongest intensity two days before the monsoon onset (their Fig. 4), while the lower-tropospheric cyclonic anomaly to the south of Pakistan reaches its peak on the onset day (their Fig. 5).This lag of the lower-tropospheric anomaly to the upper-tropospheric anomaly suggests that the former may be affected by the latter. That is, the upper-tropospheric anticyclonic anomaly over West Asia, which is associated with upstream circulation anomalies such as the summer North Atlantic Oscillation (SNAO), may induce the lower-tropospheric cyclonic anomaly, which is further favorable to enhanced rainfall in Pakistan.

    Summer rainfall in Pakistan shows strong interannual variability, particularly relative to seasonal rainfall amounts(Naheed and Rasul, 2011). The coefficient variation, i.e.,the ratio of the interannual standard deviation to the climatological average of seasonal rainfall amounts, is as high as 26% for summer rainfall in Pakistan (Hussain and Lee,2016). Therefore, it is of great importance to better understand the mechanisms responsible for the interannual variability of summer rainfall in Pakistan.

    It has been well documented that, on the interannual time scale, the upper-tropospheric anticyclonic/cyclonic anomalies over West Asia are closely related to more/less precipitation in Pakistan or India (Kripalani et al., 1997; Wu, 2002,2017; Krishnan et al., 2009; Hong and Lu, 2016; Latif et al.,2017; Yadav, 2017). The upper-tropospheric anticyclonic/cyclonic anomalies are closely associated with the circumglobal teleconnection (CGT; Ding and Wang, 2005) or Silk Road Pattern (SRP; Lu et al., 2002; Enomoto et al., 2003),which are the two teleconnection patterns similar to each other at least over the Eurasian continent (Zhou et al.,2019). Previous studies have indicated that the CGT and SNAO affect summer rainfall in North Pakistan and Northwest India (Ding and Wang, 2005; Syed et al., 2012). Syed et al. (2012) showed that the positive phases of the CGT and SNAO, characterized by upper-tropospheric anticyclonic anomalies over the Northeast Atlantic and West Asia, respectively, are related to enhanced rainfall over Pakistan during the core summer season (July?August). Kazmi et al. (2016)utilized April?June upper-tropospheric zonal wind at 200 hPa over the northeastern Atlantic as a predictor for the interannual variation of July?August Pakistan rainfall, and explained this relationship through the appearance of the CGT in July?August.

    The lower-tropospheric circulation anomalies associated with the interannual variation of Pakistan summer rainfall have been relatively less well investigated, as compared to their upper-tropospheric counterparts. Ullah and Gao(2012) suggested that the route of water vapor transport from the Arabian Sea to South Asia, which is associated with the monsoon trough, is crucial for the wet and dry conditions in Pakistan during summer. They arrived at this conclusion by analyzing two contrasting summers: 1994, an above-normal rainfall year in Pakistan, and 2002, a below-normal rainfall year. It has been well documented that the monsoon trough, which appears in the lower troposphere, is an important component of the South Asian summer monsoon,favoring active monsoon rainfall (e.g., Krishnamurti and Bhalme, 1976; Bollasina and Nigam, 2011). Therefore, it is reasonable to hypothesize that the lower-tropospheric circulation anomalies may play a crucial role in affecting the interannual variation of Pakistan summer rainfall.

    In this study, we hypothesize that the interannual variability of Pakistan summer rainfall is a result of the tropical?extratropical interaction, which is characterized by the interaction between the tropical monsoon troughs and extratropical disturbances. Actually, some previous studies have demonstrated that such tropical?extratropical interaction plays a crucial role in triggering extreme rainfall weather events in Pakistan or South Asia (e.g., Hong et al., 2011; Vellore et al., 2016). On the other hand, this interaction may also be manifested as that between the upper and lower troposphere, i.e., the interaction between the anticyclonic anomaly to the northwest of Pakistan in the upper troposphere and the stronger monsoon trough in the lower troposphere.Therefore, the aim of this study is to investigate the interactions of the circulations associated with Pakistan rainfall variability, both meridionally and vertically.

    The rest of this paper is organized as follows. The data and methods are described in section 2. The circulation anomalies associated with monthly Pakistan rainfall from June to September are shown in section 3, which indicates that these monthly circulation anomalies are similar in June,July and August. Therefore, in section 4, focusing on the June-July-August (JJA) averages, we show the upper- and lower-tropospheric circulation anomalies and compare these anomalies. Section 5 presents the main conclusions of this study and discusses some associated issues.

    2. Data and methods

    In this study, we used the monthly rainfall data of 30 weather stations, evenly distributed across Pakistan,obtained from the Pakistan Meteorological Department(PMD). The dataset has been widely used in previous studies (e.g., Iqbal and Athar, 2018). We analyzed the rainfall data during the period 1981?2017, i.e., 37 years.

    The monthly circulation reanalysis data are from the National Centers for Environmental Prediction?National Center for Atmospheric Research (Kalnay et al. 1996), with a horizontal resolution of 2.5° × 2.5°. In addition, the monthly precipitation data from the Climatic Research Unit(CRU), which are based on weather station records (version 4.04) (Harris et al. 2014), at a 0.5° × 0.5° horizontal resolution, were also applied. The time span of all the datasets is 1981?2017, consistent with that of the Pakistan station precipitation data. The anomalies stated in this study, for both precipitation and circulation, indicate the deviations relative to the climatological mean during 1981?2017.

    3. Circulation anomalies associated with the interannual variation of monthly rainfall in Pakistan

    Figure 1 shows the interannual variations of station-averaged (over all 30 stations from the PMD dataset) monthly precipitation in Pakistan from June to September. July and August are the months with the highest rainfall in Pakistan,and less rainfall occurs in the other two months (June and September). This is consistent with previous studies, which indicated that Pakistan rainfall in July and August is the strongest when compared with other months (Ding and Ke,2013; Kazmi et al., 2016; Wang et al., 2019), and dominates the interannual variability of seasonal monsoon rainfall(Ahmad et al. 2014). Furthermore, the monthly precipitation exhibits a large interannual variability. For instance,July precipitation is as little as 19.4 mm d?1in 2002, but as high as 137.6 mm d?1in 1995. Even for the months with lower precipitation, i.e., June and September, the rainfall can be quite heavy in some years. For instance, the precipitation is 83.7 in June 2007 and 102.0 mm d?1in September 2012, both being higher than the average of July precipitation (70.5 mm d?1).

    Fig. 1. Interannual variations of station-averaged monthly precipitation (unit: mm) in Pakistan for (a) June, (b) July, (c)August and (d) September.

    Fig. 2. 200-hPa horizontal wind anomalies (vectors; units: m s?1) regressed onto the standardized monthly Pakistan rainfall for (a) June, (b) July, (c) August and (d) September. Shading indicates zonal or meridional wind anomalies significant at the 0.05 level based on the Student’s t-test.

    Figure 2 shows the 200-hPa horizontal wind anomalies regressed onto monthly Pakistan rainfall each month from June. There is a significant anticyclonic anomaly to the northwest of Pakistan for all three months from June to August.This is consistent with previous studies, which reported that the upper-tropospheric anticyclonic anomaly over West Asia is related to more precipitation in Pakistan or India during summer (e.g., Kripalani et al., 1997; Wu, 2002, 2017).However, the September precipitation corresponds to weak circulation anomalies to the northwest of Pakistan (Fig. 2d),or in other words, the anticyclonic anomaly over this region appearing in the preceding months (Figs. 2a?c) is almost invisible in September (Fig. 2d). We have also analyzed the composite circulation anomalies in the excessive and deficient rainfall years, respectively, for each month from June to September, and found that there is an anticyclonic anomaly over West Asia in excessive years and a cyclonic anomaly in deficient years for each month from June to August,but the wind anomalies are weak in both the excessive and deficient years for September (not shown), consistent with the regression results.

    In addition to the anticyclonic anomaly over West Asia,an anomalous anticyclone appears over East Asia for all the months from June to August. A cyclonic anomaly appears between these two anticyclonic anomalies, but exhibits apparent changes in shape and location from month to month. Despite the monthly changes, the circulation anomalies over the Eurasian continent generally resemble the CGT or SRP, confirming the relationship between these teleconnection patterns and Pakistan rainfall. In September, there is no clear teleconnection pattern along the Asian jet (Fig. 2d), again showing a distinct change compared to the preceding months.

    Fig. 3. 700-hPa horizontal wind anomalies (vectors; units: m s?1) and the moisture flux divergence integrated from 925 to 600 hPa (shading; units: 10?6 kg m?2 s?1) regressed onto the standardized monthly Pakistan rainfall for (a)June, (b) July, (c) August and (d) September. The blue and red shading indicates anomalous convergence and divergence of the moisture flux, respectively, and the dots represent the regions significant at the 0.05 level based on the Student’s t-test. The zonal or meridional wind anomalies significant at the 0.05 level are emboldened. The solid lines delineate the topography of 3000 m.

    Figure 3 shows the lower-tropospheric wind anomalies and moisture flux divergence regressed on the standardized monthly Pakistan rainfall from June to September. There is a cyclonic anomaly to the south of Pakistan associated with increased June, July and August rainfall (Figs. 3a?c), indicating a strengthened monsoon trough and more water vapor transport into Pakistan as a result. This is further confirmed by the moisture flux divergence integrated from 925 to 600 hPa regressed onto the standardized monthly Pakistan rainfall. There are negative anomalies, which represent inward flux of moisture, over most regions of Pakistan, particularly over the southern part of the country. Weak outward moisture flux appears over the central-north region of Pakistan,which is inconsistent with heavier rainfall. This discrepancy is possibly due to the errors of moisture flux in the reanalysis data related to topographic effects, and needs more investigation. Associated with the cyclonic circulation anomaly, significant easterly anomalies appear over Pakistan.The monthly circulation anomalies are similar to the vertically integrated water vapor transport anomalies shown in Ullah and Gao (2012, their Fig. 2), who contrasted two typical years (wet year of 1994 vs dry year of 2002). Note that the anomalous moisture fluxes integrated from 925 to 600 hPa (not shown) have almost the same directions as the 700-hPa anomalous winds from June to August, suggesting that the winds play a crucial role in inducing the moisture fluxes. For September, by contrast, the cyclonic anomaly shifts westwards, the wind anomalies are much weaker over Pakistan, and the inward flux of moisture is also much weakened and appears hardly significant over the southern part of Pakistan. These suggest that the wind and moisture transport anomalies associated with monthly rainfall variability in September are distinct from the other three months,i.e., June, July, and August.

    4. Circulation anomalies associated with the interannual variation of summer rainfall

    The results reported in the preceding section indicate that both the upper- and lower-tropospheric circulation anomalies associated with monthly Pakistan rainfall are similar between June, July and August, but distinct from those for September. Therefore, considering a concise presentation of results, in the following we investigate the variability of JJA precipitation and associated circulation anomalies.

    The time series of standardized JJA-average precipitation variation (PKprec) is shown in Fig. 4a. Figure 5 shows the horizontal wind anomalies in the upper and lower troposphere regressed onto PKprec. Both the upper-tropospheric anticyclonic anomaly to the northwest of Pakistan and the lower-tropospheric cyclonic anomaly to the south of Pakistan tend to be clearer in comparison with the anomalies associated with monthly precipitation (Figs. 2 and 3).The anticyclonic anomaly in the upper troposphere is also manifested as a component of the SRP, and the cyclonic anomaly over central Asia tends to be weaker and southward-shifted in comparison with the anticyclonic anomalies over West Asia and East Asia (Fig. 5a), somewhat similar to the monthly results for June, July and August (Figs.2a?c). On the other hand, in the lower troposphere, the significant easterly anomalies at the southern side of the Tibetan Plateau and moisture transport into the southern part of Pakistan are associated with the cyclonic anomaly to the south (Fig. 5b). We have also analyzed the composite wind anomalies for the above-normal and below-normal precipitation years, respectively, and found that the results show approximately out-of-phase circulation patterns and similar circulation intensities in both the upper and lower troposphere (figures not shown, but can be partially found in Figs. 4b and c), suggesting a general linear relationship between the Pakistan summer precipitation with the upper-tropospheric anticyclonic/cyclonic anomaly to the northwest of Pakistan and the lower-tropospheric easterly/westerly anomaly over Pakistan.

    Fig. 4. (a) Standardized time series of JJA mean precipitation anomalies in Pakistan (PKprec), the ACI, and the EASTI. (b?d) Scatterplots of the (b) PKprec and ACI, (c) PKprec and EASTI, and (d) ACI and EASTI. Red symbols in (b?d) indicate the values in 2010.

    Fig. 5. Horizontal wind anomalies (vectors; units: m s?1) at (a) 200 hPa and (b) 700 hPa regressed onto the standardized Pakistan JJA rainfall index (PKprec). Shading in (a) indicates zonal or meridional wind anomalies significant at the 0.05 level based on the Student’s t-test, and in (b) indicates the moisture flux divergence integrated from 925 to 600 hPa (units: 10?6 kg m?2 s?1) regressed onto the standardized PKprec. The blue and red shading indicates anomalous convergence and divergence of the moisture flux, respectively, and dots represent regions significant at the 0.05 level based on the Student’s t-test. Zonal or meridional wind anomalies significant at the 0.05 level are emboldened in (b). Red boxes indicate the domain used to define the ACI in (a) and EASTI in (b). The solid line in (b) delineates the topography higher than 3000 m.

    In the following we use the upper-tropospheric anticyclonic anomaly and lower-tropospheric easterly anomaly to represent the typical circulation anomalies responsible for more rainfall in Pakistan, and accordingly define two indexes.The anticyclonic anomaly index (ACI) is defined as the 200-hPa vorticity anomalies averaged over the region(35°?45°N, 55°?80°E) and multiplied by minus one so positive values correspond to anticyclonic flow. In addition, the easterly anomaly index (EASTI) is defined as the 700-hPa zonal wind anomalies averaged over the region (25°?35°N,60°?80°E) and likewise multiplied by minus one so positive values correspond to easterly flow. Note that these two regions do not overlap. The time series for these indices are included in Fig. 4a, and they are both significantly correlated to the Pakistan rainfall index. The correlation coefficients between the ACI/EASTI and Pakistan rainfall index are 0.39/0.67, significant at the 0.05 and 0.01 level, respectively. The higher correlation with the lower-tropospheric circulation index implies the crucial role of moisture transport in affecting the Pakistan summer rainfall. Scatterplots between the Pakistan rainfall index with both the ACI (Fig.4b) and EASTI (Fig. 4c) also confirm the close relationship between the rainfall and upper- and lower-tropospheric circulation anomalies: Most dots (26 in Fig. 4b and 30 in Fig. 4c out of 37) are scattered in the first and third quadrants. In Fig. 4c, particularly, there is not a dot in the second quadrant, and the EASTI is weak (lower than 0.5 standard deviations) in the fourth quadrant, confirming the close relationship between the Pakistan rainfall and lower-tropospheric wind anomalies. On the other hand, the ACI and EASTI also have a significant relationship with each other (Fig.4d), with the correlation efficient being 0.49. This suggests a vertical coupling between the larger-scale upper- and lower-tropospheric anomalies, possibly via vertical motion and diabatic heating associated with the rainfall anomalies,although we cannot strictly exclude the possibility that the relationship is induced by other unknown factor(s).

    Figure 6 shows the spatial distributions of station rainfall anomalies regressed onto the standardized time series of the ACI, EASTI and Pakistan rainfall index, respectively.The ACI is positively related to precipitation over most stations (Fig. 6a). Particularly, the six stations with a significant (0.05 level) relation all have positive regression. The nine stations with negative relation are all in northern and central Pakistan, where rainfall is relatively weak. By contrast, the EASTI is more clearly related to the station rainfall and there is only one station showing a negative relationship (Fig. 6b), in agreement with the higher correlation coefficient between the EASTI and Pakistan rainfall index. The Pakistan rainfall index is positively correlated with rainfall at almost all the stations, and there are only five stations showing an insignificant relationship.

    Figure 7 shows the horizontal wind anomalies at 200 and 700 hPa regressed onto the standardized ACI. The anticyclonic anomaly to the northwest of Pakistan at 200 hPa is stronger and more significant (Fig. 7a), in comparison with that associated with the Pakistan rainfall index (Fig. 5a). Correspondingly, the SRP is clearer, manifested as the anticyclonic anomalies over West Asia and East Asia and the cyclonic anomaly over central Asia. The meridional wind anomalies associated with these anticyclonic and cyclonic anomalies exhibit alternately southerly and northerly anomalies along the Asian upper-tropospheric westerly jet.

    On the other hand, the ACI is associated with the significant easterly anomalies in the lower troposphere at the southern side of the Tibetan Plateau (Fig. 7b). The relationship between the ACI and lower-tropospheric easterly anom-alies suggests a vertical coupling of circulation anomalies over South Asia. The correlation coefficient between the ACI and EASTI is 0.49, significant at the 0.01 level, confirming the vertical coupling. However, we could not find a clear cyclonic anomaly to the south of Pakistan, which appears in the anomalies regressed onto the Pakistan rainfall index (Fig. 5b), and the anomalous moisture transport into Pakistan (Fig. 7b) seems to be weaker than that related to the Pakistan rainfall (Fig. 5b).

    Fig. 6. Rainfall anomalies (units: mm d?1) at the stations used in this study regressed onto the standardized time series of the (a) ACI, (b) EASTI and(c) Pakistan rainfall index. The definitions of these three indexes can be found in the text. Solid dots indicate anomalies significant at the 0.05 level.

    The above results suggest that there may be a vertical coupling between the upper- and lower-tropospheric circulation anomalies. In the following, we further investigate the vertical structure of circulation anomalies associated with Pakistan rainfall variability. Figure 9 shows a latitude?height cross section of zonal wind anomalies averaged over 60°?80°E, regressed onto the standardized time series of Pakistan rainfall index, ACI and EASTI, respectively. The zonal winds can depict well both the anticyclonic anomaly in the upper troposphere and easterly anomaly in the lower troposphere. Above-normal rainfall in Pakistan is associated with easterly anomalies approximately over 25°?35°N and westerly anomalies over 40°?50°N (Fig. 9a). These zonal wind anomalies are strongest at about 200 hPa, and roughly in a barotropic structure, but with a northward tilt with height. The easterly anomalies extend to the surface, confirming that the easterly anomalies along the foothills of the Himalayas favor more rainfall in Pakistan. The ACI corresponds to a similar distribution of zonal wind anomalies(Fig. 9b). These anomalies also tilt northward with height,but the westerly anomalies tend to tilt southwards in the lower troposphere. The easterly anomalies also tend to tilt southwards in the lower troposphere, but note that these anomalies may be artificial due to the plateaus over 30°?40°N. In addition, the easterly anomalies are mainly concentrated in the upper troposphere, and decrease sharply in the lower troposphere, especially when compared with the anomalies associated with the Pakistan rainfall index. The ACI-related easterly anomalies at 200 hPa are roughly two times larger than those associated with the Pakistan rainfall index, but at 700 hPa the former are smaller than the latter. The EASTI-related anomalies (Fig. 9c) are similar to those associated with the Pakistan rainfall (Fig. 9a), but they are relatively stronger throughout the entire troposphere, which is consistent with the fact that the EASTI is more highly correlated with the ACI than the Pakistan rainfall index (correlation coefficients 0.49 vs 0.39). In addition, there are westerly anomalies over 10°?20°N, and they are limited basically below the mid troposphere, suggesting the existence of significant cyclonic anomalies to the south of Pakistan. Despite the existence of differences, all the vertical structures of zonal wind anomalies shown in Fig. 9, particularly the easterly anomalies, suggest a vertical coupling, i.e., there is a close relationship in the vertical direction.

    5. Conclusions and discussion

    5.1. Conclusions

    In this study, we investigate the circulation anomalies associated with interannual variations of Pakistan rainfall during summer. It is found that wind anomalies associated with

    Fig. 8. As in Fig. 5 but for anomalies regressed onto the standardized EASTI.

    Fig. 9. Latitude?height cross section of zonal wind anomalies (contours; units: m s?1) regressed onto the standardized time series of the (a) Pakistan rainfall index, (b) ACI and (c) EASTI. These anomalies are averaged over 60°?80°E. Contour intervals are 0.5 m s?1, and zero contours are omitted. Shading indicates anomalies significant at the 0.05 level.

    Fig. 7. As in Fig. 5 but for anomalies regressed onto the standardized ACI.Pakistan monthly rainfall are similar between June, July and August, in both the upper and lower troposphere. Therefore,we focus on the summer-mean (i.e., JJA mean) results.

    The results indicate that more Pakistan summer rainfall is significantly associated with an upper-tropospheric anticyclonic anomaly to the northwest of Pakistan and lower-tropospheric easterly anomalies along the foothills of the Himalayas, including Pakistan. The lower-tropospheric easterly anomalies, in comparison with the upper-tropospheric anticyclonic anomaly, show a closer relationship with Pakistan rainfall, suggesting the important role of lower-tropospheric zonal wind anomalies in affecting rainfall through modulating water vapor transport. These upper- and lower-tropospheric circulation anomalies also show a close relationship,and the vertical section of anomalies indicates a roughly barotropic structure throughout the entire troposphere.

    On the other hand, the upper- and lower-tropospheric anomalies show separate features. The upper-tropospheric anticyclonic anomaly is closely related to the SRP along the upper-tropospheric Asian jet and the easterly anomalies along the foothills of the Himalayas, but not to the lower-tropospheric cyclonic anomaly to the south of Pakistan. By contrast, the lower-tropospheric easterly anomalies are closely related to the cyclonic anomaly to the south of Pakistan, but do not correspond well to the SRP, although showing a close relationship to the upper-tropospheric anticyclonic anomaly. Thus, attention should also be paid to different physical mechanisms responsible for these upper- and lower-tropospheric circulation anomalies associated with Pakistan rainfall.

    5.2. Discussion

    5.2.1. Comparison with flooding in 2010

    An unusually heavy rainfall event occurred during late July to early August 2010 in Pakistan, causing tremendous disasters, including widespread landslides, and a death toll of nearly 3000 in the country. The event triggered many studies (Hong et al., 2011; Galarneau et al., 2012; Lau and Kim,2012; Ullah and Gao, 2012; Martius et al., 2013;Rasmussen et al., 2015), the results of which suggest that this heavy rainfall event came about due to a combination of upper-tropospheric extratropical Rossby waves and lower-tropospheric monsoonal depressions. However, as noted by Wang et al. (2011), a cyclonic anomaly appeared to the northwest of Pakistan in the upper troposphere during the Pakistan flood period in summer 2010, in sharp contrast to the anticyclonic anomaly associated with more summer rainfall in Pakistan on the interannual time scale shown by previous studies (e.g., Syed et al., 2012; Wu, 2017) and the present study. In our analysis, the upper-tropospheric ACI has its largest negative value in 2010 in the 1981?2017 analysis period (see Fig. 4), confirming the existence of a cyclonic anomaly to the west of Pakistan in summer 2010.

    Such a discrepancy may be attributable to differences in time scales. The flooding during summer 2010 essentially was caused by synoptic events with time scales of days, while in this study we focused on the interannual variability of monthly or seasonal rainfall. The rains that caused the flooding do not contribute much to the monthly precipitation in Pakistan, as can be seen from the slightly above-normal July and August 2010 rainfall in Fig. 1. Moreover, as mentioned by Webster et al. (2011), there were some other factors resulting in the devastating flooding of summer 2010, in addition to the fact that rainfall was high, but not unexceptionally so. These factors include the precursory severe drought and resultant sparser vegetation during 2010,deforestation in northern Pakistan, and the concurrence between heavy rainfall and the mountainous region. Furthermore, the flooding during summer 2010 tended to appear over a limited spatial scope, which can be supported by the findings of Vellore et al. (2016): 34 extreme rainfall events occurred in the northwest region of India, which is the eastern neighbor of Pakistan, were selected in their study, but none of them occurred during summer 2010.

    5.2.2. Subseasonal change

    Two - Bingo, the dog who had lived in our house with its previous occupants, started to pay us unannounced and traumatic() visits. This gave Bill an opportunity to demonstrate his way with dogs.

    The similarity of circulation anomalies associated with monthly Pakistan rainfall from June to August, shown in Figs. 2 and 3, suggests that the relationship between circulation and rainfall is stable during summer. However, there are also some appreciable differences in the monthly relationship between circulation and rainfall, including that the upper-tropospheric cyclonic anomaly over southern China is strong in June and August but weak in July, and the lower-tropospheric cyclonic anomaly to the south of Pakistan is strongest in June but weakest in August. These monthly differences suggest that there would be a subseasonal change in the relationship between circulation and Pakistan rainfall.

    In the following, we further discuss the subseasonal difference based on a comparison with the results of Wang et al. (2019), who identified a negative relationship in July?August precipitation between the southeastern Tibetan Plateau and Pakistan on the interannual time scale, and suggested that the former affects the latter. Figure 10 shows that in July and August a negative rainfall anomaly appears over the southeastern Tibetan Plateau, corresponding to more rainfall in Pakistan, consistent with Wang et al.(2019), but the rainfall anomaly is weak in June. Thus, this difference in rainfall anomaly over the southeastern Tibetan Plateau may partially explain the monthly different circulation anomalies associated with Pakistan rainfall, and there would be other influencing factor(s), especially for June. In summary, the subseasonal difference in the relationship between rainfall variability over Pakistan and circulation anomalies should be further investigated.

    5.2.3. Vertical coupling

    Fig. 10. CRU rainfall anomalies (shading; units: mm d?1) regressed onto the standardized time series of the Pakistan rainfall index in (a) June, (b) July, (c) August, and (d) JJA. Marked areas indicate anomalies significant at the 0.05 level, and solid lines delineate the topography higher than 3000 m. The red parallelogram in (d) is used to obtain the domain-averaged precipitation anomalies, in order to calculate its temporal correlation coefficient with the Pakistan station rainfall.

    In this study, we suggest that both the upper- and lower-tropospheric circulation anomalies are crucial in affecting Pakistan summer rainfall. These circulation anomalies are related to each other through vertical coupling. Regarding the causality between the upper-tropospheric anticyclonic anomaly and lower-tropospheric easterly anomalies,we incline to the view that the former affects the later, considering that the former can be induced by the upstream circulation anomalies which are associated with the SRP.However, caution should be applied to any conclusion on causality. First, the upper-tropospheric anticyclonic anomaly and lower-tropospheric easterly anomalies show some unique features, suggesting that they might be affected by different factors. The upper-tropospheric anomalies tend to be more closely related to the extratropical circulation anomalies over the Eurasian continent, while the lower-tropospheric ones more closely related to the tropical or monsoonal anomalies. Second, the anomalous rainfall, as a diabatic heating forcing, may in turn have a feedback to the upper- and lower-tropospheric circulation anomalies.

    The large-scale circulation anomalies associated with Pakistan summer rainfall can be helpful for an understanding of the predictability of rainfall. On the one hand, the upper-tropospheric anomalies, which are related to the extratropical SRP, might be of low predictability (Kosaka et al., 2012). On the other hand, the more influential lower-tropospheric easterly anomalies are in the tropical regions and thus may be of greater predictability. The predictive relation between circulation factors such as the EASTI and ACI and Pakistan summer rainfall would be an interesting topic for further investigation.

    5.2.4. Other issues

    Although this study focuses on the rainfall variation in Pakistan, the results can be extended to Northwest India. On the one hand, Pakistan rainfall is closely associated with rainfall anomalies in Northwest India, in each month of summer and in terms of the summer mean (Fig. 10). On the other hand, various previous studies have indicated that the SRP or CGT can influence rainfall in Northwest India, as well as in Pakistan (e.g., Kripalani et al., 1997; Wu, 2002,2017; Ding and Wang, 2005; Syed et al., 2012).

    In addition, there is good agreement between the CRU precipitation data and Pakistan station rainfall data, which can be found in Fig. 10. This consistency between the two datasets can be further confirmed by the high correlation coefficients between them. We used the parallelogram in Fig. 10d to represent Pakistan rainfall in the CRU precipitation data and obtained the correlation coefficients between the two datasets: 0.77 for June, 0.88 for July, 0.81 for August, and 0.90 for JJA. This result suggests that both datasets are reliable for investigating the interannual variability of precipitation in Pakistan.

    Situated at the northernmost edge of the Asian monsoonal region and to the south of the Asian upper-tropospheric westerly jet, Pakistan experiences complex climate variability in summer, and further studies on Pakistan’s climate variability may be of great help towards a better understanding of tropical?extratropical interaction related to the strongest monsoon (i.e., Asian summer monsoon) in the world.

    Acknowledgements. We sincerely thank the Executive Editor-in-Chief, anonymous reviewers and Editor, who gave us insightful and constructive comments and suggestions. This research was jointly sponsored by the National Natural Science Foundation of China (Grant Nos. 41705044 and 41721004), the Joint Open Project of KLME and CIC-FEMD (Grant No.KLME20 1802), and the second Tibetan Plateau Scientific Expedition and Research (STEP) Program (Grant No. 2019QZKK0102).

    欧美另类亚洲清纯唯美| 国产精品一区二区三区四区久久| www.www免费av| 欧美中文日本在线观看视频| 免费看光身美女| 最近在线观看免费完整版| 校园春色视频在线观看| 免费看光身美女| 国产精品影院久久| 床上黄色一级片| 十八禁人妻一区二区| 女人被狂操c到高潮| 一区二区三区国产精品乱码| 999精品在线视频| 免费av毛片视频| 免费搜索国产男女视频| 精品久久蜜臀av无| 天堂av国产一区二区熟女人妻| 一级黄色大片毛片| 亚洲成人中文字幕在线播放| 性色av乱码一区二区三区2| 亚洲成人久久爱视频| 夜夜躁狠狠躁天天躁| 九九在线视频观看精品| 色噜噜av男人的天堂激情| 免费观看人在逋| 久久这里只有精品19| 精品国产乱子伦一区二区三区| 亚洲成人免费电影在线观看| 国产精品一区二区三区四区久久| 欧美成狂野欧美在线观看| 午夜激情福利司机影院| 国产精品 欧美亚洲| 色老头精品视频在线观看| 午夜福利18| 久久久久久大精品| 麻豆国产av国片精品| 最新美女视频免费是黄的| 在线观看免费视频日本深夜| 男女床上黄色一级片免费看| 毛片女人毛片| 日日摸夜夜添夜夜添小说| 我的老师免费观看完整版| 国产成年人精品一区二区| 99riav亚洲国产免费| 丰满人妻熟妇乱又伦精品不卡| 曰老女人黄片| 国产伦精品一区二区三区视频9 | 制服丝袜大香蕉在线| 三级毛片av免费| 国产1区2区3区精品| 免费在线观看影片大全网站| 成人特级av手机在线观看| 日韩有码中文字幕| 啦啦啦韩国在线观看视频| 亚洲一区二区三区色噜噜| 美女高潮喷水抽搐中文字幕| 美女 人体艺术 gogo| 国产1区2区3区精品| 欧美日韩国产亚洲二区| 每晚都被弄得嗷嗷叫到高潮| a级毛片a级免费在线| 国产成人av教育| 亚洲av中文字字幕乱码综合| 亚洲 国产 在线| 久久久久久久午夜电影| 久久久国产成人免费| 高潮久久久久久久久久久不卡| 国产极品精品免费视频能看的| 日韩欧美三级三区| 禁无遮挡网站| 日日摸夜夜添夜夜添小说| 又粗又爽又猛毛片免费看| 亚洲国产欧美人成| 可以在线观看的亚洲视频| 亚洲在线观看片| 九色国产91popny在线| 亚洲第一欧美日韩一区二区三区| 久久久久性生活片| 国产v大片淫在线免费观看| 国产在线精品亚洲第一网站| 国产又黄又爽又无遮挡在线| 99精品在免费线老司机午夜| 天天一区二区日本电影三级| 久久香蕉精品热| 18禁观看日本| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区三| 成年女人永久免费观看视频| 男女之事视频高清在线观看| 欧美成人一区二区免费高清观看 | 亚洲 欧美 日韩 在线 免费| 欧美av亚洲av综合av国产av| 老鸭窝网址在线观看| 欧美国产日韩亚洲一区| av天堂中文字幕网| 国产私拍福利视频在线观看| 国产美女午夜福利| 日韩欧美在线乱码| 亚洲精品在线观看二区| 俄罗斯特黄特色一大片| 精品久久蜜臀av无| 老司机在亚洲福利影院| 一进一出好大好爽视频| 男人的好看免费观看在线视频| 精品一区二区三区av网在线观看| 此物有八面人人有两片| 99国产综合亚洲精品| www国产在线视频色| 日韩精品中文字幕看吧| 久久精品夜夜夜夜夜久久蜜豆| 亚洲五月天丁香| 中文在线观看免费www的网站| 亚洲精品在线美女| 国产精品九九99| 两人在一起打扑克的视频| 又黄又粗又硬又大视频| 男女下面进入的视频免费午夜| 亚洲精品久久国产高清桃花| 精品久久久久久久人妻蜜臀av| 精品人妻1区二区| 香蕉丝袜av| 99热这里只有精品一区 | 观看美女的网站| 精品午夜福利视频在线观看一区| 一级黄色大片毛片| 亚洲av成人av| 亚洲精品456在线播放app | 麻豆成人午夜福利视频| 少妇的逼水好多| 人人妻人人澡欧美一区二区| 中国美女看黄片| 亚洲在线观看片| 神马国产精品三级电影在线观看| 性色av乱码一区二区三区2| tocl精华| 中文字幕高清在线视频| 亚洲中文字幕日韩| 久久精品夜夜夜夜夜久久蜜豆| 男人的好看免费观看在线视频| 美女被艹到高潮喷水动态| 国产av在哪里看| 黑人操中国人逼视频| 欧美日韩综合久久久久久 | 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美 国产精品| or卡值多少钱| 级片在线观看| 国产男靠女视频免费网站| 岛国在线观看网站| 欧美大码av| 国产高清视频在线播放一区| 精品国产乱码久久久久久男人| 美女扒开内裤让男人捅视频| 国产精品综合久久久久久久免费| 香蕉丝袜av| www.熟女人妻精品国产| 亚洲成人久久性| 99国产综合亚洲精品| 法律面前人人平等表现在哪些方面| 久久精品91蜜桃| 琪琪午夜伦伦电影理论片6080| 欧美黄色淫秽网站| 国产黄片美女视频| 免费高清视频大片| 国产日本99.免费观看| www日本在线高清视频| 久久天堂一区二区三区四区| av片东京热男人的天堂| 网址你懂的国产日韩在线| 亚洲av成人av| 中出人妻视频一区二区| 九九在线视频观看精品| 午夜福利免费观看在线| 老司机福利观看| 国产一区二区三区视频了| av在线天堂中文字幕| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| 少妇裸体淫交视频免费看高清| 日本三级黄在线观看| 12—13女人毛片做爰片一| 叶爱在线成人免费视频播放| 怎么达到女性高潮| 村上凉子中文字幕在线| a在线观看视频网站| 成人精品一区二区免费| 99re在线观看精品视频| 一本一本综合久久| 91九色精品人成在线观看| 久久婷婷人人爽人人干人人爱| 最近在线观看免费完整版| 亚洲专区字幕在线| 日韩欧美三级三区| 日本黄色片子视频| 亚洲熟女毛片儿| 中文字幕熟女人妻在线| 国产精品 欧美亚洲| 国产爱豆传媒在线观看| 最近最新中文字幕大全免费视频| 亚洲av中文字字幕乱码综合| 国产91精品成人一区二区三区| 亚洲欧美激情综合另类| 黄色成人免费大全| 亚洲av日韩精品久久久久久密| 亚洲avbb在线观看| 日韩免费av在线播放| 亚洲国产精品成人综合色| 窝窝影院91人妻| 99re在线观看精品视频| 99在线视频只有这里精品首页| 在线十欧美十亚洲十日本专区| 亚洲熟妇中文字幕五十中出| 亚洲第一电影网av| 午夜福利在线在线| 欧美成人一区二区免费高清观看 | 精品免费久久久久久久清纯| 免费观看精品视频网站| 五月伊人婷婷丁香| 国产欧美日韩一区二区精品| 三级男女做爰猛烈吃奶摸视频| 视频区欧美日本亚洲| 日本三级黄在线观看| 波多野结衣高清无吗| 国产精品久久视频播放| 亚洲av片天天在线观看| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 女同久久另类99精品国产91| 亚洲七黄色美女视频| 中文字幕久久专区| 色噜噜av男人的天堂激情| 国产成人aa在线观看| h日本视频在线播放| 欧美性猛交黑人性爽| 高清在线国产一区| 免费看a级黄色片| 给我免费播放毛片高清在线观看| 老熟妇乱子伦视频在线观看| 最近视频中文字幕2019在线8| av黄色大香蕉| 一a级毛片在线观看| 操出白浆在线播放| 欧美日韩综合久久久久久 | 美女被艹到高潮喷水动态| 精品国产三级普通话版| 99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 欧美日韩乱码在线| 窝窝影院91人妻| 久久久成人免费电影| 精品无人区乱码1区二区| 亚洲成人免费电影在线观看| 国产精品久久视频播放| 国产伦一二天堂av在线观看| 床上黄色一级片| svipshipincom国产片| 2021天堂中文幕一二区在线观| 亚洲精华国产精华精| 波多野结衣巨乳人妻| 欧美午夜高清在线| 亚洲精品456在线播放app | 亚洲欧美精品综合久久99| 99久久久亚洲精品蜜臀av| 久久天躁狠狠躁夜夜2o2o| netflix在线观看网站| 美女高潮喷水抽搐中文字幕| 中文在线观看免费www的网站| 国产精品久久久av美女十八| 可以在线观看的亚洲视频| 久久香蕉精品热| 欧美成人性av电影在线观看| 中文字幕高清在线视频| 五月玫瑰六月丁香| 日韩欧美在线二视频| 精品人妻1区二区| 欧美色视频一区免费| 国产黄片美女视频| 一个人免费在线观看电影 | 日韩av在线大香蕉| 午夜日韩欧美国产| 国产av不卡久久| 成人18禁在线播放| 一级a爱片免费观看的视频| 三级毛片av免费| 久久精品91蜜桃| 亚洲国产欧美一区二区综合| 欧美精品啪啪一区二区三区| 露出奶头的视频| 国产精品一区二区免费欧美| 在线永久观看黄色视频| 婷婷丁香在线五月| 午夜福利在线在线| 亚洲人成电影免费在线| 日韩中文字幕欧美一区二区| 18美女黄网站色大片免费观看| 欧美绝顶高潮抽搐喷水| 亚洲国产欧美人成| 欧美一级a爱片免费观看看| 久久婷婷人人爽人人干人人爱| 亚洲精品456在线播放app | www.精华液| 国产激情久久老熟女| 欧美日本亚洲视频在线播放| 亚洲人成伊人成综合网2020| 神马国产精品三级电影在线观看| 午夜激情福利司机影院| 日韩欧美 国产精品| 国产精华一区二区三区| 国产乱人伦免费视频| 亚洲欧美一区二区三区黑人| 老司机深夜福利视频在线观看| 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 精品国产美女av久久久久小说| 老司机福利观看| 无遮挡黄片免费观看| 国产成人av教育| 两个人的视频大全免费| 亚洲五月婷婷丁香| 一级黄色大片毛片| 天天躁日日操中文字幕| 超碰成人久久| 国产欧美日韩精品亚洲av| x7x7x7水蜜桃| 国产高潮美女av| 国产成人精品久久二区二区91| 国产 一区 欧美 日韩| 午夜a级毛片| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 成年女人毛片免费观看观看9| 午夜福利在线在线| 三级毛片av免费| 欧美一区二区精品小视频在线| 国模一区二区三区四区视频 | 1024手机看黄色片| 国产精品久久久久久久电影 | 在线播放国产精品三级| 国内揄拍国产精品人妻在线| 蜜桃久久精品国产亚洲av| 在线视频色国产色| 最新在线观看一区二区三区| 亚洲欧美日韩无卡精品| 国产精品亚洲av一区麻豆| 99国产极品粉嫩在线观看| 精品一区二区三区视频在线 | 日本精品一区二区三区蜜桃| 淫妇啪啪啪对白视频| 久久久久久久久久黄片| 九色国产91popny在线| 91麻豆精品激情在线观看国产| 国产精品永久免费网站| 热99在线观看视频| 亚洲国产精品久久男人天堂| 免费av毛片视频| 法律面前人人平等表现在哪些方面| 91av网一区二区| 免费看日本二区| av中文乱码字幕在线| 国产成人av教育| 成人国产一区最新在线观看| 97超级碰碰碰精品色视频在线观看| 午夜福利在线在线| xxxwww97欧美| 亚洲国产精品sss在线观看| 又爽又黄无遮挡网站| 欧美av亚洲av综合av国产av| 色av中文字幕| 欧美av亚洲av综合av国产av| 校园春色视频在线观看| 亚洲国产看品久久| 欧美日韩亚洲国产一区二区在线观看| 美女大奶头视频| 狂野欧美激情性xxxx| 一级a爱片免费观看的视频| 亚洲欧美精品综合久久99| 少妇人妻一区二区三区视频| 婷婷丁香在线五月| 国产精品久久久av美女十八| 蜜桃久久精品国产亚洲av| 精品国产三级普通话版| 国内久久婷婷六月综合欲色啪| 男女做爰动态图高潮gif福利片| 校园春色视频在线观看| 国产精品av视频在线免费观看| svipshipincom国产片| 丰满的人妻完整版| 国产人伦9x9x在线观看| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 国产精品九九99| 国产成人精品久久二区二区免费| 俄罗斯特黄特色一大片| 老司机深夜福利视频在线观看| 亚洲欧美一区二区三区黑人| 免费av不卡在线播放| 日韩高清综合在线| 国产99白浆流出| 一本一本综合久久| 国内少妇人妻偷人精品xxx网站 | 国产主播在线观看一区二区| 久久国产精品人妻蜜桃| 天堂√8在线中文| 国产欧美日韩精品一区二区| 嫩草影院精品99| 国产成+人综合+亚洲专区| 国产伦人伦偷精品视频| 国产毛片a区久久久久| a级毛片在线看网站| 亚洲人成网站在线播放欧美日韩| 麻豆成人午夜福利视频| 久久欧美精品欧美久久欧美| 精品久久蜜臀av无| 婷婷六月久久综合丁香| 男女那种视频在线观看| 久久香蕉国产精品| 国产熟女xx| 亚洲七黄色美女视频| 国产精品99久久99久久久不卡| 欧美激情在线99| 三级男女做爰猛烈吃奶摸视频| 黄片大片在线免费观看| 亚洲成人久久性| 很黄的视频免费| 国产成人系列免费观看| 国产精品 欧美亚洲| 日韩欧美三级三区| 丁香六月欧美| 91av网站免费观看| 激情在线观看视频在线高清| 国产亚洲精品久久久com| 国产一区二区三区视频了| 国产精品美女特级片免费视频播放器 | 俄罗斯特黄特色一大片| 亚洲国产精品合色在线| 18禁裸乳无遮挡免费网站照片| 午夜福利在线观看免费完整高清在 | 免费av毛片视频| 天堂网av新在线| 九九热线精品视视频播放| 欧美日韩国产亚洲二区| 亚洲av熟女| 又紧又爽又黄一区二区| 久久午夜亚洲精品久久| 国产三级中文精品| 狂野欧美激情性xxxx| 婷婷六月久久综合丁香| 人妻丰满熟妇av一区二区三区| 久久久久久久午夜电影| 老司机深夜福利视频在线观看| 天天躁日日操中文字幕| 999精品在线视频| 中文字幕精品亚洲无线码一区| 日本与韩国留学比较| 国产激情久久老熟女| 亚洲av美国av| 亚洲国产欧美网| 一个人免费在线观看电影 | 久久久久久久久中文| 露出奶头的视频| 亚洲欧美日韩东京热| a级毛片在线看网站| 精品国产亚洲在线| xxxwww97欧美| 黄色女人牲交| 亚洲 欧美一区二区三区| or卡值多少钱| 日韩精品中文字幕看吧| 久久久色成人| 久久欧美精品欧美久久欧美| 欧美中文日本在线观看视频| 18禁国产床啪视频网站| 欧美高清成人免费视频www| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 麻豆成人午夜福利视频| 美女大奶头视频| 一本精品99久久精品77| 一个人免费在线观看的高清视频| 日本熟妇午夜| 黄片大片在线免费观看| 男人舔奶头视频| 中出人妻视频一区二区| 国产精品免费一区二区三区在线| 少妇的逼水好多| 国产成人精品无人区| 国产精品98久久久久久宅男小说| 成年人黄色毛片网站| 免费看a级黄色片| 国产久久久一区二区三区| 日韩欧美国产一区二区入口| 午夜福利成人在线免费观看| 亚洲18禁久久av| 亚洲无线观看免费| 免费看十八禁软件| 神马国产精品三级电影在线观看| 日韩欧美在线二视频| 亚洲熟妇中文字幕五十中出| 免费在线观看视频国产中文字幕亚洲| 在线观看午夜福利视频| 亚洲av成人不卡在线观看播放网| 99热这里只有精品一区 | 国产精品久久电影中文字幕| www.999成人在线观看| 国产伦在线观看视频一区| 18美女黄网站色大片免费观看| 999久久久国产精品视频| 久久精品国产99精品国产亚洲性色| 精品国产美女av久久久久小说| 麻豆一二三区av精品| 两个人视频免费观看高清| 国产99白浆流出| 欧美黄色淫秽网站| 亚洲av美国av| 国产精品99久久久久久久久| 国产亚洲av嫩草精品影院| 国产三级在线视频| 午夜久久久久精精品| 欧美性猛交黑人性爽| 网址你懂的国产日韩在线| 操出白浆在线播放| 亚洲片人在线观看| 一级毛片精品| 日韩av在线大香蕉| 999久久久国产精品视频| 精品99又大又爽又粗少妇毛片 | 波多野结衣巨乳人妻| cao死你这个sao货| 夜夜看夜夜爽夜夜摸| 成年免费大片在线观看| 亚洲激情在线av| 91在线观看av| 桃红色精品国产亚洲av| 日韩欧美精品v在线| 两个人看的免费小视频| 国产精品av久久久久免费| 麻豆成人av在线观看| 免费在线观看亚洲国产| 18禁国产床啪视频网站| 无限看片的www在线观看| 两性夫妻黄色片| 国产三级在线视频| 久久久久久久久免费视频了| 日本免费a在线| 国产99白浆流出| 中出人妻视频一区二区| 久久久水蜜桃国产精品网| 99国产精品一区二区三区| 午夜精品在线福利| 国产一区二区在线观看日韩 | 国产成+人综合+亚洲专区| 毛片女人毛片| 五月伊人婷婷丁香| 国产激情偷乱视频一区二区| 国产精品久久视频播放| 亚洲真实伦在线观看| 中文字幕最新亚洲高清| 久久中文字幕人妻熟女| 精品国产亚洲在线| 婷婷丁香在线五月| 国产三级黄色录像| 小说图片视频综合网站| 女人被狂操c到高潮| 久久精品影院6| 老熟妇乱子伦视频在线观看| 国产美女午夜福利| 欧美午夜高清在线| 国产精品99久久久久久久久| 麻豆国产av国片精品| 午夜a级毛片| 亚洲中文日韩欧美视频| 一本久久中文字幕| av国产免费在线观看| 国产精品野战在线观看| 最近最新免费中文字幕在线| 最好的美女福利视频网| 99热这里只有是精品50| 欧美zozozo另类| 伊人久久大香线蕉亚洲五| 午夜久久久久精精品| 午夜亚洲福利在线播放| 一进一出抽搐gif免费好疼| 美女扒开内裤让男人捅视频| 国产一区二区三区在线臀色熟女| 欧美三级亚洲精品| 成人三级做爰电影| 久久久久久国产a免费观看| 国产亚洲av高清不卡| www国产在线视频色| 亚洲精品国产精品久久久不卡| 三级男女做爰猛烈吃奶摸视频| 国产69精品久久久久777片 | 无限看片的www在线观看| 午夜精品久久久久久毛片777| 国产精品av视频在线免费观看| 中文资源天堂在线| 亚洲色图 男人天堂 中文字幕| 日本 av在线| 亚洲成av人片在线播放无| 男插女下体视频免费在线播放| 啦啦啦免费观看视频1| 国产av麻豆久久久久久久| www.www免费av| 精华霜和精华液先用哪个| 可以在线观看毛片的网站| 99riav亚洲国产免费| 啦啦啦免费观看视频1| 午夜日韩欧美国产| tocl精华| 国产蜜桃级精品一区二区三区| 国产免费男女视频| 黄色片一级片一级黄色片| 88av欧美| 欧美不卡视频在线免费观看|