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

    Validation of the Modified Becker's Split-Window Approach for Retrieving Land Surface Temperature from AVHRR

    2015-11-21 11:23:31QUANWeijun權(quán)維俊CHENHongbin陳洪濱HANXiuzhen韓秀珍andMAZhiqiang馬志強
    Journal of Meteorological Research 2015年5期
    關(guān)鍵詞:志強

    QUAN Weijun(權(quán)維?。?,CHEN Hongbin(陳洪濱),HAN Xiuzhen(韓秀珍),and MA Zhiqiang(馬志強)

    1 Key Laboratory for Middle Atmosphere and Global Environment Observation,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    2 Environmental Meteorology Forecast Center of Beijing-Tianjin-Hebei,Beijing Meteorological Bureau,Beijing 100089

    3 National Satellite Meteorological Center,China Meteorological Administration,Beijing 100081

    Validation of the Modified Becker's Split-Window Approach for Retrieving Land Surface Temperature from AVHRR

    QUAN Weijun1,2?(權(quán)維?。?,CHEN Hongbin1(陳洪濱),HAN Xiuzhen3(韓秀珍),and MA Zhiqiang2(馬志強)

    1 Key Laboratory for Middle Atmosphere and Global Environment Observation,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    2 Environmental Meteorology Forecast Center of Beijing-Tianjin-Hebei,Beijing Meteorological Bureau,Beijing 100089

    3 National Satellite Meteorological Center,China Meteorological Administration,Beijing 100081

    To further verify the modified Becker's split-window approach for retrieving land surface temperature(LST)from long-term Advanced Very High Resolution Radiometer(AVHRR)data,a cross-validation and a radiance-based(R-based)validation are performed and examined in this paper.

    LST,split-window approach,MODTRAN,AVHRR,validation

    1.Introduction

    Land surface temperature(LST)is a key parameter in estimating the radiation budget on the land surface because it combines the results of all surfaceatmosphere interactions and energy fluxes between atmosphere and ground(Karnieli et al.,2010;Qian et al.,2013).It is closely related to soil moisture and canopy evaportraspiration and therefore,can be used for monitoring vegetation water stress,assessing surface energy balance,detecting soil drought,investigating urban climate and environment,and indicating land surface thermal conditions(Wan et al.,2004; Weng et al.,2004;Wang et al.,2007;Kustas and Anerson,2009;Weng,2009;Zhu et al.,2015).Furthermore,LST was assimilated into weather prediction models to improve the simulation of land-atmosphere exchange(Price,1982;Diak and Whipple,1993).Due to the strong heterogeneity of land surface characteristics such as vegetation,topography,and soil(Liu Yuanbo et al.,2006;Neteler,2010),and due to the fact that LST changes rapidly on both spatial and temporal scales,LST measurements with detailed spatial and temporal sampling are required in order to capture its distribution and temporal evolution.Given the complexity of surface temperature over land,ob-taining ground measurements is not practical over wide areas.With the development of remote sensing from space,satellite data offer the only possibility for measuring LST over the entire globe with sufficiently high temporal resolution and complete spatially averaged values rather than point values(Wang and Lu,2006;Li et al.,2013a).

    Supported by the Meteorological Key Technology Integration and Application Program of the China Meteorological Administration(CMAGJ2013Z20 and CMAGJ2015Z19).

    ?Corresponding author:quanquan78430@163.com.

    ?The Chinese Meteorological Society and Springer-Verlag Berlin Heidelberg 2015

    Up to now,several LST retrieval algorithms have been proposed in terms of various sensors aboard different satellites. These algorithms can be roughly grouped into three categories:the single channel algorithms(Ottl′e and Vidal-Madjar,1992;Qin et al.,2001;Jim′enez-Mu?noz and Sobrino,2003;Jim′enez-Mu?noz et al.,2009),the split-window algorithms(Becker,1987;Becker and Li,1990;Coll et al.,1994;Prata,1994;Sobrino et al.,1994;Wan and Dozier,1996;Tang et al.,2008;Quan et al.,2012),and the multi-angle algorithms(Prata,1993;Sobrino et al.,1996;Li et al.,2001),provided that the land surface emissivities(LSEs)are known a priori.If the LSEs are not known,the algorithms can be categorized into three types:stepwise retrieval method,simultaneous retrieval of LSEs and LST with known atmospheric information,and simultaneous retrieval with unknown atmospheric information(Li et al.,2013a).

    Validation of satellite derived LST product is also important because it not only provides useful information for users to understand the reliability and accuracy of LST product but also offers feedback for algorithm developers to further improve the retrieval algorithm. There are three methods,to date,which are generally used to verify LST retrieved from space:the temperature-based method(T-based),the radiance-based method(R-based),and the crossvalidation method.The T-based method fulfills verification by comparing the retrieved LST with the coincident ground measured LST(Wan et al.,2002;Coll et al.,2005;Pinker et al.,2009).Its main merit is that it provides a direct evaluation of the retrieved LST with the ground measured LST(Li et al.,2013a). However,the reliability of the T-based validation depends heavily on the accuracy of the ground LST values and how well they represent the LST on the satellite pixel scale.Because of pixel scale discrepancies between satellite and ground instruments and the heterogeneities of ground,validating retrieved LST by using the T-based method is a complex and difficult task.The R-based method is an advanced alternative method for validating space-based LST measurements(Wan and Li,2008;Coll et al.,2012).This method does not rely on ground-measured LST values but it does require the LSE spectra and the measured atmospheric profiles over the validation site at the time of the satellite overpass(Wan and Li,2008;Li et al.,2013a).The R-based method uses the satellitederived LST,the measured atmospheric profiles,and the LSEs as input parameters to an atmospheric radiative transfer model(RTM),such as the Moderate Resolution Transmittance Model(MODTRAN)(Berk et al.,2003),to simulate the radiance reaching the top of atmosphere(TOA)when the satellite overpasses.The inputting LST is not adopted until the simulated radiance matches the satellite-measured radiance.Only then is the input LST recognized as a true LST value and is further used as a criterion to evaluate the accuracy of the retrieved LST.The success of the R-based method depends on the accuracies of the atmospheric RTM,the atmospheric profiles,and the LSEs on pixel scale.The cross-validation method uses a well validated LST product,e.g.,the MYD11A1,as a reference and compares the satellite-derived LST to the reference product(Frey et al.,2012).The success of this method relies on the spatial-temporal matching accuracy and the reliability of the reference product.

    In our previous work,we have compared the Advanced Very High Resolution Radiometer(AVHRR)LST image retrieved from the modified Becker and Li's approach with the corresponding Moderate Resolution Imaging Spectroradiometer(MODIS)LST product over Beijing(Quan et al.,2012).To fully verify the performance of the modified Becker's approach,however,it is still necessary to validate the LST products of AVHRR retrieved from the modified Becker's approach by using diverse validation methods over wide regions.The current paper serves this purpose.In this paper,Section 2 reviews the modified Becker's split-window approach for retrieving LST from longterm AVHRR data.A brief description of the MODISLST/LSE products and the cross-validation procedure and results are also given in this section.Section 3 presents the procedure of the R-based method to validate the modified Becker's LST retrieval.The last section gives the conclusions and discussion.

    2.Cross-validation

    2.1 Scheme of cross-validation

    The key points of the cross-validation method depend on the matching accuracies in observing time,location,and view angle between two LST products as well as quality of the reference LST product.It is well known that the observing times for all pixels in one image are nearly the same while the local solar times for them are quite different.For instance,in a polar-orbit satellite image such as an MODIS image,the difference in the local solar time between the east and the west pixels along the scanning line can be up to 1.5 h,which means that the east pixels are exposed to solar irradiation approximately 1.5 h before the west ones if the sky is clear(Duan et al., 2014).As a result,the LST products derived from the same satellite cannot be compared if the differences in the local solar times of the pixels are significant.This phenomenon also affects LST products acquired by different satellites at different times,and significantly limits the applicability of the instantaneous LST products(Li et al.,2013a).Therefore,in this paper,a point-based comparison scheme is chosen instead of the image-based comparison to guarantee that each comparison is highly spatially-temporally matched.This minimizes potential discrepancies due to location,observation time,and view zenith angle.

    In this paper,548 weather stations are selected as validation sites(Fig.1).The land cover types include forests,shrublands,grasslands,croplands,urban areas,barren lands,desert,etc.(Fig.2).

    First,the LST pairs from the AVHRR LST and the MYD11A1 products for 2002-2008 are extracted according to the positions(longitude and latitude)of the validation sites.The LST pairs thus selected must satisfy four criteria:(1)the local solar time difference for each pixel pair should less than 10 min to minimizethe temporal match error as much as possible;(2)the quality control(QC)bit flag of the MYD11A1 must be equal to 0,which means that the MODIS LST value selected as reference has good data quality,is cloud free,has a small error in emissivity,and so on;(3)the AVHRR LST pixel must be cloud free;and(4)the view zenith angle difference between AVHRR and MODIS is less than 5°.Then,statistical analysis is performed to validate the LST pairs by using the MODIS MYD11A1 LST value as a reference to illustrate the accuracy and consistency of the AVHRR LST product.

    Fig.1.Locations of the validation sites(red dots)and radiosonde stations(blue rectangles).

    Fig.2.Land cover types(from Liu Yonghong et al.,2006).

    2.2 AVHRR LST products

    The AVHRR LST product is obtained by using the modified Becker's split-window approach(Quan et al.,2012)from long-term AVHRR data.The equation of the local split-window LST retrievalalgorithm proposed by Becker and Li(1990)is expressed as:

    where Tsis the LST;T4and T5are brightness temperatures in channels 4 and 5 of the AVHRR,respectively.The regression coefficients A0,P,and M are calculated by using the following equation: where ε=(ε4+ε5)/2 is the average emissivity,and Δε= (ε4-ε5)is the emissivity difference between channels 4 and 5.It is noted that the original regression coefficients A0,α,β,γ′,α′,and β′in Eq.(2)were calculated in terms of the spectral response functions(SRFs)of channels 4 and 5 of NOAA-9 AVHRR.The SRFs of different AVHRRs always have some discrepancies,such as spectral bandwidth and spectral sensitivity,which could result in different brightness temperatures even though the radiance into the sensors is the same.In other words,when the Becker's splitwindow algorithm is applied to the AVHRRs other than NOAA-9 AVHRR,the error of LST could be as large as 2.3 K(Czajkowski et al.,1998).To make the Becker's split-window algorithm applicable to various AVHRRs,the regression coefficients in Eq.(2)were recalculated by using the MODTRAN model simulations along with the SRFs of channels 4 and 5 for AVHRRs aboard nine NOAA satellites from NOAA-7 to NOAA-18(see Table 4 in Quan et al.,2012).

    2.3 MODIS LST and LSE products

    The MYD11A1 product of version V005 over China from 2002 to 2008 is selected as a reference to validate the AVHRR LST product. The MYD11A1 product is retrieved from the observations of MODIS aboard the Aqua satellite,whoseoverpass time is close to that of NOAA afternoon satellites. The MYD11A1 dataset is downloaded from the Land Processes Distributed Active Archive Center(LP DAAC),which is one of several discipline-specific data centers within the NASA Earth Observing System Data and Information System(EOSDIS)(https://lpdaac.usgs.gov/products/modisproducts-table/myd11a1). There are 12 scientific datasets in the MYD11A1 product.Besides the LST,quality control(QC),view time(local solar time),view angle,and clear-sky coverage,the emissivities in bands 31 and 32 are also included.The emissivities in bands 31 and 32 are estimated based on the land cover type in the pixel,which is determined by the input data of 250 m resolution land cover(MYD12Q1)and daily 1-km snow cover(MYD10-L2),using the land cover classification-based emissivity method(Snyder et al.,1998).

    The MYD11A1 product is generated from the level-2 LST product MYD11-L2,in which the LST is derived by using the generalized split-window LST algorithm(Wan and Dozier,1996):

    where Tsis the LST;T31and T32are brightness temperatures in channels 31 and 32 of the MODIS,respectively;ε= (ε31+ε32)/2 is the average emissivity,and Δε= (ε31-ε32)is the emissivity difference between channels 31 and 32 of the MODIS;C,A1,A2,A3,B1,B2,and B3are regression coefficients.

    2.4 Cross-validation results

    2.4.1 LST comparison results

    First,extracting is applied to the region of interest(ROI)along with the locations of 548 weather stations(denoted as red points in Fig.1)for the longterm AVHRR LST and MYD11A1 products for producing an LST pair file,in which each line consists of the LST pair and some auxiliary data that are extracted from the AVHRR and MYD11A1 at the same location,e.g.,the observing time and view zenith angle of AVHRR and MODIS,the emissivities in channels 4 and 5 of AVHRR,the emissivities in channels 31 and 32 of MODIS,the QC flag of MYD11A1,etc.Finally,3481 LST pairs that satisfy the criteria described in Section 2.1 are selected from this file and further analyzed.The AVHRR LSTs are higher than the corresponding MYD11A1 in most of the cases(Fig.3a). On the other hand,the AVHRR LSTs,adjusted with a linear regression equation,are close to the MYD11A1 values(Fig.3b).It can be seen from Fig.3b that a good linear relationship exists between the AVHRR LST and the MYD11A1 product.The determination coefficient of regression(R2)is equal to 0.91,the bias is -0.03 K,and the root-mean-square-deviation(RMSD)is 3.0 K.

    To further illustrate the performance of the modified Becker's split-window approach,the validation LST pairs are grouped into four categories according to the seasons. Simultaneously,in each category,six sub-categories are also generated in terms of the land cover type,i.e.,forests(evergreen needleleaf forest,evergreen broadleaf forest,deciduous needleleaf forest,deciduous broadleaf forest,and mixed forest),shrublands(closed shrublands and open shrublands),grasslands(woody savannas,savannas,grasslands,and high-cold steppe),croplands(single cropping system,double cropping system,dry and paddy cropping system,and double cropping paddy),deserts(barren lands and desert),and others(urban areas,marsh wetlands,coastal wetlands,water bodies,and snow and ice).The statistical results for the crossvalidation of each category and sub-category are listed in Table 1.

    Table 1 indicates that the biggest difference occurs on desert lands for all seasons,for which the AVHRR LST is at least 1.5 K less than the MYD11A1. This phenomenon was also found on bare soil in the previous investigation,in which the errors in the general split-window(GSW)retrieved MODIS LSTs can be greater than near surface temperature by more than 16 K(Wan and Li,2008).On the other hand,the AVHRR LSTs are greater than the MYD11A1 by 0.9 and 0.7 K on forest land in spring and winter,respectively.The LST differences on croplands are close to0 in all seasons except winter,in which the LST difference has a value of 2.0 K.However,this value is not reliable due to the small sample number.Due to the larger seasonal variations of grasslands,the LST difference on grasslands is less than-0.1 K in spring and autumn but greater than 0.7 K in summer and winter,respectively.Table 1 also shows that the RMSD between the AVHRR LST and MYD11A1 ranges from1.6 to 4.1 K,and the mean correlation coefficients(R)for all land cover types are 0.932,0.761,0.957,and 0.954 in spring,summer,autumn,and winter,respectively.The lower correlation coefficient between the AVHRR LST and the MYD11A1 for summer may stem from total water vapor amount(WV)being at its maximum during this season and the different sensitivities of the two retrieval algorithms to WV.

    Fig.3.Scatter plots of MYD11A1 versus AVHRR LST derived from the modified Becker's split-window approach.(a)Before regression and(b)after regression.

    Table 1.Statistical results of the cross comparison between the MYD11A1 and AVHRR LST products for different land cover types

    2.4.2 LSE comparison results

    It is necessary to explain the spectral characters of the thermal-infrared(TIR)channels of AVHRR and MODIS before comparing their LSE products.The SRFs of Aqua/MODIS and NOAA-18/AVHRR(represents AVHRR sensors)are plotted in Fig.4.

    Figure 4 indicates that the spectral width of the MODIS TIR channel is narrower than that of the corresponding AVHRR TIR channel.Additionally,there is almost no overlap between the two adjacent MODIS TIR channels while somewhat overlap exists between channels 4 and 5 of AVHRR.On the other hand,the central wavelength of channel 31 of MODIS tends towards a longer wavelength relative to that of channel 4 of AVHRR.

    It should be noted that the LSE estimation methods used in the MYD11A1 and AVHRR LST products are quite different.The land surface emissivities adopted in retrieving MODIS LST are estimated with the classification-based emissivity method(CBEM),which assumes that similarly classified land covers exhibit similar LSEs.On contrast,the land surface emissivities used in deriving AVHRR LST are calculated based on the normalized difference vegetation index(NDVI)value,with a so-called NDVI-based emissivity method(NBEM).The NBEM takes advantage of a statistical relationship between the NDVI calculated from the visible-near infrared(VNIR)data and the LSE in the TIR channels(Van De Griend and Owe,1993;Valor and Caselles,1996;Sobrino and Raissouni,2000).

    The reliability of CBEM depends on accuracy of land cover classification,prior knowledge for each class,and precision of the component(endmember)weight within a satellite pixel with coarse resolution(Li et al.,2013b).Snyder et al.(1998)pointed out that the accuracy of CBEM would degrade as a result of the uncertainties in the soil moisture,the annual biophysical cycle of vegetation,and the appearance of snow and ice.In some cases,there is not a direct correlation between the land cover based on the VNIR classification and the LSE in the TIR channel.

    Since different formulae are adopted in the calculation of LSEs for soil,mixed soil,vegetation,and dense vegetation types,the main deficiency of the NBEM is the lack of continuity in the LSE values during the land cover type transitioning from soil to vegetation(Sobrino et al.,2008).Difficulty in obtaining a priori knowledge of the emissivities of soil and vegetation is another souce of error in this method(Jim′enez-Mu?noz and Sobrino,2006).

    Owing to the different LSE estimation methodsand the SRF discrepancies in TIR channels between MODIS and AVHRR,it seems unsound to compare the LSEs applied in the MYD11A1 and AVHRR LST products.In this paper,scatter points of the difference in LSE for 3481 validation LSE pairs are plotted in Fig.5 to help users roughly understand the relationships between them.It can be seen from Fig. 5 that the points are symmetrically distributed along the average lines,though the point distribution patterns are dispersive.Meanwhile,the general trends of the LSE differences are also plotted in Fig.5 with red lines,in which the data points are smoothed by using a moving average.Generally speaking,the LSE difference between AVHRR channel 4 and MODIS channel 31 is larger than that between AVHRR channel 5 and MODIS channel 32,and the average LSE difference(standard deviation)is-0.0092(0.0104)and-0.0038(0.0058),respectively.

    Fig.4.Aqua/MODIS and NOAA-18/AVHRR spectral response functions.

    3.Radiance-based validation

    3.1 Scheme of R-based validation

    Wan and Li(2008)presented a radiance-based validation method to verify the MODIS LST V5 product.The core steps of this method contain three parts:1)calculate the TOA radiance in MODIS channel 31 by using the MODTRAN 4 code with the assumption that the surface emissivity spectrum and the atmospheric profiles have been precisely measured;2)simulate the radiance in MODIS channel 31 using the MODTRAN model by adjusting the input LST value. If the simulated radiance in MODIS channel 31 is equal to that measured by the MODIS sensor,the last input LST is considered the true LST value;and 3)comparethe retrieved MODIS LST with the last input LST in the second step.This method is suitable for validating the LST product in both day and night conditions and does not rely on ground-measured LST values.

    Fig.5.Scatter plots of LSE differences(a)between AVHRR channel 4 and MODIS channel 31 and(b)between AVHRR channel 5 and MODIS channel 32.Blue points denote the differences in LSE between AVHRR and MODIS,and red lines are the results of those points after being smoothed by a moving average.Mean and Stddev are the mean value and standard devation calculated from the smoothed values,respectively.

    The flow chart of the R-based validation method employed in this study is shown in Fig.6.First,the AVHRR LST product is derived from the AVHRR L1B data by using the modified Becker's split-window algorithm. Second,the TOA brightness temperature in AVHRR channel 4 and retrieved LST at four radiosonde stations(Beijing,Chifeng,Qingdao,and Zhangqiu)are extracted by using the ROI extracting method along with the positions of the stations from the AVHRR L1B and the AVHRR LST product,respectively.Third,atmospheric temperature and humidity profiles for the whole atmoshpere at the radiosonde stations are rebuilt by assembling the radiosonde balloon measurements,whose maximum detecting altitudes are lower than 30 km in most cases,and the standard atmospheric profiles with the altitude higher than 30 km.Fourth,the land surface emissivity estimated with the NDVI by using the NBEM method,the atmospheric profiles,and the initial surface temperature(Ts)are input into the MODTRAN 4.1 model to calculate the radiance reaching the TOA. The radiance is then further converted to the brightness temperature by using the inverse of Planck's radiation equation(Robel,2007,2009).This step does not stop until the simulated brightness temperature matches that observed by the AVHRR at the stations. Then,the input Tsis recognized as the true value of LST,i.e.,the inverted LST.Finally,the statistical results are given on the basis of the comparison between the inverted and the retrieved LST data.

    3.2 Atmospheric profiles

    Accurate temperature profile and total water vapor amount are two key parameters in the R-based LST validation method. It is worth noting that the observation times of routine radiosonde balloons are usually fixed at 0000 and 1200 UTC,which differs from the overpass times of the afternoon NOAA satellites.Fortunately,two additional observations at 0200 and 1400 BT(Beijing Time)during summer were performed at some radiosonde stations in China.In this paper,21 atmospheric profiles at 4radiosonde stations,observing at 1400 BT under clear sky conditions,are selected as input parameters into the MODTRAN and plotted in Fig.7.The four radiosonde stations(denoted with blue rectangles in Fig. 1)include Beijing(54511),Chifeng(54218),Qingdao(54857),and Zhangqiu(54727)stations.The differ-ences between the radiosonde balloon observation time and the satellite overpassing time are limited to 30 min in order to guarantee that these atmospheric profiles can represent the atmospheric conditions when the satellites are overpassing.These radiosonde balloon data come from the mandatory pressure level dataset of China Radiosonde Archive version 1.0(CRA-M1.0),produced by the National Meteorological Information Center of the China Meteorological Administration.At the same time,the WV values are also calculated by integrating the humidity profiles for the whole atmosphere(Yang and Qiu,2002),and are listed in Table 2.

    Fig.6.Flow chart of the R-based method to validate AVHRR LST product.

    Fig.7.Atmospheric temperature profiles used in the R-based validation.(a)Beijing,(b)Chifeng,(c)Qingdao,and(d)Zhangqiu.ST-ID is the station ID,and the observing times for these profiles are also denoted,e.g.,2007072714 is 1400 BT 27 July 2007.

    Figure 7 shows all atmospheric profiles that were obtained in the summer from 2006 to 2012.The maximum sounding altitudes in these cases are below 30 km except the one observed at Chifeng on 26 August 2007(Fig.7b),which was 31.4 km.Additionally,there were no clouds or temperature inversions occurring in these cases.

    Table 2.Results of the R-based LST validations

    3.3 Estimation of surface emissivity

    Surface emissivity is another important input to the MODTRAN model.In this paper,the LSEs at the radiosonde stations are estimated by using the NBEM method with the NDVI values.For details,refer to Section 4.1 in Quan et al.(2012).

    3.4 Results

    The overall results of the R-based validation are listed in Table 2.It can be seen that the satellite overpassing times in all cases are close to the radiosonde time(1400 BT).The WV values range from 7.1 to 44.7 mm,i.e.,the atmosphere is relatively moist.Additionally,the atmospheric transmittance,which ranges from 0.486 to 0.928,generally decreases with the increase of WV.The brightness temperature in channel 4 of AVHRR at the validation sites is extracted from the AVHRR channel 4 images.The LSE in channel 4 of AVHRR,which is estimated from the NDVI value by using the NBEM method,is also listed.

    Table 2 indicates that the retrieved LSTs from the AVHRR data by using the modified Becker's splitwindow approach are generally less than the inverted LSTs,especially in cases when the WV is greater than 15 mm.To further understand the water vapor impact on the retrieved LST,the average ΔT is calculated forhigh and low WV values.The result shows that the average ΔT is-1.3 K when the WV is less than 20 mm and up to-2.5 K when the WV is greater than 20 mm. Besides the WV,the LSEs and the view angles also influence the accuracy of retrieved LST.Therefore,the atmosphere-independent coefficients in the modified Becker's split-window algorithm are only suitable for the scenarios in which the atmosphere is not too wet(WV<20 mm).To make the split-window algorithm applicable to more general atmospheric conditions,the regression coefficients in the modified Becker's algorithm need to be considered as functions of WV,view zenith angle,etc.

    4.Conclusions and discussion

    Validation of the LST retrieved from the satellite observations is an important step towards improving the accuracy of the retrieval algorithm and implemeting better use of the LST product.However,it is a hard task due to many variables such as the true surface temperature,the accurate pixel matching,and the suitable atmospheric profiles.In this paper,we take advantage of two different methods,i.e.,the crossvalidation method and the R-based validation method,to verify the LST product retrieved from long-term AVHRR data by using the modified Becker's splitwindow algorithm(Quan et al.,2012).

    The cross-validations are carried out via comparison between 3481 LST validation pairs. Using data from 548 weather stations in China from 2002 to 2008,these LST validation pairs are extracted from the AVHRR LST and the corresponding MYD11A1 products. The results show that there is a strong linear relationship between the AVHRR LST and the MYD11A1,with a bias of 0.03 K and a RMSD of 3.0 K.The cross-validation method takes advantage of the QC bit flag in the MYD11A1 to choose the high quality MYD11A1 LST as a reference.However,the MYD11A1 product is also retrieved surface temperature rather than a true value of the LST.Thus,this comparison can only be used to indicate the similarity between the AVHRR LST product and the reference MODIS LST product.

    On the other hand,the R-based LST validation method does not need in-situ LST measurements,so it is suitable for the validation of both daytime and nighttime LSTs. To apply the R-based validation approach,accurate atmospheric temperature and water vapor profiles around the satellite overpass time are needed as input parameters into the MODTRAN model in order to simulate the satellite sensor's measurements. Fortunately,a few Chinese radiosonde stations provided additional observations at 1400 BT during summer,which made the accurate atmospheric profiles around the satellite overpass time available. In this paper,21 accurate atmospheric temperature and humidity profiles at 4 radiosonde stations are used along with the estimated LSE as the input parameters into the MODTRAN 4.1 model to fulfill the R-based validation.The result shows that the AVHRR LST retrieved from the modified Becker's split-window approach underestimates the LST in most of the cases and with an average LST difference of-1.3 K when the WV is less than 20 mm,and-2.5 K when the WV is larger than 20 mm.

    In the future,considerable work needs to be conducted to improve the accuracy and consistency of the LST retrieval algorithm.For example,a split-window LST retrieval algorithm can be developed to derive the LST from long time series AVHRR data,taking into consideration the view zenith angle,LSE,WV,etc. Of course,the LST validation method and validation scheme are still necessary to guarantee that the LST retrieval algorithm is reliable.

    Acknowledgments.We would like to thank the anonymous reviewers and editors for their valuable and stimulating comments that have greatly improved the paper.We also appreciate the Jet Propulsion Laboratory,California Institute of Technology of USA for providing the ASTER spectral library,which contains the John Hopkins University spectral emissivity data used in this paper.We are also very grateful to Mrs. Hu Tianjie and Mr.Ma Jingjin from the Beijing Meteorological Bureau for providing the radiosonde data.

    Becker,F(xiàn).,1987:The impact of spectral emissivity on the measurement of land surface temperature froma satellite.Int.J.Remote Sens.,8,1509-1522.

    Becker,F(xiàn).,and Li Zhaoliang,1990:Towards a local split window method over land surfaces.Int.J.Remote Sens.,11,369-393.

    Berk,A.,G.P.Anderson,P.K.Acharya,et al.,2003:MODTRAN4 Version 3 Revision 1 User's Manual,Air Force Research Laborary,North Andover,MA,USA,10-35.

    Coll,C.,V.Caselles,J.A.Sobrino,et al.,1994:On the atmospheric dependence of the split-window equation for land surface temperature.Int.J.Remote Sens.,15,105-122.

    Coll,C.,V.Caselles,J.M.Galve,et al.,2005:Ground measurements for the validation of land surface temperatures derived from AATSR and MODIS data. Remote Sens.Environ.,97,288-300.

    Coll,C.,E.Valor,J.M.Galve,et al.,2012:Long-term accuracy assessment of land surface temperatures derived from the Advanced Along-Track Scanning Radiometer.Remote Sens.Environ.,116,211-225.

    Czajkowski,K.P.,S.N.Goward,and H.Ouaidrari,1998:Impact of AVHRR filter functions on surface temperature estimation from the split window approach. Int.J.Remote Sens.,19,2007-2012.

    Diak,G.R.,and M.S.Whipple,1993:Improvements to models and methods for evaluating the land-surface energy balance and‘effective'roughness using radiosonde reports and satellite-measured‘skin'temperature data.Agric.Forest Meteor.,63,189-218.

    Duan,S.B.,Z.L.Li,B.H.Tang,et al.,2014:Generation of a time-consistent land surface temperature product from MODIS data.Remote Sens.Environ.,140,339-349.

    Frey,C.M.,C.Kuenzer,and S.Dech,2012:Quantitative comparison of the operational NOAA-AVHRR LST product of DLR and the MODIS LST product V005. Int.J.Remote Sens.,33,7165-7183.

    Jim′enez-Mu?noz,J.C.,and J.A.Sobrino,2003:A generalized single-channel method for retrieving land surface temperature from remote sensing data.J. Geophys.Res.,108,doi:10.1029/2003JD003480.

    Jim′enez-Mu?noz,J.C.,and J.A.Sobrino,2006:Error sources on the land surface temperature retrieved from thermal infrared single channel remote sensing data.Int.J.Remote Sens.,27,999-1014.

    Jim′enez-Mu?noz,J.C.,J.Cristobal,J.A.Sobrino,et al.,2009:Revision of the single-channel algorithm for land surface temperature retrieval from Landsat thermal-infrared data.IEEE Trans.Geosci.Remote Sens.,47,339-349.

    Karnieli,A.,N.Agam,R.T.Pinker,et al.,2010:Use of NDVI and land surface temperature for drought assessment:Merits and limitations.J.Climate,23,618-633.

    Kustas,W.,and M.Anderson,2009:Advances in thermal infrared remote sensing for land surface modeling. Agric.Forest Meteor.,149,2071-2081.

    Li Zhaoliang,M.P.Stoll,Zhang Renhua,et al.,2001:On the separate retrieval of soil and vegetation temperatures from ATSR data.Sci.China(Ser.D),44,97-111.

    Z.H.Li,B.H.Tang,H.Wu,et al.,2013a:Satellitederived land surface temperature:Current status and perspectives.Remote Sens.Environ.,131,14-37.

    Li,Z.L.,H.Wu,N.Wang,et al.,2013b:Land surface emissivity retrieval from satellite data.Int.J.Remote Sens.,34,3084-3127.

    Liu Yuanbo,T.Hiyama,Y.Yamaguchi,2006:Scaling of land surface temperature using satellite data:A case examination on ASTER and MODIS products over a heterogeneous terrain area. Remote Sens. Environ.,105,115-128.

    Liu Yonghong,Niu Zheng,Xu Yongming,et al.,2006:Design of land cover classification system for China and its application research based on MODIS data. Trans.CSAE,22,99-104.(in Chinese)

    Neteler,M.,2010:Estimating daily land surface temperatures in mountainous environments by reconstructed MODIS LST data.Remote Sens.,2,333-351.

    Ottl′e,C.,and D.Vidal-Madjar,1992:Estimation of land surface temperature with NOAA9 data.Remote Sens.Environ.,40,27-41.

    Pinker,R.T.,D.L.Sun,M.P.Hung,et al.,2009:Evaluation of satellite estimates of land surface temperature from GOES over the United States.J.Appl. Meteor.Climatol.,48,167-180.

    Prata,A.J.,1993:Land surface temperatures derived from the advanced very high resolution radiometer and the along-track scanning radiometer.1:Theory. J.Geophys.Res.,98,16689-16702.

    Prata,A.J.,1994:Land surface temperature determination from satellites.Adv.Spac.Res.,14,15-26.

    Price,J.C.,1982:On the use of satellite data to infer surface fluxes at meteorological scales.J.Appl.Meteor.,21,1111-1122.

    Qian,Y.G.,Z.L.Li,and F.Nerry,2013:Evaluation of land surface temperature and emissivities retrieved from MSG/SEVIRI data with MODISland surfacetemperatureand emissivity products.Int.J.Remote Sens.,34,3140-3152,doi:10.1080/01431161.2012.716538.

    Qin,Z.,A.Karnieli,and P.Berliner,2001:A monowindow algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region.Int.J.Remote Sens.,22,3719-3746.

    Quan Weijun,Chen Hongbin,Han Xiuzhen,et al.,2012:A modified Becker's split-window approach for retrieving land surface temperature from AVHRR and VIRR.Acta Meteor.Sinica,26,229-240,doi:10.1007/s13351-012-0208-y.

    Robel,J.,2007:AVHRR Level 1b data base.NOAA polar orbiter data user's guide,U.S.Department of Commerce,Asheville,USA,22-28.

    Robel,J.,2009:Calibration of NOAA KLM instruments.NOAA KLM user's guide with NOAA-N,-N'Supplement,National Oceanic and Atmospheric Administration,Asheville,USA,1-11.

    Snyder,W.C.,Z.Wan,Y.Zhang,et al.,1998:Classification-based emissivity for land surface temperature measurement from space.Int.J.Remote Sens.,19,2753-2774.

    Sobrino,J.A.,Z.L.Li,M.P.Stoll,et al.,1994:Improvements in the split-window technique for land surface temperature determination.IEEE Trans.Geosci. Remote Sens.,32,243-253.

    Sobrino,J.A.,Z.L.Li,M.P.Stoll,et al.,1996:Multichannel and multi-angle algorithms for estimating sea and land surface temperature with ATSR data. Int.J.Remote Sens.,17,2089-2114.

    Sobrino,J.A.,and N.Raissouni,2000:Toward remote sensing methods for land cover dynamic monitoring:Application to Morocco.Int.J.Remote Sens.,21,353-366.

    Sobrino,J.A.,J.A.Jim′enez-Mu?noz,C.S`oria,et al.,2008:Land surface emissivity retrieval from different VNIR and TIR sensors.IEEE Trans.Geosci. Remote Sens.,46,316-327.

    Tang Bohui,Bi Yuyun,Li Zhaoliang,et al.,2008:Generalized split-window algorithm for estimate of land surface temperature from Chinese geostationary FengYun meteorologicalsatellite(FY-2C)data. Sensors,8,933-951.

    Valor,E.,and V.Caselles,1996:Mapping land surface emissivity from NDVI:Application to European,African and South American areas.Remote Sens. Environ.,57,167-184.

    Van De Griend,A.A.,and M.Owe,1993:On the relationship between thermal emissivity and the normalized difference vegetation index for natural surfaces. Int.J.Remote Sens.,14,1119-1131.

    Wan,Z.M.and J.Dozier,1996:A generalized splitwindow algorithm for retrieving land-surface temperature from space.IEEE Trans.Geosci.Remote Sens.,34,892-905.

    Wan,Z.M.,Y.L.Zhang,Q.C.Zhang,et al.,2002:Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data.Remote Sens.Environ.,83,163-180.

    Wan,Z.M.,P.Wang,and X.Li,2004:Using MODIS land surface temperature and normalized difference vegetation index products for monitoring drought in the southern Great Plains,USA.Int.J.Remote Sens.,25,61-72.

    Wan,Z.M.,and Z.L.Li,2008:Radiance-based validation of the V5 MODIS land-surface temperature product.Int.J.Remote Sens.,29,5373-5395.

    Wang,K.C.,P.C.Wang,Z.Q.Li,et al.,2007:A simple method to estimate actual evapotranspiration from a combination of net radiation,vegetation index,and temperature.J.Geophys.Res.,112,D15107,doi:10.1029/2006JD008351.

    Wang Minyan and L¨u Daren,2006:Diurnal and seasonal variation of clear-sky land surface temperature of several representative land surface types in China retrieved by GMS-5.Acta Meteor.Sinica.,20,475-488.

    Weng Qihao,2009:Thermal infrared remote sensing for urban climate and environmental studies:Methods,applications,and trends.ISPRS J.Photogramm. Remote Sens.,64,335-344.

    Weng,Q.H.,D.S.Lu,and J.Schubring,2004:Estimation of land surface temperature-vegetation abundance relationship for urban heat island studies. Remote Sens.Environ.,89,467-483.

    Yang Jingmei and Qiu Jinhuan,2002:A method for estimating precipitable water and effective water vapor content from ground humidity parameters.Chinese J.Atmos.Sci.,26,9-22.(in Chinese)

    Zhu Xiaochen,Qiu Xinfa,Zeng Yan,et al.,2015:A remote sensing model to estimate sunshine duration in the Ningxia Hui autonomous region,China.J. Meteor.Res.,29,144-154,doi:10.1007/s13351-015-4059-1.

    Quan Weijun,Chen Hongbin,Han Xiuzhen,et al.,2015:Validation of the modified Becker's splitwindow approach for retrieving land surface temperature from AVHRR.J.Meteor.Res.,29(5),823-836,

    10.1007/s13351-015-4110-2.

    (Received November 6,2014;in final form July 18,2015)

    In the cross-validation,3481 LST data pairs are extracted from the AVHRR LST product retrieved with the modified Becker's approach and compared with the Moderate Resolution Imaging Spectroradiometer(MODIS)LST product(MYD11A1)for the period 2002-2008,relative to the positions of 548 weather stations in China.The results show that in most cases,the AVHRR LST values are higher than the MYD11A1. When the AVHRR LSTs are adjusted with a linear regression,the values are close to the MYD11A1,showing a good linear relationship between the two datasets(R2=0.91).

    In the R-based validation,comparison is made between AVHRR LST retrieved from the modified Becker's approach and the inversed LST from the Moderate Resolution Transmittance Model(MODTRAN)consolidated with observed temperature and humidity profiles at four radiosonde stations.The results show that the retrieved AVHRR LST deviates from the MODTRAN inversed LST by-1.3(-2.5)K when the total water vapor amount is less(larger)than 20 mm.This provides useful hints for further improvement of the LST retrieval algorithms'accuracy and consistency.

    猜你喜歡
    志強
    NFT與絕對主義
    趙志強書法作品
    學習“集合”,學什么
    李志強·書法作品稱賞
    袁志強 始終奮戰(zhàn)在防疫第一線
    盧志強 用心于畫外
    海峽姐妹(2019年4期)2019-06-18 10:39:00
    送別張公志強
    寶藏(2018年12期)2019-01-29 01:50:50
    ON ENTIRE SOLUTIONS OF SOME TYPE OF NONLINEAR DIFFERENCE EQUATIONS?
    Analysis of Tibetan Plateau Vortex Activities Using ERA-Interim Data for the Period 1979-2013
    志強的石
    中華奇石(2014年12期)2014-07-09 18:30:22
    日本 欧美在线| 国产精品综合久久久久久久免费| www.精华液| av超薄肉色丝袜交足视频| 欧美日韩黄片免| 精品国产乱子伦一区二区三区| 国产欧美日韩一区二区精品| 国产成人欧美在线观看| 在线免费观看的www视频| 国产三级在线视频| 99久久综合精品五月天人人| 夜夜爽天天搞| 免费高清视频大片| 少妇粗大呻吟视频| 免费在线观看视频国产中文字幕亚洲| 午夜精品一区二区三区免费看| 欧美黄色淫秽网站| 亚洲国产欧美一区二区综合| 婷婷六月久久综合丁香| 国产亚洲精品久久久久久毛片| 亚洲精品av麻豆狂野| 91麻豆av在线| 别揉我奶头~嗯~啊~动态视频| 亚洲人成77777在线视频| 两个人免费观看高清视频| 好男人电影高清在线观看| 妹子高潮喷水视频| 怎么达到女性高潮| 久久午夜亚洲精品久久| 少妇裸体淫交视频免费看高清 | 一a级毛片在线观看| 午夜免费观看网址| 黑人欧美特级aaaaaa片| 国产三级黄色录像| 亚洲专区字幕在线| 免费看美女性在线毛片视频| 最近在线观看免费完整版| 国产欧美日韩一区二区精品| 级片在线观看| 免费无遮挡裸体视频| ponron亚洲| 2021天堂中文幕一二区在线观| 草草在线视频免费看| 九九热线精品视视频播放| 最近最新中文字幕大全免费视频| 操出白浆在线播放| 男男h啪啪无遮挡| 又紧又爽又黄一区二区| 国内揄拍国产精品人妻在线| 欧美在线一区亚洲| 91麻豆av在线| 国产精品乱码一区二三区的特点| 亚洲成人精品中文字幕电影| 亚洲人成77777在线视频| 后天国语完整版免费观看| 蜜桃久久精品国产亚洲av| 亚洲男人的天堂狠狠| 18禁黄网站禁片免费观看直播| 亚洲午夜精品一区,二区,三区| 亚洲黑人精品在线| 日本熟妇午夜| 亚洲美女黄片视频| 九色国产91popny在线| 国产一区二区三区视频了| 特大巨黑吊av在线直播| 亚洲精品国产精品久久久不卡| 国产精品99久久99久久久不卡| 女人被狂操c到高潮| 一个人免费在线观看的高清视频| 日韩欧美国产一区二区入口| 欧美日本视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美一区二区三区黑人| 成人av在线播放网站| 国产高清videossex| 又黄又爽又免费观看的视频| 国产蜜桃级精品一区二区三区| 成人午夜高清在线视频| 亚洲精品国产一区二区精华液| 欧美 亚洲 国产 日韩一| 神马国产精品三级电影在线观看 | 亚洲成a人片在线一区二区| 99精品在免费线老司机午夜| 每晚都被弄得嗷嗷叫到高潮| 国产精品日韩av在线免费观看| 国产成人影院久久av| 亚洲男人的天堂狠狠| 欧美激情久久久久久爽电影| 99久久精品热视频| 高清毛片免费观看视频网站| 高潮久久久久久久久久久不卡| 国产视频一区二区在线看| 久久精品aⅴ一区二区三区四区| 18禁美女被吸乳视频| 国产日本99.免费观看| 欧美乱妇无乱码| 91在线观看av| 亚洲va日本ⅴa欧美va伊人久久| 国产精华一区二区三区| 亚洲激情在线av| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲精品av在线| 免费在线观看成人毛片| 99精品欧美一区二区三区四区| 国产成人精品无人区| 狂野欧美激情性xxxx| 国产亚洲精品一区二区www| 午夜福利在线在线| 黄频高清免费视频| 国产精品永久免费网站| 亚洲av美国av| 老司机午夜十八禁免费视频| 在线观看午夜福利视频| 午夜福利在线观看吧| 久久久久久久午夜电影| 亚洲国产精品sss在线观看| 亚洲国产精品合色在线| 国内久久婷婷六月综合欲色啪| 青草久久国产| 国内毛片毛片毛片毛片毛片| 99久久久亚洲精品蜜臀av| www.熟女人妻精品国产| 亚洲欧美精品综合一区二区三区| 天堂√8在线中文| 真人一进一出gif抽搐免费| 成人18禁高潮啪啪吃奶动态图| 又紧又爽又黄一区二区| 中文字幕人妻丝袜一区二区| 热99re8久久精品国产| 看片在线看免费视频| 免费在线观看影片大全网站| 国产精品自产拍在线观看55亚洲| 欧洲精品卡2卡3卡4卡5卡区| 精品国产亚洲在线| 欧美日韩中文字幕国产精品一区二区三区| 成年女人毛片免费观看观看9| 欧美乱色亚洲激情| 精品国产美女av久久久久小说| 色精品久久人妻99蜜桃| 亚洲国产中文字幕在线视频| 亚洲 欧美 日韩 在线 免费| 一a级毛片在线观看| 久久精品综合一区二区三区| 午夜精品久久久久久毛片777| 欧美丝袜亚洲另类 | 99久久99久久久精品蜜桃| 日韩欧美国产一区二区入口| 国产精品 国内视频| 欧美久久黑人一区二区| 免费看a级黄色片| 成人午夜高清在线视频| 淫妇啪啪啪对白视频| 亚洲专区中文字幕在线| av视频在线观看入口| 欧美中文综合在线视频| 午夜福利高清视频| 最近最新免费中文字幕在线| 精品国内亚洲2022精品成人| 18禁黄网站禁片免费观看直播| 欧美成人一区二区免费高清观看 | 我的老师免费观看完整版| 国产三级黄色录像| √禁漫天堂资源中文www| 欧美三级亚洲精品| 久久久国产精品麻豆| 99久久综合精品五月天人人| 老司机靠b影院| 我要搜黄色片| 久久国产精品人妻蜜桃| 午夜福利18| 久久精品国产亚洲av高清一级| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 久久久精品欧美日韩精品| 老司机福利观看| 一边摸一边做爽爽视频免费| 国产亚洲精品久久久久5区| 日本三级黄在线观看| 亚洲电影在线观看av| 色综合站精品国产| 久久伊人香网站| 欧美一区二区精品小视频在线| 欧美日韩亚洲综合一区二区三区_| 国产成人啪精品午夜网站| 亚洲精品在线美女| 亚洲av熟女| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩 | 久久天堂一区二区三区四区| 免费在线观看成人毛片| 三级男女做爰猛烈吃奶摸视频| 国产精品九九99| 国产精品亚洲av一区麻豆| 一二三四在线观看免费中文在| 国产欧美日韩精品亚洲av| av片东京热男人的天堂| 亚洲电影在线观看av| 天天添夜夜摸| 亚洲精品国产一区二区精华液| 成人av一区二区三区在线看| 听说在线观看完整版免费高清| 欧美日韩乱码在线| 免费在线观看影片大全网站| 久久精品人妻少妇| 一进一出好大好爽视频| 国产成人精品无人区| 久久久久国产一级毛片高清牌| 韩国av一区二区三区四区| 97碰自拍视频| 国产一区二区激情短视频| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 熟女电影av网| 在线永久观看黄色视频| 99久久无色码亚洲精品果冻| 国产精品久久电影中文字幕| 亚洲国产精品合色在线| 99re在线观看精品视频| 男女之事视频高清在线观看| 亚洲av中文字字幕乱码综合| 久久久久久久久久黄片| 欧美另类亚洲清纯唯美| 午夜福利在线在线| 99精品久久久久人妻精品| 亚洲午夜精品一区,二区,三区| 国产av不卡久久| 亚洲男人的天堂狠狠| 亚洲全国av大片| 69av精品久久久久久| 看片在线看免费视频| 久久这里只有精品中国| 精品无人区乱码1区二区| 成年版毛片免费区| 在线观看美女被高潮喷水网站 | 日本a在线网址| 国产精品综合久久久久久久免费| 亚洲美女视频黄频| 国产黄色小视频在线观看| 国产激情久久老熟女| 成人av一区二区三区在线看| 一级片免费观看大全| 国产黄a三级三级三级人| 国产亚洲欧美98| 精品国产超薄肉色丝袜足j| 国产精品自产拍在线观看55亚洲| 男女那种视频在线观看| 日韩成人在线观看一区二区三区| 毛片女人毛片| 亚洲精品一卡2卡三卡4卡5卡| 视频区欧美日本亚洲| 法律面前人人平等表现在哪些方面| 亚洲精品在线观看二区| 久9热在线精品视频| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 99精品在免费线老司机午夜| 美女高潮喷水抽搐中文字幕| 伦理电影免费视频| 亚洲国产看品久久| xxxwww97欧美| 亚洲 欧美 日韩 在线 免费| 狠狠狠狠99中文字幕| 日日夜夜操网爽| 久久这里只有精品中国| 人妻久久中文字幕网| 欧美+亚洲+日韩+国产| 男插女下体视频免费在线播放| 亚洲最大成人中文| 听说在线观看完整版免费高清| 欧美一级a爱片免费观看看 | 精品久久久久久成人av| 亚洲乱码一区二区免费版| 香蕉国产在线看| 午夜激情av网站| 嫩草影院精品99| 久久久久久久久免费视频了| 国产精品久久久人人做人人爽| 国产成人欧美在线观看| 久久久久久人人人人人| 两性夫妻黄色片| 又黄又爽又免费观看的视频| а√天堂www在线а√下载| 人人妻人人看人人澡| 亚洲,欧美精品.| 国产激情久久老熟女| 国产亚洲精品第一综合不卡| 亚洲欧美激情综合另类| 午夜视频精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲电影在线观看av| 男女视频在线观看网站免费 | 久久天躁狠狠躁夜夜2o2o| 久久精品亚洲精品国产色婷小说| 国产精品精品国产色婷婷| 亚洲欧美日韩高清在线视频| 国内揄拍国产精品人妻在线| 首页视频小说图片口味搜索| 久久欧美精品欧美久久欧美| 制服诱惑二区| 亚洲一区二区三区色噜噜| 免费看十八禁软件| 一二三四在线观看免费中文在| 午夜亚洲福利在线播放| 搞女人的毛片| 男女之事视频高清在线观看| 超碰成人久久| 国内毛片毛片毛片毛片毛片| 亚洲中文日韩欧美视频| 长腿黑丝高跟| 成人18禁在线播放| 久久精品aⅴ一区二区三区四区| 免费看日本二区| 好男人电影高清在线观看| 黄片大片在线免费观看| 黑人巨大精品欧美一区二区mp4| 国产伦人伦偷精品视频| 99国产精品一区二区三区| 成在线人永久免费视频| 亚洲国产欧美网| 黄色视频不卡| 成人欧美大片| 男女床上黄色一级片免费看| 国产成人av激情在线播放| 1024香蕉在线观看| 最新美女视频免费是黄的| 亚洲国产看品久久| 欧美日韩亚洲国产一区二区在线观看| 欧美av亚洲av综合av国产av| 首页视频小说图片口味搜索| 亚洲国产精品999在线| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区视频了| 国产成人影院久久av| 欧美不卡视频在线免费观看 | 国产精品98久久久久久宅男小说| 两个人视频免费观看高清| 中文资源天堂在线| 嫩草影院精品99| 国产精品野战在线观看| 久久久国产精品麻豆| 日韩大尺度精品在线看网址| www.999成人在线观看| 正在播放国产对白刺激| 婷婷六月久久综合丁香| 99热这里只有精品一区 | 日韩精品免费视频一区二区三区| 99久久精品国产亚洲精品| 999精品在线视频| 十八禁人妻一区二区| 欧美不卡视频在线免费观看 | 精品国产超薄肉色丝袜足j| 中亚洲国语对白在线视频| 久久久久国产精品人妻aⅴ院| 美女黄网站色视频| 真人一进一出gif抽搐免费| 露出奶头的视频| 欧美色视频一区免费| 真人做人爱边吃奶动态| 身体一侧抽搐| videosex国产| 三级毛片av免费| 老司机深夜福利视频在线观看| 国产精品久久久人人做人人爽| 午夜福利视频1000在线观看| 91在线观看av| 国产精品久久久久久精品电影| 麻豆成人午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 久久香蕉激情| 岛国在线观看网站| 淫秽高清视频在线观看| 99久久综合精品五月天人人| 999精品在线视频| а√天堂www在线а√下载| 一级毛片女人18水好多| 99国产极品粉嫩在线观看| 人成视频在线观看免费观看| 日韩大尺度精品在线看网址| 美女午夜性视频免费| 亚洲五月婷婷丁香| 亚洲中文字幕日韩| 给我免费播放毛片高清在线观看| 亚洲精品中文字幕在线视频| 给我免费播放毛片高清在线观看| 免费无遮挡裸体视频| 亚洲精华国产精华精| www日本在线高清视频| 久久久久久久精品吃奶| 欧美精品亚洲一区二区| 男人的好看免费观看在线视频 | 亚洲人成电影免费在线| 又粗又爽又猛毛片免费看| 日韩有码中文字幕| 国产精品久久久久久亚洲av鲁大| 久久精品国产99精品国产亚洲性色| 久久国产乱子伦精品免费另类| 国产成人欧美在线观看| 久久欧美精品欧美久久欧美| 精品免费久久久久久久清纯| 一边摸一边抽搐一进一小说| 亚洲男人天堂网一区| 岛国在线免费视频观看| 欧美久久黑人一区二区| 亚洲 国产 在线| 国产精品久久视频播放| 欧美高清成人免费视频www| 久久久久久久久久黄片| 欧美激情久久久久久爽电影| av福利片在线观看| 老熟妇乱子伦视频在线观看| 在线永久观看黄色视频| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产看品久久| 亚洲五月婷婷丁香| 国产精品一区二区精品视频观看| 成在线人永久免费视频| 高潮久久久久久久久久久不卡| 一a级毛片在线观看| bbb黄色大片| 亚洲人成网站在线播放欧美日韩| 最近视频中文字幕2019在线8| 国产又黄又爽又无遮挡在线| 色噜噜av男人的天堂激情| 亚洲成a人片在线一区二区| 少妇粗大呻吟视频| 男插女下体视频免费在线播放| 欧美性猛交╳xxx乱大交人| 制服诱惑二区| 欧美不卡视频在线免费观看 | 日日摸夜夜添夜夜添小说| 亚洲精品色激情综合| 精品欧美一区二区三区在线| 12—13女人毛片做爰片一| 欧美成人一区二区免费高清观看 | 午夜激情av网站| 欧美日韩乱码在线| 成人欧美大片| 免费av毛片视频| 亚洲欧美精品综合久久99| 国产三级中文精品| 久久午夜亚洲精品久久| 少妇粗大呻吟视频| 久久久久久亚洲精品国产蜜桃av| av在线播放免费不卡| 亚洲第一电影网av| 国产高清视频在线观看网站| 丁香六月欧美| 叶爱在线成人免费视频播放| 曰老女人黄片| 一边摸一边抽搐一进一小说| 悠悠久久av| 国产视频内射| 少妇人妻一区二区三区视频| 99热这里只有精品一区 | 午夜福利免费观看在线| 国产熟女xx| 午夜亚洲福利在线播放| 亚洲精品一卡2卡三卡4卡5卡| 久久久久九九精品影院| 日韩精品青青久久久久久| 亚洲第一欧美日韩一区二区三区| 午夜a级毛片| 99热这里只有是精品50| 一个人免费在线观看的高清视频| 一本大道久久a久久精品| 色播亚洲综合网| 全区人妻精品视频| 女警被强在线播放| 国产成人系列免费观看| 观看免费一级毛片| videosex国产| 成年版毛片免费区| 制服诱惑二区| 日韩欧美精品v在线| 人妻久久中文字幕网| 黄片小视频在线播放| 一级毛片精品| 精品免费久久久久久久清纯| 一边摸一边做爽爽视频免费| 日本一区二区免费在线视频| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久久久久久久| 日本精品一区二区三区蜜桃| 欧美高清成人免费视频www| 亚洲一码二码三码区别大吗| 亚洲电影在线观看av| x7x7x7水蜜桃| 欧美在线黄色| 国内精品久久久久精免费| 亚洲av电影在线进入| 国产伦一二天堂av在线观看| 最近最新免费中文字幕在线| 99久久久亚洲精品蜜臀av| 又爽又黄无遮挡网站| 在线永久观看黄色视频| 丁香六月欧美| 69av精品久久久久久| 亚洲电影在线观看av| 无人区码免费观看不卡| 午夜福利欧美成人| 亚洲在线自拍视频| 久久中文字幕人妻熟女| 又黄又粗又硬又大视频| a级毛片a级免费在线| 午夜福利在线观看吧| 免费观看精品视频网站| 色综合亚洲欧美另类图片| 少妇的丰满在线观看| 日本免费一区二区三区高清不卡| 99在线视频只有这里精品首页| 一区二区三区激情视频| 在线观看午夜福利视频| 白带黄色成豆腐渣| 两个人视频免费观看高清| 无限看片的www在线观看| 俄罗斯特黄特色一大片| 久久 成人 亚洲| 长腿黑丝高跟| 50天的宝宝边吃奶边哭怎么回事| 一进一出好大好爽视频| 一本精品99久久精品77| 午夜免费观看网址| 久久人妻av系列| 丰满人妻一区二区三区视频av | 国产精品免费视频内射| 欧美成人性av电影在线观看| 亚洲av片天天在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美+亚洲+日韩+国产| 熟女少妇亚洲综合色aaa.| 波多野结衣高清作品| 亚洲av成人一区二区三| www.www免费av| 午夜福利高清视频| √禁漫天堂资源中文www| 黄色视频不卡| 午夜激情av网站| 亚洲人成电影免费在线| 一区福利在线观看| 国产精品野战在线观看| 日本 欧美在线| 天天躁夜夜躁狠狠躁躁| 美女扒开内裤让男人捅视频| 亚洲精品粉嫩美女一区| 午夜激情av网站| 在线观看免费午夜福利视频| 国产精品av视频在线免费观看| 成人欧美大片| 日本 欧美在线| 亚洲专区中文字幕在线| www.熟女人妻精品国产| 亚洲av电影在线进入| 欧美性长视频在线观看| 亚洲黑人精品在线| 午夜激情av网站| 精品国产亚洲在线| 禁无遮挡网站| 国产精品久久久久久久电影 | 国语自产精品视频在线第100页| 老鸭窝网址在线观看| 激情在线观看视频在线高清| 中国美女看黄片| 亚洲熟女毛片儿| 国模一区二区三区四区视频 | 精品一区二区三区视频在线观看免费| 国产熟女午夜一区二区三区| 老汉色∧v一级毛片| 午夜老司机福利片| 欧美日本亚洲视频在线播放| 国产精品亚洲av一区麻豆| 无遮挡黄片免费观看| 国产激情久久老熟女| 床上黄色一级片| 欧美不卡视频在线免费观看 | 国产一区在线观看成人免费| 精品国内亚洲2022精品成人| 99久久无色码亚洲精品果冻| 999久久久国产精品视频| 在线播放国产精品三级| or卡值多少钱| 欧美日韩精品网址| 美女 人体艺术 gogo| 国产精品久久久人人做人人爽| 99久久99久久久精品蜜桃| 午夜免费激情av| 亚洲美女黄片视频| 国产99久久九九免费精品| 国产成人精品无人区| 91字幕亚洲| 亚洲欧洲精品一区二区精品久久久| 变态另类丝袜制服| 久久国产精品影院| 久久这里只有精品19| 宅男免费午夜| 国产伦一二天堂av在线观看| 欧美乱妇无乱码| 亚洲片人在线观看| 日日摸夜夜添夜夜添小说| 精品久久久久久久末码| 国产亚洲精品第一综合不卡| 日本一本二区三区精品| 日日干狠狠操夜夜爽| 久久婷婷成人综合色麻豆| 亚洲欧美精品综合久久99| 91成年电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久亚洲精品国产蜜桃av| 88av欧美| 国产精品一区二区三区四区免费观看 | 美女免费视频网站| 日日夜夜操网爽| 麻豆av在线久日| 亚洲美女黄片视频| 午夜精品一区二区三区免费看|