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

    日本-黄色视频高清免费观看| 国产色婷婷99| 久久草成人影院| 又黄又爽又刺激的免费视频.| 欧美人与善性xxx| 国产爱豆传媒在线观看| 久久精品久久久久久噜噜老黄 | 国产熟女欧美一区二区| 一级a爱片免费观看的视频| 99在线人妻在线中文字幕| 午夜老司机福利剧场| 99久国产av精品| 青春草视频在线免费观看| 免费av不卡在线播放| av国产免费在线观看| 激情 狠狠 欧美| 亚洲欧美日韩高清专用| 国产精品一区二区免费欧美| av.在线天堂| 亚洲欧美成人综合另类久久久 | 国产一区二区在线av高清观看| 亚洲欧美清纯卡通| 99热这里只有精品一区| 欧美最新免费一区二区三区| 精品国产三级普通话版| 亚洲五月天丁香| 特大巨黑吊av在线直播| 特级一级黄色大片| 久久精品国产自在天天线| 男女之事视频高清在线观看| 晚上一个人看的免费电影| 国产精品亚洲美女久久久| 男人舔奶头视频| 精品久久国产蜜桃| 少妇人妻精品综合一区二区 | 久久精品国产亚洲av涩爱 | 男女之事视频高清在线观看| 欧美高清性xxxxhd video| 精品少妇黑人巨大在线播放 | 悠悠久久av| 亚洲熟妇中文字幕五十中出| 一级黄片播放器| 99久久精品一区二区三区| 国产亚洲精品久久久久久毛片| 国内揄拍国产精品人妻在线| 午夜老司机福利剧场| 国产大屁股一区二区在线视频| 三级经典国产精品| 亚州av有码| 国产精品久久视频播放| 亚洲天堂国产精品一区在线| 永久网站在线| 国产av不卡久久| 高清毛片免费观看视频网站| 中文字幕久久专区| 丰满人妻一区二区三区视频av| 99九九线精品视频在线观看视频| 特大巨黑吊av在线直播| 国产亚洲精品久久久久久毛片| 久久精品国产清高在天天线| 禁无遮挡网站| 波多野结衣巨乳人妻| 不卡一级毛片| 国产一区二区激情短视频| 色播亚洲综合网| 在现免费观看毛片| 免费人成在线观看视频色| 91在线观看av| 亚洲精品在线观看二区| 欧美zozozo另类| 女的被弄到高潮叫床怎么办| 欧美区成人在线视频| 久久中文看片网| 国产成人影院久久av| 精品久久久久久久久久久久久| 变态另类成人亚洲欧美熟女| 麻豆国产97在线/欧美| av女优亚洲男人天堂| 人妻丰满熟妇av一区二区三区| 国产女主播在线喷水免费视频网站 | 欧美成人一区二区免费高清观看| 精品人妻一区二区三区麻豆 | 国产黄色小视频在线观看| 欧美丝袜亚洲另类| 搞女人的毛片| 日本 av在线| 少妇猛男粗大的猛烈进出视频 | 在线观看免费视频日本深夜| 国产成人91sexporn| 久久久久国内视频| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| 麻豆国产97在线/欧美| 美女 人体艺术 gogo| 18+在线观看网站| av中文乱码字幕在线| 国产精品女同一区二区软件| 欧洲精品卡2卡3卡4卡5卡区| 波野结衣二区三区在线| 五月玫瑰六月丁香| 国产日本99.免费观看| 91在线观看av| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清专用| 国产精品亚洲美女久久久| 黄片wwwwww| 少妇猛男粗大的猛烈进出视频 | 国内精品美女久久久久久| 国产欧美日韩精品亚洲av| 人妻少妇偷人精品九色| 国产精品乱码一区二三区的特点| 久久久久久久久久黄片| 麻豆乱淫一区二区| 伦理电影大哥的女人| 亚洲美女视频黄频| 久久国内精品自在自线图片| 中国美白少妇内射xxxbb| 成人三级黄色视频| 最近手机中文字幕大全| 国产亚洲精品av在线| 亚洲自拍偷在线| 一区二区三区四区激情视频 | 成人亚洲欧美一区二区av| 五月玫瑰六月丁香| 成年免费大片在线观看| 婷婷精品国产亚洲av| 精品一区二区免费观看| 亚洲va在线va天堂va国产| 亚洲在线观看片| 老熟妇仑乱视频hdxx| 国产精品爽爽va在线观看网站| 一级av片app| 婷婷亚洲欧美| 亚洲av成人av| 亚洲精品456在线播放app| 色综合色国产| 亚洲18禁久久av| 精品人妻熟女av久视频| 婷婷亚洲欧美| 97在线视频观看| 国产精品一二三区在线看| 伊人久久精品亚洲午夜| 在线免费观看的www视频| 久99久视频精品免费| 欧美+日韩+精品| 成人国产麻豆网| 日本一二三区视频观看| 亚洲成av人片在线播放无| 日本色播在线视频| 1024手机看黄色片| 人妻制服诱惑在线中文字幕| 亚洲国产欧洲综合997久久,| 亚洲av不卡在线观看| 插逼视频在线观看| 99九九线精品视频在线观看视频| 国产高清三级在线| 无遮挡黄片免费观看| 亚洲国产精品成人综合色| 丝袜喷水一区| 国产女主播在线喷水免费视频网站 | 国产精品久久久久久亚洲av鲁大| 国产亚洲精品久久久com| 乱系列少妇在线播放| 啦啦啦韩国在线观看视频| 日本五十路高清| 亚洲欧美日韩卡通动漫| 国产探花在线观看一区二区| 亚洲国产精品sss在线观看| 最后的刺客免费高清国语| 欧美最黄视频在线播放免费| 日本与韩国留学比较| 夜夜爽天天搞| 午夜激情欧美在线| 熟妇人妻久久中文字幕3abv| 国产亚洲91精品色在线| 亚洲av电影不卡..在线观看| 久久精品91蜜桃| 露出奶头的视频| 成人亚洲精品av一区二区| 中文亚洲av片在线观看爽| 欧美一区二区国产精品久久精品| 日本三级黄在线观看| 国产精品一区二区性色av| 嫩草影院精品99| 亚洲久久久久久中文字幕| 欧美激情久久久久久爽电影| 日本 av在线| 最近在线观看免费完整版| 天堂影院成人在线观看| 国产午夜福利久久久久久| 91久久精品国产一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲乱码一区二区免费版| 一区二区三区免费毛片| 日韩在线高清观看一区二区三区| 黑人高潮一二区| 国产人妻一区二区三区在| 黄色配什么色好看| 欧美性猛交黑人性爽| 亚洲久久久久久中文字幕| 国产爱豆传媒在线观看| 最近手机中文字幕大全| 欧美不卡视频在线免费观看| 丰满乱子伦码专区| 一本精品99久久精品77| 美女免费视频网站| aaaaa片日本免费| av天堂中文字幕网| 亚洲av一区综合| 国产aⅴ精品一区二区三区波| av在线观看视频网站免费| 亚洲婷婷狠狠爱综合网| 亚洲,欧美,日韩| 免费电影在线观看免费观看| 极品教师在线视频| 久久精品人妻少妇| 精品久久久久久久末码| 少妇的逼水好多| 久久精品夜夜夜夜夜久久蜜豆| 国产美女午夜福利| 亚洲第一区二区三区不卡| or卡值多少钱| 国产欧美日韩精品一区二区| 精品一区二区免费观看| 成熟少妇高潮喷水视频| 欧美日本亚洲视频在线播放| 国产美女午夜福利| 国产精品一区二区性色av| 女人十人毛片免费观看3o分钟| 青春草视频在线免费观看| 在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 在线a可以看的网站| 日韩强制内射视频| 国产 一区 欧美 日韩| 人人妻人人澡欧美一区二区| 国产成人aa在线观看| 看十八女毛片水多多多| 搡老妇女老女人老熟妇| 亚洲中文日韩欧美视频| 国产高清视频在线观看网站| 日韩欧美国产在线观看| 少妇被粗大猛烈的视频| a级毛色黄片| 亚洲电影在线观看av| 两个人视频免费观看高清| 看片在线看免费视频| 国产伦精品一区二区三区视频9| 热99re8久久精品国产| 变态另类成人亚洲欧美熟女| 色播亚洲综合网| 在线播放无遮挡| 丰满人妻一区二区三区视频av| 免费看日本二区| 国内少妇人妻偷人精品xxx网站| 一边摸一边抽搐一进一小说| 国产白丝娇喘喷水9色精品| 日韩强制内射视频| 麻豆一二三区av精品| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲国产欧洲综合997久久,| 亚洲人成网站在线播| 嫩草影视91久久| 亚洲欧美清纯卡通| 国产老妇女一区| 69人妻影院| 亚洲中文日韩欧美视频| 亚洲国产欧美人成| 黄色欧美视频在线观看| 一进一出好大好爽视频| 欧美国产日韩亚洲一区| 久久精品夜夜夜夜夜久久蜜豆| 久久久久性生活片| 一进一出抽搐gif免费好疼| 久久天躁狠狠躁夜夜2o2o| 日日摸夜夜添夜夜添av毛片| 综合色av麻豆| 国产在视频线在精品| 3wmmmm亚洲av在线观看| 久久99热6这里只有精品| 日韩欧美 国产精品| 在线观看午夜福利视频| 色综合色国产| 丝袜美腿在线中文| 久久精品国产亚洲av涩爱 | 韩国av在线不卡| 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看| 免费av不卡在线播放| 日日干狠狠操夜夜爽| 日日摸夜夜添夜夜添av毛片| 中国国产av一级| 99精品在免费线老司机午夜| 久久99热这里只有精品18| 久久草成人影院| 国产男靠女视频免费网站| 成人性生交大片免费视频hd| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| 三级经典国产精品| 久久韩国三级中文字幕| 日本熟妇午夜| 男女之事视频高清在线观看| 搞女人的毛片| 超碰av人人做人人爽久久| 国产欧美日韩精品一区二区| 国产一区二区激情短视频| 国产成人a区在线观看| 久久精品国产亚洲av香蕉五月| 性欧美人与动物交配| 97碰自拍视频| 国产黄色视频一区二区在线观看 | 久久久久国产精品人妻aⅴ院| 亚洲最大成人av| 人妻丰满熟妇av一区二区三区| 国产私拍福利视频在线观看| 亚洲国产精品久久男人天堂| 国产一区二区在线观看日韩| 亚洲人成网站在线播| 97人妻精品一区二区三区麻豆| 日本撒尿小便嘘嘘汇集6| 国产高清视频在线播放一区| 国产成人aa在线观看| 精品一区二区三区人妻视频| 精品一区二区三区视频在线| 精品一区二区三区人妻视频| 国产淫片久久久久久久久| 99视频精品全部免费 在线| 欧美激情在线99| 国产伦在线观看视频一区| 欧美一区二区国产精品久久精品| a级一级毛片免费在线观看| 日韩精品有码人妻一区| 国产69精品久久久久777片| 国产女主播在线喷水免费视频网站 | 免费一级毛片在线播放高清视频| 精品乱码久久久久久99久播| 十八禁国产超污无遮挡网站| 久久久久久伊人网av| 91在线精品国自产拍蜜月| 亚洲av成人精品一区久久| 国产精品免费一区二区三区在线| 日本一本二区三区精品| 亚洲精品色激情综合| 国产 一区 欧美 日韩| 日日撸夜夜添| 淫妇啪啪啪对白视频| 搡女人真爽免费视频火全软件 | 精品人妻偷拍中文字幕| 久99久视频精品免费| 悠悠久久av| 欧美色视频一区免费| 不卡一级毛片| 成人性生交大片免费视频hd| 亚洲av中文av极速乱| 亚洲欧美中文字幕日韩二区| 久久久久久九九精品二区国产| 亚洲国产色片| 91久久精品电影网| 亚洲人与动物交配视频| 久久精品人妻少妇| 国产亚洲欧美98| 亚洲乱码一区二区免费版| 成年免费大片在线观看| 精华霜和精华液先用哪个| 色视频www国产| 午夜影院日韩av| 在线看三级毛片| 神马国产精品三级电影在线观看| aaaaa片日本免费| av中文乱码字幕在线| 夜夜夜夜夜久久久久| 人人妻,人人澡人人爽秒播| 99热6这里只有精品| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲网站| 午夜影院日韩av| 成人二区视频| 久久精品久久久久久噜噜老黄 | 男人狂女人下面高潮的视频| 色噜噜av男人的天堂激情| 国产91av在线免费观看| 日韩精品有码人妻一区| 中文在线观看免费www的网站| 精品一区二区免费观看| 精品一区二区三区av网在线观看| 女生性感内裤真人,穿戴方法视频| 看片在线看免费视频| 欧美性猛交黑人性爽| 97超碰精品成人国产| 露出奶头的视频| 91狼人影院| 日本爱情动作片www.在线观看 | 麻豆国产av国片精品| 91久久精品电影网| 国产视频内射| 色噜噜av男人的天堂激情| 亚洲美女搞黄在线观看 | 亚洲丝袜综合中文字幕| 欧美一区二区国产精品久久精品| 亚洲性久久影院| 超碰av人人做人人爽久久| 亚洲最大成人中文| 亚洲欧美日韩卡通动漫| 九九在线视频观看精品| 午夜精品在线福利| 国产高清视频在线播放一区| 久久99热6这里只有精品| 成人三级黄色视频| 亚洲不卡免费看| 国产蜜桃级精品一区二区三区| 直男gayav资源| av免费在线看不卡| 日韩大尺度精品在线看网址| 色哟哟·www| 在线播放无遮挡| 精品免费久久久久久久清纯| 一级黄色大片毛片| 在现免费观看毛片| 蜜臀久久99精品久久宅男| 午夜激情福利司机影院| 婷婷色综合大香蕉| 日本在线视频免费播放| 久久精品国产清高在天天线| 国产精品电影一区二区三区| 深爱激情五月婷婷| 特级一级黄色大片| 国模一区二区三区四区视频| 一进一出抽搐gif免费好疼| 一级黄片播放器| 91av网一区二区| 亚洲国产精品久久男人天堂| 国产欧美日韩精品亚洲av| 在线国产一区二区在线| 69人妻影院| 国产高清不卡午夜福利| 国产人妻一区二区三区在| av中文乱码字幕在线| 特大巨黑吊av在线直播| 亚洲一级一片aⅴ在线观看| 国产色婷婷99| 偷拍熟女少妇极品色| 99久久九九国产精品国产免费| 六月丁香七月| 身体一侧抽搐| 少妇熟女欧美另类| 亚洲欧美成人综合另类久久久 | 在线播放国产精品三级| 十八禁网站免费在线| 18+在线观看网站| 91在线精品国自产拍蜜月| 变态另类成人亚洲欧美熟女| 日韩亚洲欧美综合| 亚洲美女黄片视频| 又爽又黄a免费视频| 精品午夜福利在线看| 午夜免费男女啪啪视频观看 | 女的被弄到高潮叫床怎么办| 国产精品福利在线免费观看| 国内久久婷婷六月综合欲色啪| 黑人高潮一二区| 看黄色毛片网站| 日日干狠狠操夜夜爽| 九九在线视频观看精品| 少妇人妻精品综合一区二区 | 免费高清视频大片| 欧美高清成人免费视频www| 亚洲av熟女| 亚洲国产高清在线一区二区三| 最近2019中文字幕mv第一页| 午夜福利在线在线| 亚洲不卡免费看| 精品乱码久久久久久99久播| 97超碰精品成人国产| 日韩av在线大香蕉| 国产一级毛片七仙女欲春2| 99热6这里只有精品| 少妇丰满av| 日本一本二区三区精品| 深爱激情五月婷婷| 精品不卡国产一区二区三区| 国产真实乱freesex| 深夜a级毛片| 国产综合懂色| 黄色视频,在线免费观看| 国产精品国产高清国产av| 日韩高清综合在线| 此物有八面人人有两片| 国产单亲对白刺激| 久久婷婷人人爽人人干人人爱| 一级毛片aaaaaa免费看小| 亚洲国产欧美人成| 老熟妇乱子伦视频在线观看| 在线免费十八禁| 日韩人妻高清精品专区| 女的被弄到高潮叫床怎么办| 国产精品一区二区三区四区久久| 日韩在线高清观看一区二区三区| 午夜福利在线观看吧| 日韩一本色道免费dvd| 午夜福利成人在线免费观看| 国产成年人精品一区二区| 中国美白少妇内射xxxbb| 中国美女看黄片| 久久久久久九九精品二区国产| 亚洲欧美精品综合久久99| 国产精品人妻久久久影院| 日韩在线高清观看一区二区三区| 一区二区三区免费毛片| 神马国产精品三级电影在线观看| 国产av一区在线观看免费| 亚洲内射少妇av| 免费看光身美女| 国产成人影院久久av| av中文乱码字幕在线| 天堂影院成人在线观看| 亚洲性夜色夜夜综合| 人人妻人人澡人人爽人人夜夜 | 精品一区二区三区视频在线观看免费| 日本黄色片子视频| 免费不卡的大黄色大毛片视频在线观看 | 国产免费一级a男人的天堂| 欧美丝袜亚洲另类| 日日撸夜夜添| 久久久久免费精品人妻一区二区| 日韩高清综合在线| 欧美日本亚洲视频在线播放| 国产av一区在线观看免费| 狂野欧美激情性xxxx在线观看| 欧美最新免费一区二区三区| 韩国av在线不卡| 久久久久久伊人网av| 在线播放无遮挡| 一个人观看的视频www高清免费观看| 日韩欧美 国产精品| 欧美另类亚洲清纯唯美| 免费不卡的大黄色大毛片视频在线观看 | 人人妻人人澡欧美一区二区| 99久国产av精品| 免费不卡的大黄色大毛片视频在线观看 | 午夜a级毛片| 国产乱人偷精品视频| 日本精品一区二区三区蜜桃| 国产伦在线观看视频一区| 一级毛片我不卡| 免费观看在线日韩| 最近在线观看免费完整版| 欧美日韩在线观看h| 国产成人a区在线观看| 精品免费久久久久久久清纯| 12—13女人毛片做爰片一| 中出人妻视频一区二区| 床上黄色一级片| 国产精品国产高清国产av| 人妻少妇偷人精品九色| 男女做爰动态图高潮gif福利片| 日日摸夜夜添夜夜添小说| 69av精品久久久久久| 我的老师免费观看完整版| 中文字幕免费在线视频6| 国产国拍精品亚洲av在线观看| 国产免费男女视频| 国产片特级美女逼逼视频| 色噜噜av男人的天堂激情| 午夜视频国产福利| 亚洲高清免费不卡视频| 久久久午夜欧美精品| 老女人水多毛片| 超碰av人人做人人爽久久| 久久精品91蜜桃| 国产精品一二三区在线看| 国产精品,欧美在线| 免费一级毛片在线播放高清视频| 午夜福利在线观看吧| 午夜影院日韩av| 日韩大尺度精品在线看网址| 夜夜爽天天搞| 免费看a级黄色片| 久久久精品大字幕| 一个人免费在线观看电影| 日本黄大片高清| 成年av动漫网址| 亚洲国产欧美人成| 99视频精品全部免费 在线| 久久精品国产99精品国产亚洲性色| 国产精品人妻久久久影院| 国产精品一区二区三区四区免费观看 | .国产精品久久| 亚洲成a人片在线一区二区| 久久婷婷人人爽人人干人人爱| 免费av观看视频| 最近手机中文字幕大全| 最好的美女福利视频网| 又爽又黄a免费视频| 校园人妻丝袜中文字幕| 精品人妻偷拍中文字幕| 久久久久久久久久久丰满| 日韩欧美免费精品| 老熟妇仑乱视频hdxx| 又爽又黄a免费视频| 国产一区二区激情短视频| 亚洲最大成人中文| 日韩高清综合在线| 亚洲激情五月婷婷啪啪| 久久久国产成人精品二区| 国产极品精品免费视频能看的| 欧美性猛交黑人性爽| 日韩 亚洲 欧美在线| 黄片wwwwww| 又黄又爽又刺激的免费视频.| 亚洲精品粉嫩美女一区| 国产视频内射| 国内揄拍国产精品人妻在线| 一级毛片aaaaaa免费看小| 搡女人真爽免费视频火全软件 |