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

    Evaluation of Reanalysis Products with in situ GPS Sounding Observations in the Eastern Himalayas

    2014-12-08 07:33:50
    關(guān)鍵詞:件產(chǎn)品不合格品

    Department of Lower Atmosphere Observation Research (LAOR), Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029, China

    1 Introduction

    Located in the southern region of the Tibetan Plateau(TP), the Himalayans are the highest mountain ranges in the world, with an average height of 6 km and a length of 2400 km. Because these mountains border the TP and insulate the warm, moist air in the south from the cold,dry air in the north, they form a boundary of the climate systems of the TP and those of the adjacent Indian subcontinent (Ye and Gao, 1979; Boos and Kuang, 2010).The Yarlung Tsangpo River valley, which is located in the Eastern Himalayas, is a key water vapor channel for the TP (Gao et al., 1985; Yang et al., 1987; Tian et al., 2001).In summer, the southerlies from the Bay of Bengal bring large amounts of water vapor into the hinterland of TP through this channel, which profoundly influences the synoptic and climate systems of the TP as well as the downstream area (Simmonds et al., 1999; Zhou and Yu,2005).

    Due to the high altitude, severe weather, and tough environment, few in situ observations have been conducted in the TP, particularly in the high Himalayan region (Zhang et al., 2012). Therefore, reanalysis products, which include full space-time coverage, are widely used in atmospheric research of the TP. However, significant discrepancies exist among reanalysis products due to differences in physical structures, parameterization schemes,resolutions, and assimilation processes in these model systems.

    Many studies have evaluated the reliabilities of the reanalysis products in the TP, and most focus on surface variables such as surface temperature, precipitation, and surface heat fluxes (Zhao and Fu, 2006a, b; Ma et al.,2008, 2009; Mao et al., 2010; Wang and Zeng, 2012; Zhu et al., 2012). Such studies have demonstrated that reanalysis products can capture the climatology of most variables in the TP but are unable to reproduce their climatic variations. In addition, the observation data that have been assimilated into reanalysis systems are sparse and unevenly distributed, which contributes to the uncertainties in these products. Moreover, because the performances of reanalysis products in the TP are regionally dependent, no single reanalysis project is appropriate for all of the statistical quantities for surface variables (Wang and Zeng, 2012). Few studies have evaluated reanalysis data at high levels in the TP because of the difficulty in obtaining sounding data in this region. Bao and Zhang(2013) compared the datasets from Climate Forecast System Reanalysis (CFSR), National Center for Environmental Prediction and the Department of Energy Reanalysis 1 (NCEP-R1), Interim European Centre for Medium Range Weather Forecasts Reanalysis (ERAInterim), and ERA 40-year Reanalysis (ERA-40) with sounding observations from 11 stations in the TP from May to August 1998. Their results showed that although the mean structures among the products were consistent with the observation, large diurnal variations were present in their performances in different regions of the TP and at various levels. Thus, reanalysis products should be carefully evaluated for specific application prior to conducting meteorological research in the TP.

    In this study, we compared six widely used reanalysis products including NCEP Reanalysis 2 (NCEP-R2), NCEP Final Analysis (FNL), Japanese 25-year Reanalysis (JRA-25), ERA-Interim, Modern Era Retrospective Analysis for Research and Applications (MERRA), and CFSR with the independent sounding observations recorded in the Yarlung Tsangpo River valley in the Eastern Himalayas in June 2010. We evaluated air temperature, specific humidity, u-wind, and v-wind of these reanalysis datasets from 100 hPa to 650 hPa. This research provides a reference for studies in the Eastern Himalayas that consider reanalysis data.

    2 Data and method

    The Global Positioning System (GPS) sounding data were obtained from an observation campaign conducted in the Yarlung Tsangpo River valley (29.448°N, 94.691°E,2930 m) in the Eastern Himalayas. The Vaisala GPS radiosonde was operated twice daily, at 0000 UTC and 1200 UTC, from 1 June to 30 June 2010. The sounding data were collected with a time interval of 2 s; variables included height, barometric pressure, air temperature, relative humidity, wind speed, and wind direction. The observation data were linearly interpolated between 100 hPa and 700 hPa at intervals of 50 hPa.

    Six reanalysis products including CFSR (Saha et al.,2010), ERA-Interim (Simmons et al., 2006), FNL(NCEP/NWS/NOAA, 2000), JRA-25 (Onogi et al., 2007),MERRA (Rienecker et al., 2011), and NCEP-R2 (Kanamitsu et al., 2002) data were evaluated in this study. It should be noted that although FNL is not a reanalysis product, it was regarded as such in this study for convenience. The horizontal resolutions are 0.5° × 0.5° for CFSR, 0.75° × 0.75° for ERA-Interim, 1.0° × 1.0° for FNL,1.25° × 1.25° for JRA-25, 0.67° × 0.50° for MERRA, and 2.5° × 2.5° for NCEP-R2. At each pressure level, the gridded data of 0000 UTC and 1200 UTC were linearly interpolated to the observation site. The surface pressures of these analyses, which were between 652 hPa and 677 hPa, were lower than the 710 hPa observation value. In this study, variables from 100 hPa to 650 hPa were analyzed. All of the products had a vertical resolution of 50 hPa between 100 hPa and 300 hPa; CFSR, ERA-Interim,MERRA, and FNL had a resolution of 50 hPa between 300 hPa and 650 hPa; and JRA-25 and NCEP-R2 had a resolution of 100 hPa between 300 hPa and 600 hPa.

    Temperature (T), specific humidity (q), u-wind, and v-wind were evaluated by comparing the statistical quantities of the average, mean bias (MB), root-mean-square difference (RMSD), and correlation coefficient (R) of the reanalysis datasets with those of the observations.

    3 Results

    3.1 Temperature

    Figure 1 a shows that all of the products effectively captured the average profile of the temperature. Large errors were present at low levels below 500 hPa, and all products except for NCEP-R2 underestimated the temperature below 550 hPa. MERRA had the largest negative MB of –3.3°C at 650 hPa (Fig. 1b). NCEP-R2 had the largest positive MB of 0.7°C at 600 hPa and the largest negative MB above 400 hPa. Although JRA-25 had the smallest MB of –0.2°C at 600 hPa, this product had a large MB of 0.8°C at 400 hPa. Throughout the column,CFSR had the smallest MB, followed by ERA-Interim.

    Figure 1 Vertical profiles of the temperature statistical quantities between the six reanalysis datasets examined in this study and Global Positioning System (GPS) sounding in the Eastern Himalayas in June 2010 including (a) average, (b) mean bias (MB), (c) root-mean-square difference (RMSD), and (d) correlation coefficient (R). The reanalysis projects include Climate Forecast System Reanalysis (CFSR), Interim European Centre for Medium Range Weather Forecasts Reanalysis(ERA-Interim), Japanese 25-year Reanalysis (JRA-25), Modern Era Retrospective Analysis for Research and Applications (MERRA), National Center for Environmental Prediction and the Department of Energy Reanalysis 2 (NCEP-R2), and NCEP Final Analysis (FNL).

    The reanalysis products underestimated the temperature at lower and upper levels, which is consistent with previous studies of the surface temperature (Ma et al.,2008; Wang and Zeng, 2012); however they overestimated the temperature at middle levels, which differs from that reported by Bao and Zhang (2013). In previous studies, NCEP-R2 temperature showed the largest negative MB (Wang et al., 2012; Ma et al., 2008; Zhao and Fu,2006a; Bao and Zhang, 2013). However, our result shows that NCEP-R2 overestimated the temperature between 500 hPa and 600 hPa in the Eastern Himalayas in June 2010.

    The RMSD of temperature was small between 150 hPa and 500 hPa and large at low levels (Fig. 1c). MERRA had the largest RMSD at 650 hPa. NCEP-R2 had a larger RMSD at both lower and higher levels. Generally, the RMSDs of JRA-25 and ERA-Interim were the smallest throughout the column. The correlation coefficients for all products except CFSR and MERRA were larger than 0.8 below 250 hPa (Fig. 1d). The correlations increased with height and reached their maxima at 250 hPa and suddenly decreased to their minima at 200 hPa. This correlation coefficient profile differs from that reported by Bao and Zhang (2013), in which the temperature correlation coefficients between observation and CFSR and ERA-Interim increased with height and reached their maxima at 100 hPa.

    3.2 Humidity

    The mean profile of specific humidity was captured by the reanalysis products (Fig. 2a). The water vapor was overestimated at lower levels below 600 hPa and upper levels above 350 hPa and was underestimated at middle levels (Figs. 2a and 2b). At 650 hPa, CFSR had the largest MB at 0.67 g kg–1, and ERA-Interim had the smallest MB at 0.22 g kg–1. Negative MBs were larger at 500 hPa in NCEP-R2, JRA-25, and FNL; that of NCEP-R2 was the largest at –1.14 g kg–1. The negative MB of CFSR and ERA-Interim were the smallest. Above 300 hPa, the MBs were very small. The MB of reanalysis data for the Eastern Himalayas differs significantly from that reported by Bao and Zhang (2013), in which the RH was overestimated below 250 hPa and underestimated above 250 hPa.Throughout the column, ERA-Interim had the smallest MB for specific humidity, and CFSR has the second-smallest with the exception of a large negative MB at 450 hPa.

    The profiles of RMSDs were large and variable with height below 400 hPa and decreased with height above 400 hPa (Fig. 3c). At 650 hPa, CFSR has the largest RMSD, while FNL has the smallest. Throughout the column, CFSR and ERA-Interim had the smallest RMSD whereas NCEP-R2 had the largest; the RMSDs of JRA-25 and FNL were very similar.

    <1),且各件產(chǎn)品是否為不合格品相互獨(dú)立.

    The correlation coefficients of specific humidity decreased with height from 600 hPa. This result differs from the correlation of RH in other regions of the TP, where the correlation increases with height (Bao and Zhang, 2013).The correlation of specific humidity was smaller than that of the temperature, which indicates the reanalysis products were unable to capture the variation of the water vapor as well as the temperature. All analyses had the largest correlation coefficient at 600 hPa. NCEP-R2 had the weakest relationship with the observation at nearly all levels, even showing a negative value of –0.10 at 200 hPa. The correlation coefficients of other datasets were strong below 300 hPa; above this height, the correlation of ERA-Interim, MERRA, and JRA-25 were larger.Throughout the column, ERA-Interim showed the strongest correlation. JRA-25 also had a large correlation coefficient but it had no data at 650 hPa due to its coarse vertical resolution.

    Figure 2 Vertical profiles of the specific humidity between the six reanalysis datasets examined in this study (same as Fig. 1) and GPS sounding in the Eastern Himalayas in June 2010 including (a) average,(b) MB, (c) RMSD, and (d) R.

    Figure 3 Vertical profiles of the u-wind between the six reanalysis datasets examined in this study (same as Fig. 1) and GPS sounding in the Eastern Himalayas in June 2010 including (a) average, (b) MB, (c)RMSD, and (d) R.

    3.3 Wind

    The averaged observed u-wind in June 2010 was easterly below 550 hPa and shifted to westerly above (Fig.3a). The west wind increased with height and reached its maximum at 200 hPa, then decreased with height. ERAInterim, CFSR, and FNL reproduced both the low level easterly and the high level westerly; the u-wind of ERAInterim was closest to observation. All reanalysis products captured the variations of the u-wind above 600 hPa but showed large MBs near 550 hPa (Figs. 3a and 3b). All products overestimated the westerly between 600 and 500 hPa. NCEP-R2 had the largest positive MB of 7.6 m s–1at 600 hPa and the largest negative MB of –1.6 m s–1at 200 hPa. Throughout the column, the MB of ERA-Interim was the smallest, within 1.0 m s–1. For the u-wind RMSD,ERA-Interim and NCEP-R2 the largest and smallest, respectively (Fig. 3c). The RMSDs of CFSR and FNL were close and were small below 400 hPa and large above.

    The correlation coefficients of u-wind increased with height, showing small values of less than 0.5 below 500 hPa and a maximum of larger than 0.95 at 200 hPa. All reanalysis products except ERA-Interim and CFSR had negative correlation coefficients below 550 hPa; JRA-25 had the largest negative value at –0.27. The weak correlation below 500 hPa indicates that the reanalysis products were not able to capture the daily variation of the u-wind at lower levels in the Eastern Himalayas. Throughout the vertical column, the u-wind of ERA-Interim and CFSR showed the strongest correlation with the observation.

    The averaged observed v-wind increased with height from 650 hPa, reached maximum at 550 hPa, then decreased with height and shifted to northerly, showing a negative value at 200 hPa, and reached a maximum at 150 hPa (Fig. 4a). CFSR, ERA-Interim, FNL, and MERRA captured the southerly wind maximum at 550 hPa, because their vertical resolutions are higher than those of JRA-25 and NCEP-R2. All reanalysis products overestimated the v-wind below 550 hPa and underestimate it between 500 hPa and 200 hPa (Fig. 4b). MERRA had the largest positive MB of 5.2 m s–1at 600 hPa, and CFSR had the largest negative MB of –3.2 m s–1at 500 hPa.Throughout the column, ERA-Interim had the smallest bias.

    The v-wind RMSDs were between approximately 3.0 m s–1and 4.5 m s–1throughout the column (Fig. 4c). In general, ERA-Interim had the smallest RMSD for v-wind,and CFSR had the largest. The correlation coefficient of v-wind increased with height from 650 hPa, with small values of less than 0.5 appearing below 550 hPa (Fig. 4c).ERA-Interim, CFSR, and FNL showed negative correlation at 650 hPa. Throughout the column, ERA-Interim and JRA-25 had the closest relationship with the observation, whereas CFSR and FNL had the smallest correlation coefficients.

    4 Discussion and conclusions

    In this study, six reanalysis products were evaluated with independent sounding observations obtained in the Yarlung Tsangpo River valley in the Eastern Himalayas in June 2010. Our results show that the averaged vertical profiles of temperature, specific humidity, u-wind, and v-wind were captured by all reanalysis products; How-ever, large biases appeared at low levels, and each product showed different performances for different variables at different heights.

    Figure 4 Vertical profiles of the v-wind between the six reanalysis datasets examined in this study (same as Fig. 1) and GPS sounding in the Eastern Himalayas in June 2010 including (a) average, (b) MB, (c)RMSD, and (d) R.

    For the temperature, all products except NCEP-R2 had cold biases at lower and upper levels and warm biases at middle levels. NCEP-R2 showed an overestimation below 400 hPa and an underestimation above, which differs from that in other regions in the TP (Zhao and Fu, 2006a;Ma et al., 2008; Wang and Zeng, 2012; Bao and Zhang,2013). For specific humidity, all products had an overestimation at lower and upper levels and underestimation at middle levels, and their correlations with observations decreased with height from 600 hPa. These result differed from the RH in other region of the TP (Bao and Zhang,2013), in which RH was overestimated at the lower levels and underestimated at higher levels, and the correlations with observations increased with height. For u- and v-wind, all products had an overestimation at lower levels and underestimation at middle levels, and their correlations increased with height. In the study of Bao and Zhang (2013), CFSR and ERA-interim underestimated u-wind at all levels, whereas the MBs of v-wind were small.

    Throughout the column, ERA-Interim was closest to the observations for all the four variables in the Eastern Himalayas, whereas NCEP-R2 had the largest uncertainties. In addition, CFSR showed good performance in temperature and specific humidity. Therefore, although horizontal resolution has a relationship with the performances of reanalysis products for the Eastern Himalayas, it is not a decisive factor. Many factors such as model physical parameterizations and assimilation schemes can affect the accuracy of the reanalysis products. We found that in the Eastern Himalayas, although the resolutions of CFSR and MERRA were higher than those of other products, their surface pressures had the largest negative biases. This result may reflect the difficulties in topography processing in the Himalayas and may relate to their performances in this region.

    Therefore, no reanalysis product has been determined as superior to others for all variables at all heights.Therefore, users should select the appropriate datasets for their research purposes according to the performances of specific variables and related statistical quantities. Because large biases appeared at lower levels, cautions must be heeded when using these products as the initial and boundary fields for model simulations, particularly simulations of near-surface processes in the Himalayas.

    Acknowledgements.This study was supported by the Special Fund for Meteorological Research in the Public Interest(GYHY2012 06041), the National Natural Science Foundation of China (Grant No. 40905067), and the Ministry of Science and Technology of the People’s Republic of China (2009CB421403).The authors would like to acknowledge the data providers including NCEP for NCEP-R2, FNL, and CFSR data; European Centre for Medium-Range Weather Forecasts (ECMWF) for ERA-Interim data; the Japan Meteorological Agency (JMA) and the Central Research Institute of Electric Power Industry (CRIEPI) for JRA-25 data; and the Global Modeling and Assimilation Office (GMAO) for MERRA data.

    Bao, X., and F. Zhang, 2013: Evaluation of NCEP-CFSR, NCEPNCAR, ERA-Interim, and ERA-40 reanalysis datasets against independent sounding observations over the Tibetan Plateau, J.Climate, 26, 206–214.

    Boos, W. R., and Z. Kuang, 2010: Dominant control of the South Asian monsoon by orographic insulation versus plateau heating,Nature, 463, 218–222.

    Gao, D., H. Zou, and W. Wang, 1985: Influence of water vapor pass along the Yarlungzangbo River on precipitation, Mountain Res.(in Chinese), 3(4), 239–249.

    Kanamitsu, M., W. Ebisuzaki, J. Woollen, et al., 2002: NCEP-DOE AMIP-II Reanalysis (R-2), Bull. Amer. Meteor. Soc., 83, 1631–1643.

    Ma, L., T. Zhang, O. W. Frauenfeld, et al., 2009: Evaluation of precipitation from the ERA-40, NCEP-1, and NCEP-R2 Reanalyses and CMAP-1, CMAP-2, and GPCP-2 with ground-based measurements in China, J. Geophys. Res., 114, doi:10.1029/2008JD 011178.

    Ma, L., T. Zhang, Q. Li, et al., 2008: Evaluation of ERA-40,NCEP-1, and NCEP-R2 reanalysis air temperatures with groundbased measurements in China, J. Geophys. Res., 113, doi:10.1029/2007JD009549.

    Mao, J., X. Shi, L. Ma, et al., 2010: Assessment of reanalysis daily extreme temperatures with China’s homogenized historical dataset during 1979–2001 using probability density functions, J.Climate, 23, 6605–6623.

    NCEP/NWS/NOAA/U.S. Department of Commerce, 2000: NCEP FNL operational model global tropospheric analyses, continuing from July 1999, available at http://rda.ucar.edu/datasets/ds083.2.

    Onogi, K., and J. Tsutsui, H. Koide, et al., 2007: The JRA-25 Reanalysis, J. Meteor. Soc. Japan, 85, 369–432.

    Rienecker, M. R., M. J. Suarez, R. Gelaro, et al., 2011: MERRANASA’s Modern-ERA retrospective analysis for research and applications, J. Climate, 24, 3624–3648.

    Saha, S., S. Moorthi, H.-L. Pan, et al., 2010: The NCEP climate forecast system reanalysis, Bull. Amer. Meteor. Soc., 91, 1015–1057.

    Simmonds, I., D. Bi, and P. Hope, 1999: Atmospheric water vapor flux and its association with rainfall over China in summer, J.Climate, 12, 1353–1367.

    Simmons, A., S. Uppala, D. Dee, et al., 2006: ERA-Interim: NewECMWF reanalysis products from 1989 onwards, ECMWF News-lett., 110, 25–35.

    Tian, L., V. Masson-Delmotte, M. Stievenard, et al., 2001: Tibetan Plateau summer monsoon northward extent revealed by measurements of water stable isotopes, J. Geophys. Res.,106(D22), 28081–28088.

    Wang, A., and X. Zeng, 2012: Evaluation of multireanalysis products with in situ observations over the Tibetan Plateau, J.Geophys. Res., 117(D05102), doi:10.1029/2011JD016553

    Yang, Y., D. Gao, and B. Li, 1987: Primremary studying water vapor pass along the downstream of Yarlung Zangbo, Sci. China Ser.B-Chem. (in Chinese), 8, 893–902.

    Ye, D., and Y. Gao, 1979: Meteorology of the Qinghai-Xizang Plateau (in Chinese), CAS Qinghai-Tibet Plateau Compr. Sci.Invest. Team (Eds.), Science Press, Beijing, 279pp.

    Zhang, R., T. Koike, X. Xu, et al., 2012: A China-Japan cooperative JICA atmospheric observing network over the Tibetan Plateau(JICA/Tibet project): An overviews, J. Meteor. Soc. Japan, 90,1–16.

    Zhao, T., and C. Fu, 2006a: Preliminary comparison and analysis between ERA-40, NCEP-2 reanalysis and observations over China, Climate Environ. Res. (in Chinese), 11, 14–32.

    Zhao, T., and C. Fu, 2006b: Comparison of products from ERA-40,NCEP-2, and CRU station data for summer precipitation over China, Adv. Atmos. Sci., 23(4), 593–604.

    Zhou, T.-J., and R.-C. Yu, 2005: Atmospheric water vapor transport associated with typical anomalous summer rainfall patterns in China, J. Geophys. Res., 110, doi:10.1029/2004JD005413.

    Zhu, X., Y. Liu, and G. Wu, 2012: An assessment of summer sensible heat flux on the Tibetan Plateau from eight data sets, Sci.China Ser. D-Earth Sci., 55, 779–786.

    猜你喜歡
    件產(chǎn)品不合格品
    Dose-Dependent,Frequency-Dependent,and Cumulative Effects on Cardiomyocyte Injury and Autophagy of 2.856 GHz and 1.5 GHz Microwave in Wistar Rats
    不合格品管理中的誤區(qū)辨析和應(yīng)把握的幾個(gè)問題
    2020年1-11月軟件業(yè)務(wù)收入同比增12.5%
    Bone marrow mesenchymal stem cells induce M2 microglia polarization through PDGF-AA/MANF signaling
    軍工產(chǎn)品不合格品審理程序的分析與探究
    新目標(biāo)英語七年級(上)
    鋰電池企業(yè)不合格品管理
    科技視界(2014年22期)2014-09-02 00:01:00
    武器裝備研制不合格品管理方法研究
    av播播在线观看一区| 2022亚洲国产成人精品| 热re99久久精品国产66热6| 中文精品一卡2卡3卡4更新| 日韩伦理黄色片| 亚洲怡红院男人天堂| 日韩精品有码人妻一区| 日本-黄色视频高清免费观看| 一级毛片黄色毛片免费观看视频| 国产一级毛片在线| 国产精品伦人一区二区| 久久午夜福利片| xxx大片免费视频| 伊人亚洲综合成人网| 午夜免费观看性视频| tube8黄色片| 久久 成人 亚洲| av天堂久久9| 美女xxoo啪啪120秒动态图| 久久人人爽人人爽人人片va| 大码成人一级视频| 午夜日本视频在线| 国产欧美日韩精品一区二区| av在线老鸭窝| 国产精品国产三级国产专区5o| 少妇 在线观看| 老司机亚洲免费影院| 99九九线精品视频在线观看视频| 一区在线观看完整版| 日本与韩国留学比较| 只有这里有精品99| 国内揄拍国产精品人妻在线| 久久久久久久久大av| 一级a做视频免费观看| 天堂8中文在线网| 男人狂女人下面高潮的视频| 中文资源天堂在线| 亚洲精品日韩在线中文字幕| 美女xxoo啪啪120秒动态图| 麻豆乱淫一区二区| 人人妻人人爽人人添夜夜欢视频 | 啦啦啦视频在线资源免费观看| 成人18禁高潮啪啪吃奶动态图 | av在线app专区| 亚洲精品久久午夜乱码| 国产高清有码在线观看视频| 99久久精品国产国产毛片| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久久性| 久久久久久人妻| a级一级毛片免费在线观看| 人妻夜夜爽99麻豆av| 久久久久久伊人网av| 国产精品久久久久久av不卡| 久久久久久久国产电影| 九九在线视频观看精品| 这个男人来自地球电影免费观看 | 狂野欧美激情性bbbbbb| 一级毛片久久久久久久久女| av卡一久久| 国产精品久久久久成人av| 国产老妇伦熟女老妇高清| 高清毛片免费看| 97超碰精品成人国产| 成人免费观看视频高清| 黄片无遮挡物在线观看| 久久ye,这里只有精品| 一本久久精品| 香蕉精品网在线| 亚洲av成人精品一区久久| 熟女av电影| 久久久久久久久久成人| 亚洲丝袜综合中文字幕| 亚洲国产精品999| 亚洲成人av在线免费| 亚洲精品视频女| 国内揄拍国产精品人妻在线| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品一区二区性色av| 爱豆传媒免费全集在线观看| 国产精品一区www在线观看| 久久久久久久亚洲中文字幕| 91午夜精品亚洲一区二区三区| 午夜久久久在线观看| 男人爽女人下面视频在线观看| 中文字幕久久专区| 日日啪夜夜撸| 亚洲情色 制服丝袜| 2021少妇久久久久久久久久久| 一区二区三区乱码不卡18| 男人添女人高潮全过程视频| 如何舔出高潮| 天堂俺去俺来也www色官网| 中文字幕亚洲精品专区| 欧美精品高潮呻吟av久久| 高清视频免费观看一区二区| 亚洲欧美一区二区三区黑人 | 各种免费的搞黄视频| 春色校园在线视频观看| 欧美变态另类bdsm刘玥| 日日啪夜夜撸| 日日爽夜夜爽网站| 新久久久久国产一级毛片| 精品熟女少妇av免费看| 成人影院久久| 在线观看www视频免费| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 能在线免费看毛片的网站| 精品久久久噜噜| 欧美精品一区二区大全| 天美传媒精品一区二区| 国产日韩欧美亚洲二区| 蜜臀久久99精品久久宅男| 18禁在线播放成人免费| 女的被弄到高潮叫床怎么办| 久久人妻熟女aⅴ| 午夜久久久在线观看| 人妻系列 视频| 亚洲在久久综合| 国产欧美日韩综合在线一区二区 | 日日摸夜夜添夜夜添av毛片| av播播在线观看一区| 免费观看无遮挡的男女| 又黄又爽又刺激的免费视频.| 国产成人精品无人区| 一级a做视频免费观看| 人妻系列 视频| 黄色配什么色好看| 亚洲av男天堂| 久久99热这里只频精品6学生| 国产乱来视频区| 草草在线视频免费看| 菩萨蛮人人尽说江南好唐韦庄| 3wmmmm亚洲av在线观看| 日产精品乱码卡一卡2卡三| 菩萨蛮人人尽说江南好唐韦庄| 国产精品偷伦视频观看了| 亚洲欧美成人综合另类久久久| 精品少妇久久久久久888优播| 国产中年淑女户外野战色| 又黄又爽又刺激的免费视频.| 黄色日韩在线| 国内揄拍国产精品人妻在线| 一级二级三级毛片免费看| 简卡轻食公司| 国产日韩欧美在线精品| 国产一区二区三区综合在线观看 | 黄色毛片三级朝国网站 | 久久人妻熟女aⅴ| 久久久久精品久久久久真实原创| 久久午夜综合久久蜜桃| 2021少妇久久久久久久久久久| 韩国av在线不卡| 高清欧美精品videossex| 久久精品久久久久久久性| 精品一区二区三卡| 国产在线视频一区二区| 欧美最新免费一区二区三区| 男女啪啪激烈高潮av片| 日韩中字成人| 午夜免费鲁丝| 午夜激情福利司机影院| 成人综合一区亚洲| 成年美女黄网站色视频大全免费 | 国产白丝娇喘喷水9色精品| 欧美日韩视频精品一区| 亚洲高清免费不卡视频| 国产91av在线免费观看| 少妇精品久久久久久久| 国产精品久久久久久av不卡| 啦啦啦视频在线资源免费观看| 最近2019中文字幕mv第一页| 久久久久久久久久久丰满| 女性生殖器流出的白浆| 国产 一区精品| 亚洲av福利一区| 午夜视频国产福利| 国产亚洲一区二区精品| 午夜影院在线不卡| 欧美日韩一区二区视频在线观看视频在线| 久久99蜜桃精品久久| 欧美另类一区| 亚洲欧美中文字幕日韩二区| 免费少妇av软件| av天堂中文字幕网| 国产极品粉嫩免费观看在线 | 亚洲国产av新网站| 亚洲欧美日韩另类电影网站| 日韩,欧美,国产一区二区三区| 中文天堂在线官网| 亚洲美女视频黄频| 国产亚洲最大av| 国语对白做爰xxxⅹ性视频网站| 欧美xxxx性猛交bbbb| 极品教师在线视频| 纯流量卡能插随身wifi吗| 精品一品国产午夜福利视频| 国产欧美日韩精品一区二区| 欧美精品人与动牲交sv欧美| 久久久久久久久大av| 日日爽夜夜爽网站| 亚洲欧美日韩另类电影网站| 中文在线观看免费www的网站| 内射极品少妇av片p| 午夜日本视频在线| 欧美国产精品一级二级三级 | 色5月婷婷丁香| 亚洲国产精品一区二区三区在线| 熟妇人妻不卡中文字幕| 男人爽女人下面视频在线观看| 亚洲国产欧美日韩在线播放 | 日韩不卡一区二区三区视频在线| 自拍欧美九色日韩亚洲蝌蚪91 | 日日啪夜夜爽| 精品亚洲乱码少妇综合久久| 一个人看视频在线观看www免费| 欧美三级亚洲精品| 精品久久久久久电影网| 国产成人精品婷婷| 亚洲国产精品专区欧美| 精品国产乱码久久久久久小说| 国产高清国产精品国产三级| 国产精品久久久久久精品电影小说| 国产色婷婷99| 男女免费视频国产| 有码 亚洲区| 国产精品一区www在线观看| 久久久久国产网址| 我要看日韩黄色一级片| 男女无遮挡免费网站观看| 日韩精品有码人妻一区| 免费看av在线观看网站| 精品人妻熟女毛片av久久网站| av播播在线观看一区| 日韩强制内射视频| a 毛片基地| 久久狼人影院| 亚洲欧美一区二区三区黑人 | 欧美97在线视频| 国产成人一区二区在线| 国产成人精品无人区| 我的女老师完整版在线观看| 午夜老司机福利剧场| 免费黄频网站在线观看国产| 国产一区二区三区综合在线观看 | 免费大片18禁| 久久韩国三级中文字幕| 下体分泌物呈黄色| 纯流量卡能插随身wifi吗| 黄色配什么色好看| 九九久久精品国产亚洲av麻豆| 国产免费一区二区三区四区乱码| 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 欧美少妇被猛烈插入视频| 久久人人爽人人片av| 三级经典国产精品| 少妇裸体淫交视频免费看高清| 国产成人精品一,二区| 亚洲av欧美aⅴ国产| 一边亲一边摸免费视频| 99久久人妻综合| 免费人妻精品一区二区三区视频| 一区二区三区乱码不卡18| 精品亚洲成a人片在线观看| 国产精品一区二区性色av| 精品午夜福利在线看| 久久人人爽av亚洲精品天堂| 超碰97精品在线观看| 少妇的逼水好多| 亚洲国产精品一区二区三区在线| 久久国产乱子免费精品| 亚洲av福利一区| 老司机影院成人| av在线app专区| 日本免费在线观看一区| 国产成人freesex在线| 成年人免费黄色播放视频 | 精华霜和精华液先用哪个| 菩萨蛮人人尽说江南好唐韦庄| 日韩av免费高清视频| 欧美精品亚洲一区二区| 亚洲精品一区蜜桃| 在线精品无人区一区二区三| 22中文网久久字幕| 日韩av免费高清视频| 国产成人aa在线观看| 久久久久久久久久成人| 亚洲欧美成人综合另类久久久| 亚洲第一区二区三区不卡| 日韩伦理黄色片| 18禁在线无遮挡免费观看视频| 日日啪夜夜爽| 欧美xxⅹ黑人| 男女边吃奶边做爰视频| 欧美3d第一页| .国产精品久久| 日韩av免费高清视频| 麻豆乱淫一区二区| 久久精品国产亚洲av天美| 国产中年淑女户外野战色| 久久人人爽人人爽人人片va| 成人无遮挡网站| 精品亚洲成a人片在线观看| 天天操日日干夜夜撸| 激情五月婷婷亚洲| 久热久热在线精品观看| 99视频精品全部免费 在线| 男女边吃奶边做爰视频| 久久久久久久亚洲中文字幕| 亚洲精品日韩av片在线观看| 插阴视频在线观看视频| 在线观看国产h片| 久久精品夜色国产| 欧美高清成人免费视频www| 一区二区三区精品91| 午夜福利视频精品| 国产精品麻豆人妻色哟哟久久| 91午夜精品亚洲一区二区三区| 国精品久久久久久国模美| 观看免费一级毛片| 精品熟女少妇av免费看| 久久久国产一区二区| 午夜激情福利司机影院| 伊人久久精品亚洲午夜| 各种免费的搞黄视频| 五月天丁香电影| 国产在线男女| 亚洲av日韩在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 久久精品国产鲁丝片午夜精品| 成年人免费黄色播放视频 | 国国产精品蜜臀av免费| 99热6这里只有精品| 精品少妇黑人巨大在线播放| 少妇丰满av| 亚洲av免费高清在线观看| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久| 国产精品嫩草影院av在线观看| 亚洲精品一二三| 少妇人妻 视频| a级毛色黄片| 美女大奶头黄色视频| 亚洲av免费高清在线观看| 国产av码专区亚洲av| 美女内射精品一级片tv| 国产高清不卡午夜福利| 中文字幕精品免费在线观看视频 | 精品久久久噜噜| 深夜a级毛片| 亚洲综合精品二区| 男女无遮挡免费网站观看| 国内少妇人妻偷人精品xxx网站| 2018国产大陆天天弄谢| 国产精品一区www在线观看| 日韩精品免费视频一区二区三区 | 91午夜精品亚洲一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 久久韩国三级中文字幕| 国产一区亚洲一区在线观看| 天美传媒精品一区二区| 亚洲av在线观看美女高潮| 中文天堂在线官网| 一二三四中文在线观看免费高清| 免费黄频网站在线观看国产| 国产极品天堂在线| 一二三四中文在线观看免费高清| 亚洲成人一二三区av| 欧美xxⅹ黑人| 亚洲国产色片| 欧美精品一区二区大全| 狂野欧美白嫩少妇大欣赏| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 日产精品乱码卡一卡2卡三| 韩国高清视频一区二区三区| 国产精品偷伦视频观看了| 丰满饥渴人妻一区二区三| 亚洲一区二区三区欧美精品| 免费在线观看成人毛片| 91aial.com中文字幕在线观看| 麻豆精品久久久久久蜜桃| 日产精品乱码卡一卡2卡三| 久久人人爽人人片av| 国产极品天堂在线| 免费看日本二区| 国产精品无大码| 丝袜在线中文字幕| 日韩不卡一区二区三区视频在线| av在线app专区| 亚洲一级一片aⅴ在线观看| 丝袜喷水一区| 国产精品免费大片| 久久久欧美国产精品| av天堂久久9| 国产伦在线观看视频一区| 在线观看国产h片| 国产伦在线观看视频一区| 天美传媒精品一区二区| 日韩熟女老妇一区二区性免费视频| 99久久综合免费| 99精国产麻豆久久婷婷| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| 99九九在线精品视频 | 91久久精品国产一区二区三区| 久久青草综合色| 观看美女的网站| 日日啪夜夜爽| 国产精品一区二区性色av| 日韩欧美 国产精品| 日本av手机在线免费观看| 高清在线视频一区二区三区| 亚洲第一区二区三区不卡| 亚洲内射少妇av| 激情五月婷婷亚洲| 狂野欧美白嫩少妇大欣赏| 国产成人午夜福利电影在线观看| 26uuu在线亚洲综合色| 国产淫片久久久久久久久| 不卡视频在线观看欧美| 成人国产av品久久久| 丝袜脚勾引网站| 伊人久久精品亚洲午夜| 欧美少妇被猛烈插入视频| 丁香六月天网| 视频区图区小说| 久久精品久久久久久久性| 午夜福利网站1000一区二区三区| 能在线免费看毛片的网站| www.av在线官网国产| 激情五月婷婷亚洲| 大又大粗又爽又黄少妇毛片口| 亚洲精品日韩av片在线观看| 亚洲精品日本国产第一区| 成年人免费黄色播放视频 | 午夜激情久久久久久久| av女优亚洲男人天堂| 伦精品一区二区三区| 午夜日本视频在线| 在线观看免费高清a一片| 免费看av在线观看网站| 伊人亚洲综合成人网| 韩国av在线不卡| 成人黄色视频免费在线看| 精品99又大又爽又粗少妇毛片| 国产成人91sexporn| 精品久久久噜噜| 激情五月婷婷亚洲| 久久鲁丝午夜福利片| 18禁在线无遮挡免费观看视频| 亚洲第一区二区三区不卡| 精品人妻熟女av久视频| 国产 一区精品| 老熟女久久久| 国产亚洲欧美精品永久| 国产伦在线观看视频一区| 日日摸夜夜添夜夜爱| 熟女电影av网| 最近的中文字幕免费完整| 国产精品福利在线免费观看| 高清毛片免费看| 国产精品国产三级国产av玫瑰| 久久婷婷青草| 色视频在线一区二区三区| 亚洲成色77777| 一本一本综合久久| 免费人妻精品一区二区三区视频| 中文字幕制服av| av免费在线看不卡| 国产欧美日韩精品一区二区| 一级二级三级毛片免费看| 成人亚洲精品一区在线观看| 日韩亚洲欧美综合| 2022亚洲国产成人精品| 免费黄网站久久成人精品| 视频区图区小说| 在线亚洲精品国产二区图片欧美 | 欧美精品一区二区大全| 在线观看免费高清a一片| kizo精华| 十分钟在线观看高清视频www | 大又大粗又爽又黄少妇毛片口| 亚洲欧洲国产日韩| 国产淫片久久久久久久久| 国产午夜精品一二区理论片| 国产伦精品一区二区三区视频9| 免费看av在线观看网站| 国产永久视频网站| 亚洲欧美清纯卡通| 国产精品.久久久| 精品少妇久久久久久888优播| 97超碰精品成人国产| 一级毛片电影观看| 人妻 亚洲 视频| 高清av免费在线| 夜夜爽夜夜爽视频| 精品午夜福利在线看| 91精品一卡2卡3卡4卡| 在线观看免费高清a一片| 日本-黄色视频高清免费观看| 日韩 亚洲 欧美在线| 美女主播在线视频| 精品一区二区三卡| 亚洲国产精品一区三区| 成人亚洲欧美一区二区av| 黄色视频在线播放观看不卡| 99热国产这里只有精品6| 免费少妇av软件| 国产精品女同一区二区软件| 日日摸夜夜添夜夜添av毛片| 国产亚洲最大av| 国产伦精品一区二区三区四那| 爱豆传媒免费全集在线观看| 精品亚洲乱码少妇综合久久| 精华霜和精华液先用哪个| 久久精品国产亚洲网站| 最新的欧美精品一区二区| 一级黄片播放器| 汤姆久久久久久久影院中文字幕| 午夜激情福利司机影院| freevideosex欧美| 国产视频内射| 九九久久精品国产亚洲av麻豆| 在线观看免费视频网站a站| av国产精品久久久久影院| 日本黄大片高清| 超碰97精品在线观看| 亚洲一级一片aⅴ在线观看| a级片在线免费高清观看视频| 91精品伊人久久大香线蕉| 国产免费一级a男人的天堂| 人妻系列 视频| 色婷婷av一区二区三区视频| 亚洲精华国产精华液的使用体验| 精品国产国语对白av| 亚洲精品亚洲一区二区| 精品人妻一区二区三区麻豆| 我要看日韩黄色一级片| 国产av精品麻豆| 国产乱来视频区| 99久久中文字幕三级久久日本| 亚洲欧美成人精品一区二区| 亚洲婷婷狠狠爱综合网| 国产黄色免费在线视频| 97在线人人人人妻| 激情五月婷婷亚洲| 久久久久久久久久人人人人人人| 最近2019中文字幕mv第一页| 亚洲国产av新网站| 国产成人午夜福利电影在线观看| 人妻人人澡人人爽人人| 亚洲欧美精品专区久久| 日本91视频免费播放| 免费大片18禁| 高清黄色对白视频在线免费看 | 少妇人妻一区二区三区视频| 黄色欧美视频在线观看| 我要看日韩黄色一级片| 全区人妻精品视频| 超碰97精品在线观看| 精品午夜福利在线看| 国产一区二区三区综合在线观看 | 午夜福利影视在线免费观看| 亚洲欧美精品专区久久| 亚洲国产毛片av蜜桃av| 国产男人的电影天堂91| 人人妻人人看人人澡| 久久久精品94久久精品| 精品视频人人做人人爽| 亚洲内射少妇av| 亚洲丝袜综合中文字幕| 国产国拍精品亚洲av在线观看| 大又大粗又爽又黄少妇毛片口| 精品午夜福利在线看| tube8黄色片| 亚洲av福利一区| 欧美精品高潮呻吟av久久| 色吧在线观看| 精品人妻一区二区三区麻豆| 亚洲精品色激情综合| 一个人免费看片子| 国产欧美另类精品又又久久亚洲欧美| 涩涩av久久男人的天堂| 久久久久视频综合| 伊人久久国产一区二区| 啦啦啦在线观看免费高清www| av专区在线播放| 伦理电影免费视频| 伊人久久精品亚洲午夜| 九九爱精品视频在线观看| 免费播放大片免费观看视频在线观看| 成人美女网站在线观看视频| av免费观看日本| 中文字幕人妻熟人妻熟丝袜美| 久久久国产一区二区| 中国美白少妇内射xxxbb| 国产综合精华液| 人体艺术视频欧美日本| 亚洲欧美日韩另类电影网站| 亚洲av国产av综合av卡| 狠狠精品人妻久久久久久综合| 色婷婷av一区二区三区视频| av在线老鸭窝| 欧美国产精品一级二级三级 | 七月丁香在线播放| 赤兔流量卡办理| 最近中文字幕高清免费大全6| 人体艺术视频欧美日本| 成年人午夜在线观看视频| 亚洲美女黄色视频免费看| 国产中年淑女户外野战色| 人人妻人人澡人人看|