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

    Study on Maize-water Model for Supplemental Irrigation in Loess Plateau

    2015-02-05 12:08:07XingLIMangmangGOU
    Agricultural Science & Technology 2015年5期

    Xing LI,Mangmang GOU

    1.Inner Mongolia Engineering Research Center for Water-saving Agriculture,Inner Mongolia Normal University,Hohhot 010022, China

    2.Institute of Water Conservancy Science,Inner Mongolia,Hohhot 010020,China

    Study on Maize-water Model for Supplemental Irrigation in Loess Plateau

    Xing LI1*,Mangmang GOU2

    1.Inner Mongolia Engineering Research Center for Water-saving Agriculture,Inner Mongolia Normal University,Hohhot 010022, China

    2.Institute of Water Conservancy Science,Inner Mongolia,Hohhot 010020,China

    The Loess Plateau has a typical semi-arid climate,and the area suffers from very harsh ecological environment,severe soil erosion and water runoff,and uneven distributed precipitation.Due to the relatively low holding capacity,current rainwater-collecting and conservation facilities can only supplement a maximum of 18 mm of water for crop production in each irrigation.In this study,mathematical models were constructed to identify the water requirement critical period of maize crop by evaluating response of each individual developmental stage to supplemental irrigation with harvested rainwater.In the transformed Jensen model,ETmin/Eta was used as the index of relative evapotranspiration.The use of relative yield and relative crop evapotranspiration was able to eliminate influences from unintended environmental factors.A BP neural network crop-water model for extreme water deficit condition was constructed using the index of relative evapotranspiration as the input and the index of relative yield as the output after iterative training and adjustment of weight values.Comparison of measured maize yields to those predicted by the two models confirmed that the BP neural network crop-water model is more accurate than the transformed Jensen model in predicting the sensitivity index to waterdeficit at various growth stages and maize yield when provided with supplemental irrigation with harvested rainwater.

    BP neural network;Model of crop response to water;The transformed Jensen model;Supplemental irrigation with harvested rainwater

    A s water resources crisis worsens on the Loess Plateau in China,rain-fed agriculture has become a major activity in agricultural development in this arid and semi-arid region.Rainwater harvesting projects collect rainwater across multiple spatial and temporal scales,thus providing an effective approach to using natural precipitation[1-7].But due to small water storage capacity,current rain water collection facilities can only harvest a limited amount of water.The relationship between plants and water is very complex[8-10].Current models for water use as a function of crop production are constructed using the addition or multiplication mathematical principles[11-16].None of those models can accurately simulate the relationship between water and plants.For instance,when Jensen model was used to fit the experimental data under water deficit conditions,the resultant non-linear model predicted a negative sensitivity index at certain growth stages.To resolve this issue,Guo et al.proposed that the relationship between accumulation function of sensitivity index and the plant growth periods should fit to a sine curve[17-18].The logistic growth curve was also used to fit the accumulation function of sensitive index.Both models were able to produce reasonable results.However, the curve model for the accumulation function of sensitivity index is a nonphysical model and it is liable to cause some errors,although the scale of the errors and the subsequent impact on simulation accuracy is not yet known[19].Jensen model is the most commonly used simulation method in under-sufficient irrigation experiments where the well-irrigated treatments are used as the control samples[20-21].Various correction coefficients from a combination of different factors are also integrated into Jensen model[22],and the construction of a transformed Jensen model that can predict the crop-water relationship for the system of supplemental irrigation with harvested rainwater is critically important for the Loess Plateau region.

    Based on the record at the Rainwater Harvest Experimental Station in Junggar,Qiandu town,Inner Magnolia,this area has very limited water resources.Under this circumstance,this study designed the experiment where supplementing rainwater at each of the maize growth stages was the treated whereas no-irrigation across the whole growth period was the control. The relationship between final yield and the water status at each growth stage was treated as a complex nonlinear approximation problem.A BP neural network maize-water model under extremely water deficit condition was constructed.This model was used to predict the optimal timing of supplemental irrigation with rainwater to produce the highest yield.This study has provided theoretically bases for developing a crop production system using harvested rain water resources to produce high yield and high water use efficiency in arid and semi-arid regions.

    Materials and Methods

    Field experiment

    A field experiment was performed from 2003 to 2007 on a Water Collection and Supplemental Irrigation Experimental Station in northern Junggar,Ordos,Inner Mongolia.This typical semi-arid is on the Loess Plateau facing serious water erosion problem. Over 30%of the land is drenches and gullies connected with each other.Elevation ranges from 1030 m in the deepest gullies to 1 320 m at the top of mountains.The scarce rainfall is concentrated in a very short time period.

    In the spring and fall,this area experiences large daily temperature fluctuation.The multiple-year average temperature is 7.2℃,the annual sunshine duration hours is 3 117.3 h,the accumulated temperature at above≥10℃is 3 500℃,the annual wind speed is 2.3 m/s,and the average non-frost season is 145 d.The multiple year average precipitation is 392.1 mm.Rainfall is distributed unevenly across spaces and time,it occurs mostly in summer and fall.Plants receive 60%of the life-time water from July to September.The multiple year evaporation is 2 093 mm,which is five times of total precipitation.

    In the field experiment,each treatment group consisted of three replicate plots of 34 m2each,and 2 m wide barrier crop zone was set up between treatments.Maize cultivar‘Zhedan No.7’was sown at a rate of 37.5 kg/hm2and spacing of 30 cm x 50 cm.Crop density in the field was approximately 5.25×10 plants/hm.

    Supplemental irrigation was applied using a ground level fixed plastic hose pressure-controlled drip irrigation system controlled through a water meter.Based on the amount of water in the reservoirs and the capacity of collecting rain water and the efficiency of the irrigation system,supplemental irrigation of 18 mm at each grow stage was scheduled.The experimental design in 2006 was described in Table 1.

    The transformed Jensen model under extreme water deficit condition

    A mathematical model of watercrop production function simulates the relationship between crop yield and periodic water consumption,and it also provides the basic knowledge for designing water saving irrigation and management plans.Jensen model is the most frequently used periodic water-production function model.In this model crop evapotranspiration is the variable to be used to predict effects of different water deficit levels at each growth stage on final crop yield.The equation for Jensen model is written as follows:

    where n—the number of crop growth stages;

    i—code of crop developmental stage;

    λi—the sensitivity indexes at growth stage i;

    ETai—the actual evapotranspiration at growth stage i (mm);

    ETmi—the evapotranspiration potential at growth stage i corresponding to Ym(mm);

    Ya—the actual yield(kg);

    Ym—the highest yield potential under well-irrigated condition(kg).

    In previous studies,sufficient irrigation was used as the control for the non-sufficient irrigation treatments, results from those experiment represent the effect of water deficit at each growth stage on final yield.Due to the natural environment and climate conditions in the Loess Plateau and the small amount of precipitation collected and would be available for supplemental irrigation,in this study non-irrigation across the whole growth period was used as the control for treatments of irrigation at each growth stage.Therefore,results from this study would represent sensitivity of yield to supplemental water at each growth stage.The equation of Jensen model was modified as follows:

    where:ETmin—the amount of actual evapotranspiration corresponding to Ymin(mm);

    Ymin—the lowest crop yield(kg); ETmin/ETa—index of relative evapotranspiration;

    Ymin/Ya—index of relative yield.

    BP neural network model

    Artificial neural network(ANN) started in the 1940s,and it gained momentum in the 1980s.During the last 20 years,a second wave of interest in exploring this modeling system has formed.ANN is a powerful method that simulates human brain basic structure and function to solve non-linear high dimension and high-order data information questions.The data processing method is characterized by strong self-organizing,selfadaptable learning and memory association ability.Over 80%of ANN models use the BP reverse path forwarding transmission network modeling scheme.The BP modeling can process input and output identification of nonlinear dynamical systems,and the weight values are adjusted by using the back propagation learning algorithm for training the network[23].

    Table 1 The supplemental irrigation schedule of maize crop in 2006

    Determination of the number of layers in the network model

    BP network model has multiple layers which fit to the nonlinear differentiable function by performing weight training, therefore complex non-linear function approximation problems can be resolved using the BP network containing single hidden layer.In this study,a three layer network structure was used.

    The activation function for each layer of the network in the BP model

    The BP network’s capacity in nonlinear approximation is embedded in the sigmoid activation function, therefore hidden layers normally use such function.The sigmoid or linear transfer function can be used at the output layer.The sigmoid transfer function reaches a non-linear approximation faster,but the output is restricted in the range of 0-1 or-1-1.After comprehensive analysis,the tan-sigmoid transfer function was used at the hidden layers using the following program:f(x)=(1-e-2x)/(1+e-2x).The output layer transfer function was purelin [f(x)=x],which can take on any values.

    Determination of the node number per layer in the BP network model

    Selection of units in the hidden layers is the key for determining the node number per layer in the BP network. The accuracy of the model can be comprised if too few nodes are selected.But using more nodes in the network will increase the amount of training time.Normally the trial-and-error method is used to determine the optimum number of hidden nodes.Selection of input and output layer depends directly on the input and output variables for solving the real problem, therefore these parameters are easy to determine.In this study,according to physiological properties of maize crop,the whole growth cycle was divided into five stages:Emergence-jointing stage,Jointing-tassel,Heading-early milk,Milk-nearly mature,and Matureharvest.

    In the model,the input layer had five 5 neurons each corresponding to the relative evapotranspiration at one of the five growth stages.One neuron was selected for maize yield at the output layer.The learning rate was derived by changing the training parameters for all the items,it was adjusted with the error-correction rule during the training process.The initial learning rate had very minimal effect on the network training output.After several iterations,the learning rate reached 0.05 when the fitting error was 0.0001 at eight nodes.The final network structure was 5-8-1.

    Results and Analysis

    Field experiment

    The collected data include soil water content,the amount of supplemental irrigation,and precipitation. The evapotranspiration level at each growth stage was calculated using the following equation:

    where

    ETa—the measured evapotranspiration(mm);

    Woand Wt—water content of soil in the top 1 cm layer at the begging and the end of each treatment period, respectively(mm);

    Po—the effective precipitation during each growth stage(mm);

    K—the amount of underground water fed to plants during each treatment period.In the study area,underground water level was below 50 m and no underground water was available to the crop,therefore K value was 0(mm);

    M—the amount of irrigation for each time period(mm);

    I—the amount of precipitation infiltrated into deep soil layers.Water content was measured to be 0,therefore I was 0(mm);

    U,R—the amount of lateral water feed and seepage.The experimental site was located on the land,therefore the U and R value were both set at 0(mm).

    Results from the 2006 field experiment and the dataset derived using equation(3)were presented in Table 2.

    Training the BP neural network model

    In the BP neural network cropwater model,the index of relative evapotranspiration was the input and the corresponding yield indices from different treatments at the same growth stage were the output values. This would have eliminated the unintended influences from climate factors. In the network training process if the weight values are set too big,with the increased weight gain from all the neurons the input of the activation functions would fall in the saturation region of the sigmoid transfer function,then it will lead to small derivative values thus disabling the adjustment procedure of the corresponding network weights. For input data value in the 0-1 range, a preventive strategy is to place the input data with bigger weight values inthe sigmoid activation function region with large gradient at the hidden layer.

    Table 2 Maize yield under supplemental irrigation of collected rainwater at different growth stages

    The BP neural network training program was written using the MatLab language.Experimental data in Table 2 were converted to the indices of relative evapotranspiration and relative yield using the matrix computation.The link weight values of the layers after training are as follows:

    The trained BP neural network model was then used to simulate maize yield from supplemental irrigation treatments(Table 3).The relative evapotranspiration index(ETmin/ETa)was used as the input to predict the relative yield index(Ymin/Ya) from each irrigation treatment,which was then used as the index of yield sensitivity to supplemental irrigation at each growth stage.

    According to the constructed BP crop-water model for supplemental irrigation with harvested rainwater,decreasing water supply led to reduced yields but higher relative yield indices,which can be explained by the general crop-water relationship(Table 3).Based on the sensitivity index,supplemental irrigation at the Heading-early milk growth stage had the greatest effect on increasing yield,and the smallest effect was found at the Emergence-jointing stage.On the other hand,the simulation results also confirmed the beneficial effect from water withholding(drought acclimation)at seedling stage on crop yield.Results in Table 3also showed that sensitivity of maize crop to water deficit status decreased in the following sequence:Heading-Early milk>Jointing-tassel>Milk-mature>Emergence-jointing stages.

    Table 3The indices of relative yield(Ymin/Ya)simulated using the BP Model

    Table 4 The indices of relative yield(Ymin/Ya)of maize in 2007

    Table 5 The indices of relative yields(Ymin/Ya)of measured maize yield and simulated data

    Table 6 The water deficit sensitivity index simulated using the Jensen model

    Furthermore,the accuracy of the constructed BP neural network models was validated by comparing the predicted value to the actual maize yield. To simplify the process,the treatments in 2007 were drip irrigation at growth stages 1,3 and 4.The actual yield index at the corresponding stages was ETmin/ETa=0.887 3,ETmin/ETa= 0.523 9,and ETmin/ETa=0.691 7,and the predicted value for the respective counterpart was ETmin/ETa=0.9, ETmin/ETa=0.5,and ETmin/ETa=0.7 (Table 4).Error between the predicted and actual yields was less than 3% (Table 5),which confirmed that the predicted yield is very close to the actual measurement.Therefore,the BP neural network crop-water model for supplemental irrigation with harvested rainwater can be used to simulate the relationship between supplemental irrigation at various growth stage of maize and the final yield under extreme water deficit condition.

    The transformed Jensen model

    To test the applicability of the transformed Jensen model,a multivariate regression program was written in MatLab language to obtain the sensitivity indices at all the growth stages(Table 6).Five schemes were designed to predict maize yield at different levels and timing of supplemental irrigation.

    The sensitivity index is 0.027 1 at Emergence-jointing,0.269 0 at Jointing-tassel,0.337 3 at Heading-early milk,and 0.255 9 at Milk-nearly mature stage.Then five supplemental irrigation schemes consisting of different amounts of water and at different crop growth stages were designed to predict yields under each of the scenarios(Table 7).The measured maize yield in 2007 and the simulated yield using the transformed Jensen models are presented in Table 8.The relative error between the two sets of yield data was below 5%which is higher than the 3%error level using the BP neural network model.In general,the predicted maize yields using the two models were very close to each otheras well as to the actual yield,and the BP neural network model is slightly better than the transformed Jensen model.

    Table 7 The predicted relative yield index using the transformed Jensen-model(Ymin/Ya)

    Table 8 The measured and simulated relative yield index(Ymin/Ya)of maize

    According to the sensitivity index from the transformed Jensen Model, heading to milking stage is the most sensitive to water deficit(Table 6). The second most critical stage is the elongation-heading stage.These results agree to the simulated water deficit sensitivity index using the BP crop-water model for supplemental irrigation with harvested rainwater under extremely water deficit condition. Such a high level of correlation confirmed that the BP neural network crop-water model can be used to accurately simulate maize yield under various water deficit conditions.

    Conclusions

    (1)Supplemental irrigation at the flag leaf to grain filling stages has the most pronounced effect on improving maize yield.Moderate water constraint at seedling emergence-stem elongation stages will not cause obvious yield reduction.

    (2)The BP neural network cropwater model for supplemental irrigation with harvested rainwater was able to simulate the generic relationship between yield and water status at various growth stages of maize crop.Both the BP and the transformed Jensen model predicted similar sensitive patterns to water deficit conditions at each of the maize growth stages.

    (3)Compared to the transformed Jensen model,the BP neural network crop-water model was more accurate in simulating maize yield on the Loess Plateau.Thus the BP model system is preferred under such situation.But the experimental data for constructing the model were collected for a short period of time and from a small number of locations,the constructed model should be used only in areas with similar environmental conditions.

    Furthermore,artificial neural network program has its inherent deficiencies,for instance the number of hidden layers and the hidden nodes have to be determined using the trial and error method.Additionally,parameters in the neural network model have not a physical meaning,they cannot represent the connections between component variables in the network.Future research needs to be focused on how to improve the BP neural network crop water-model and the network training methods.

    [1]HERRMANN T,SCHMIDA U.Rainwater utilization in Germany:efficiency,dimensioning,hydraulic and environmental aspects[J].Urban Water,2000,1: 307-316.

    [2]ROCKSTROM J,BARRON J,Fox P. Rainwater management for increased productivity among small-holder farmer in drought prone environments[J]. Physics and Chemistry of the Earth, 2002,27:949-959.

    [3]SAZAKLI E,ALEXOPOULOS A,Leotsinidis M.Rainwater harvesting,quality assessment and utilization in Kefalonia Island,Greece[J].Water Research, 2007,41:2039-2047.

    [4]PANIGRAHI B,PANDA S N,MAL B C. Rainwater conservation and recycling by optimal size on-farm reservoir[J]. Resources,2007,50:459-474.

    [5]KAHINDA JM,LILLIE E S B,TAIGBENU A E,et al.Developing suitability maps for rainwater harvesting in South Africa[J].Physics and Chemistry of the Earth,2008,33:788-799.

    [6]LEE JY,YANG JS,HAN M,et al.Comparison of the microbiological and chemical characterization of harvested rainwater and reservoir water as alternative water resources[J].Science of the Total Environment,2010,408:896-905.

    [7]LI Z,BOYLE F,REYNOLDS A.Rainwater harvesting and grey water treatment systems for domestic application in Ireland[J].Desalination,2010,260: 1-8.

    [8]LUO JJ,YANG FK,GAo SM.Suitable irrigation program for water-harvesting areas in semi-arid loess plateau[J]. Journal Irrigation and Drainage,2003, 22:25-28.

    [9]WANG Y.Effect of the supplemental irrigation of harvested rainwater on dryland winter wheat yield and water use [J].Research of Soil and Water Conservation,2003,10:104-107.

    [10]ZHANG ZS,WEI XH,LI XR,et al.Analysis on investment and benefit of harvested rainwater utilization in the northwest loess Plateau[J].Advances in Water Science,2004,15:813-818.

    [11]LUIS R.SALADO-NAVARRO,THOMAS R.Sinclair.Crop rotations in Argentina:Analysis of water balance and yield using crop models[J].Agricultural Systems,2009,102:11-16.

    [12]ALI REZA SEPASKHAH,SAGHAR FAHANDEZH-SAADI,SHAHROKH ZAND-PARSA.Logistic model application for prediction of maize yield under water and nitrogen management[J]. Agricultural Water Management,2011, 99:51-57.

    [13]PATRICIO GRASSINI,HAISHUN YANG,SUAT IRMAK,JOHN THORBURN,CHARLES BURR,KENNETH G.CASSMAN.High-yield irrigated maize in the Western U.S.Corn Belt: II.Irrigation management and crop water productivity[J].Field Crops Research,2011,120:133-141.

    [14]DOM NGUEZ A,MART NEZ R S, JUAN J A DE,MART NEZ-ROMERO A,TARJUELO J M.Simulation of maize crop behavior under deficit irrigation using MOPECO model in a semi-arid environment[J].Agricultural Water Management,2012,107:42-53.

    [15]M.ABEDINPOUR,A.SARANGI,T.B. S.RAJPUT,MAN SINGH,H. PATHAK,T.AHMAD.Performance evaluation of Aqua Crop model for maize crop in a semi-arid environment [J].Agricultural Water Management,2012,110:55-56.

    [16]MR SALAZAR,JE HOOK,A GARCIA Y GARCIA,JO PAZ,B.CHAVES,G. HOOGENBOOM.Estimating irrigation water use for maize in the Southeastern USA:A modeling approach[J].A-gricultural Water Management,2012. 107:104-111.

    [17]GUO QS,LEI ZD,YANG SX.Senstive Index of Jensen’s Water Production Function Model of Wheat[J].Advances in Water Science,1996,7:20-25.

    [18]WANG YR,LEI ZD,YANG SX.Cumulative function of sensitive index for winter wheat[J].Journal of Hydraulic Engineering,1997,5:28-35.

    [19]JIAO XY,PENG SZ.Reasons and Solutions of Negative Values of Sensitive Indexes of Jensen’s Model[J].Journal of Shen Yang Agriculture University, 2004,35:439-442.

    [20]ZHOU ZW,SHANG SH,LEI ZD. Jensen model and ANN model for water—fertilizer production function of winter wheat[J].Advances in Water Science,2003,14:280-284.

    [21]CUI YL,MA CX,SHEN XZ,et al.Predicting reference evaportranspiration based on artificial neural network with genic arithmetic[J].Advances in Water Science,2005,16:76-81.

    [22]GUO XP,ZHU CL.Model of crop response to water with influence of after—effect of water stress[J].Advances in Water Science,2004,15: 463-466.

    [23]WEI ZM,CHEN YX,SHI HB,et al. Preliminary study on spring wheat response to water with BP neural network method[J].Journal Irrigation and Drainage,2002,21:12-16.

    Responsible editor:Xiaoxue WANG

    Responsible proofreader:Xiaoyan WU

    Supported by Inner Mongolia water conservancy"Twelfth five-year"Major Science and Technology Demonstration Project-scientific Support Project for New Water-saving Irrigation Area of Four ten Million mu in Inner Mongolia in China(20121036)and the National Natural Science Foundation of China(No.51469026,2012MS0621)and the Guided Reward Fund for Scientific and Technological Innovation,Inner Mongolia,China.

    *Corresponding author.E-mail:lixinggmm@163.com

    Received:March 3,2015 Accepted:April 17,2015

    十八禁人妻一区二区| 亚洲自拍偷在线| 日韩欧美国产一区二区入口| 亚洲熟妇中文字幕五十中出| 俺也久久电影网| 中文字幕精品亚洲无线码一区| 两性夫妻黄色片| 又黄又粗又硬又大视频| 91麻豆av在线| 欧美色欧美亚洲另类二区| 日本一区二区免费在线视频| 麻豆国产av国片精品| 久久天躁狠狠躁夜夜2o2o| 日日干狠狠操夜夜爽| 国产精品久久久久久人妻精品电影| 亚洲精华国产精华精| 久久精品成人免费网站| 欧美在线一区亚洲| a级毛片在线看网站| 高清毛片免费观看视频网站| 可以在线观看的亚洲视频| 级片在线观看| 少妇熟女aⅴ在线视频| www.自偷自拍.com| 欧美日韩亚洲综合一区二区三区_| 日本一区二区免费在线视频| 欧美乱码精品一区二区三区| 亚洲 国产 在线| 长腿黑丝高跟| 午夜精品在线福利| 窝窝影院91人妻| 亚洲一区二区三区色噜噜| 香蕉av资源在线| 18美女黄网站色大片免费观看| 亚洲精品国产一区二区精华液| 超碰成人久久| 亚洲自拍偷在线| 久久人妻福利社区极品人妻图片| 亚洲精品国产一区二区精华液| 757午夜福利合集在线观看| 亚洲欧美日韩无卡精品| 亚洲专区国产一区二区| 波多野结衣高清无吗| 国产一区二区激情短视频| 麻豆成人午夜福利视频| 日韩大尺度精品在线看网址| 成人国产一区最新在线观看| 久久人妻福利社区极品人妻图片| 亚洲avbb在线观看| 国产精品乱码一区二三区的特点| 日本 av在线| 久久久久久久精品吃奶| 午夜久久久久精精品| 变态另类成人亚洲欧美熟女| 黄片小视频在线播放| 日韩高清综合在线| 亚洲 欧美一区二区三区| 老汉色∧v一级毛片| 男女床上黄色一级片免费看| 一进一出抽搐动态| 欧美日韩国产亚洲二区| 午夜福利欧美成人| 黄片大片在线免费观看| 丁香六月欧美| 嫩草影视91久久| 久久久久久大精品| 国产单亲对白刺激| 在线观看66精品国产| 丰满的人妻完整版| 麻豆国产av国片精品| 窝窝影院91人妻| 久久精品人妻少妇| 精品国产乱码久久久久久男人| 国产区一区二久久| av免费在线观看网站| 国内揄拍国产精品人妻在线| 国模一区二区三区四区视频 | x7x7x7水蜜桃| 黑人操中国人逼视频| 日韩精品青青久久久久久| www.999成人在线观看| 好男人在线观看高清免费视频| 欧美一级a爱片免费观看看 | 亚洲av电影不卡..在线观看| 欧美成人午夜精品| 亚洲精品国产精品久久久不卡| 免费在线观看黄色视频的| 久久 成人 亚洲| 好男人电影高清在线观看| 90打野战视频偷拍视频| 欧美一级毛片孕妇| 女同久久另类99精品国产91| 亚洲精品粉嫩美女一区| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 亚洲人成77777在线视频| 国产亚洲欧美98| 亚洲va日本ⅴa欧美va伊人久久| 成人三级黄色视频| cao死你这个sao货| 色尼玛亚洲综合影院| 亚洲av成人一区二区三| 一夜夜www| a级毛片在线看网站| 国产成人一区二区三区免费视频网站| 免费在线观看黄色视频的| 欧美中文日本在线观看视频| 亚洲第一电影网av| 国产成人精品久久二区二区免费| 国产精品一区二区三区四区免费观看 | 午夜久久久久精精品| 99久久99久久久精品蜜桃| 三级男女做爰猛烈吃奶摸视频| 男女午夜视频在线观看| 国内毛片毛片毛片毛片毛片| 在线播放国产精品三级| 999久久久国产精品视频| 午夜福利免费观看在线| 熟女少妇亚洲综合色aaa.| 国产黄片美女视频| 亚洲国产欧美人成| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 88av欧美| 午夜老司机福利片| 一本大道久久a久久精品| 12—13女人毛片做爰片一| 老汉色∧v一级毛片| 热99re8久久精品国产| 久99久视频精品免费| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区av网在线观看| 激情在线观看视频在线高清| 日韩精品青青久久久久久| 丰满的人妻完整版| 一夜夜www| 真人一进一出gif抽搐免费| 精品欧美一区二区三区在线| 国产激情偷乱视频一区二区| av免费在线观看网站| 午夜影院日韩av| 99在线人妻在线中文字幕| 亚洲五月天丁香| 99久久精品国产亚洲精品| 脱女人内裤的视频| 久久 成人 亚洲| 很黄的视频免费| 女人被狂操c到高潮| 久久精品亚洲精品国产色婷小说| 18禁黄网站禁片免费观看直播| 国内毛片毛片毛片毛片毛片| 久久草成人影院| 老熟妇乱子伦视频在线观看| 欧美黑人精品巨大| 国产精品久久久久久人妻精品电影| 亚洲成a人片在线一区二区| 免费一级毛片在线播放高清视频| 亚洲精品在线美女| 亚洲熟妇熟女久久| 亚洲欧美一区二区三区黑人| 无人区码免费观看不卡| 亚洲精品av麻豆狂野| 欧美三级亚洲精品| 精品午夜福利视频在线观看一区| 首页视频小说图片口味搜索| 在线十欧美十亚洲十日本专区| 美女大奶头视频| 性色av乱码一区二区三区2| 精品日产1卡2卡| 91九色精品人成在线观看| 亚洲乱码一区二区免费版| 欧美日韩亚洲综合一区二区三区_| 热99re8久久精品国产| 免费一级毛片在线播放高清视频| 亚洲熟女毛片儿| 天堂动漫精品| 国产成人av教育| 日韩欧美在线乱码| 三级国产精品欧美在线观看 | 激情在线观看视频在线高清| 12—13女人毛片做爰片一| 国产精品久久久久久精品电影| av在线天堂中文字幕| 欧美色视频一区免费| 校园春色视频在线观看| 亚洲专区字幕在线| 男人舔女人的私密视频| 男女床上黄色一级片免费看| 丝袜人妻中文字幕| 午夜福利成人在线免费观看| 在线观看舔阴道视频| 午夜福利高清视频| 亚洲激情在线av| 久99久视频精品免费| 亚洲最大成人中文| 久久久久久国产a免费观看| 特大巨黑吊av在线直播| 国产爱豆传媒在线观看 | 亚洲中文字幕日韩| 亚洲九九香蕉| 久久香蕉精品热| 免费看a级黄色片| 精品久久久久久久末码| 无限看片的www在线观看| 午夜免费观看网址| 老司机福利观看| 妹子高潮喷水视频| 亚洲国产精品成人综合色| 欧美最黄视频在线播放免费| 18禁黄网站禁片免费观看直播| 成人手机av| 欧美乱色亚洲激情| 长腿黑丝高跟| 成人三级做爰电影| av中文乱码字幕在线| a级毛片在线看网站| 一级毛片高清免费大全| 99热只有精品国产| 99热6这里只有精品| 亚洲午夜理论影院| 亚洲av片天天在线观看| 日本撒尿小便嘘嘘汇集6| 看免费av毛片| 一边摸一边做爽爽视频免费| 国产又色又爽无遮挡免费看| 最近视频中文字幕2019在线8| 国产精品美女特级片免费视频播放器 | 免费看美女性在线毛片视频| 国产成人系列免费观看| 国内毛片毛片毛片毛片毛片| 99riav亚洲国产免费| 国产片内射在线| www日本在线高清视频| 亚洲真实伦在线观看| 中国美女看黄片| 亚洲国产日韩欧美精品在线观看 | 国产精品1区2区在线观看.| 国产亚洲欧美在线一区二区| 亚洲欧洲精品一区二区精品久久久| 日韩精品青青久久久久久| 国模一区二区三区四区视频 | 在线视频色国产色| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲人成77777在线视频| 欧美中文综合在线视频| 亚洲一区中文字幕在线| 亚洲国产精品合色在线| 精品一区二区三区四区五区乱码| 欧美另类亚洲清纯唯美| 国产精品一区二区三区四区免费观看 | 色哟哟哟哟哟哟| 亚洲,欧美精品.| 亚洲五月婷婷丁香| 久久国产精品人妻蜜桃| 久久香蕉精品热| 亚洲av美国av| 午夜精品久久久久久毛片777| 淫妇啪啪啪对白视频| 大型黄色视频在线免费观看| 在线观看美女被高潮喷水网站 | 日本 av在线| 亚洲欧美激情综合另类| 国产精品久久久久久久电影 | 中文字幕熟女人妻在线| 国产99久久九九免费精品| 中文字幕人成人乱码亚洲影| 国产97色在线日韩免费| 999久久久精品免费观看国产| 亚洲av成人av| www国产在线视频色| 亚洲专区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲电影在线观看av| 成人午夜高清在线视频| 欧美性长视频在线观看| 国产亚洲精品久久久久5区| 亚洲精品av麻豆狂野| 国产欧美日韩一区二区三| 国产单亲对白刺激| 91字幕亚洲| 亚洲五月婷婷丁香| 亚洲熟女毛片儿| 国产成人精品无人区| 亚洲国产日韩欧美精品在线观看 | 日韩有码中文字幕| 999精品在线视频| 色播亚洲综合网| 欧美日本视频| 国产亚洲欧美在线一区二区| 成人一区二区视频在线观看| 亚洲av美国av| 中文在线观看免费www的网站 | 人妻丰满熟妇av一区二区三区| 久久婷婷人人爽人人干人人爱| 两人在一起打扑克的视频| 999久久久国产精品视频| 男插女下体视频免费在线播放| 激情在线观看视频在线高清| 亚洲一码二码三码区别大吗| 亚洲欧美日韩东京热| 久久久国产精品麻豆| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 两人在一起打扑克的视频| 天天躁夜夜躁狠狠躁躁| 99在线视频只有这里精品首页| 一边摸一边做爽爽视频免费| 久久久久国内视频| 午夜激情av网站| 亚洲七黄色美女视频| 9191精品国产免费久久| 久久久久久国产a免费观看| 91老司机精品| 岛国视频午夜一区免费看| 亚洲国产中文字幕在线视频| 亚洲成人中文字幕在线播放| 首页视频小说图片口味搜索| 日韩成人在线观看一区二区三区| 99久久国产精品久久久| 久久欧美精品欧美久久欧美| 99热6这里只有精品| 99久久久亚洲精品蜜臀av| 亚洲五月婷婷丁香| 国产99白浆流出| 嫁个100分男人电影在线观看| 777久久人妻少妇嫩草av网站| 国产高清有码在线观看视频 | 国产一区二区激情短视频| 91老司机精品| 国产高清视频在线观看网站| 精品久久久久久久久久免费视频| 亚洲精品色激情综合| 色尼玛亚洲综合影院| 日韩精品免费视频一区二区三区| 欧美中文综合在线视频| 日韩欧美一区二区三区在线观看| 国产片内射在线| 男女做爰动态图高潮gif福利片| 成熟少妇高潮喷水视频| 一区福利在线观看| 国产久久久一区二区三区| 美女午夜性视频免费| 欧美中文日本在线观看视频| 村上凉子中文字幕在线| 久久久国产成人精品二区| 亚洲精品一区av在线观看| 久久这里只有精品19| 国产高清视频在线播放一区| 一进一出抽搐动态| 97碰自拍视频| 国产日本99.免费观看| www.自偷自拍.com| 两人在一起打扑克的视频| 午夜福利在线在线| 国产区一区二久久| 国产成人系列免费观看| 欧美极品一区二区三区四区| 两个人免费观看高清视频| 人妻久久中文字幕网| 亚洲国产精品sss在线观看| e午夜精品久久久久久久| 国产精品影院久久| 欧美中文日本在线观看视频| 国产一区在线观看成人免费| 国产单亲对白刺激| 三级国产精品欧美在线观看 | 中文字幕久久专区| ponron亚洲| 美女午夜性视频免费| 欧美成人一区二区免费高清观看 | 人妻丰满熟妇av一区二区三区| 欧美av亚洲av综合av国产av| 久久久久国内视频| 亚洲专区字幕在线| АⅤ资源中文在线天堂| 午夜福利成人在线免费观看| 久久精品人妻少妇| 午夜福利在线观看吧| 亚洲av成人一区二区三| 国产三级中文精品| 亚洲精品中文字幕在线视频| 日韩欧美免费精品| 男女那种视频在线观看| 美女午夜性视频免费| 国产精品 欧美亚洲| 一二三四在线观看免费中文在| 亚洲无线在线观看| 午夜成年电影在线免费观看| 午夜视频精品福利| 夜夜爽天天搞| 波多野结衣巨乳人妻| 久久这里只有精品中国| 国产69精品久久久久777片 | 亚洲 欧美 日韩 在线 免费| 国产高清有码在线观看视频 | 亚洲aⅴ乱码一区二区在线播放 | 欧美黄色片欧美黄色片| 最新在线观看一区二区三区| 一边摸一边做爽爽视频免费| 久久中文字幕人妻熟女| 亚洲av成人精品一区久久| 亚洲乱码一区二区免费版| 又大又爽又粗| 超碰成人久久| av有码第一页| 精品第一国产精品| 国产三级中文精品| 久久九九热精品免费| 成人国产综合亚洲| 亚洲av电影在线进入| 99热这里只有是精品50| 久久精品夜夜夜夜夜久久蜜豆 | 欧美日韩国产亚洲二区| 九色成人免费人妻av| 亚洲欧洲精品一区二区精品久久久| 村上凉子中文字幕在线| 精品欧美国产一区二区三| 白带黄色成豆腐渣| 99国产综合亚洲精品| 国产精品电影一区二区三区| 色综合亚洲欧美另类图片| 久久久久久亚洲精品国产蜜桃av| 大型黄色视频在线免费观看| 两性夫妻黄色片| 亚洲色图av天堂| 午夜福利在线观看吧| 精品国产亚洲在线| 久久久精品国产亚洲av高清涩受| 无遮挡黄片免费观看| 亚洲avbb在线观看| 欧美黑人欧美精品刺激| 成人精品一区二区免费| 欧美日韩乱码在线| 国产精品98久久久久久宅男小说| 亚洲五月婷婷丁香| 日韩精品青青久久久久久| 性欧美人与动物交配| avwww免费| 国产片内射在线| 五月玫瑰六月丁香| 亚洲av成人一区二区三| 午夜精品一区二区三区免费看| 男女下面进入的视频免费午夜| 国内久久婷婷六月综合欲色啪| 国产又色又爽无遮挡免费看| 两个人看的免费小视频| 亚洲熟妇熟女久久| 又黄又粗又硬又大视频| 精品久久久久久久末码| 成人av一区二区三区在线看| 久久久久久亚洲精品国产蜜桃av| 99久久精品热视频| 怎么达到女性高潮| 精品国产乱码久久久久久男人| 白带黄色成豆腐渣| 国产欧美日韩一区二区精品| 亚洲美女黄片视频| 国产av一区在线观看免费| 长腿黑丝高跟| 特级一级黄色大片| a在线观看视频网站| 午夜福利高清视频| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩一区二区三| 欧美最黄视频在线播放免费| 午夜亚洲福利在线播放| 午夜两性在线视频| 久久这里只有精品19| av国产免费在线观看| 久久午夜综合久久蜜桃| 露出奶头的视频| 少妇熟女aⅴ在线视频| 十八禁人妻一区二区| 亚洲专区中文字幕在线| 色av中文字幕| 亚洲精品美女久久久久99蜜臀| 亚洲国产欧美网| 99国产精品一区二区蜜桃av| 亚洲精品中文字幕在线视频| 久久久久久亚洲精品国产蜜桃av| 久久天堂一区二区三区四区| 中文字幕av在线有码专区| 男女之事视频高清在线观看| 一区二区三区激情视频| 国内精品久久久久久久电影| av欧美777| 国产成人av激情在线播放| 欧美色视频一区免费| 午夜影院日韩av| 国内久久婷婷六月综合欲色啪| 亚洲熟妇熟女久久| 日本熟妇午夜| 男人的好看免费观看在线视频 | 国产v大片淫在线免费观看| 香蕉av资源在线| 成人精品一区二区免费| 欧美黑人巨大hd| 亚洲一区二区三区色噜噜| 成年免费大片在线观看| 国产精品一及| 91在线观看av| 又大又爽又粗| av天堂在线播放| 国产精品 欧美亚洲| 麻豆国产97在线/欧美 | 淫秽高清视频在线观看| 成人国产综合亚洲| 久久香蕉激情| 巨乳人妻的诱惑在线观看| 亚洲人成网站高清观看| 少妇裸体淫交视频免费看高清 | 啪啪无遮挡十八禁网站| 两个人的视频大全免费| 中出人妻视频一区二区| 免费高清视频大片| 在线十欧美十亚洲十日本专区| 欧美乱妇无乱码| 人人妻人人澡欧美一区二区| 久久午夜综合久久蜜桃| 亚洲精品美女久久久久99蜜臀| 人妻丰满熟妇av一区二区三区| АⅤ资源中文在线天堂| 国产真实乱freesex| 人妻丰满熟妇av一区二区三区| 91麻豆精品激情在线观看国产| svipshipincom国产片| 日韩欧美一区二区三区在线观看| 黑人欧美特级aaaaaa片| 丝袜美腿诱惑在线| 久久精品91蜜桃| 女人被狂操c到高潮| 88av欧美| 国产精品乱码一区二三区的特点| tocl精华| 中文字幕人成人乱码亚洲影| 日日爽夜夜爽网站| 母亲3免费完整高清在线观看| 一区二区三区高清视频在线| 亚洲欧美精品综合久久99| 亚洲熟妇中文字幕五十中出| 欧美日韩瑟瑟在线播放| 国产真实乱freesex| 男女下面进入的视频免费午夜| 可以在线观看毛片的网站| 两个人视频免费观看高清| 悠悠久久av| 国产精品久久电影中文字幕| 亚洲av日韩精品久久久久久密| 国产精品98久久久久久宅男小说| 99热6这里只有精品| 国产欧美日韩一区二区三| 中文资源天堂在线| 国产亚洲av嫩草精品影院| 亚洲国产精品sss在线观看| 国内精品久久久久精免费| 日日干狠狠操夜夜爽| 一进一出抽搐gif免费好疼| 国产精品久久电影中文字幕| 淫秽高清视频在线观看| 熟女电影av网| 狂野欧美激情性xxxx| 国产av不卡久久| 欧美日韩精品网址| 久久久久性生活片| 欧美黄色片欧美黄色片| 少妇熟女aⅴ在线视频| 欧美日韩福利视频一区二区| 亚洲精品久久成人aⅴ小说| 人人妻人人看人人澡| 18禁裸乳无遮挡免费网站照片| 亚洲自拍偷在线| 免费观看精品视频网站| 国产成人精品久久二区二区91| 日韩免费av在线播放| 国产高清videossex| 日韩欧美免费精品| 少妇粗大呻吟视频| 在线播放国产精品三级| 99国产综合亚洲精品| 亚洲精品中文字幕在线视频| 国产精品久久电影中文字幕| 嫩草影院精品99| 亚洲中文字幕日韩| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久| 精品久久久久久久人妻蜜臀av| 天堂√8在线中文| 欧美在线一区亚洲| 岛国在线观看网站| 丝袜人妻中文字幕| 三级国产精品欧美在线观看 | 亚洲国产欧美一区二区综合| 全区人妻精品视频| 日韩有码中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| www.精华液| 亚洲18禁久久av| 中文资源天堂在线| 99久久无色码亚洲精品果冻| 亚洲一区二区三区色噜噜| x7x7x7水蜜桃| 久久精品夜夜夜夜夜久久蜜豆 | 999久久久国产精品视频| 又爽又黄无遮挡网站| 欧美午夜高清在线| 亚洲精品av麻豆狂野| 黄频高清免费视频| 在线a可以看的网站| 麻豆国产av国片精品| 午夜福利18| 级片在线观看| 九九热线精品视视频播放| 亚洲国产精品合色在线| 国产精华一区二区三区|