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

    An Observational and Modeling Study of Extratropical Transition of Hurricane Sandy in 2012

    2015-03-15 01:43:38FUDanLIPengyuanandFUGang
    Journal of Ocean University of China 2015年5期
    關(guān)鍵詞:氣象站普及氣象

    FU Dan, LI Pengyuan, and FU Gang

    Department of Meteorology, Ocean University of China, Qingdao 266100, P. R. China

    An Observational and Modeling Study of Extratropical Transition of Hurricane Sandy in 2012

    FU Dan, LI Pengyuan, and FU Gang*

    Department of Meteorology, Ocean University of China, Qingdao 266100, P. R. China

    ? Ocean University of China, Science Press and Springer-Verlag Berlin Heidelberg 2015

    Around 30 October 2012, Hurricane Sandy made landfall along the New Jersey shoreline after its completion of extratropical transition and transformation into an extratropical cyclone. The strong gale induced a catastrophic storm surge, and caused 72 death and damage of more than $50 billion. In this paper, the evolutionary process and spatial structure of the Hurricane Sandy during its extratropical transition were investigated by using Weather Research and Forecasting (WRF) version 3.3.1 modeling results and National Center for Environmental Prediction (NCEP) Coupled Forecast System model version 2 reanalysis datasets (CFSv2). It is found that during the upper-level trough interaction on 29 October, Sandy gradually fused with a pre-existing mid-latitude low-pressure system, and finished the re-intensification. WRF modeling results showed that the second peak occurred mainly due to the enhanced vertical motion, reduced vertical wind shear as well as the supplement of potential vorticity resulting from trough interaction over the southeast of Great Lakes. The cold continental air from the back of trough was encircled within the warm core system cyclonically, forming the characteristic of warm seclusion.

    Hurricane Sandy; extratropical transition; mid latitude trough; WRF modeling

    1 Introduction

    In meteorological community, Hurricane Sandy from 22 October to 31 October 2012 has widely been regarded to be a classic late-season tropical cyclone (TC) over the North Atlantic, as well as the deadliest and most destructive hurricane in the 2012 Atlantic hurricane season. During its movement crossing seven countries, at least 286 people were killed. On 29 October 2012, a devastating storm surge induced by the Hurricane Sandy attacked the United States stretching from the southern New Jersey to Rhode Island, causing 72 deaths and approximately $50 billion damage (Blake et al., 2013).

    Previous studies indicated that some of TCs might be transformed into extratropical cyclones at the end of their lives, usually occurring in mid-latitudes, where there is sufficient forcing from the upper-level troughs (Wang et al., 2012; DiMego and Bosart, 1982; Hart and Evans, 2001; Klein et al., 2000). This process in which TCs lose their tropical characteristics and become extratropical cyclones is usually defined as extratropical transition (ET) (Sekioka, 1970, 1972; Klein et al., 2000; Ritchie and Elsberry, 2001; Jones et al., 2003). During ET, the cyclone gradually loses its warm core structure and is distorted in shape, becoming less symmetric in circulation,convection, and humidity, which results in strong winds, extremely heavy rainfall as well as other potential secondary disasters like flooding and mudslide (Evans and Hart, 2008). Although in hurricane seasonextratropical transition of TC has been documented over the Northwest Pacific, Northeast Pacific, North Atlantic, and Southeast Indian Ocean, the systematic investigation of ET in later hurricane season is rare.

    Hurricane Sandy underwent ET process and became an extratropical cyclone around 21 UTC 29 October 2012 (Blake et al., 2013). During its ET process, Sandy produced torrential rainfalls that mainly fell into the south and west of Sandy’s moving track, which is consistent with previous studies (Bergeron, 1954; Palmén, 1958; Carr and Bosart, 1978; Foley and Hanstrum, 1994). The extreme precipitation was reported in Bellevue, Maryland with a peak amount of 325.88 mm, and an average of 127–178 mm of rainfall in the east of Maryland and Virginia, the south of Delaware and New Jersey. Due to the deep convection, even a snowfall occurred in an area stretching from North Carolina to the southwestern Pennsylvania, and the amount being over 90 cm in local accumulation. Additionally, the ET process of Sandy also was accompanied with a catastrophic storm surge along the northeastern United States coastline driven by the sharp pressure gradient and strong surface easterly flow.

    Several previous studies have been made to investigate the ET process of hurricane. Evans and Hart (2003) indicated that although cold core is the predominant factor tojudge the completion of ET for some cases, sometimes, the cold core is not so obvious even in the extratropical cyclone situation. Thorncroft and Jones (2000) found that a cold core develops at low levels while a weak warm core is always maintained above 600 hPa during the ET of Hurricane Felix in 1995. Browning et al. (1998) investigated the ET of Hurricane Lili in 1996, and illustrated how Lili developed into a system with warm core seclusion structure, which was similar to the final phase of extratropical frontal cyclone studied by Shapiro and Keyser (1990). Browning (2004) also indicated that the TCs with warm seclusion phase might easily produce damage near sea surface.

    For Hurricane Sandy, it reached a second peak during its ET process associated with warm seclusion stage, and the upper-level warm core was still maintained even when being a mature extratropical cyclone. The present paper aims to use the observational data as well as the simulation results from the Weather Research and Forecasting (WRF) Model to make a systematic study about the evolution process of Hurricane Sandy’s extratropical transition, including how Sandy loses its tropical structure and fuses with a pre-existing mid-latitude upper-level trough, as well as how Sandy is re-intensified during the warm seclusion phase. The paper is organized as follows. Section 2 displays an overview of Hurricane Sandy’s life cycle. Section 3 introduces the data and methodology. Section 4 observationally and numerically depicts the interaction between Sandy’s warm core and the midlatitude low-pressure system. Finally, Section 5 presents the conclusions.

    2 Overview of Hurricane Sandy

    Hurricane Sandy initially formed over the southern part of Kingston, Jamaica (13.3?N, 78.2?W) from an easterly wave which moved westward on 22 October, 2012. At that time, it was first identified as a tropical depression (TD) by National Hurricane Center (NHC) (Blake et al., 2013). After then, Sandy moved westward and developed into a tropical storm around 1800 UTC 22 October, being contributed by the slight vertical wind shear and warm SST. This system strengthened into a category-3 (Simpson) hurricane with about 954 hPa minimum central sea level pressure and 50 m s-1surface wind speed slightly before it went northward and made landfall in San Diego, Cuba around 0525 UTC 25 October. Meanwhile, a well- defined‘eye’ structure with deep eyewallcould be observed from the satellite image (Figs.1a and 1c). Sandy moved northwestward and was downgraded to a tropical storm soonafter crossing Cuba around 0000 UTC 27 October, due to the upper-level cold and dry air mass accompanied with strong southwesterly wind shear (Michael, 2012), with the convection becoming asymmetric (Fig.1a).

    Fig.1 (a), (b) Hurricane Sandy track from HURDAT2 overlaid by CIRA-RAMMB 4-km remapped enhanced infrared satellite imagery (shaded; unit: ℃) at 0625 UTC 25, 2345 UTC 25, 1825 UTC 28 and 0015 UTC 30 October respectively. The Sandy centers are marked by symbols of ‘L’,‘L’ for tropical depression, tropical storm, hurricane an d e xtratropical low pressure system according to HURDAT2. (c) Time series of HURDAT2 minimum central seapressure (SLP;solid line with dots; unit: hPa), maximum sustained wind speed (solid line with circles; unit:This figure isadapted from the Fig.1 of Galarneau et al. (2013).

    As Sandy turned northeastward, a well defined vortex structure was reconstructed and the convection center was mainly located on Sandy’s northwest flank by 0000 UTC 28 October (Fig.1b). During the day-time of October 29, Sandy approached a second pre-existing upper-level trough, and moved to northwestward towards the northeastern United States shoreline due to the midtropospheric blocking high pressure pattern. In the meantime, the diameter of whole Sandy system reached 1850 km, and the storm was re-intensified into a category-2 hurricane with 940 hPa minimum central sea level pressure and 45 m s-1wind speed at 1800 UTC 29 October (Fig.1c). The rapid intensification of wind speed and sea level pressure gradient drove a devastating storm surge that caused damage of tens of billion dollars. On 5 November 2012, NHC ranked the Hurricane Sandy as the second costliest US hurricane since 1900, just behind Hurricane Katrina in 2005.

    3 Data and Modeling Results

    3.1 Observational and Reanalysis Data

    In this paper, we employ the following data to study the ET process of Hurricane Sandy systematically.

    1) Hurricane data 2 (HURDAT2) offered by NHC, which provides the hurricane’s best track position, minimum sea level pressure, and maximum sustained wind speed. This data is available at 0000, 0600, 1200, 1800 UTC every day. Available from: http://www.nhc.noaa.gov/ data/ hurdat/hurdat2-1851-2013-052714.txt.

    2) 4 km×4 km resolution Remapped Color Enhanced Infrared Imagery from the Cooperative Institute for Research in the Atmosphere (CIRA, Colorado State University)– Regional and Mesoscale Meteorology Branch (RAMMB). Available from: http://rammb.cira.colostate.edu/products/ tc_realtime/.

    3) Multiplatform Satellite Surface Wind Analysis product, which combines the information from different data sources including Advanced Microwave Sounding Unit (AMSU) winds, Cloud-drift/IR/WV winds, IR-proxy winds and Scatterometer winds. Available from: http://rammb.cira. colostate.edu/products/tc_realtime/.

    4) Sounding data of station 72501 downloaded from University of Wyoming, which is available 2 times per day. Available from: http://weather.uwyo.edu/upperair/sounding. html.

    5) Climate Forecast System version 2 (CFSv2) data from NCEP (National Center for Environmental Prediction), the spatial resolution is 0.5?×0.5? and 37 layers in the vertical direction (1000, 975, 950, 925, 900, 875, 850, 825, 800, 775, 750, 700, 650, 600, 550, 500, 450, 400, 350, 300, 250, 225, 200, 175, 150, 125, 100, 70, 50, 30, 20, 10, 7, 5, 3, 2 , 1 hPa), the time interval is 6 hours. Available from: http://rda.ucar.edu/datasets/ds094.0/.

    6) Real-time, global, sea surface temperature (RTG_ SST_HR) data from NCEP, the spatial resolution is 0.083?× 0.083?, and it is a daily mean product. Available from: ftp://polar.ncep.noaa.gov/pub/history/sst.

    3.2 Modeling Results

    The numerical simulation of ET process of Hurricane Sandy was performed by WRFv3.3.1 model. Two domains were used in the simulation with an outer domain of 21 km×21 km horizontal resolution which covers the whole east of United States and west of Atlantic basin, and an inner nest domain of 7 km×7 km (Fig.2a). The total integration times for outer domain and inner domain were 60 hours (initialized at 0000 UTC 28 October 2012) and 48 hours (initialized at 1200 UTC 28 October 2012), respectively. Both two domains had 40 layers vertically, up to 50 hPa. For the configuration, Yonsei University Planenary Boundary Layer (PBL) scheme (YSU; Noh et al., 2003; Hong et al., 2006), Noah-MP Land surface physics (Chen and Dudhia, 2001), Dudhia (1989) short wave radiation scheme, Rapid Radiative Transfer Model (RRTM) long wave radiation scheme (Mlawer et al., 1997), Betts-Miller-Janjic cumulus parameterization (Betts, 1986; Betts and Miller, 1986; Janjic, 1994) and WDM 6-class microphysics scheme (Lin et al., 1983) were used in the simulation (Table 1).

    The simulation results of inner domain during the ET period of Hurricane Sandy are shown in Fig.2. Generally, the track simulation is relative accurate, with less than 30 km deviation compared with the position from HURDAT2 data. Basically, the nested simulation of WRFv3.3.1 is able to capture the interaction between Sandy and second upper-level trough as well as Sandy’s northwestward steering toward New Jersey shoreline. However, the intensity simulation is not as good as track simulation. This is perhaps due to the impotency of initial condition, neither simulated wind speed nor central sea level pressure was robust enough at first 3-h running. After a period of adjustment, the WRF model successfully captured Sandy’s second peak and rapid dissipation both in pressure and wind speed field. But the simulated wind speed at 10-m height was roughly 2–7 m s-1smaller than the observations (Fig.2b). The comparison of the 10-m wind speed between the satellite data and WRF simulation at 1800 UTC 29 October is shown in Figs.2c and 2d. WRF model could reconstruct the Sandy’s circulation as well as the‘eye’ structure during the its secondary peak, but the area with more than 50 knot wind speed seemed to be too large both at the south and the north of Sandy’s center, while the wind speed at the south of 35?N was too weak.

    As the initial condition was merely derived from CFSv2 reanalysis data, which may have some biases with observation, the simulated intensity is less accurate compared with the track simulation. Anyhow, the WRFv3.3.1 model is able to reproduce the whole ET process of Sandy successfully. It is hoped in future that data assimilation and ensemble simulation will improve the simulation accuracy.

    Fig.2 (a) Comparison of Hurricane Sandy track between WRF simulation (red solid line with dot) and HURDAT2 (black solid line with dot). The nested domain used for the simulation is also shown (red dash box). (b) Time series of HURDAT2 minimum central sea level pressure (SLP; black solid line with dots; unit: hPa), intensity (Vmax; black solid line with circles; unit: m s-1) as well as WRF simulated minimum central sea level pressure (red line; unit: hPa) and 10 m wind speed (green line; unit: m s-1) of Hurricane Sandy. (c) Multiplatform Satellite Surface Wind Analysis at 1800 UTC 29 October (unit: knot). (d) WRF simulated 10 m wind speed (unit: knot) and barb at 1800 UTC 29 October.

    Table 1 Model parameterization for Hurricane Sandy (2012) simulation

    4 Analyses of Sandy’s ET Process

    4.1 Interaction with Mid-Latitude Trough

    When Sandy moved poleward to mid-latitude, the environmental conditions changed a lot, characterized by the enhancement of Coriolis force, decreasing of sea surface temperature as well as increased atmospheric baroclinicity. A pre-existing upper-level trough deepened due to the building of mid-tropospheric trough over the central Atlantic basin and located over the central America from 28 to 30 October. Additionally, the mid-latitudetrough dove southeastward from central America and formed a negative tilt late on 28 October. Under this condition, northwestward moving Sandy gradually interacted with this upper-tropospheric trough. As Sandy was located downstream, the upper-level trough not only provided great divergence but also decreased the vertical wind shear, especially at the west flank of Hurricane Sandy (Figs.3a and 4a). Increased upper-level divergence in conjunction with decreased vertical wind shear further enhanced vertical motion near the Sandy’s center, making Sandy more convective. At this time, the Sandy’s central air mass mean temperature was distinctly higher than the surrounding.

    Around 0600 UTC 29 October, Sandy moved closer to the upper-level trough. At the 500 hPa level, Sandy’s 5600 gpm isobar had been fused with mid-latitude trough, and deep convection near Sandy’s warm core center and northwest flank wasmore concentrated (Figs.3b and 4b). Late on 29 October, a new but smaller low-pressure system was drawn out from negative tilt trough at 200 hPa and approached Sandy (Fig.3c). Simultaneously, the cold continental air from the back of mid-latitude tough intruded southeastward and entered Sandy’s warm core cyclonically prior to its landfall along New Jersey shoreline (Figs.3c and 3d). At mid-level, Sandy’s isobars had fused with mid-latitude low-pressure system, forming a new and big low-pressure center and even being deepened.

    Using sea surface and satellite observations, we may review this period more clearly. At 1500 UTC 29 October, an expanded low-pressure system was observed at surface, with a stationary front and an occlusion located on the western and northern sides of Sandy, respectively (Fig.5a). A dense overcast near cyclone center could still be found from the satellite image at this time, but disappeared in the next 6 hours (Fig.5b). By 2100 UTC 29 October, the latter cold front caught up with the former warm front largely, and occlusion was even wrapped into the cyclone center. Moreover, Sandy’s warm core was lifted by the latter cold air and further secluded, which made Sandy acquire the characteristics of a warm seclusion type cyclone (Galarneau et al., 2013). As Sandy no longer maintained its deep convection and organized frontal structure, it could be confirmed that Sandy had become an extratropical cyclone after finishing its ET process by the end of 29 October.

    Twelve hours later, the upward vertical motion near the Sandy’s center became very weak, and Sandy had totally fused with pre-existing upper-level trough both at 200 hPa and 500 hPa levels. In addition, the cold air mass at 850 hPa was embedded into warm core structure totally, but the 1000–500 hPa thickness field formed a ‘tongue-like’structure (Figs.3d and 4d) near the low-pressure center. This ‘tongue-like’ geopotential thickness field suggested the incomplete vanishing of warm core, as the 500–1000 hPa thickness value is used to define ‘bulk’ air mass mean temperature.

    and 1000–500 hPathick ness (blue contour, un it: gpm) at (a) 1800 UTC 28, (b) 0600 UTC 29, (c) 0000 UTC 30, and (d) 1200 UTC 30 October2012. The cyclone centers are indicated with symbols ‘’ for hurricane or ‘L’ for extratropicallow pressure system.

    Fig.4 500 hPa geopotential height (black contour, unit: gpm), 850 hPa potential temperature (shaded, lower than 291 K, unit: K), 200 hPa minus 850 hPa wind shear (arrows, unit: m s-1) and 600–400 hPa mean ascent (red contour, beginning from -0.5 Pa s-1, unit: Pa s-1) at (a) 1800 UTC 28, (b) 0600 UTC 29, (c) 0000 UTC 30, and (d) 1200 UTC 30 October 2012.The cyclone centers are indicated with symbols ‘’ for hurricane or ‘L’ for extratropical low pressure system.

    Fig.5 (a) Surface observations at 1500 UTC 29 October and GOES-E visible satellite image at 1445 UTC 29 October overlapped by fronts and isobars (4 hPa interval) analysis that provided by NHC at 1500 UTC 29 October. (b) Surface observations at 2100 UTC 29 October and GOES-E visible satellite image at 2045 UTC 29 October overlapped by fronts and isobars (4 hPa interval) analysis that provided by NHC at 2100 UTC 29 October (from Fig.22 and Fig.23 of Blake et al., 2013).

    4.2 Intensification During Sandy’s ET Process

    From both observational data and reanalysis data we can find that Hurricane Sandy experienced the ET process from late 28 October to late 29 October, and it became an extratropical cyclone prior to its landfall along the New Jersey shoreline around 2330 UTC October 29. During that period, Sandy reintensified from a weak category-1 hurricane with a wind speed around 35 m s-1to a category-2 hurricane with a wind speed around 45 m s-1near surface. Moreover, the minimum central sea level pressure even dropped down to 940 hPa at 1800 UTC 29 October simultaneously. This sharply increased wind speed accompanied by a rapidly decreased sea level pressure caused a devastating storm surge–the averaged water level rise at New York City came up to 3.96 m above normal tide level, and an easterly fetch that affected the entire eastern coast of United States (Blake et al., 2013).In order to investigate this extraordinary reintensification deeply, we will use more observational data as well as the WRF simulation results to analyze this process.

    We now analyse the vertical sections which respectively crossing the Sandy center at four different times (1800 UTC 28, 0900 UTC 29, 2100 UTC 29 and 1200 UTC 30 October) to illustrate the evolutionary processes of Sandy’s secondary peak from the perspectives of geopotential height, potential temperature, three dimensional wind speed, as well as potential vorticity. Each vertical section is designed in 20? longitude by 20? latitude and center crossed Sandy’s inner core (Fig.6).

    Fig.6 Sandy track simulated by WRF from 1200 UTC 28 to 1200 UTC 30 October 2012, and the orientations of cross sections. Each section is 20? longitude by 20? latitude, centering with Sandy inner core. Number 1–4 indicate the time at 1800 UTC 28, 0900 UTC 29, 2100 UTC 29 and 1200 UTC 31 October 2012, respectively.

    TC is a rapidly rotating vortex system characterized by a low-pressure center, surrounded by great pressure gradient confined to the inner core and decreasing with height (Frank, 1977). By 1800 UTC 28 October, the deep negative geopotential height anomaly could reach to 300 hPa, while the 0 isoline had a great northwestward tilting over 400 hPa, resulting from the approaching of upstream upper-level trough. The whole system was really deep, with obvious warm core inside (Figs.7a, 7d). From both geopotential height anomaly and wind speed fields, we could find a well-defined eyewall structure, while deep convection was only limited in a narrow band located in the west and north of Sandy (Figs.7b, 7e). Additionally, a upper-level jet stream stayed in the south of Sandy near 200 hPa height, and the horizontal wind speed in the south was apparently greater than that of the north of Sandy center. The high PV value was located near the kernel structure with high PV gradient confined to the inner core (Figs.7c, 7f), just like a mature TC potential vorticity distribution (Kofron et al., 2010). Hurricane Sandy had a typical mature TC characteristic at this time.

    Early on 29 October, Sandy made its forth turning, and advanced to the northeastern coastline of the United States. At 0900 UTC 29 October, due to the interaction with upper-level trough, the negative geopotential height anomaly was not as robust as 15 hours ago, and the system got a prominent trend to westward tilt above 600 hPa (Fig.8a). The upward motion in the western side of Sandy was strengthened, which was maybe caused by the increased upper-level divergence (Figs.3b, 8b). Despite the upward motion, the wind speed also had a significant increase on both of the northern and southern sides from bottom to over 400 hPa (Figs.8b, 8e). The weakened vertical wind shear further contributed to Sandy’s intensification (Fig.4b). What is more, besides the inner core high PV axis, a upper-level high PV was center located in west of low-level PV center (Fig.8c), which would strength the vortex through a series of adjustments (Hoskins et al., 1985). After that, Sandy moved across the main axis of Gulf Stream near 1200 UTC 29 October (Fig.9a), where there is a great enhancement in SST. Higher surface temperature is theoretically favorable to supply more energy to TC (Evans, 1993; Schade, 2000). The sounding profile observed at the station 72501 (40.86?N, 72.86?W, marked in Fig.2a) at 1200 UTC 29 October is shown in Fig.9b. As the dew-point deficit was almost zero under 500 hPa, it suggested that environmental condition was really humid due to the high SST. Additionally, the hodograph and wind barbs suggested that vertical wind shear was also really small both in magnitude and direction. The increased upper-level divergence and reduced vertical wind shear together with the higher SST caused Sandy’s reintensification.

    Sandy reached its second peak at 1800 UTC 29 Octo-ber. By 2100 UTC 29 October, Sandy still maintained its strong westward tilting structure, but cold air mass had occupied the south and west of Sandy, and invaded the warm core from upper-level to bottom in the east of inner center, i.e., cold continental air was encircled inside Sandy cyclonically (Figs.10a, 10d). Although both vertical and horizontal windswere reinforced remarkably, Sandy began to lose its symmetric wind field gradually (Figs.10b, 10e). At the same time, potential vorticity displayed an interesting vertical distribution: although the major high PV axis slightly tilted to the west from top to bottom, another 500 hPa high PV center traveled ahead of the low level center (Figs.10c, 10f). This unusual PV configuration indicated the maturity of Sandy’s secondary peak as well as the upcoming dissipation (Hoskins et al., 1985).

    After Sandy was reintensified, it moved over much colder water area, in which SST dropped from 25℃ to 17℃. Sandy made its landfall to the United States around 2330 UTC 29 October. Surface friction and cold air mass in conjunction with lowerSST weakened this system and hastened its loss of tropical structure. By 1200 UTC 30 October, Sandy had been completely integrated intopreexisting mid-latitude low-pressure system. The inner core of Sandy no longermaintained its structure as tightly as previous, though the central pressure was still lower than the surrounding and pressure gradient weakened significantly (Figs.11a, 11d). Due to the invasion of cold conti-nental air mass, the warm core structure almost vanished below 700 hPa but was still activated above 500 hPa. The wind field weakened and became asymmetric (Figs.11b, 11e). The potential vorticity near the inner core also reduced dramatically and became less concentrated (Fig.11c), which suggests that Sandy became an extratropical cyclone and dissipated much.

    Fig.7 Vertical cross section analyses at 1800 UTC 28 October 2012. (a)–(c) W1–E1 and (d)–(f) N1–S1 orientations are shown in Fig.6. (a), (d): geopotential height anomaly (contours, zero line is marked by blue color; unit: gpm) and potential temperature anomaly (shaded; unit: K), each mean field is determined from the longitude/latitude that shown above; (b),(e): horizontal windspeed (contours, red contours mark the wind speed greater than and upward mo-tion (shaded; unit: (c), (f): potential vorticity (unit: PVU). The cyclone center is denoted by the symbol ‘’.

    Fig.8 As in Fig.7, but at 0900 UTC 29 October 2012. (a)–(c) W2–E2 and (d)–(f) S2–N2 orientations are shown in Fig.6.The cyclone center is denoted by the symbol ‘’.

    Fig.9 (a) Threedays (28, 29, 30 October) mean sea surface temperature (SST, unit: ℃). The cyclone centers are marked bythe symbol of ‘’ and ‘’. (b) Sounding data and hodograph of the station 72501 (Fig.2a) at 1200 UTC 29 October 2012. Solid line indicates air temperature (unit: ℃), and dashed line indicates dew-point temperature (unit: ℃).

    5 Conclusions

    In this paper, we made a systematic investigation regarding the ET process of Hurricane Sandy observationally and numerically. WRF3.3.1 model successfully reproduces the ET of Sandy that began on 29 October 2012 when Sandy interacted with a pre-existing upper-level negative tilt trough. The negative tilt trough enhanced upper-level divergence and decreased vertical wind shear significantly, which further contributed to deep convection. In addition to higher SST near Gulf Stream, Hurricane Sandy began to reintensified results from the cyclonic penetration of cold continental air that from back of upper-level trough and wrapped into the Sandy’s warm core. This strong continental cold air mass further caught up the warm front which located at the east of Sandy’s center, making Sandy get the characteristic of warm seclusion cyclone. This warm seclusion structure during ET process also provided a positive environmental condition for Sandy’s reintensification. Sandy reached its second peak around 1800 UTC 29 October with 940 hPa central sea level pressure and 45 m s-1wind speed. Late on 29 October, the cold air mass invaded the warm core from upper-level to low-level in the east of Sandy. Surface fric-tion, lower SST and involvement of cold continental air mass accelerated Sandy’s dissipation. Early on 30 October, both the decreased inner core pressure and PV gradient associated with asymmetric wind field suggested Sandy’s loss of tropical structure. At last, Sandy completely fused with a small low-pressure system from midlatitude negative tilt trough and became an extratropical cyclone.

    Fig.10 As in Fig.7, but at 2100 UTC 29 October 2012. (a)–(c) W3–E3 and (d)–(f) S3–N3 orientationsare shown in Fig.6.The cyclone center is denoted by the symbol ’.

    Fig.11 As in Fig.7, but at 1200 UTC 30 October 2012. (a)–(c) W4–E4 and (d)–(f) S4–N4 orientations are shown in Fig.6. The cyclone center is denoted by the symbol ‘L’.

    Acknowledgements

    References

    Bergeron, T., 1954. Review of modern meteorology – 12. The problem of tropical hurricanes. Quarterly Journal of the Royal Meteorological Society, 80: 131-164.

    Betts, A. K., 1986. A new convective adjustment scheme. Part I: Observational and theoretical basis. Quarterly Journal of the Royal Meteorological Society, 112: 677-691.

    Betts, A. K., and Miller, M. J., 1986. A new convective adjustment scheme. Part II: Single column tests using GATE wave, BOMEX, and arctic air-mass data sets. Quarterly Journal of the Royal Meteorological Society, 112: 693-709.

    Blake, E. S., Kimberlain, T. B., Berg, R. J., Cangialosi, J. P., and Beven II, J. L., 2013. Tropical Cyclone Report Hurricane Sandy (AL182012) 22–29 October 2012. National Hurricane Center.

    Browning, K. A., 2004. The sting at the end of the tail: Damaging winds associated with extratropical cyclones. Quarterly Journal of the Royal Meteorological Society, 130: 375-399.

    Browning, K. A., Vaughan, G., and Panagi, P., 1998. Analysis of an ex-tropical cyclone after reintensifying as a warm-core extratropical cyclone. Quarterly Journal of the Royal Meteorological Society, 124: 2329-2356.

    Carr, F. H., and Bosart, L. F., 1978. A diagnostic evaluation of rainfall predictability for tropical storm Agnes, June 1972. Monthly Weather Review, 106: 363-374.

    科技的發(fā)展使得現(xiàn)階段我國(guó)農(nóng)業(yè)生產(chǎn)得到不斷的提高,但是農(nóng)業(yè)由于受到氣象因素的影響較大,因此良好的氣象監(jiān)測(cè)是實(shí)現(xiàn)農(nóng)業(yè)進(jìn)一步發(fā)展的基礎(chǔ)。我國(guó)現(xiàn)階段自動(dòng)氣象站的使用已經(jīng)得到普及,基于自動(dòng)氣象站的各項(xiàng)技術(shù)也日趨成熟,得到大范圍的使用。

    Chen, F., and Dudhia, J., 2001. Coupling an advanced landsurface/hydrology model with the Penn State/NCAR MM5 modeling system. Part I: Model description and implementation. Monthly Weather Review, 129: 569-585.

    DiMego, G. J., and Bosart, L. F., 1982. The transformation of Tropical Storm Agnes into an extratropical cyclone. Part I: The observed fields and vertical motion computations. Monthly Weather Review, 110: 385-411.

    Dudhia, J., 1989. Numerical study of convection observed during the winter monsoon experiment using a mesoscale twodimensional model. Journal of the Atmospheric Sciences, 46: 3077-3107.

    Evans, C., and Hart, R. E., 2008. Analysis of the wind field evolution associated with the extratropical transition of Bonnie (1998). Monthly Weather Review, 136: 2047-2065.

    Evans, J. L., 1993. Sensitivity of tropical cyclones to sea surface temperature. Journal of Climate, 6: 1133-1140.

    Evans, J. L., and Hart, R. E., 2003. Objective indicators of the life cycle evolution of extratropical transition for Atlantic tropical cyclones. Monthly Weather Review, 131: 909-925.

    Foley, G. R., and Hanstrum, B. N., 1994. The capture of tropical cyclones by cold fronts off the west coast of Australia. Weather and Forecasting, 9: 577-592.

    Frank, W. M., 1977. The structure and energetics of the tropical cyclone I. Storm structure. Monthly Weather Review, 105: 1119-1135.

    Galarneau, T. J., Davis, C. A., and Shapiro, M. A., 2013. Intensification of Hurricane Sandy (2012) through extratropical warm core seclusion. Monthly Weather Review, 141: 4296-4321.

    Hart, R. E., and Evans, J. L., 2001. A climatology of the extratropical transition of Atlantic tropical cyclones. Journal of Climate, 14: 546-564.

    Hong, S.-Y., Noh, Y., and Dudhia, J., 2006. A new vertical diffusion package with an explicit treatment of entrainment processes. Monthly Weather Review, 134: 2318-2341.

    Hoskins, B. J., McIntyre, M. E., and Robertson, A. W., 1985. On the use and significance of isentropic potential vorticity maps. Quarterly Journal of the Royal Meteorological Society, 111: 877-946.

    Janjic, Z. I., 1994. The step-mountain eta coordinate model: Further developments of the convection, viscous sublayer, and turbulence closure schemes. Monthly Weather Review, 122: 927-945.

    Jones, S. C., Harr, P. A., Abraham, J., Bosart, L. F., Bowyer, P. J., Evans, J. L., Hanley, D. E., Hanstrum, B. N., Hart, R. E., Lalaurette, F., Sinclair, M. R., Smith, R. K., and Thorncroft, C., 2003. The extratropical transition of tropical cyclones: Forecast challenges, current understanding, and future directions. Weather and Forecasting, 18: 1052-1092.

    Klein, P. M., Harr, P. A., and Elsberry, R. L., 2000. Extratropical transition of western North Pacific tropical cyclones: An overview and conceptual model of the transformation stage. Weather and Forecasting, 15: 373-396.

    Kofron, D. E., Ritchie, E. A., and Tyo, J. S., 2010. Determination of a consistent time for the extratropical transition of tropical cyclones. Part II: Potential vorticity metrics. Monthly Weather Review, 138: 4344-4361.

    Lin, Y.-L., Farley, R. D., and Orville, H. D., 1983. Bulk parameterization of the snow field in a cloud model. Journal of Climate and Applied Meteorology, 22: 1065-1092.

    Michael, B., 2012. Hurricane Sandy Discussion NuhPaer 14 (Report). National Hurricane Center.

    Mlawer, E. J., Taubman, S. J., Brown, P. D., Iacono, M. J., and Clough, S. A., 1997. Radiative transfer for inhomogeneous atmosphere: RRTM, a validated correlated-k model for the longwave. Journal of Geophysical Research, 102 (D14): 16663-16682.

    Noh, Y., Cheon, W. G., Hong, S. Y., and Raasch, S., 2003. Improvement of the K-profile model for the planetary boundary layer based on large eddy simulation data. Boundary-Layer Meteorology, 107: 401-427.

    Palmén, E., 1958. Vertical circulation and release of kinetic energy during the development of hurricane Hazel into an extratropical storm. Tellus, 10: 1-23.

    Ritchie, E. A., and Elsberry, R. L., 2001. Simulations of the transformation stage of the extratropical transition of tropical cyclones. Monthly Weather Review, 129: 1462-1480.

    Schade, L. R., 2000. Tropical cyclone intensity and sea surface temperature. Journal of the Atmospheric Sciences, 57: 3122-3130.

    Sekioka, M., 1970. On the behaviour of cloud pattern as seen on satellite photographs in the transformation of a typhoon into an extratropical cyclone. Journal of the Meteorological Society of Japan, 48: 224-233.

    Sekioka, M., 1972. Note on the extratropical transformation of a typhoon in relation with old outbreaks. Archiv für Meteorologie, Geophysik und Bioklimatologie, A21: 413-418.

    Shapiro, M. A., and Keyser, D., 1990. Fronts, jet streams, and the tropopause. Extratropical Cyclones: The Erik Palmén Memorial Volume. Newton, C. W., and Holopaninen, E. O., eds., American Meteorological Society, 167-191.

    Thorncroft, C. D., and Jones, S. C., 2000. The extratropical transition of Hurricanes Felix and Iris in 1995. Monthly Weather Review, 128: 947-972.

    Wang, Q., Q. Li, and G. Fu, 2012: Determining the Extratropical Transition Onset and Completion Times of Typhoons Mindulle (2004) and Yagi (2006) Using Four Methods. Wea. Forecasting, 27, 1394–1412.

    (Edited by Xie Jun)

    (Received September 30, 2014; revised November 28, 2014; accepted April 13, 2015)

    J. Ocean Univ. China (Oceanic and Coastal Sea Research)

    DOI 10.1007/s11802-015-2770-2

    ISSN 1672-5182, 2015 14 (5): 783-794

    http://www.ouc.edu.cn/xbywb/

    E-mail:xbywb@ouc.edu.cn

    * Corresponding author. Tel: 0086-532-66781920

    E-mail: fugang@ouc.edu.cn

    猜你喜歡
    氣象站普及氣象
    氣象
    氣象樹(shù)
    珠峰上架起世界最高氣象站
    《內(nèi)蒙古氣象》征稿簡(jiǎn)則
    5G手機(jī)打響“普及戰(zhàn)”
    海峽姐妹(2020年7期)2020-08-13 07:49:28
    心靈氣象站
    大國(guó)氣象
    自動(dòng)氣象站應(yīng)該注意的一些防雷問(wèn)題
    自動(dòng)氣象站常見(jiàn)故障判斷與維護(hù)
    河南科技(2014年12期)2014-02-27 14:10:40
    天文知識(shí)普及
    視野(2012年2期)2012-07-26 02:50:20
    三上悠亚av全集在线观看| 天天躁日日躁夜夜躁夜夜| 精品一区在线观看国产| 美女扒开内裤让男人捅视频| 精品一区在线观看国产| 久久综合国产亚洲精品| 免费在线观看日本一区| 国产片内射在线| 夫妻午夜视频| 久久久久久久久免费视频了| 国产男女内射视频| 777米奇影视久久| 国产精品偷伦视频观看了| 亚洲精品久久久久久婷婷小说| 1024视频免费在线观看| 亚洲av日韩在线播放| 亚洲性夜色夜夜综合| 五月开心婷婷网| 亚洲精品粉嫩美女一区| 久久 成人 亚洲| 久久这里只有精品19| 久久久久久久久久久久大奶| 欧美在线一区亚洲| 麻豆av在线久日| 国产男女超爽视频在线观看| 19禁男女啪啪无遮挡网站| 久久久欧美国产精品| 老汉色av国产亚洲站长工具| 一二三四社区在线视频社区8| 深夜精品福利| 欧美激情极品国产一区二区三区| 99国产精品免费福利视频| 99精国产麻豆久久婷婷| 人人妻人人澡人人爽人人夜夜| 欧美激情 高清一区二区三区| 啦啦啦 在线观看视频| 香蕉丝袜av| 国产福利在线免费观看视频| 亚洲美女黄色视频免费看| 两人在一起打扑克的视频| 伦理电影免费视频| 国产亚洲午夜精品一区二区久久| 亚洲免费av在线视频| 日本av免费视频播放| 国产成人免费观看mmmm| 18禁裸乳无遮挡动漫免费视频| 亚洲国产中文字幕在线视频| 亚洲国产精品一区三区| 国产成人精品久久二区二区免费| 午夜福利视频精品| 国产在视频线精品| 久久青草综合色| 免费在线观看视频国产中文字幕亚洲 | 国产精品久久久久成人av| 国产亚洲av高清不卡| 另类精品久久| 亚洲av成人不卡在线观看播放网 | 国产一区二区三区av在线| 桃花免费在线播放| 伦理电影免费视频| 99国产精品99久久久久| 色婷婷久久久亚洲欧美| 自线自在国产av| 欧美老熟妇乱子伦牲交| 国产欧美亚洲国产| 日韩,欧美,国产一区二区三区| 久久精品国产亚洲av高清一级| av福利片在线| 亚洲第一欧美日韩一区二区三区 | a在线观看视频网站| 亚洲精品乱久久久久久| 精品第一国产精品| 国产免费一区二区三区四区乱码| 久久国产亚洲av麻豆专区| 久久久久国产精品人妻一区二区| 久久人人97超碰香蕉20202| 国产精品麻豆人妻色哟哟久久| 午夜福利乱码中文字幕| 成年人免费黄色播放视频| 国产精品 国内视频| 一进一出抽搐动态| 99热国产这里只有精品6| 久久香蕉激情| 国产精品亚洲av一区麻豆| 亚洲国产av新网站| 日韩制服骚丝袜av| 日本猛色少妇xxxxx猛交久久| 亚洲情色 制服丝袜| 国产精品久久久久成人av| 狠狠精品人妻久久久久久综合| 男女午夜视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 精品国产一区二区三区四区第35| 五月开心婷婷网| 在线观看免费日韩欧美大片| 亚洲美女黄色视频免费看| 久久久久视频综合| 99久久99久久久精品蜜桃| 少妇粗大呻吟视频| 欧美黑人精品巨大| 国产不卡av网站在线观看| 亚洲av美国av| 99国产精品一区二区蜜桃av | 亚洲熟女毛片儿| 乱人伦中国视频| 国产精品二区激情视频| a级片在线免费高清观看视频| 一本一本久久a久久精品综合妖精| 亚洲av电影在线进入| 国产色视频综合| 婷婷丁香在线五月| 欧美大码av| 日日摸夜夜添夜夜添小说| 精品少妇一区二区三区视频日本电影| av在线老鸭窝| 90打野战视频偷拍视频| 男女下面插进去视频免费观看| 人人妻人人澡人人看| 99热网站在线观看| 狠狠婷婷综合久久久久久88av| 一个人免费在线观看的高清视频 | 成人国语在线视频| 一级毛片精品| 午夜两性在线视频| 亚洲欧美色中文字幕在线| 国产亚洲精品第一综合不卡| 欧美黑人欧美精品刺激| 欧美xxⅹ黑人| 国产成人av教育| 精品国产乱子伦一区二区三区 | 性色av乱码一区二区三区2| 国产黄频视频在线观看| 国产成人免费观看mmmm| www.熟女人妻精品国产| 欧美日韩av久久| 亚洲欧美色中文字幕在线| 女性生殖器流出的白浆| 在线观看免费日韩欧美大片| 菩萨蛮人人尽说江南好唐韦庄| 淫妇啪啪啪对白视频 | 精品视频人人做人人爽| 日韩熟女老妇一区二区性免费视频| 麻豆国产av国片精品| www.精华液| 久久精品国产综合久久久| 欧美日韩黄片免| 午夜福利一区二区在线看| 高清黄色对白视频在线免费看| 精品久久久久久电影网| 窝窝影院91人妻| 国产1区2区3区精品| 天天躁日日躁夜夜躁夜夜| 99久久人妻综合| 色94色欧美一区二区| 日韩有码中文字幕| 国产亚洲午夜精品一区二区久久| 精品国产乱子伦一区二区三区 | 亚洲人成77777在线视频| 在线观看免费午夜福利视频| 丰满人妻熟妇乱又伦精品不卡| 国产免费一区二区三区四区乱码| 黄网站色视频无遮挡免费观看| 久久精品国产a三级三级三级| 国产av一区二区精品久久| 最近最新中文字幕大全免费视频| 久久久精品国产亚洲av高清涩受| 法律面前人人平等表现在哪些方面 | 丝袜喷水一区| 香蕉丝袜av| 可以免费在线观看a视频的电影网站| 亚洲成人手机| 动漫黄色视频在线观看| 亚洲成av片中文字幕在线观看| 日韩大码丰满熟妇| 亚洲成国产人片在线观看| 国产精品久久久久成人av| 国产欧美日韩一区二区三区在线| 国产三级黄色录像| 99热国产这里只有精品6| 丁香六月欧美| 秋霞在线观看毛片| 日韩,欧美,国产一区二区三区| av有码第一页| e午夜精品久久久久久久| 精品国产国语对白av| 国产精品影院久久| 国产一区二区三区综合在线观看| 亚洲五月婷婷丁香| 国产精品一区二区免费欧美 | 制服诱惑二区| 脱女人内裤的视频| 久久人人爽人人片av| 新久久久久国产一级毛片| 国产三级黄色录像| 热99re8久久精品国产| 黑人操中国人逼视频| 亚洲精品中文字幕在线视频| 亚洲av美国av| 亚洲精品自拍成人| 2018国产大陆天天弄谢| 99热全是精品| 91麻豆精品激情在线观看国产 | 一级片免费观看大全| √禁漫天堂资源中文www| 中文字幕av电影在线播放| 亚洲午夜精品一区,二区,三区| 老司机影院成人| 精品亚洲成国产av| 午夜福利一区二区在线看| 欧美激情高清一区二区三区| 在线观看免费午夜福利视频| 精品第一国产精品| av免费在线观看网站| 自线自在国产av| av有码第一页| 国产熟女午夜一区二区三区| 丰满迷人的少妇在线观看| 欧美性长视频在线观看| 久久精品国产a三级三级三级| 亚洲国产日韩一区二区| 国产av又大| 免费av中文字幕在线| 精品福利观看| 午夜福利乱码中文字幕| 亚洲欧美精品自产自拍| 一级毛片电影观看| 亚洲精品一卡2卡三卡4卡5卡 | 一本色道久久久久久精品综合| 精品国产超薄肉色丝袜足j| 少妇裸体淫交视频免费看高清 | avwww免费| 精品熟女少妇八av免费久了| 水蜜桃什么品种好| 国产精品自产拍在线观看55亚洲 | 99热网站在线观看| 婷婷成人精品国产| 最新的欧美精品一区二区| 国产在视频线精品| 久久久国产精品麻豆| 久久久精品94久久精品| 免费黄频网站在线观看国产| 精品国产一区二区久久| 韩国高清视频一区二区三区| 午夜福利在线观看吧| 欧美乱码精品一区二区三区| 日韩三级视频一区二区三区| 久久人人爽av亚洲精品天堂| 久久精品aⅴ一区二区三区四区| 18禁黄网站禁片午夜丰满| 黄色视频在线播放观看不卡| 精品国产一区二区三区久久久樱花| 91老司机精品| 国产av一区二区精品久久| e午夜精品久久久久久久| 国产无遮挡羞羞视频在线观看| 国产亚洲精品久久久久5区| 亚洲 欧美一区二区三区| h视频一区二区三区| 亚洲欧美成人综合另类久久久| 狂野欧美激情性bbbbbb| 纵有疾风起免费观看全集完整版| 一本综合久久免费| videosex国产| 免费一级毛片在线播放高清视频 | 日本av免费视频播放| 9热在线视频观看99| 大陆偷拍与自拍| 青青草视频在线视频观看| 天天影视国产精品| av片东京热男人的天堂| 亚洲专区中文字幕在线| 成人黄色视频免费在线看| 亚洲国产精品成人久久小说| 少妇裸体淫交视频免费看高清 | 国产欧美日韩精品亚洲av| 熟女少妇亚洲综合色aaa.| 极品少妇高潮喷水抽搐| bbb黄色大片| 天堂俺去俺来也www色官网| cao死你这个sao货| 岛国在线观看网站| 亚洲av日韩在线播放| 99久久99久久久精品蜜桃| 精品少妇黑人巨大在线播放| 91老司机精品| 一本一本久久a久久精品综合妖精| 欧美日韩亚洲国产一区二区在线观看 | 日韩大片免费观看网站| 欧美+亚洲+日韩+国产| 香蕉国产在线看| 美女国产高潮福利片在线看| 在线看a的网站| 热re99久久国产66热| 十八禁网站网址无遮挡| 欧美久久黑人一区二区| 在线观看免费视频网站a站| 成人国语在线视频| 日本91视频免费播放| 日本一区二区免费在线视频| 久久久精品区二区三区| 亚洲激情五月婷婷啪啪| 欧美精品啪啪一区二区三区 | 啦啦啦中文免费视频观看日本| kizo精华| www.熟女人妻精品国产| 国产精品久久久人人做人人爽| 各种免费的搞黄视频| 18禁观看日本| 美女主播在线视频| 欧美亚洲日本最大视频资源| 国产黄频视频在线观看| 久久国产精品影院| 久久久久国内视频| 日韩一区二区三区影片| 这个男人来自地球电影免费观看| 黑人巨大精品欧美一区二区mp4| 日本a在线网址| 国产成人精品久久二区二区免费| 人人妻人人澡人人爽人人夜夜| 久久99热这里只频精品6学生| 午夜福利乱码中文字幕| 免费av中文字幕在线| 日韩视频在线欧美| 大码成人一级视频| 国产精品久久久av美女十八| 亚洲av男天堂| 一个人免费看片子| 老司机午夜福利在线观看视频 | 久久天堂一区二区三区四区| 亚洲欧洲日产国产| 国产主播在线观看一区二区| 多毛熟女@视频| 男女无遮挡免费网站观看| 少妇人妻久久综合中文| 99国产极品粉嫩在线观看| 久久香蕉激情| 久久人人97超碰香蕉20202| 日本av免费视频播放| 建设人人有责人人尽责人人享有的| 久久天躁狠狠躁夜夜2o2o| 手机成人av网站| 国产伦人伦偷精品视频| 亚洲专区字幕在线| 精品久久久久久久毛片微露脸 | 一本一本久久a久久精品综合妖精| 国产高清videossex| 精品人妻在线不人妻| 亚洲国产中文字幕在线视频| 久热这里只有精品99| 国产精品免费大片| 久久精品aⅴ一区二区三区四区| 午夜激情久久久久久久| 欧美黑人精品巨大| 亚洲精品国产色婷婷电影| 国产av一区二区精品久久| 午夜久久久在线观看| av欧美777| 亚洲avbb在线观看| 在线永久观看黄色视频| 久久ye,这里只有精品| 在线永久观看黄色视频| 两个人看的免费小视频| 一本一本久久a久久精品综合妖精| 91大片在线观看| 宅男免费午夜| 久久久久久久久免费视频了| 欧美av亚洲av综合av国产av| 久久免费观看电影| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 亚洲av片天天在线观看| 国产精品二区激情视频| 婷婷丁香在线五月| 亚洲国产欧美日韩在线播放| 中亚洲国语对白在线视频| 久久 成人 亚洲| 日韩欧美一区二区三区在线观看 | 精品国产乱码久久久久久小说| 一级a爱视频在线免费观看| 久久精品人人爽人人爽视色| 18禁裸乳无遮挡动漫免费视频| 一二三四在线观看免费中文在| 啦啦啦 在线观看视频| 亚洲精品乱久久久久久| 亚洲欧美日韩另类电影网站| 99国产精品免费福利视频| 日韩大片免费观看网站| 亚洲自偷自拍图片 自拍| 少妇精品久久久久久久| 亚洲国产欧美一区二区综合| 真人做人爱边吃奶动态| 国产精品香港三级国产av潘金莲| 国产成人精品久久二区二区免费| 十八禁高潮呻吟视频| 麻豆av在线久日| 国产精品99久久99久久久不卡| 免费观看人在逋| 国产成人免费无遮挡视频| 男人添女人高潮全过程视频| 亚洲美女黄色视频免费看| 美国免费a级毛片| 成在线人永久免费视频| 91九色精品人成在线观看| 黄色视频,在线免费观看| 9191精品国产免费久久| 性色av乱码一区二区三区2| 男人爽女人下面视频在线观看| cao死你这个sao货| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 亚洲精品国产av蜜桃| 精品久久久久久电影网| 国产精品麻豆人妻色哟哟久久| 一级片'在线观看视频| 国产av一区二区精品久久| av天堂久久9| 亚洲精品美女久久久久99蜜臀| 精品国内亚洲2022精品成人 | 国产高清视频在线播放一区 | 美女中出高潮动态图| 一级毛片电影观看| 国产在线视频一区二区| 丝袜美足系列| 天堂中文最新版在线下载| 亚洲专区字幕在线| 亚洲精华国产精华精| 精品少妇一区二区三区视频日本电影| 久9热在线精品视频| 三上悠亚av全集在线观看| 免费不卡黄色视频| 美女午夜性视频免费| 高清av免费在线| 亚洲精品一二三| 50天的宝宝边吃奶边哭怎么回事| 国产精品.久久久| 9191精品国产免费久久| 国产精品久久久久久人妻精品电影 | 国产亚洲av高清不卡| 波多野结衣一区麻豆| 老熟女久久久| 欧美一级毛片孕妇| 久久毛片免费看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 日韩视频一区二区在线观看| 男人添女人高潮全过程视频| 黑人猛操日本美女一级片| 亚洲男人天堂网一区| 国产免费视频播放在线视频| 精品国产乱码久久久久久男人| 曰老女人黄片| 欧美日韩视频精品一区| videos熟女内射| 午夜久久久在线观看| 日本精品一区二区三区蜜桃| 在线观看舔阴道视频| 黄频高清免费视频| 久久国产精品影院| 日本一区二区免费在线视频| 国产又爽黄色视频| 两个人看的免费小视频| 老司机影院毛片| 美女国产高潮福利片在线看| 亚洲美女黄色视频免费看| 久久精品国产a三级三级三级| 91麻豆av在线| 欧美老熟妇乱子伦牲交| 午夜福利视频精品| cao死你这个sao货| 老汉色∧v一级毛片| 色婷婷av一区二区三区视频| 午夜福利乱码中文字幕| 日日夜夜操网爽| 国产亚洲精品久久久久5区| 狂野欧美激情性bbbbbb| 欧美日韩福利视频一区二区| 国产日韩欧美在线精品| 中国国产av一级| 性高湖久久久久久久久免费观看| www日本在线高清视频| 亚洲一区二区三区欧美精品| 久久人人97超碰香蕉20202| 久久久久国产一级毛片高清牌| kizo精华| 精品国产国语对白av| av网站免费在线观看视频| 亚洲精品国产区一区二| 99热网站在线观看| 国产又色又爽无遮挡免| 99久久综合免费| 在线看a的网站| 热re99久久精品国产66热6| 人人妻人人澡人人看| 亚洲国产看品久久| av国产精品久久久久影院| 狠狠婷婷综合久久久久久88av| 黑人巨大精品欧美一区二区mp4| 免费在线观看黄色视频的| tube8黄色片| 亚洲国产精品一区三区| 午夜福利在线免费观看网站| 午夜日韩欧美国产| 国产精品.久久久| av天堂久久9| 女性被躁到高潮视频| 成年美女黄网站色视频大全免费| 亚洲自偷自拍图片 自拍| 香蕉丝袜av| 欧美激情久久久久久爽电影 | 亚洲专区中文字幕在线| 久久久精品国产亚洲av高清涩受| 国产亚洲av片在线观看秒播厂| 人妻久久中文字幕网| 亚洲精品一二三| 一区二区三区四区激情视频| 日韩视频一区二区在线观看| 免费在线观看视频国产中文字幕亚洲 | 中文精品一卡2卡3卡4更新| 国产淫语在线视频| 成人18禁高潮啪啪吃奶动态图| 日韩,欧美,国产一区二区三区| 国产免费一区二区三区四区乱码| 日本a在线网址| 亚洲国产精品成人久久小说| 亚洲精品日韩在线中文字幕| 侵犯人妻中文字幕一二三四区| 国产黄频视频在线观看| 午夜精品国产一区二区电影| 久久这里只有精品19| 欧美97在线视频| √禁漫天堂资源中文www| 一区福利在线观看| 国产人伦9x9x在线观看| 免费在线观看完整版高清| 亚洲 欧美一区二区三区| 性高湖久久久久久久久免费观看| 黄片播放在线免费| 青春草视频在线免费观看| 法律面前人人平等表现在哪些方面 | 精品人妻一区二区三区麻豆| 99久久精品国产亚洲精品| 亚洲伊人色综图| 精品一区在线观看国产| 大陆偷拍与自拍| 一区二区三区激情视频| 少妇的丰满在线观看| 国产伦理片在线播放av一区| 国产精品久久久久成人av| 后天国语完整版免费观看| 欧美97在线视频| 国产成人免费无遮挡视频| 一级片免费观看大全| av线在线观看网站| 精品一区在线观看国产| 大陆偷拍与自拍| 韩国高清视频一区二区三区| 精品久久蜜臀av无| 一二三四社区在线视频社区8| 色婷婷久久久亚洲欧美| 国产欧美日韩综合在线一区二区| av欧美777| 亚洲成人手机| 亚洲五月婷婷丁香| 国产又爽黄色视频| 国产1区2区3区精品| 亚洲精品自拍成人| 国产亚洲av片在线观看秒播厂| 黄色a级毛片大全视频| 天堂中文最新版在线下载| 欧美 日韩 精品 国产| 法律面前人人平等表现在哪些方面 | 中国国产av一级| 黄频高清免费视频| 亚洲国产成人一精品久久久| svipshipincom国产片| 99国产综合亚洲精品| 蜜桃国产av成人99| www.精华液| 亚洲性夜色夜夜综合| 超碰97精品在线观看| 亚洲avbb在线观看| 国产精品自产拍在线观看55亚洲 | 97精品久久久久久久久久精品| 两人在一起打扑克的视频| 99re6热这里在线精品视频| 免费日韩欧美在线观看| 国产极品粉嫩免费观看在线| 欧美黑人精品巨大| 精品久久久精品久久久| 免费观看a级毛片全部| 一区二区日韩欧美中文字幕| 女人爽到高潮嗷嗷叫在线视频| 国产熟女午夜一区二区三区| 午夜福利乱码中文字幕| 91九色精品人成在线观看| 成人黄色视频免费在线看| 欧美黑人精品巨大| 午夜福利影视在线免费观看| 最黄视频免费看| 蜜桃在线观看..| 色婷婷久久久亚洲欧美| 久久这里只有精品19| 男女床上黄色一级片免费看| 97人妻天天添夜夜摸| e午夜精品久久久久久久| 日韩,欧美,国产一区二区三区| 一级a爱视频在线免费观看| 国产精品亚洲av一区麻豆| 大香蕉久久网| 亚洲三区欧美一区| 熟女少妇亚洲综合色aaa.| 99精国产麻豆久久婷婷| 欧美 日韩 精品 国产| 搡老乐熟女国产| 51午夜福利影视在线观看| 国产深夜福利视频在线观看|