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

    Quantitative retrieval of soil salt content based on measured spectral data

    2016-12-19 06:34:52HanChenDuanTaoWangXianXueCuiHuaHuangChangZhenYan
    Sciences in Cold and Arid Regions 2016年6期

    HanChen Duan,Tao Wang,Xian Xue,CuiHua Huang,ChangZhen Yan

    Key Laboratory of Desert and Desertification,Northwest Institute of Eco-Environment and Resources,Chinese Academy of Sciences,Lanzhou,Gansu 730000,China

    Quantitative retrieval of soil salt content based on measured spectral data

    HanChen Duan,Tao Wang,Xian Xue*,CuiHua Huang,ChangZhen Yan

    Key Laboratory of Desert and Desertification,Northwest Institute of Eco-Environment and Resources,Chinese Academy of Sciences,Lanzhou,Gansu 730000,China

    Choosing the Minqin Oasis,located downstream of the Shiyang River in Northwest China,as the study area,we used field-measured hyperspectral data and laboratory-measured soil salt content data to analyze the characteristics of saline soil spectral reflectance and its transformation in the area,and elucidated the relations between the soil spectral reflectance,reflectance transformation,and soil salt content.In addition,we screened sensitive wavebands.Then,a multiple linear regression model was established to predict the soil salt content based on the measured spectral data,and the accuracy of the model was verified using field-measured salinity data.The results showed that the overall shapes of the spectral curves of soils with different degrees of salinity were consistent,and the reflectance in visible and near-infrared bands for salinized soil was higher than that for non-salinized soil.After differential transformation,the correlation coefficient between the spectral reflectance and soil salt content was obviously improved.The first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance produced the highest accuracy and stability in the bands at 462 and 636 nm; the determination coefficient was 0.603,and the root mean square error was 5.407.Thus,the proposed model provides a good reference for the quantitative extraction and monitoring of regional soil salinization.

    spectral reflectance; soil salt content; salinization; multiple linear regression; Minqin Oasis

    1 Introduction

    Soil salinization,which is one of the main types of land desertification and degradation,has led to global resource,environmental,and ecological problems (Li and Wu,2002; Hussain et al.,2006).Soil salinization can weaken and even destroy land productivity,reduce crop yields,and threaten the ecological environment and biosphere.It has become a major barrier to agricultural development in arid and semi-arid regions and threatens the security and stability of oases (Dehaan and Taylor,2002; Metternicht and Zinck,2003).Therefore,understanding the nature,scope,and extent of saline soils can assist in the effective governance of soil salinization,facilitating local agricultural production and sustainable development of regional economies.

    At present,the application of remote sensing technology in mapping and monitoring soil salinization has shown great accuracy and efficiency,and multispectral remote sensing technology is widely used to monitor soil salinization (Khan and Sato,2001; Khan et al.,2005).However,because of a relatively low spectral resolution,it is difficult to distin-guish the mixed spectral characteristics of complex surfaces,leading to loss of some spectral details and making it impossible to accurately estimate the soil salt content (Pu et al.,2012).In contrast to multispectral technology,hyperspectral technology can generate continuous spectral information of ground surface features on the nanoscale and distinguish subtle spectral differences among objects.Thus,hyperspectral technology makes possible the identification of surface features based on diagnostic spectral absorption characteristics,remote sensing quantitative analysis,and chemical composition research (Zhang et al.,2012b; Xu et al.,2014).Given the advantages of hyperspectral technology,in recent years it has been widely used for the quantitative extraction of soil salinization characteristics.Based on measured spectral data and EO-1 hyperspectral images,Weng et al.(2010) accurately estimated the soil salt content using the salinization index.By analyzing the spectral characteristics of several typical objects in the delta oasis of the Weigan and Kuqa rivers of the Tarim Basin,Zhang et al.(2010) established a ground spectral database.Ding et al.(2012) used saline soil,vegetation,and saline soil spectral reflectance to establish a hyperspectral model for monitoring salinization based on a measured comprehensive spectral index.Zhao et al.(2014) analyzed the spectral features of saline soil by combining spectral reflectance data with hyperspectral satellite images and extracted salinization information using the spectral angle classification method.

    In this study,using the Minqin Oasis in Northwest China as the study area,the spectral reflectance of soils with different degrees of salinity was obtained by field measurements,and the relation between the reflectance and soil salt content was analyzed.Finally,sensitive wavebands were screened.Using multiple linear regression analysis,we established a model to quantitatively estimate soil salinity based on measured data.This model facilitates reliable and efficient monitoring and evaluation of soil salinization on a large scale.

    2 Study area

    The Minqin Oasis is located in the lower reaches of the Shiyang River in Gansu Province at latitudes of 38°05′N-39°26′N and longitudes of 101°59'E-104°12'E (Figure 1).The Minqin Oasis is bordered on the east by the Tengger Desert and on the northwest by the Badain Jaran Desert,which is surrounded by the Gobi Desert on three sides.The total area of Minqin County is 1.6×104km2,wherein the deserts,eroded mountains,and salinas account for 91%,and the oasis accounts for only 9% (Gao et al.,2006; Pang et al.,2014).The study area is arid and has a typical temperate continental climate; its annual average temperature is 7.6 °C,annual average rainfall is 110 mm,and annual average evaporation is 2,644 mm (about 24 times the amount of precipitation) (Chen et al.,2014).

    The soil types in the study area include sandy,gray-brown desert,meadow,and meadow-boggy soils.The natural vegetation,which is divided into desert and meadow types,is of the typical extremely arid desert steppe type (Yue et al.,2011; Zhang et al.,2014).Minqin is a typical oasis agricultural county; however,the influence of drought,shortage of water resources,and strong evaporation result in salt accumulation in low-lying land.This soil salinization seriously restricts the sustainable development of oasis agriculture and the social economy.

    Figure 1 Study area and distribution of sampling sites

    3 Methodology

    3.1Collection of soil samples and measurement of salt content

    In this study,field sampling was conducted on October 2-10,2014 (after the autumn harvest and before autumn plowing and winter irrigation).We collected 42 soil sampling units from the surface (0~10 cm) of typical soil salinization regions.The dimensions of a sampling unit were 100m×100m,and three samples (forming a triangular shape) were obtained from each unit,resulting in a total of 126 soil samples.The coordinates of the center position of each sampling unit were simultaneously recorded via GPS.After natural drying and grinding,the soil samples were filtered using a 1-mm aperture sieve and mixed with deionized water (5:1 by mass) to extract the salt content.The mean salt content of the three soil samples in each sampling unit was considered as the measured value of a sampling unit.Errors and abnormal values caused by human operation and instruments were eliminated.In total,we obtained 39 data points from the sampling units (117 soil samples),29 of which were used to construct the model,and the remaining 10 were used to test the model's accuracy.

    According to the classification standard of soil salinization of Minqin County (Chen et al.,2013),and the distribution of the field survey sampling points,we classified the soil salinization degree by choosing the soil salt content as the main factor (Table 1).It can be seen from Table 1 that the degree of soil salinization was divided into five types in the study area:non-salinization,slight salinization,moderate salinization,severe salinization,and saline.The sampling points were mainly distributed in abandoned land,unused land within the oasis,and the border region between the oasis and desert,all of which have a certain representativeness.

    Table 1 Classification of soil salinization

    3.2Field spectral data acquisition

    The spectral data were acquired using an ASD Field Spec.Hand Held 2 spectrometer (Analytical Spectral Devices,Boulder,CO,USA) with a spectral range of 325~1,075 nm,spectral resolution of <3 nm,sampling interval of 1.4 nm,and wavelength accuracy of ±1 nm.

    The field spectral measurements and soil sampling were simultaneously performed on the same ground.Field spectra were obtained on clear,cloudless days with low wind speeds and sufficient,stable sunlight intensity.To reduce the influence of changes in the solar altitude angle on the spectral measurement results,the measurement period was 10:00-14:00 local time,because during this period the field of view was directly illuminated by the sun and the atmospheric transmittance was the highest.The view angle of the probe (25°) was perpendicular to the surface (20 cm) during the measurements.Before measurement,the influence of dark current was removed,and the instrument was optimized based on a white reference panel.The spectral data were collected 10 different times at each measurement point.After removing any abnormal curves,the arithmetic average value was determined as the spectral reflectance value of a sampling point.The average value of the spectral reflectance of the three sampling points in each unit was determined as the spectral reflectance of a sampling unit to reduce the difference among sampling points in the same area.

    3.3Spectral data processing

    Because of the difference in energy responses among spectrometer bands,a spectral curve has "burr" noise; thus,optimizing the spectral curve is necessary.In this study,the nine-point weighted moving average method was used to smooth spectral reflectance data (Peng et al.,2013).

    Previous studies (e.g.,He et al.,2006; Ma,2014) demonstrated that the appropriate mathematical transformation of a spectral reflectance can effectively limit the influence of low-frequency noise on a target spectrum,eliminate the correlation among bands,and enhance the spectral difference among features.It is easy to identify a high-correlation band with soil salinity,which can improve the inversion accuracy.Inaddition to the analysis of the original spectral reflectance (R),various mathematical transformations have been employed:root mean squarelogarithm (lgR),reciprocal (1/R),logarithmic reciprocal (1/lgR),reciprocal logarithm (lg(1/R)),first-order differential (R′),first-order differential of root mean squarefirst-order differential of logarithm ((lgR)′),first-order differential of reciprocal ((1/R)′),first-order differential of logarithmic reciprocal ((1/lgR)′),and first-order differential of reciprocal logarithm ((lg(1/R))′).

    Correlation analysis helps identify bands that are sensitive to soil properties and establish a regression model for soil properties (Lei et al.,2014).Based on the measured spectral reflectance and its transformation forms as feature vectors,the correlation coefficient ribetween each band's feature vector and soil salt content was obtained,and the characteristic bands of the salt content in the soil were screened.A prediction model was then established using a curve regression method based on the following formula:

    where riis the correlation coefficient between spectral reflectance,its transformation forms,and soil salt content; n is the total number of sampling units; j is the band number; Rijis the reflectance value of a spectrum and its transformation forms in the j band of the ithsampling unit;is the average value of the spectral reflectance and transformation forms of n units in the j band; SCiis the soil salt content in the ithunit; andis the average value of the soil salt content in n units.

    3.4Model construction and testing

    In this study,we randomly selected 29 soil samples from a total of 39 samples to establish a curve regression model; the other 10 samples were used for model validation.The total root mean square error (RMSE) was used to evaluate the accuracy of the estimation model; a relatively smaller RMSE value indicates a relatively higher model precision and more reliable simulation results.RMSE is defined as follows:

    where n is the total number of sampling units,is the soil salt content predicted by the model,and yiis the measured soil salt content.

    4 Results

    4.1Analysis of spectral characteristics of saline soil

    Figure 2 shows the reflectance curves of soils with different degrees of salinity.Although some differences were observed,the overall shapes of the reflectance curves of these soils were consistent.The spectral reflectance increased with wavelength; the reflectance increased sharply in the wavelength range of 400~600 nm and gently in the range of 600~900 nm.With the exception of a weak absorption band in the vicinity of 751 nm,no obvious peak and trough were observed in the spectral curves.In the wavelength range of 350~1,000 nm,the spectral reflectance did not increase with the soil salt content.In contrast,as the degree of soil salinization increased,the reflectivity gradually decreased,and the spectral reflectance of the slightly salinized soil was the highest.This result was attributed to the serious effect of soil moisture on the spectral curves.In the study area,the soil water content is high,and salt accumulates because of strong evaporation,resulting in soil salinization.

    Figure 2 Spectral curves of soil reflectance with different degrees of salinity

    The correlation analysis in this study indicated a significantly positive correlation between soil salinity and soil water content (Figure 3); the correlation coefficient was 0.852,which was significant at the 0.01 level.These results indicate that the higher the soil moisture,the greater the soil salinity,and vice versa.Previous studies showed that the soil in the study area contains a large amount of high-moisture-absorbing salt (MgCl2) (Qi et al.,2010; Zhang and Qi,2013).Thus,under natural conditions,the water in the surrounding environment can be easily absorbed,increasing soil moisture and resulting in a negative correlation between the spectral reflectance and soil saltcontent.In general,the spectral reflectance of a salinized soil in the visible and near-infrared bands is higher than that of a non-salinized soil.

    Figure 3 Relation between soil salinity and soil water

    4.2Screening of sensitive bands of salt content in soil

    There are significant differences between the spectral reflectance and soil salinity,and the sensitive bands of soil salinity can be screened from measured soil spectral curves.Thus,it is necessary to quantitatively analyze the relation between the soil spectral reflectance and soil salinity.To identify the sensitive bands of soil salinity,we calculated the correlation coefficient between the soil spectral reflectance,its transformation,and soil salt content (Figure 4).Since the reflectance signal had more noise at wavelengths below 400 nm and above 900 nm,we selected a wavelength range of 400~900 nm to analyze the correlation between the soil spectral reflectance and soil salinity.As shown in Figure 4,the correlation coefficient changed with wavelength,and both positive and negative correlations were observed between spectral reflectance,its transformation,and soil salt content.Soil salt content with spectral reflectance (R),root mean squareand logarithm (lgR) was negatively correlated,whereas that with reciprocal (1/R),logarithmic reciprocal (1/lgR),and reciprocal logarithm (lg(1/R)) was positively correlated.Soil salt content with first-order differential (R′),first-order differential of root mean squarefirst-order differential of logarithm ((lgR)′),first-order differential of reciprocal ((1/R)′),first-order differential of logarithmic reciprocal ((1/lgR)′),and first-order differential of reciprocal logarithm ((lg(1/R))′) were both positively and negatively correlated.After differential transformation,the correlation coefficient showed a fluctuating trend around the 0 axis,and the fluctuation range of the correlation coefficient was significantly greater than that before differential transformation.

    From the correlation coefficient peak values and their corresponding bands between the spectral reflectance,reflectance transformation,and the soil salt content,we can see that (Table 2) the correlation coefficients of the peak band between the spectral reflectance (R),root mean squarelogarithm (lgR),first-order differential of reciprocal ((1/R)′),first-order differential of reciprocal logarithm ((lg(1/R))′),and the soil salt content was negative; the correlation coefficients of the first three variables was significantly correlated at the 0.05 level,and that of the last two variables were significantly correlated at the 0.01 level.

    However,not all transformation forms were able to improve the sensitivity of the spectral reflectance to the soil salt content.On the contrary,some transformations reduced the correlation between the reflectance and the soil salt content,such as the root mean squarelogarithm (lgR),reciprocal (1/R),and reciprocal logarithm (lg(1/R)); the correlation coefficient between them and the soil salt content was smaller than that between the original spectral reflectance and the soil salt content.The correlation coefficients of the peak band between the reciprocal (1/R),logarithmic reciprocal (1/lgR),reciprocal logarithm (lg(1/R)),first-order differential of logarithm ((lgR)′),first-order differential of logarithmic reciprocal ((1/lgR)′) and soil salt content were positive; the correlation coefficients of the peak band of reciprocal (1/R) and reciprocal logarithm (lg(1/R)) show significantly correlated at the 0.05 level,whereas the rest were significantly correlated at the 0.01 level.Among them,the soil salt content and the first-order differential of logarithm ((lgR)′) had the most close correlations; the corresponding correlation coefficients were 0.453 and 0.516,respectively.The correlation coefficients of the peak band between the first-order differential (R′),first-order differential of root mean squareand the soil salt content were both positive and negative,but all were significantly correlated at the 0.01 level.These results indicate a linear relation between soil salt content and the spectral indices of the corresponding points.

    After the differential transformation of the spectral reflectance,the subtle information contained in the original spectral data was amplified,significantly enhancing the correlation coefficient between the salt content and spectral reflectance.On this basis,we determined that the corresponding bands of the correlation coefficient's peak value of the spectral different transformations were the sensitive bands to the soil salt content (Table 2).Based on the correlation analysis between the spectral reflectance and the soil salt content,the sensitive bands of the soil salt content were located between the visible and near-infrared wave bands.This conclusion is consistent with the results of previous research (Gu et al.,2011; Zhang et al.,2011).

    Figure 4 Correlation coefficients of spectral reflectance and its different transforms with soil salinity

    Table 2 Peak values and the corresponding bands of correlation coefficients between spectral different transformations and soil salinity

    4.3Construction of the model to predict soil salt content and verification of its accuracy

    According to the correlation coefficients between spectral reflectance,its different transformations,and soil salinity,the corresponding bands of the correlation coefficient peak values were selected.We applied the multiple linear regression method to construct an inversion model of soil salinity based on different independent variables (Table 3).The accuracy of the model was verified using the field-measured soil salinity data,and the optimal inversion model of soil salinity was determined based on the measured spectral data.

    By comparison and analysis the soil salt content inversion models (Table 3),the sensitive bands of thesoil salt content of R,,lgR,1/R,1/lgR,and lg(1/R) were mostly located around 897 and 898 nm,and other forms after differential transform were mostly located around 462 and 636 nm.The first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance had the highest determination coefficient (R2= 0.603),the lowest RMSE (5.407),and the most reliable stability.Followed by the first-order differential of the logarithm model,its determination coefficient R2and RMSE were 0.578 and 5.605,respectively.The logarithmic reciprocal model had the lowest determination coefficient (R2= 0.395),the highest RMSE (6.406),and the worst reliability.Therefore,the first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance ((lg(1/R))′) was used to retrieve and extract the soil salinity in the study area.

    We considered the 10 soil sample units that were not involved in constructing the model as true values in order to verify the stability and prediction ability of the model.The determination coefficient (R2) was used to test the stability of the model,and the total RMSE was used to test the model's predictive capability; the predicted values of soil salinity were obtained from the first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance,and then the regressive relation was established between the predicted and measured values.The scatter diagram of the predicted and measured values is shown in Figure 5.It can be seen that the determination coefficient (R2) of the test samples was 0.806,and the total RMSE was 3.662,satisfying the accuracy requirement of this study.Therefore,compared with the other models in Table 3,the first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance had a better accuracy and higher stability; thus,it can better achieve the quantitative extraction and retrieval of the soil salt content in the study area.

    Figure 5 Scatterplot of the measured and estimated values of the curve regression model

    Table 3 Comparison of multiple linear regression models with different dependent variables

    5 Discussion

    Previous research has shown that the correlation between soil salt content and soil spectral reflectance under artificial light is significantly higher than that under natural light.The spectral reflectance increases with soil salinity (Rao et al.,1995; Farifteh et al.,2008; Zhao et al.,2012).The quantitative inversion model established in this paper is mainly based on the field-measured hyperspectral data,the purpose of which is to provide a scientific basis for quantitative monitoring of soil salinization by using hyperspectralremote sensing.However,we found that,due to the influence of the sampling number,local climate,soil type,soil moisture content,soil texture,and many other factors,there are some differences between the soil spectral curves (Elvidge et al.,1990; Zhang et al.,2012a).Although the accuracy of the soil salt content estimation model based on the reflectance of artificial light is higher than that based on natural light (Zhang et al.,2012c),it is not the actual reflection of the land surface.It is difficult to establish the corresponding relationship between the estimation model and the remote sensing images.So,for monitoring the soil salinization quantitatively,spectrum analysis for modeling under natural light is more applicable to connect the satellite images,which helps achieve rapid and accurate monitoring on the regional scale.Therefore,compared with the reflectance of artificial light in the laboratory,the soil salt content estimation model based on the spectral reflectance of natural light is more practical.

    During our research,we found that there were some wavelengths that had significant correlations with the soil salt content.However,due to the existence of collinearity effects between bands,the high-autocorrelation bands were removed to avoid the impact of collinearity effects on the results of the study.Only two sensitive bands that had higher correlation and lower autocorrelation with the soil salt content were selected to establish the soil salinity estimation model.Moreover,the linear regression analysis method used in this study has some limitations:it requires the independent and dependent variables to have a good linear relation.If there is no linear relation between them,this method will not be able to accurately estimate the salinity of the soil.Therefore,we need to develop a method based on a combination of linear and nonlinear relations to model soil salinity and to make the quantitative monitoring methods more accurate.Also,the number of samples we used for model establishment and verification was too limited,which may have induced adverse influences to the results.The relatively small sampling number may have lessened the simulation accuracy and created an error in the result,thereby restricting the model extent to only a large regional scale.In future research we will increase the number of soil samples as much as possible,using different soil types and textures,so as to improve the accuracy and applicability of the model.

    6 Conclusions

    In this paper,based on field-measured soil spectral reflectance and their 11 transformation forms,the coupling relations between the soil salt content and the spectral reflectance,as well as its transformation forms,were analyzed.Based on the results,we constructed a quantitative model to predict the soil salt content and verified the model using measured soil salt content data.The main conclusions are as follows:

    1) The spectral curves of soils with different degrees of salinity were consistent as a whole,and the spectral reflectance increased with wavelength.Spectral reflectance increased sharply in the wavelength range of 400~600 nm and gently in the range between 600 and 900 nm.Due to the effects of soil moisture,the soil spectral reflectance decreased with increasing degree of soil salinization; thus,spectral reflectance and soil salt content were significantly negatively correlated.By analyzing the spectral characteristics of the saline soil,we found that the spectral reflectance in the visible and near-infrared regions was higher for salinized soil than that for non-salinized soil.

    2) Both positive and negative correlations were found between the soil spectral reflectance,its transformation forms,and the soil salt content.After the differential transformation of spectral reflectance,the correlation coefficient between soil salt content and spectral reflectance was significantly enhanced.The sensitive bands of the soil salt content of the different spectral transformation forms were determined based on the corresponding bands of the correlation coefficient peak values,which were all located between the visible and near-infrared wavelengths.

    3) By comparing and analyzing the inversion models of the soil salt content based on different independent variables,we found that the optimal sensitive bands of the first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance were the bands at 462 and 636 nm; the correlation coefficient was maximized for these bands,and the RMSE was minimized.The inversion model was the optimum model for the estimation of soil salt content.We used the measured soil salinity data to test the stability and predictive ability of the model; the results showed that the determination coefficient (R2) of the test samples was 0.806,and the RMSE was 3.662.Thus,compared with other regression models,the first-order differential transformation model based on the logarithm of the reciprocal of saline soil spectral reflectance had a better accuracy and higher stability,and this model can better predict the soil salt content in the study area.

    Acknowledgments:

    This study was financially supported by the National Natural Science Foundation of China (No.41401109),the Foundation for Excellent Youth Scholars of CAREERI,CAS (No.Y551D21001),and the Open Fund Project of the Key Laboratory of Desert and Desertification,CAS (No.Y452J71001).The authors would like to thank Enago (www.enago.cn) for the English language review.

    References:

    Chen LJ,Feng Q,Cheng AF,2013.Spatial distribution of soil water and salt contents and reasons of saline soils' development in the Minqin Oasis.Journal of Arid Land Resources and Environment,27(11):99-105.DOI:10.3969/ j.issn.1003-7578.2013.11.017.(in Chinese)

    Chen XS,Gao BB,Wang XJ,et al.,2014.Changes and status of desertification land in Minqin,Gansu,China.Journal of Desert Research,34(4):970-974.DOI:10.7522/ j.issn.1000-694X.2014.00077.(in Chinese)

    Dehaan RL,Taylor GR,2002.Field-derived spectra of salinized soils and vegetation as indicators of irrigation induced soil salinization.Remote Sensing of Environment,80(3):406-417.DOI:10.1016/S0034-4257(01)00321-2.

    Ding JL,Wu MC,Liu HX,et al.,2012.Study on the soil salinization monitoring based on synthetical hyperspectral index.Spectroscopy and Spectral Analysis,32(7):1918-1922.DOI:10.3964/j.issn.1000-0593(2012)07-1918-05.(in Chinese)

    Elvidge CD,1990.Visible and near-infrared reflectance characteristics of dry plant materials.International Journal of Remote Sensing,11(10):1775-1795.DOI:10.1080/01431169008955129.

    Farifteh J,Meer F,Meijde M,et al.,2008.Spectral characteristics of salt-affected soils:A laboratory experiment.Geoderma,145(3-4):196-206.DOI:10.1016/j.geoderma.2008.03.011.

    Gao ZH,Li ZY,Wei HD,et al.,2006.Quantitative monitoring of vegetation cover change by using remotely sensed data over Minqin Oasis,Gansu.Geographical Research,25(4):587-595.DOI:10.3321/j.issn:1000-0585.2006.04.004.(in Chinese)

    Gu XH,Wang K,Pan YC,et al.,2011.Remote sensing quantitative retrieval of farmland organic matter based on HJ1A-HIS image.Geography and Geo-Information Science,27(6):69-73.(in Chinese)

    He T,Wang J,Lin ZJ,et al.,2006.Spectral features of soil organic matter.Geomatics and Information Science of Wuhan University,31(11):975-979.(in Chinese)

    Hussain N,Al-Rawahy SA,Rabee J,et al.,2006.Causes,origin,genesis and extent of soil salinity in the Sultanate of Oman.Pakistan Journal of Agricultural Sciences,43(1/2):1-6.

    Khan NM,Rastoskuev VV,Sato Y,2005.Assessment of hydrosaline land degradation by using a simple approach of remote sensing indicators.Agricultural Water Management,77(1-3):96-109.DOI:10.1016/j.agwat.2004.09.038.

    Khan NM,Sato Y,2001.Monitoring hydro-salinity status and its impact in irrigated semi-arid areas using IRS-1B LISS-II data.Asian Journal of Geoinform,1(3):63-73.

    Lei L,Tiyip T,Ding JL,et al.,2014.Study on the soil salinization monitoring based on measured hyperspectral and HSI data.Spectroscopy and Spectral Analysis,34(7):1948-1953.DOI:10.3964/j.issn.1000-0593(2014)07-1948-06.(in Chinese)

    Li FQ,Wu ZR,2002.A study on soil salinization early-warning in semi-arid area:As a sample in west Jilin province.Bulletin of Soil and Water Conservation,22(1):57-59.(in Chinese)

    Ma C,2014.Research on soil salinization using remote sensing of HJ-1A hyper spectral images.Journal of Arid Land Resources and Environment,28(2):180-184.(in Chinese)

    Metternicht GI,Zinck JA,2003.Remote sensing of soil salinity:potentials and constraints.Remote Sensing of Environment,85(1):1-20.DOI:10.1016/S0034-4257(02)00188-8.

    Pang GJ,Wang T,Sun JH,et al.,2014.Quantitative analysis of soil salinity content with hyperspectra data in Minqin,Gansu,China.Journal of Desert Research,34(4):1073-1079.DOI:10.7522/j.issn.1000-694X.2014.00045.(in Chinese)

    Peng J,Xiang HY,Wang JQ,et al.,2013.Inversion models of soil water content using hyperspectral measurements in fields of the arid region farmland.Agricultural Research in the Arid Areas,31(2):241-246.(in Chinese)

    Pu Z,Yu RD,Yin CY,et al.,2012.Optimal hyperspectral indices for soil salt content estimation on typical saline soil in arid areas.Bulletin of Soil and Water Conservation,32(6):129-133.(in Chinese)

    Qi HP,Weng YL,Zhao DY,et al.,2010.A study of salinity characteristics and spectral characteristics of salt affected soil in Caka-Gonghe basin.Remote Sensing for Land & Resources,(Supp.1):4-8.(in Chinese).

    Rao BRM,Sankar TR,Dwivedi RS,et al.,1995.Spectral behavior of salt-affected soils.International Journal of Remote Sensing,16(12):2125-2136.DOI:org/10.1080/01431169508954546.

    Weng YL,Gong P,Zhu ZL,2010.A spectral index for estimating soil salinity in the yellow river delta region of China using EO-1 Hyperion data.Pedosphere,20(3):378-388.DOI:10.1016/S1002-0160(10)60027-6.

    Xu YJ,Ma HC,Peng SY,2014.Study on identification of altered rock in hyperspectral imagery using spectrum of field object.Ore Geology Reviews,56:584-595.DOI:org/10.1016/j.oregeorev.2013.07.004.

    Yue DX,Du J,Gong J,et al.,2011.Dynamic analysis of farmland ecosystem service value and multiple regression analysis of the influence factors in Minqin Oasis.Acta Ecological Sinica,31(9):2567-2575.(in Chinese)

    Zhang F,Tiyip T,Ding JL,et al.,2010.The reflectance spectroscopy of geo-objects in Xinjiang on delta oasis of Weigan and Kuqa Rivers.Journal of Infrared and Millimeter Wave,29(3):190-195.(in Chinese)

    Zhang F,Tiyip T,Ding JL,et al.,2012a.Spectral reflectance properties of major objects in desert oasis:a case study of the Weigan-Kuqa River Delta Oasis in Xinjiang,China.Environmental Monitoring and Assessment,184(8):5105-5119.DOI:10.1007/s10661-011-2326-x.

    Zhang F,Tiyip T,Ding JL,et al.,2012b.The study on the reflectance spectral characteristic of salt-affected soil in typical oasis of the middle reaches of Tarim River.Progress in Geography,31(7):921-932.(in Chinese)

    Zhang F,Xiong FM,Luan FM,et al.,2011.Characteristics of field measured spectral response to alkalinization soil.Journal of Infrared and Millimeter Wave,30(1):55-60.(in Chinese)

    Zhang F,Xiong HG,Ding JL,et al.,2012c.Characteristics of laboratory-field measured spectra responding to alkalinized soil and conversion.Transactions of the Chinese Society of Agricultural Engineering,28(5):101-107.DOI:10.3969/j.issn.1002-6819.2012.05.017.(in Chinese)

    Zhang JM,Qi WW,2013.Soil salt components and their spectral characteristics in Minqin Oasis of Northwest China.Chinese Journal of Ecology,32(10):2620-2626.(in Chinese)

    Zhang M,Gao TP,Zhu YR,et al.,2014.Spatial distribution characteristics of soil salinity in habitat soil of Halogeton glomeratus in Minqin Oasis.Soils,46(4):756-760.DOI:10.13758/j.cnki.tr.2014.04.028.(in Chinese)

    Zhao ZL,Tiyip T,Sun Q,et al.,2014.Soil spectrum characteristics and information extraction of salinization:A case study in Weigan-Kuqa Oasis in Xinjiang.Progress in Geography,33(2):280-288.DOI:10.11820/dlkxjz.2014.02.014.(in Chinese)

    Zhao ZL,Tiyip T,Zhang F,et al.,2012.Spectral characteristics of soil salt content in typical oasis of Tarim River's middle reaches.Journal of Natural Disasters,21(5):72-78.(in Chinese)

    Duan HC,Wang T,Xue X,et al.,2016.Quantitative retrieval of soil salt content based on measured spectral data.Sciences in Cold and Arid Regions,8(6):0507-0515.

    10.3724/SP.J.1226.2016.00507.

    *Correspondence to:Ph.D.,Xian Xue,Professor of Northwest Institute of Eco-Environment and Resources,Chinese Academy of Sciences.No.320,West Donggang Road,Lanzhou,Gansu 730000,China.E-mail:xianxue@lzb.ac.cn

    April 28,2016Accepted:July 29,2016

    各种免费的搞黄视频| 五月玫瑰六月丁香| 69人妻影院| 国产美女午夜福利| 成人欧美大片| 亚洲av中文av极速乱| 亚洲精品色激情综合| 国产成人免费无遮挡视频| 少妇人妻精品综合一区二区| 夜夜爽夜夜爽视频| av网站免费在线观看视频| 久久久成人免费电影| 亚洲内射少妇av| 丰满少妇做爰视频| av线在线观看网站| 国产熟女欧美一区二区| 国产成人91sexporn| 少妇猛男粗大的猛烈进出视频 | 美女脱内裤让男人舔精品视频| 五月伊人婷婷丁香| 色视频在线一区二区三区| 亚洲成人久久爱视频| 在线观看一区二区三区激情| 欧美最新免费一区二区三区| 日本wwww免费看| 免费观看性生交大片5| 欧美丝袜亚洲另类| 亚洲欧美成人精品一区二区| 亚洲精品,欧美精品| 亚洲精品国产成人久久av| 国产老妇伦熟女老妇高清| 国产 一区 欧美 日韩| 夫妻性生交免费视频一级片| 我的女老师完整版在线观看| 五月天丁香电影| 真实男女啪啪啪动态图| 永久网站在线| 欧美国产精品一级二级三级 | 国产免费一区二区三区四区乱码| 涩涩av久久男人的天堂| 身体一侧抽搐| 国产黄a三级三级三级人| 一本久久精品| 天堂俺去俺来也www色官网| 久久久久久久国产电影| 精品一区二区三卡| 久久久久久九九精品二区国产| 国产成年人精品一区二区| 春色校园在线视频观看| 国产一区有黄有色的免费视频| 国产老妇女一区| 草草在线视频免费看| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久久免| 免费av不卡在线播放| 国产精品秋霞免费鲁丝片| 欧美zozozo另类| 欧美日韩亚洲高清精品| av天堂中文字幕网| av又黄又爽大尺度在线免费看| 国产亚洲精品久久久com| 99热这里只有是精品50| 男女边摸边吃奶| 欧美另类一区| 国产v大片淫在线免费观看| 亚洲在久久综合| 国产真实伦视频高清在线观看| 成人漫画全彩无遮挡| 免费在线观看成人毛片| 免费看光身美女| 欧美日韩视频精品一区| 亚洲精品影视一区二区三区av| 高清毛片免费看| 亚洲成色77777| 免费av毛片视频| 一区二区三区乱码不卡18| 老师上课跳d突然被开到最大视频| 在线a可以看的网站| 国产精品久久久久久精品电影| 一区二区三区乱码不卡18| 真实男女啪啪啪动态图| 婷婷色麻豆天堂久久| 久久影院123| 久久99热6这里只有精品| av在线亚洲专区| 亚洲成人精品中文字幕电影| 国产av国产精品国产| 内地一区二区视频在线| 大话2 男鬼变身卡| 国产亚洲最大av| 三级国产精品片| 欧美+日韩+精品| 国产一区二区在线观看日韩| 久久午夜福利片| 涩涩av久久男人的天堂| 男人舔奶头视频| 99九九线精品视频在线观看视频| av线在线观看网站| 国产亚洲5aaaaa淫片| 91久久精品国产一区二区三区| 久久人人爽人人片av| 欧美日本视频| 亚洲伊人久久精品综合| 亚洲精品国产av成人精品| 精品久久久噜噜| 国产永久视频网站| 久久精品久久精品一区二区三区| 免费黄网站久久成人精品| 欧美最新免费一区二区三区| 欧美少妇被猛烈插入视频| 久久精品国产亚洲网站| 成人亚洲精品一区在线观看 | 国产一区二区亚洲精品在线观看| 晚上一个人看的免费电影| 99久久九九国产精品国产免费| 美女国产视频在线观看| 精品一区二区三区视频在线| 高清视频免费观看一区二区| www.色视频.com| 禁无遮挡网站| 久热这里只有精品99| 欧美亚洲 丝袜 人妻 在线| 观看美女的网站| 国产欧美日韩一区二区三区在线 | 国产精品99久久99久久久不卡 | 能在线免费看毛片的网站| 在线观看一区二区三区激情| 欧美日韩视频精品一区| 成人无遮挡网站| 国产探花在线观看一区二区| av国产免费在线观看| 一个人观看的视频www高清免费观看| 亚洲精品中文字幕在线视频 | 国产精品久久久久久精品古装| 网址你懂的国产日韩在线| 色吧在线观看| 精品人妻熟女av久视频| 国产欧美亚洲国产| 2018国产大陆天天弄谢| 午夜精品国产一区二区电影 | 日日摸夜夜添夜夜爱| 九九爱精品视频在线观看| 国产精品一区二区在线观看99| 一级毛片 在线播放| 婷婷色综合www| 成人午夜精彩视频在线观看| 又大又黄又爽视频免费| 最近2019中文字幕mv第一页| 亚洲在久久综合| 街头女战士在线观看网站| 大片电影免费在线观看免费| 亚洲色图综合在线观看| 久久久久九九精品影院| 欧美xxxx性猛交bbbb| 又黄又爽又刺激的免费视频.| 亚洲,欧美,日韩| 亚洲精品乱码久久久久久按摩| a级毛片免费高清观看在线播放| 超碰97精品在线观看| 老司机影院毛片| 久久人人爽人人爽人人片va| 亚洲av不卡在线观看| 99视频精品全部免费 在线| 高清av免费在线| 成年免费大片在线观看| 99久久九九国产精品国产免费| 亚洲精品自拍成人| 精品久久久久久久久av| 日韩亚洲欧美综合| 赤兔流量卡办理| 日本午夜av视频| 一级毛片我不卡| 久久久欧美国产精品| 看非洲黑人一级黄片| 99久久精品国产国产毛片| 国产午夜精品一二区理论片| 99久久中文字幕三级久久日本| 狂野欧美白嫩少妇大欣赏| 九草在线视频观看| 成人特级av手机在线观看| 国产精品蜜桃在线观看| 91久久精品电影网| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品| 久久99精品国语久久久| 亚洲国产最新在线播放| 麻豆成人av视频| 日日撸夜夜添| 又爽又黄无遮挡网站| 你懂的网址亚洲精品在线观看| 男女国产视频网站| 观看免费一级毛片| 哪个播放器可以免费观看大片| 性色avwww在线观看| 国产精品国产av在线观看| 26uuu在线亚洲综合色| 色婷婷久久久亚洲欧美| h日本视频在线播放| 九色成人免费人妻av| 如何舔出高潮| 久久精品人妻少妇| 亚洲va在线va天堂va国产| 小蜜桃在线观看免费完整版高清| 免费不卡的大黄色大毛片视频在线观看| 99热全是精品| 天堂中文最新版在线下载 | 青春草视频在线免费观看| 91aial.com中文字幕在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲av中文字字幕乱码综合| 一级毛片 在线播放| a级毛色黄片| 在线天堂最新版资源| 亚洲av电影在线观看一区二区三区 | 又爽又黄a免费视频| 久久久午夜欧美精品| 国产免费一区二区三区四区乱码| 一个人看视频在线观看www免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av.在线天堂| 久久影院123| tube8黄色片| 亚洲成人精品中文字幕电影| 乱系列少妇在线播放| 国产精品.久久久| 国产亚洲午夜精品一区二区久久 | 一边亲一边摸免费视频| 久久精品国产鲁丝片午夜精品| 国产爱豆传媒在线观看| 久久ye,这里只有精品| 青春草国产在线视频| 看黄色毛片网站| 国产av码专区亚洲av| 亚洲美女视频黄频| 国产淫片久久久久久久久| tube8黄色片| 边亲边吃奶的免费视频| 国产男女内射视频| 丝袜脚勾引网站| 秋霞伦理黄片| 我的老师免费观看完整版| 麻豆国产97在线/欧美| 日本欧美国产在线视频| 亚洲精品国产av蜜桃| 国产精品.久久久| 成人国产av品久久久| 精品国产露脸久久av麻豆| 3wmmmm亚洲av在线观看| 精品少妇久久久久久888优播| 亚洲精品中文字幕在线视频 | 听说在线观看完整版免费高清| 最近的中文字幕免费完整| 日韩在线高清观看一区二区三区| freevideosex欧美| 深爱激情五月婷婷| 亚洲精品中文字幕在线视频 | 偷拍熟女少妇极品色| 中文字幕亚洲精品专区| 99久久九九国产精品国产免费| 亚洲国产精品专区欧美| 性插视频无遮挡在线免费观看| 色吧在线观看| 亚洲av一区综合| 国产v大片淫在线免费观看| 99久久九九国产精品国产免费| 成人国产麻豆网| 好男人视频免费观看在线| videos熟女内射| 国产精品一区二区性色av| 一级毛片电影观看| 欧美极品一区二区三区四区| 欧美区成人在线视频| 五月玫瑰六月丁香| 国产精品99久久99久久久不卡 | 久久精品国产a三级三级三级| 精品熟女少妇av免费看| 国产av不卡久久| 亚洲av不卡在线观看| 蜜臀久久99精品久久宅男| 久久国内精品自在自线图片| 久久久久久久久久人人人人人人| 久久6这里有精品| 日韩成人伦理影院| 男女国产视频网站| 国产在线男女| 日韩av免费高清视频| videos熟女内射| 99热国产这里只有精品6| 青春草国产在线视频| 欧美精品一区二区大全| 一本一本综合久久| 美女脱内裤让男人舔精品视频| 啦啦啦中文免费视频观看日本| kizo精华| 舔av片在线| 亚洲丝袜综合中文字幕| 少妇的逼好多水| 97在线人人人人妻| 日韩不卡一区二区三区视频在线| 中文字幕免费在线视频6| 亚洲精品久久久久久婷婷小说| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| 18禁在线无遮挡免费观看视频| 日韩国内少妇激情av| 一级毛片电影观看| 国产午夜精品一二区理论片| 啦啦啦在线观看免费高清www| 国产在线一区二区三区精| 婷婷色综合大香蕉| 成人亚洲精品av一区二区| 黄片无遮挡物在线观看| 青春草国产在线视频| 日韩人妻高清精品专区| 亚洲人与动物交配视频| 男女啪啪激烈高潮av片| 成年人午夜在线观看视频| 99精国产麻豆久久婷婷| 欧美成人a在线观看| 看免费成人av毛片| 丝袜美腿在线中文| 精品酒店卫生间| 亚洲精品久久午夜乱码| 免费看光身美女| 久久综合国产亚洲精品| 欧美+日韩+精品| 青春草亚洲视频在线观看| 国产淫语在线视频| 亚洲国产欧美人成| 日韩电影二区| 麻豆国产97在线/欧美| 黄色欧美视频在线观看| 日韩av免费高清视频| 成人美女网站在线观看视频| 国产精品人妻久久久影院| 国产大屁股一区二区在线视频| 久久精品久久久久久噜噜老黄| 嘟嘟电影网在线观看| 只有这里有精品99| 日韩一区二区三区影片| 亚洲成人精品中文字幕电影| 蜜桃亚洲精品一区二区三区| 亚洲成人精品中文字幕电影| 国产一级毛片在线| 黄片wwwwww| 国产精品蜜桃在线观看| 视频区图区小说| 国产欧美日韩一区二区三区在线 | 亚洲av成人精品一区久久| 免费av毛片视频| 久久久久久久久久久丰满| 成人国产麻豆网| 精品99又大又爽又粗少妇毛片| av国产免费在线观看| 亚洲精品中文字幕在线视频 | 久久久久国产网址| 在线观看免费高清a一片| 亚洲国产日韩一区二区| 国产成人a∨麻豆精品| 国产精品.久久久| 亚洲aⅴ乱码一区二区在线播放| 十八禁网站网址无遮挡 | 我要看日韩黄色一级片| 特级一级黄色大片| 国产高清有码在线观看视频| 伊人久久精品亚洲午夜| 男男h啪啪无遮挡| 成人亚洲精品一区在线观看 | 日韩 亚洲 欧美在线| 边亲边吃奶的免费视频| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 色视频www国产| 国产精品福利在线免费观看| 一级av片app| 精品人妻视频免费看| 久久99热6这里只有精品| 国产精品偷伦视频观看了| 日韩一区二区三区影片| 97人妻精品一区二区三区麻豆| 亚洲精品乱码久久久v下载方式| av国产免费在线观看| 欧美日韩亚洲高清精品| 日韩制服骚丝袜av| 在线观看一区二区三区激情| 男女边吃奶边做爰视频| 国产伦理片在线播放av一区| av在线播放精品| 岛国毛片在线播放| 午夜老司机福利剧场| 插逼视频在线观看| 色婷婷久久久亚洲欧美| 男女边摸边吃奶| 久久久久性生活片| 3wmmmm亚洲av在线观看| 精品亚洲乱码少妇综合久久| 黑人高潮一二区| 亚洲在线观看片| 亚洲av福利一区| 男人添女人高潮全过程视频| 久久ye,这里只有精品| 国产伦理片在线播放av一区| 亚洲,一卡二卡三卡| 国产白丝娇喘喷水9色精品| 六月丁香七月| 大话2 男鬼变身卡| 韩国av在线不卡| 亚洲自拍偷在线| 亚洲国产成人一精品久久久| 久久99热这里只有精品18| 国产乱人偷精品视频| 高清日韩中文字幕在线| 一级毛片aaaaaa免费看小| 激情五月婷婷亚洲| 精品少妇黑人巨大在线播放| 少妇猛男粗大的猛烈进出视频 | 狂野欧美激情性bbbbbb| 国产精品福利在线免费观看| 熟妇人妻不卡中文字幕| 99热网站在线观看| 亚洲三级黄色毛片| 在线观看国产h片| 欧美+日韩+精品| 日日摸夜夜添夜夜爱| 夫妻性生交免费视频一级片| 中文乱码字字幕精品一区二区三区| 国产永久视频网站| 在线免费十八禁| 你懂的网址亚洲精品在线观看| 久久6这里有精品| 波多野结衣巨乳人妻| 七月丁香在线播放| 国产片特级美女逼逼视频| 夫妻性生交免费视频一级片| 午夜免费鲁丝| 久久久精品免费免费高清| 亚洲精品色激情综合| 99久久人妻综合| 少妇人妻一区二区三区视频| 免费黄色在线免费观看| 制服丝袜香蕉在线| 一区二区三区免费毛片| 国产欧美亚洲国产| 超碰av人人做人人爽久久| 午夜免费男女啪啪视频观看| 中文精品一卡2卡3卡4更新| 国产日韩欧美亚洲二区| 日韩不卡一区二区三区视频在线| 国产精品蜜桃在线观看| 97超视频在线观看视频| 又爽又黄a免费视频| 一级片'在线观看视频| 久久久精品免费免费高清| 美女xxoo啪啪120秒动态图| 日韩,欧美,国产一区二区三区| av国产久精品久网站免费入址| 日本黄大片高清| 亚洲欧美日韩另类电影网站 | 日韩大片免费观看网站| 国产伦精品一区二区三区四那| 免费大片18禁| 欧美日韩综合久久久久久| 亚洲国产色片| 国产一区二区三区av在线| 别揉我奶头 嗯啊视频| 国产成年人精品一区二区| 看黄色毛片网站| 黄色怎么调成土黄色| av天堂中文字幕网| 国产美女午夜福利| 国产一区二区在线观看日韩| 国产成人精品一,二区| 日韩强制内射视频| 国产精品.久久久| 99久国产av精品国产电影| 日日啪夜夜撸| 直男gayav资源| 2018国产大陆天天弄谢| 一个人看视频在线观看www免费| 亚洲国产最新在线播放| 中文字幕av成人在线电影| 国产精品久久久久久久久免| 亚洲精品乱久久久久久| 久久亚洲国产成人精品v| 国产精品久久久久久av不卡| 欧美精品一区二区大全| 香蕉精品网在线| 久久久久久伊人网av| av专区在线播放| 黄色怎么调成土黄色| 精品国产露脸久久av麻豆| 最近的中文字幕免费完整| 欧美日韩精品成人综合77777| 免费黄频网站在线观看国产| 久久影院123| 人人妻人人澡人人爽人人夜夜| 国产片特级美女逼逼视频| 中文欧美无线码| 内地一区二区视频在线| 一级a做视频免费观看| 日本熟妇午夜| 婷婷色综合www| 汤姆久久久久久久影院中文字幕| 别揉我奶头 嗯啊视频| 亚洲内射少妇av| 久久鲁丝午夜福利片| 99久国产av精品国产电影| 日本与韩国留学比较| av线在线观看网站| av一本久久久久| 国产精品.久久久| av女优亚洲男人天堂| 哪个播放器可以免费观看大片| 日本一本二区三区精品| 青青草视频在线视频观看| 欧美区成人在线视频| 亚洲精品成人av观看孕妇| 精品国产乱码久久久久久小说| 天堂网av新在线| 视频中文字幕在线观看| 国产精品久久久久久久电影| 日本黄色片子视频| 丰满人妻一区二区三区视频av| 97超碰精品成人国产| 97热精品久久久久久| 亚洲一区二区三区欧美精品 | 国产亚洲午夜精品一区二区久久 | 国产淫片久久久久久久久| 97人妻精品一区二区三区麻豆| 国产成人a区在线观看| 亚洲av一区综合| 精品熟女少妇av免费看| 国产精品伦人一区二区| 国产精品人妻久久久影院| 99re6热这里在线精品视频| 国产亚洲av嫩草精品影院| 只有这里有精品99| 国产高清三级在线| 嫩草影院新地址| 我的老师免费观看完整版| 18禁动态无遮挡网站| 大片电影免费在线观看免费| 国产毛片a区久久久久| 国产老妇伦熟女老妇高清| 3wmmmm亚洲av在线观看| 可以在线观看毛片的网站| 色播亚洲综合网| 亚洲激情五月婷婷啪啪| 97在线人人人人妻| 青青草视频在线视频观看| 亚洲欧美一区二区三区黑人 | 久久精品久久久久久久性| 97超视频在线观看视频| 精品少妇黑人巨大在线播放| 国产片特级美女逼逼视频| 欧美bdsm另类| 80岁老熟妇乱子伦牲交| 男男h啪啪无遮挡| 亚洲av中文av极速乱| 亚洲天堂国产精品一区在线| 青春草亚洲视频在线观看| 肉色欧美久久久久久久蜜桃 | 我的女老师完整版在线观看| 久久久久久久久久成人| 成年人午夜在线观看视频| 又大又黄又爽视频免费| 国产精品国产av在线观看| 插阴视频在线观看视频| 久久国产乱子免费精品| 午夜免费观看性视频| 少妇猛男粗大的猛烈进出视频 | 亚洲图色成人| 成人漫画全彩无遮挡| 成年女人看的毛片在线观看| 特大巨黑吊av在线直播| 国产成人福利小说| 国产精品一区二区在线观看99| 三级国产精品片| av又黄又爽大尺度在线免费看| 欧美激情久久久久久爽电影| 51国产日韩欧美| 亚洲综合色惰| 久久精品国产亚洲av涩爱| 下体分泌物呈黄色| 美女国产视频在线观看| 精品少妇黑人巨大在线播放| 乱系列少妇在线播放| 少妇人妻久久综合中文| 99re6热这里在线精品视频| 亚洲欧美中文字幕日韩二区| 亚洲国产av新网站| 麻豆成人av视频| 丝袜脚勾引网站| 最近的中文字幕免费完整| 99久久中文字幕三级久久日本| 一个人看视频在线观看www免费| 国产精品一区www在线观看| 少妇人妻久久综合中文| 内射极品少妇av片p| 水蜜桃什么品种好| 最近最新中文字幕免费大全7| 久久鲁丝午夜福利片| 午夜激情福利司机影院| 国产成年人精品一区二区| 国产白丝娇喘喷水9色精品| 国产精品一二三区在线看| 国产欧美另类精品又又久久亚洲欧美| 精品酒店卫生间| 韩国高清视频一区二区三区| 高清av免费在线| 欧美区成人在线视频| 国产精品蜜桃在线观看| 国产高清国产精品国产三级 | 国产一区二区亚洲精品在线观看| 国产熟女欧美一区二区|