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

    Establishment of an Objective Standard for the Definition of Binary Tropical Cyclones in the Western North Pacific

    2020-10-15 10:09:30FuminRENYanjunXIEBiwenYINMingyangWANGandGuopingLI
    Advances in Atmospheric Sciences 2020年11期

    Fumin REN, Yanjun XIE, Biwen YIN, Mingyang WANG, and Guoping LI

    1State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China

    2School of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu 610225, China

    3Zhuzhou Meteorological Observatory, Zhuzhou 412000, China

    4Qinhuangdao Meteorological Observatory, Qinhuangdao 066000, China

    ABSTRACT

    To develop an objective standard for defining binary tropical cyclones (BTCs) in the western North Pacific (WNP),two best-track datasets, from the China Meteorological Administration and the Joint Typhoon Warning Center, were adopted for statistical analyses on two important characteristics of BTCs—two TCs approaching each other, and counterclockwise spinning. Based on the high consistency between the two datasets, we established an objective standard,which includes a main standard for defining BTCs and a secondary standard for identifying typical/atypical BTCs. The main standard includes two requirements: two coexisting TCs are a pair of BTCs if (i) the separation distance is ≤ 1800 km,and (ii) this separation maintains for at least 12 h. Meanwhile, the secondary standard defines a typical BTC as one for which there is at least one observation when the two TCs approach each other and spin counterclockwise simultaneously.Under the standard, the ratio of typical BTCs increases as the BTC duration increases or the minimum distance between the two TCs decreases. Then, using the JTWC dataset, it was found that there are 505 pairs of BTCs during the period 1951?2014, including 328 typical BTCs and 177 atypical BTCs, accounting for 65.0% and 35.0% of the total, respectively.In addition, a study of two extreme phenomena—the maximum approaching speed and the maximum counterclockwise angular velocity in typical BTCs—shows that the configuration of the circulation conditions and the distribution of the BTCs favor the formation of these extreme phenomena.

    Key words: objective standard, binary tropical cyclones, Western North Pacific

    1. Introduction

    A binary vortex interaction, in which two cyclonic vortices at a relatively close range simultaneously perform a mutual counterclockwise spin and move closer to each other, is referred to as a binary tropical cyclone (BTC) and the Fujiwhara effect in honor of the pioneering work of Fujiwhara (1921, 1923, 1931). Chen and Ding (1979) pointed out that, when two TCs are enough close to each other, for either one, along the line between their centers, the pressure gradient force decreases suddenly on the side close to the other TC while the pressure gradient force remains the same on the opposite side. This situation is conducive to the two TCs approaching each other. Haurwitz (1951) studied the motion of BTCs and reported a method to compute the angular velocity of the rotation of the binary axis.

    This typical motion of BTCs (the Fujiwhara effect)does not always occur when two TCs are close to each other(Hoover, 1961; Brand, 1970; Dong and Neumann, 1983;Lander and Holland, 1993; Carr et al., 1997; Carr and Elsberry, 1998; Wu et al., 2011). Hoover (1961) reported that a separation distance of 600 nautical miles is the critical distance determining whether hurricane pairs in the Atlantic Ocean have the characteristics of BTCs, whereas Brand(1970) suggested that 750 nautical miles is the threshold separation distance in the western North Pacific (WNP).Lander and Holland (1993) showed that most interacting cyclones show four characteristic features: approach and capture, mutual orbit, merger, and escape. Rapid merger only occurs when BTCs approach within a critical separation distance, such as 150?300 km (Ritchie and Holland, 1993).Carr et al. (1997) proposed three modes for BTCs: direct interaction, semi-direct interaction, and indirect interaction. Carr and Elsberry (1998) carried out an objective diagnosis of the interactions of BTCs and presented objective criteria for distinguishing these three modes. BTCs have different manifestations because the relative motion of TC pairs is a function of both the Fujiwhara effect and differences in the large-scale environmental steering forces acting on each system (Brand, 1970; Dong and Neumann, 1983; Wang and Fu, 1983; Kuo et al., 2000; Wu et al., 2011). The significance of the Fujiwara effect in a TC pair depends on the distance between the two centers, the size, structure, and intensity of the two TCs, as well as the environmental steering currents (Brand, 1970).

    Many different models have been proposed to explain the interactions of BTCs in terms of dynamics and thermodynamics (Chang, 1983; Holland and Dietachmayer, 1993;Falkovich et al., 1995; Wang and Holland, 1995; Khain et al., 2000; Prieto et al., 2003; Wu et al., 2003; Yang et al.,2008; Xu et al., 2011, 2013; Wu et al., 2012; Baumann et al., 2015; Liu and Tan, 2016; Xian and Chen, 2019).However, few studies have investigated the definition of BTCs and the criteria determining their formation.

    There are three existing definitions of BTCs in the WNP. Brand (1970) used the following criteria to determine a valid binary system: (1) the TCs approach within 700 nautical miles (~1300 km) of each other at one point in time; (2) both storms occur over the open ocean; and (3) the storms must reach the intensity of a typhoon at one point in time. Moreover, Dong and Neumann (1983) defined a BTC as two TCs that: (1) have coexisted for at least 48 h; (2)were separated at some time by < 1334 km; and (3) attained at least tropical storm (≥ 18 m s?1) status. Meanwhile, a BTC has been defined by Wu et al. (2011) and Jang and Chun (2016) as two TCs that: (1) had attained at least tropical storm status; (2) were separated by a minimum distance < 1600 km; and (3) had coexisted for at least 48 h.However, although all these definitions use criteria based on the separation distance, TC intensity and the coexistence time, they are different from each other, and especially there does not exist any explanation or reason about the origins of those definitions. This means that the definitions are all mainly subjective and a generally accepted standard does not exist. It is not difficult to understand that this situation is very unfavorable for further in-depth research and operational application in this area. Therefore, a scientific issue arises as to whether we can establish an objective standard to define BTCs.

    The motivation for this study is to develop an objective standard for defining BTCs in the WNP through carefully examining the observed interaction behavior between two coexisting TCs. We discuss the available data and our methods in section 2. Section 3 describes in detail how our standard is established. Using this standard, the main characteristics of BTCs are analyzed and presented in section 4. A summary and concluding remarks are given in the final section.

    2. Data and methods

    2.1. Data

    Different datasets show clear differences in the intensity and location of TCs in the WNP basin (Yu et al., 2007;Knapp and Kruk, 2009; Song et al., 2010; Ren et al., 2011).We therefore used two different best-track datasets—one from the China Meteorological Administration (CMA) and the other from the Joint Typhoon Warning Center(JTWC)—to develop an objective standard for defining BTCs. Both datasets contain information on the six-hourly track and intensity of TCs in the period 1951?2014. We also used the monthly mean National Centers for Environmental Prediction reanalysis geopotential height dataset with a spatial resolution of 2.5° × 2.5° for the same period.

    2.2. Methods

    The basic features of BTCs are “mutual counterclockwise spinning and moving closer to each other” (Fujiwhara,1921, 1923, 1931). We therefore designed a statistical analysis method to analyze the observations. Three parameters—the distance (d) between the centers of two coexisting TCs, the change in this distanceand the mutual angular velocity (r)—need to be calculated for any two coexisting TCs.

    The value of d can be approximately calculated according to the formula for the distance between two points on a plane. Δ d can be calculated as:

    where t is the number of records for two coexisting TCs and dtand dt?1represent the distances at time t and time t?1,respectively.< 0 means that the two TCs are approaching each other, whereas> 0 means that they are moving apart.

    Figure 1 shows the schematic for calculating the mutual angular velocity between the two coexisting TCs. The value of r can be calculated as:

    Equations (1) and (2) show that calculations can only be performed when t ? 1 ≥ 1 (i.e. t ≥ 2); that is, the duration T ≥ 12 h (12 = 2 × 6). Therefore, the duration of coexistence T should be at least 12 h and the three parameters only make sense when t ≥ 2.

    Fig. 1. Schematic for calculating the mutual angular velocity between two coexisting TCs.

    3. An objective standard for BTCs

    Considering that the basic features of BTCs are“mutual counterclockwise spinning and moving closer to each other”, we focus on the behavior of interaction between two coexisting TCs. Figure 2 presents frequency?distance distributions of two-TC coexistence for the CMA and the JTWC TC datasets. It can be seen that the frequencies in the two datasets both show a unimodal distribution with a peak at 2100 km, while the frequency in the JTWC dataset shows a weak secondary peak at 1900 km. In addition, as the CMA dataset contains more independent TCs than the JTWC dataset (Ren et al., 2011), the frequency of two coexisting TCs in the CMA dataset is greater than that in the JTWC dataset.

    Figure 3 shows the ratio?distance distributions of TCs that are moving apart and approaching each other for the two datasets. The CMA dataset (Fig. 3a) shows that the ratio for two TCs approaching each other increases as the distance decreases, with a ratio of 0.387 at 3000 km and 0.701 at 500 km. By contrast, the ratio for two TCs that are moving apart decreases as the distance decreases. The two ratio lines intersect over a range from 1800 to 2300 km, in which the ratios for both TCs approaching each other and TCs that are moving apart oscillate around 0.5, which means that the probabilities of occurrences for approaching and escaping are equivalent. The JTWC dataset (Fig. 3b) shows similar features with the same range of intersection. The range of intersection means that within a separation distance of 1800 km,the ratio for TCs approaching each other is > 0.5, which means that the probability of occurrence for approaching becomes dominant and increases as the distance decreases.

    Fig. 2. Frequency?distance distributions of two-TC coexistence for the CMA and the JTWC TC datasets. Here, “distance” is the amount of space between two TC centers, while“frequency” means the number of TC pairs. The distance interval for the statistics is 100 km,with the maximum value representing the interval, e.g., 1500 for (1400, 1500].

    Fig. 3. Approaching and escaping ratio?distance distributions of two TCs coexisting in the CMA and the JTWC datasets. Here, “ratio” means the ratio of number of coexisting TC pairs approaching or escaping to the total number of coexisting TC pairs. The distance interval for the statistics is 100 km, with the maximum value representing the interval, e.g., 1500 for (1400,1500]. (a) CMA dataset. (b) JTWC dataset.

    Fig. 4. As in Fig. 3 but for counterclockwise and clockwise ratio?distance distributions. Here, “ratio” indicates the ratio of number of coexisting TC pairs with counterclockwise rotation or clockwise rotation to the total number of coexisting TC pairs. (a) CMA dataset. (b) JTWC dataset.

    Figure 4 depicts the counterclockwise and clockwise ratio?distance distributions and can be used to examine the mutual angular velocity between two coexisting TCs. The JTWC dataset (Fig. 4b) shows that the counterclockwise ratio increases as the distance decreases, with a ratio of 0.416 at 3000 km and 0.846 at 500 km. By contrast, the clockwise ratio decreases as the distance decreases, with a ratio of 0.558 at 3000 km and 0.154 at 500 km. The counterclockwise and clockwise ratios intersect over a range from 1800 to 2000 km, in which both ratios oscillate around 0.5, which means that the probabilities of occurrences for counterclockwise and clockwise are equivalent. The CMA dataset (Fig.4a) also displays highly consistent features with the same range of intersection. The intersection range also means that within a separation distance of 1800 km, the counterclockwise ratio is > 0.5, which means that the probability of occurrence for counterclockwise becomes dominant and increases as the distance decreases.

    Based on these analyses, the two datasets both show that 1800 km is a key distance in defining BTCs: within a separation distance of 1800 km, the approaching ratio and counterclockwise ratio are both > 0.5 and increase as the distance decreases. As a result, 1800 km was selected as the threshold distance to define a BTC in Eq. (3). Taking into consideration the requirement for the duration T in Eq. (4), an objective standard for a BTC in the WNP can be established. If the distance between the centers of two coexisting TCs is ≤ d0(1800 km) and the duration is at least 12 h, then they can be defined as a BTC:

    The distance threshold, 1800 km, which is larger than those of the three existing definitions of BTCs, is a key parameter in the objective standard. The distance between the centers of the two TC centers defines 900 km as an important TC size. According to Chavas et al. (2016), the median and mean sizes of the outer radius of a TC in the WNP are 957.6 and 993.5 km, respectively. These three TC sizes are clearly of the same order of magnitude. This suggests that the distance threshold, 1800 km, can be understood to be twice the mean outer radius within which two TCs easily interact with each other and show clear characteristics of BTCs.

    Our results are consistent with previous studies (Brand,1970; Dong, 1980, 1981; Dong and Neumann, 1983; Kuo et al., 2000; Wu et al., 2011) in which the significance of the Fujiwara effect in a pair of TCs is dependent on the distance between the centers of two TCs. Other factors include the size, structure, and the intensity of the two TCs, as well as the environmental steering currents.

    ?

    4. Main characteristics of BTCs under the objective standard

    Under the objective standard defined by Eqs. (3) and(4) and the secondary standard of condition 2, the JTWC dataset was adopted to analyze the characteristics of BTCs. As shown in Table 1, there are 505 pairs of BTCs in the WNP during the period 1951?2014, with 328 pairs of typical BTCs and 177 pairs of atypical BTCs accounting for 65.0%and 35.0% of the total, respectively.

    We considered two statistical variables in detail: the duration of the BTC and the minimum distance between the centers of the two TCs during the duration of the BTC. The duration varied from 12 to 264 h and Fig. 5a shows the variations in the ratios for the duration of the two types of BTC.The ratio for a typical BTC is smaller than that for an atypical BTC only when the duration is < 24 h. The ratio is 0.28 at 12 h, increases rapidly between 12 and 84 h, and then varies in the range 0.80?1.0, with the value always being 1.0 when the duration of the BTC is ≥ 144 h. The minimum distance changes from 133.5 to 1800 km and Fig. 5b presents the variation in the ratios at the minimum distance for the two types of BTC. The ratio for a typical BTC increases as the minimum distance decreases and is smaller than that of an atypical BTC when the minimum distance is between 1700 and 1800 km. The ratio is 0.37 at 1800 km and is larger than that of an atypical BTC when the minimum distance is < 1700 km, with a value of 1.0 at ≤ 300 km.

    Fig. 5. Variations of rates of different variables for the two types of BTC over the WNP during the period 1951?2014: (a) duration of BTCs; (b) minimum distance between the two TCs during the duration of BTCs.

    Fig. 6. Annual variations of frequencies for the two types of BTC over the WNP during the period 1951?2014.

    Figure 6 depicts annual variations of frequencies for the two types of BTC over the WNP during the period 1951?2014. During 1951?2014, the frequency of typical BTC mainly shows an obvious interdecadal variation, with large values during the 1990s and the maximum 16 in 1994,but values of zero in 1952 and 2014. Meanwhile, the frequency of atypical BTCs also indicates a similar interdecadal variation. The above interdecadal variations with large values during the 1990s in the two types of BTC are consistent with the variations in TC number and intensity in the JTWC dataset during that period (Ren et al., 2011). Ren et al. (2011) pointed out that the TC number and intensity in the JTWC dataset seems to be overestimated since 1988, particularly from 1993 to 2003, with the possible factors being as follows: (i) no TC intensity references with high accuracy for operational centers after aircraft reconnaissance terminated in 1987; (ii) obvious differences in applications of the Dvorak technique in different operational typhoon/hurricane centers; and (iii) differences in time intervals for the maximum sustained wind speed (MSWS) at different cen-ters, especially when in-situ observation was available.

    Figure 7 presents the seasonal variations of frequencies for the total and the two types of BTCs. It clearly shows that about 93.5% occur during June to November, with the most active season between July and September, but especially August, in which 67.3% occur. Meanwhile, BTCs rarely appear in May, December and January, while no BTC occurs from February to April.

    For each BTC pair, there might exist two (different) numbers of times for approaching each other or counterclockwise rotation. Among each of the two numbers, two extreme phenomena are of particular interest: the maximum approaching speed (extreme phenomenon 1) and the maximum counterclockwise angular velocity (extreme phenomenon 2) between the two TCs. To study the characteristics of typical BTCs with these two extreme phenomena during the active season from July to September, a statistic was derived for the direction of the line segment between the two TCs. It is easy to understand thatθ, as defined in section 2, satisfies the relation 180° >θ≥ 0°. Six directions are applied in the statistics: direction 1 (30° >θ≥ 0°), direction 2 (60° >θ≥ 30°), direction 3 (90° >θ≥ 60°), direction 4 (120° >θ≥ 90°), direction 5 (150° >θ≥ 120°), and direction 6 (180° >θ≥ 150°). The results of the above statistics show that, with extreme phenomenon 1, direction 1(35.0%) and direction 6 (28.3%) are the two most frequent directions, with a total frequency ratio of 63.3%. By contrast, with extreme phenomenon 2, direction 1 (37.2%) and direction 2 (17.9%) are the two most frequent directions,with a total frequency ratio of 55.1%.

    Figure 8a shows the distributions of the BTC locations in the two most frequent directions, their mean position with extreme phenomenon 1, and the average position of subtropical highs during the period July?September. The mean position of BTCs shows a west?east direction, with the center of the western TC at (21.8°N, 127.0°E) and the center of the eastern TC at (22.0°N, 142.1°E), with a separation distance of 1557.3 km. The average position of the subtropical high clearly benefits the formation of extreme phenomenon 1. With these phenomena, the BTC is mainly controlled by an easterly stream along the southern edge of the subtropical high and the eastern TC moves rapidly toward the western TC under the guiding easterly airflow.

    Figure 8b depicts similar characteristics to Fig. 8a, but for extreme phenomenon 2. The mean position of the BTC shows movement from the northeast to southwest, with the center of the western TC locating at (21.3°N, 127.1°E) and the center of the eastern TC locating at (25.6°N, 139.0°E),with a separation distance of 1303.9 km. The location of the subtropical high clearly benefits the formation of extreme phenomenon 2. At the moment, the BTC is mainly distributed in a northeast to southwest direction, with the eastern TC locating at the southwestern edge of the subtropical high and being controlled by the southeasterly stream, thus favoring counterclockwise spinning between the BTCs.

    In order to compare the differences of the large-scale weather systems associated with the two type BTCs, representative samples for typical and atypical BTCs were chosen from the BTCs with the intensity of the stronger TC ≥32.7 m s?1and that of the weaker TC < 32.7 m s?1, as follows: (i) typical BTCs: the first moment when BTCs simultaneously show approaching and spinning counterclockwise for the typical BTCs in the two most frequent directions during July?September; and (ii) atypical BTCs: the moment when BTCs are at the minimum distance during the duration of the atypical BTCs in the two most frequent directions during July?September. Figure 9 shows the distributions of the locations of the representative samples for typical and atypical BTCs. Figure 9a shows that the mean position of representative typical BTCs shows a slightly northeast?southwest direction, with the center of the western TC at (20.9°N, 127.8°E), the center of the eastern TC at(21.5°N, 142.5°E), and a separation distance of 1521.3 km.Based on the analyses of Fig. 8 it is easy to understand that the configuration of the subtropical high and the representative typical BTCs clearly benefits the occurrence of BTCs sim-ultaneously approaching and spinning counterclockwise.Meanwhile, Fig. 9b indicates that the mean position of representative atypical BTCs shows an obvious northwest?southeast direction, with the western TC at (27.7°N, 128.1°E),the center of the eastern TC at (20.2°N, 139.8°E), and a separation distance of 1452.2 km. Taking into consideration that the western TC generally locates in the monsoon trough (figure omitted), which is mainly west of the mean position of representative atypical BTCs, the configuration of the subtropical high, the monsoon trough and the representative atypical BTCs does not benefit the occurrence of BTCs approaching or spinning counterclockwise.

    Fig. 7. Seasonal variations of frequencies for the total and the two types of BTCs over the WNP during the period 1951?2014.

    Fig. 8. Positions of typical BTCs (two blue dots connected by gray line segments) in the two most frequent directions and their mean position (the two red dots connected by a red line segment) with two different extreme phenomena during the duration of BTCs and the average position of the subtropical high during the period of frequent BTCs (July?September) over the WNP from 1951 to 2014 using the monthly mean geopotential height from the NCEP reanalysis dataset: (a) maximum speed of approaching toward each other; (b) maximum counterclockwise angular velocity.

    5. Summary and concluding remarks

    The following conclusions can be drawn based on our analyses.

    An objective standard to define BTCs has been established based on statistical observations of TCs in the WNP.The objective standard includes a main standard for defining BTCs and a secondary standard for defining typical/atypical BTCs. The main standard has two requirements: two coexisting TCs are a pair of BTCs if (i) the separation distance is ≤ 1800 km, and (ii) this separation maintains for at least 12 h. The secondary standard defines a typical BTC as one for which there is at least one observation when the two TCs approach each other and spin counterclockwise simultaneously; otherwise they are atypical BTCs.

    Fig. 9. As in Fig. 8 but for representative samples for (a) typical and (b) atypical BTCs in the two most frequent directions and their mean positions. Typical BTCs are selected at the first moment when BTCs simultaneously show approaching and spinning counterclockwise during July?September, while atypical BTCs are chosen at the moment when BTCs are at the minimum distance during the same period. The subtropical high is the average position of subtropical highs at corresponding moments.

    There are 505 pairs of BTCs during the period 1951?2014, including 328 typical BTCs and 177 atypical BTCs, accounting for 65.0% and 35.0% of the total, respectively. The duration of the BTCs varies from 12 to 264 h and the minimum distance is in the range 133.5?1800 km. The ratio of typical BTCs increases as the duration increases,with the value always being 1.0 when the duration is ≥ 144 h. The ratio increases as the minimum distance decreases,with the value always being 1.0 at ≤ 300 km.

    A study of the two extreme phenomena for typical BTCs shows that the configurations of the circulation conditions and the distribution of the BTC favor the formation of the extreme phenomena. At the time of the maximum approaching speed, when the BTCs are mainly located along the southern edge of the subtropical high, the movement of the eastern TC toward the western TC can be increased under the easterly steering airflow. At the time of the maximum counterclockwise angular velocity, the BTCs are mainly distributed in a northeast?southwest direction,with the eastern TC locating at the southwestern edge of the subtropical high and being controlled by the southeasterly stream. This favors anticlockwise spinning between the two TCs.

    In order to understand more clearly the merits of our BTC-identification standard compared to the previous BTC definitions, the differences can be summarized in two aspects: First, our standard is objective, while the previous BTC definitions are mainly subjective. Our standard is based on statistical analyses of two best-track datasets,CMA and JTWC, with the separation distance of ≤ 1800 km between the centers of two TCs equating to two times the average size of the TC outer radius in the WNP. Meanwhile, the previous BTC definitions do not provide any explanation or reason regarding the origins of the definitions. Second, our standard has a secondary standard for identifying typical/atypical BTCs, whereas the previous BTC definitions do not. Under our standard, a smaller minimum distance between the two TCs or a longer duration of interaction gives a larger proportion of typical BTCs. Meanwhile, it is notable that about 35.0% (= 177 / 505) of total BTCs are atypical BTCs that do not have any moment where they are simultaneously approaching each other and spinning counterclockwise, especially at short distances with a certain ratio of BTCs moving apart and rotating clockwise. These situations might indicate the environmental effects on BTC behavior.

    A number of new issues arise from this study. For example, are there any differences between our standard in the WNP and standards in other TC basins in the Southern Hemisphere? BTCs often produce extreme rainfall events resulting in serious damage—for example, Typhoon Morakot in Taiwan in 2009 and Typhoon Fitow in Zhejiang in 2013 in China were both the result of BTCs. However, we do not have a clear picture of the characteristics of BTC interactions in China’s offshore area. These issues require further in-depth studies.

    Acknowledgements. The authors would like to express their sincere thanks to the three anonymous reviewers and the editor for their helpful suggestions and comments to improve the manuscript. This work was supported by the National Natural Science Foundation of China (Grant No. 41675042) and the Jiangsu Collaborative Innovation Center for Climate Change. We acknowledge the use of the best-track datasets from the CMA and JTWC,and the monthly mean geopotential height from the NCEP reanalysis dataset.

    97精品久久久久久久久久精品| 国产精品99久久99久久久不卡| 老司机靠b影院| 丁香六月天网| 亚洲中文av在线| 大香蕉久久成人网| 两个人看的免费小视频| 国产黄色视频一区二区在线观看| 另类亚洲欧美激情| 一级片免费观看大全| 久久久久视频综合| 亚洲精品日本国产第一区| 久久精品国产a三级三级三级| 午夜免费观看性视频| 久久人妻熟女aⅴ| 国产成人91sexporn| 久久 成人 亚洲| 老汉色∧v一级毛片| 尾随美女入室| 久久久亚洲精品成人影院| 777米奇影视久久| 国产在视频线精品| 精品国产一区二区三区四区第35| 日韩一本色道免费dvd| 亚洲精品国产一区二区精华液| 亚洲自偷自拍图片 自拍| 亚洲国产精品国产精品| 肉色欧美久久久久久久蜜桃| 亚洲免费av在线视频| 色播在线永久视频| tube8黄色片| 精品国产乱码久久久久久男人| 亚洲国产av影院在线观看| 黄片播放在线免费| 久久久久久人人人人人| 人妻一区二区av| 18在线观看网站| 激情五月婷婷亚洲| 亚洲精品在线美女| 亚洲精品美女久久av网站| 嫁个100分男人电影在线观看 | 久久99热这里只频精品6学生| 亚洲伊人久久精品综合| 日本av免费视频播放| 制服诱惑二区| 国产成人欧美在线观看 | 国产一区亚洲一区在线观看| 国产视频首页在线观看| 国产极品粉嫩免费观看在线| 又大又黄又爽视频免费| 婷婷丁香在线五月| 黄色a级毛片大全视频| 一级毛片我不卡| 久热爱精品视频在线9| 三上悠亚av全集在线观看| 亚洲国产毛片av蜜桃av| 亚洲第一青青草原| 一级毛片黄色毛片免费观看视频| 性色av乱码一区二区三区2| 免费观看a级毛片全部| 久久久久久久久免费视频了| 91字幕亚洲| 90打野战视频偷拍视频| 国产精品人妻久久久影院| 国产精品久久久久成人av| 国产爽快片一区二区三区| 精品一区二区三卡| 免费日韩欧美在线观看| 亚洲欧洲精品一区二区精品久久久| 欧美精品av麻豆av| 男女床上黄色一级片免费看| 国产欧美日韩精品亚洲av| 女性生殖器流出的白浆| 成人黄色视频免费在线看| 欧美日韩视频高清一区二区三区二| 久久久久久免费高清国产稀缺| www.999成人在线观看| 午夜av观看不卡| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲高清精品| 天天添夜夜摸| 日本wwww免费看| 国产1区2区3区精品| 午夜精品国产一区二区电影| 国产有黄有色有爽视频| 国产精品三级大全| 亚洲专区国产一区二区| 国产精品久久久人人做人人爽| 亚洲激情五月婷婷啪啪| 欧美大码av| 亚洲成人免费电影在线观看 | 2018国产大陆天天弄谢| 亚洲三区欧美一区| 性高湖久久久久久久久免费观看| 18禁裸乳无遮挡动漫免费视频| 无限看片的www在线观看| 少妇精品久久久久久久| 男女之事视频高清在线观看 | 中文字幕制服av| 国产精品九九99| 精品国产一区二区久久| 一本久久精品| 人妻 亚洲 视频| 久久精品成人免费网站| 高清不卡的av网站| 极品少妇高潮喷水抽搐| 亚洲黑人精品在线| 国产一区二区 视频在线| 国产无遮挡羞羞视频在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲国产av影院在线观看| 夜夜骑夜夜射夜夜干| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 高潮久久久久久久久久久不卡| 国产高清视频在线播放一区 | 少妇人妻 视频| 天天躁日日躁夜夜躁夜夜| 国产在线观看jvid| 久久99精品国语久久久| 精品高清国产在线一区| av欧美777| 久久久亚洲精品成人影院| 十八禁高潮呻吟视频| 国产欧美日韩一区二区三 | 日本一区二区免费在线视频| 一本久久精品| 欧美激情 高清一区二区三区| 后天国语完整版免费观看| 亚洲 欧美一区二区三区| 亚洲av成人不卡在线观看播放网 | 首页视频小说图片口味搜索 | 啦啦啦 在线观看视频| 精品国产一区二区三区久久久樱花| 国产精品国产三级国产专区5o| 黄色视频在线播放观看不卡| 日韩视频在线欧美| 成在线人永久免费视频| 中文欧美无线码| 2018国产大陆天天弄谢| 国产99久久九九免费精品| 亚洲自偷自拍图片 自拍| 人人妻,人人澡人人爽秒播 | 精品人妻熟女毛片av久久网站| 母亲3免费完整高清在线观看| 中国国产av一级| 飞空精品影院首页| 亚洲七黄色美女视频| 国产精品 欧美亚洲| 亚洲精品国产区一区二| 成年美女黄网站色视频大全免费| 两个人看的免费小视频| 自线自在国产av| 你懂的网址亚洲精品在线观看| 国产精品三级大全| 欧美+亚洲+日韩+国产| 国产伦人伦偷精品视频| 国产成人av激情在线播放| 男人舔女人的私密视频| 电影成人av| 人人妻人人添人人爽欧美一区卜| 波多野结衣av一区二区av| 天天操日日干夜夜撸| 中国国产av一级| 久久人人97超碰香蕉20202| 免费观看人在逋| 欧美精品av麻豆av| tube8黄色片| 日韩大码丰满熟妇| 咕卡用的链子| 亚洲成人国产一区在线观看 | 亚洲欧美精品综合一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产精品亚洲av一区麻豆| 久久国产精品影院| 大陆偷拍与自拍| www.精华液| 久久久精品区二区三区| 亚洲欧美精品自产自拍| 久久 成人 亚洲| 国产真人三级小视频在线观看| 一区二区三区精品91| 精品少妇内射三级| 1024香蕉在线观看| 国产一区二区在线观看av| 欧美日韩精品网址| 性高湖久久久久久久久免费观看| av天堂久久9| 伊人久久大香线蕉亚洲五| 51午夜福利影视在线观看| 久久ye,这里只有精品| www.自偷自拍.com| 免费人妻精品一区二区三区视频| 热99国产精品久久久久久7| 黄片小视频在线播放| 色婷婷av一区二区三区视频| 久久久久久免费高清国产稀缺| 亚洲九九香蕉| 亚洲精品av麻豆狂野| 巨乳人妻的诱惑在线观看| 日本av免费视频播放| 又黄又粗又硬又大视频| 岛国毛片在线播放| 国产日韩一区二区三区精品不卡| 日本欧美视频一区| 免费在线观看黄色视频的| 美女国产高潮福利片在线看| 男女床上黄色一级片免费看| 午夜久久久在线观看| 飞空精品影院首页| 十八禁人妻一区二区| 99九九在线精品视频| 免费看不卡的av| 极品人妻少妇av视频| 日韩 亚洲 欧美在线| 国产一区二区激情短视频 | 狂野欧美激情性bbbbbb| 久久女婷五月综合色啪小说| 50天的宝宝边吃奶边哭怎么回事| 精品一品国产午夜福利视频| 99热网站在线观看| 国产欧美日韩一区二区三区在线| 欧美xxⅹ黑人| 99久久综合免费| 久久精品熟女亚洲av麻豆精品| 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲 | 久久久久久久大尺度免费视频| 天堂中文最新版在线下载| 50天的宝宝边吃奶边哭怎么回事| 国产午夜精品一二区理论片| 久久久久视频综合| av线在线观看网站| 中文字幕人妻熟女乱码| 亚洲精品第二区| 亚洲精品美女久久久久99蜜臀 | 人妻一区二区av| av又黄又爽大尺度在线免费看| 黄色 视频免费看| 交换朋友夫妻互换小说| 中文字幕精品免费在线观看视频| 人人妻人人爽人人添夜夜欢视频| 欧美另类一区| 亚洲国产精品国产精品| 婷婷色综合大香蕉| 亚洲三区欧美一区| 色精品久久人妻99蜜桃| 看免费av毛片| 亚洲伊人色综图| 国产精品.久久久| 免费高清在线观看视频在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美日韩国产mv在线观看视频| 亚洲精品一区蜜桃| 美女脱内裤让男人舔精品视频| 18禁裸乳无遮挡动漫免费视频| 免费日韩欧美在线观看| 久久狼人影院| 天天操日日干夜夜撸| 激情五月婷婷亚洲| 亚洲三区欧美一区| 亚洲精品美女久久av网站| 久久久精品国产亚洲av高清涩受| 亚洲av片天天在线观看| av线在线观看网站| 精品免费久久久久久久清纯 | 老司机影院成人| 性高湖久久久久久久久免费观看| 久久人人爽av亚洲精品天堂| 国产一卡二卡三卡精品| 91麻豆av在线| av在线老鸭窝| 亚洲av男天堂| 成年av动漫网址| 国产黄频视频在线观看| 亚洲欧美一区二区三区国产| 在线看a的网站| 热re99久久国产66热| 亚洲av美国av| 国产亚洲午夜精品一区二区久久| 精品一区二区三卡| 日韩制服骚丝袜av| 中文字幕制服av| 嫁个100分男人电影在线观看 | 国产精品一国产av| 亚洲中文字幕日韩| 婷婷色麻豆天堂久久| av不卡在线播放| 五月开心婷婷网| 9色porny在线观看| 国产伦理片在线播放av一区| 制服人妻中文乱码| 国产精品久久久人人做人人爽| 中文字幕色久视频| 久久久国产欧美日韩av| 日韩一卡2卡3卡4卡2021年| 飞空精品影院首页| 日本色播在线视频| 国产精品久久久久久精品古装| 亚洲欧美一区二区三区黑人| 无限看片的www在线观看| 亚洲中文av在线| 我的亚洲天堂| 午夜福利,免费看| 夫妻午夜视频| 一区二区日韩欧美中文字幕| 午夜久久久在线观看| 欧美av亚洲av综合av国产av| 首页视频小说图片口味搜索 | 国产精品免费大片| 日韩精品免费视频一区二区三区| 又粗又硬又长又爽又黄的视频| 97在线人人人人妻| av片东京热男人的天堂| 亚洲国产精品国产精品| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 久久精品亚洲av国产电影网| 亚洲中文字幕日韩| 免费看十八禁软件| 久久人人爽人人片av| 男的添女的下面高潮视频| 婷婷色综合www| 精品久久久久久久毛片微露脸 | 91字幕亚洲| 亚洲国产欧美网| 在线观看一区二区三区激情| 如日韩欧美国产精品一区二区三区| 久久久久网色| 视频区欧美日本亚洲| 久久午夜综合久久蜜桃| 一个人免费看片子| 后天国语完整版免费观看| 狂野欧美激情性bbbbbb| 啦啦啦 在线观看视频| 精品一区二区三区四区五区乱码 | 好男人电影高清在线观看| 国产免费福利视频在线观看| 中文字幕最新亚洲高清| www.自偷自拍.com| 母亲3免费完整高清在线观看| 脱女人内裤的视频| 精品一区二区三卡| 中文字幕人妻熟女乱码| 五月天丁香电影| 国产精品二区激情视频| 亚洲免费av在线视频| 亚洲中文字幕日韩| 大陆偷拍与自拍| 欧美大码av| 丁香六月天网| 亚洲五月婷婷丁香| 一本综合久久免费| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久午夜乱码| 久久精品久久久久久久性| 免费少妇av软件| 亚洲精品国产av蜜桃| 欧美精品一区二区大全| 丰满少妇做爰视频| 国产欧美亚洲国产| 国产视频首页在线观看| 天天影视国产精品| 两性夫妻黄色片| 欧美日韩av久久| 一级毛片黄色毛片免费观看视频| 黄色毛片三级朝国网站| 久久久欧美国产精品| 一级片'在线观看视频| 亚洲少妇的诱惑av| 99香蕉大伊视频| 妹子高潮喷水视频| 国产福利在线免费观看视频| 免费一级毛片在线播放高清视频 | 国产国语露脸激情在线看| 免费看av在线观看网站| 男女床上黄色一级片免费看| 老司机午夜十八禁免费视频| 欧美日韩黄片免| 丝袜美腿诱惑在线| 国产精品一区二区免费欧美 | 日韩视频在线欧美| 成人国语在线视频| 亚洲av成人不卡在线观看播放网 | 久久免费观看电影| 我要看黄色一级片免费的| 亚洲成色77777| 我的亚洲天堂| 久久久欧美国产精品| 国产精品免费大片| av在线app专区| 最黄视频免费看| 亚洲中文字幕日韩| 国产有黄有色有爽视频| 免费一级毛片在线播放高清视频 | av一本久久久久| 亚洲第一av免费看| 免费在线观看黄色视频的| 秋霞在线观看毛片| 热99久久久久精品小说推荐| 久久av网站| 国产麻豆69| 亚洲精品第二区| 亚洲欧美精品自产自拍| 国产人伦9x9x在线观看| 亚洲少妇的诱惑av| 波野结衣二区三区在线| 人成视频在线观看免费观看| 国产在视频线精品| 国产成人免费无遮挡视频| 午夜免费观看性视频| 在线av久久热| 亚洲专区中文字幕在线| 免费看十八禁软件| 免费在线观看视频国产中文字幕亚洲 | 激情五月婷婷亚洲| 在线观看一区二区三区激情| 又大又黄又爽视频免费| 性色av一级| 黄频高清免费视频| 国产1区2区3区精品| 日日夜夜操网爽| 久久久久久免费高清国产稀缺| 欧美在线一区亚洲| 中文乱码字字幕精品一区二区三区| 欧美日韩综合久久久久久| 亚洲精品久久午夜乱码| 男人添女人高潮全过程视频| 国产精品秋霞免费鲁丝片| 热re99久久精品国产66热6| 一二三四社区在线视频社区8| 精品免费久久久久久久清纯 | 啦啦啦在线观看免费高清www| 久久久久精品国产欧美久久久 | 日日爽夜夜爽网站| avwww免费| 性少妇av在线| 少妇粗大呻吟视频| 国产黄色免费在线视频| 中文乱码字字幕精品一区二区三区| 亚洲成色77777| 脱女人内裤的视频| 亚洲成人国产一区在线观看 | avwww免费| 亚洲国产欧美在线一区| 丝瓜视频免费看黄片| 国产精品一国产av| 国产视频一区二区在线看| 免费女性裸体啪啪无遮挡网站| 精品卡一卡二卡四卡免费| 日本黄色日本黄色录像| 制服人妻中文乱码| 男人舔女人的私密视频| 精品人妻一区二区三区麻豆| 亚洲五月婷婷丁香| 男人爽女人下面视频在线观看| 欧美激情高清一区二区三区| 日本av手机在线免费观看| 国产一卡二卡三卡精品| 国产精品三级大全| 日韩视频在线欧美| 777米奇影视久久| 欧美日韩视频高清一区二区三区二| 在线观看免费视频网站a站| 日韩中文字幕视频在线看片| 精品人妻1区二区| 日韩一卡2卡3卡4卡2021年| 老熟女久久久| 两人在一起打扑克的视频| 99久久人妻综合| 精品久久久久久电影网| 成年女人毛片免费观看观看9 | 国产片特级美女逼逼视频| 丝袜美腿诱惑在线| 色播在线永久视频| 亚洲成av片中文字幕在线观看| 久久ye,这里只有精品| 欧美变态另类bdsm刘玥| 看免费av毛片| 久久人妻熟女aⅴ| 97人妻天天添夜夜摸| 欧美成狂野欧美在线观看| 日韩制服丝袜自拍偷拍| 日韩av免费高清视频| 国产伦理片在线播放av一区| 好男人视频免费观看在线| 天天操日日干夜夜撸| 另类精品久久| 一级毛片 在线播放| 考比视频在线观看| 少妇猛男粗大的猛烈进出视频| 久久人妻熟女aⅴ| 国产精品国产三级国产专区5o| 亚洲欧洲国产日韩| 男女无遮挡免费网站观看| 欧美久久黑人一区二区| 夫妻性生交免费视频一级片| 欧美精品亚洲一区二区| 中文字幕人妻熟女乱码| 男人添女人高潮全过程视频| 亚洲欧美日韩高清在线视频 | 大香蕉久久网| 欧美国产精品va在线观看不卡| 欧美日韩成人在线一区二区| 亚洲少妇的诱惑av| 久久ye,这里只有精品| 午夜福利免费观看在线| 黄色视频在线播放观看不卡| 免费在线观看黄色视频的| 国产真人三级小视频在线观看| 天堂8中文在线网| 精品国产乱码久久久久久男人| 99热全是精品| 日本vs欧美在线观看视频| h视频一区二区三区| 久久青草综合色| 免费女性裸体啪啪无遮挡网站| 亚洲成人免费电影在线观看 | 久热这里只有精品99| 中文字幕av电影在线播放| 久久精品国产亚洲av高清一级| 激情视频va一区二区三区| 七月丁香在线播放| 免费在线观看完整版高清| 一本久久精品| 国产精品偷伦视频观看了| 女人高潮潮喷娇喘18禁视频| 在线观看国产h片| 国产免费现黄频在线看| 国产精品一区二区精品视频观看| 国产国语露脸激情在线看| 欧美黑人精品巨大| 99国产精品99久久久久| 久久天躁狠狠躁夜夜2o2o | 国产精品免费大片| 国产精品.久久久| 欧美精品亚洲一区二区| 国产精品秋霞免费鲁丝片| 午夜福利一区二区在线看| 亚洲中文日韩欧美视频| 久久精品成人免费网站| av在线app专区| 亚洲国产精品一区三区| 精品国产乱码久久久久久男人| 亚洲中文字幕日韩| 一本色道久久久久久精品综合| 黄片小视频在线播放| 午夜福利影视在线免费观看| 在线观看免费日韩欧美大片| 久久久国产欧美日韩av| 校园人妻丝袜中文字幕| 天天操日日干夜夜撸| 好男人电影高清在线观看| 男女高潮啪啪啪动态图| 最近中文字幕2019免费版| 日韩一区二区三区影片| 亚洲av男天堂| 中文字幕人妻丝袜一区二区| 免费日韩欧美在线观看| 欧美日韩福利视频一区二区| 80岁老熟妇乱子伦牲交| 久久精品国产综合久久久| 日韩大码丰满熟妇| 波多野结衣av一区二区av| 色视频在线一区二区三区| 久久精品国产亚洲av高清一级| 欧美日韩亚洲高清精品| 多毛熟女@视频| 黑丝袜美女国产一区| 久久久精品区二区三区| 国产男人的电影天堂91| 五月开心婷婷网| 亚洲五月色婷婷综合| 免费不卡黄色视频| 亚洲成av片中文字幕在线观看| 亚洲精品日本国产第一区| 欧美日韩福利视频一区二区| 欧美变态另类bdsm刘玥| 日韩免费高清中文字幕av| 男人操女人黄网站| 在线观看免费视频网站a站| 午夜日韩欧美国产| 国产精品成人在线| 亚洲国产看品久久| 亚洲av片天天在线观看| 青草久久国产| 久久国产精品男人的天堂亚洲| 国产成人免费无遮挡视频| 又大又爽又粗| 在线观看一区二区三区激情| 黄色毛片三级朝国网站| 婷婷色麻豆天堂久久| 亚洲国产毛片av蜜桃av| av国产精品久久久久影院| 老鸭窝网址在线观看| 99久久人妻综合| 国产成人av激情在线播放| 亚洲av欧美aⅴ国产| 日本五十路高清| 黄色怎么调成土黄色| 伊人亚洲综合成人网| 男女午夜视频在线观看| 黄色怎么调成土黄色| 久久久亚洲精品成人影院| 日本91视频免费播放| 一级毛片女人18水好多 | 久久久久久人人人人人| 国产视频首页在线观看| 日本a在线网址| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产精品一区三区| 纵有疾风起免费观看全集完整版| 99久久人妻综合| 校园人妻丝袜中文字幕| 日韩av在线免费看完整版不卡|