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

    A distributed eco-hydrological model and its application

    2017-02-01 08:49:12ZongxueXuLeiLiJieZhao
    Water Science and Engineering 2017年4期

    Zong-xue Xu*,Lei Li,Jie Zhao

    aCollege of Water Sciences,Beijing Normal University,Beijing 100875,China

    bBeijing Key Laboratory of Urban Hydrological Cycle and Sponge City Technology,Beijing 100875,China

    cBureau of Hydrology,Ministry of Water Resources of China,Beijing 100053,China

    1.Introduction

    Eco-hydrology was described as an independent discipline for thefirst time at the International Conference on Water and the Environment in Dublin,Ireland,held from January 26 to 31,1992.Studies on eco-hydrological processes in arid regions,whose spatial heterogeneity leads to complex and little described eco-hydrological processes,have received broad attention in thefields of hydrology and water resources recently(Vertessy et al.,1996;Moret et al.,2007).The ecological and hydrological characteristics of arid and semiarid areas have changed due to fragile ecological environments,water resources shortages,and intense disturbance by human activities.

    However,the interplay between vegetation and water is insufficient in hydrological models,and most previous studies have only emphasized the importance ofhydrological processes or ecosystems independently,while the connection between ecological systems and hydrological processes has been neglected.In the traditional hydrological models,the effects of vegetation ecosystems on hydrological processes are usually generalized as an empirical parameter,the value of which isfixed at spatial and temporal scales.In the current ecological models,hydrological processes are mainlyconceptualized in terms of vertical processes of the soilvegetation-atmosphere continuum(SPAC)and water movement in the soil,and the water interaction between plant roots and soil zone interface cannot be described in detail(Xu and Zhao,2016;Zhao et al.,2015).For soil water simulation,the unsaturated zone is generally conceptualized as a reservoir,and water is generally simulated according to the principle of the water balance.In fact,hydrological processes in nature are closely related to the ecological processes of the land,while most hydrological models or ecological models are used to simulate hydrological and ecological processes independently.

    With the development of computer science,remote sensing,and geographic information systems,and especially the improvement of earth observation system technology and digital simulation technology,eco-hydrological models have received more attention recently.At present,eco-hydrological models can be divided into two categories:loosely coupled models with ecological and hydrological simulation results,and one-way coupling models of the hydrological or ecological system.In order to reflect the relationship between water and ecology as a whole,two-way coupling of the water cycle and plant growth was realized in this study using a distributed eco-hydrological model,the ecology module for a grid-based integrated surface and groundwater model(Eco-GISMOD),which was developed to simulate eco-hydrological processes by considering the interaction between vegetation and water in different soil layers.Ecological and hydrological processes are synchronized in Eco-GISMOD through a series of processes of exchanges of water and energy between vegetation and soil,providing a new method for investigating eco-hydrological processes at the river basin scale.

    2.Model description

    Based on Eco-GISMOD,a new distributed ecohydrological simulation system with a three-soil layer structure was developed using modular technology and distributed data management.It includesfive modules:preprocessing,interpolation,evapotranspiration estimation,runoff simulation,and eco-simulation(including natural plants and crops).The model can be executed independently without another platform.Basic data of runoff generation,such asflow direction,simulation sequence,and potential evapotranspiration,are identified in the preprocessing and interpolation modules first.Then,the water allocation and movement of each layer are estimated on the basis of plant transpiration,ground evaporation,and precipitation intercepted by plants.Correspondingly,the dynamic characteristics of leaf area index for different natural vegetation types are computed according to the water requirements and water consumption for each type of vegetation in Eco-GISMOD.

    2.1.Preprocessing

    The pre-processing module is used to delineate the river basin and to extract the drainage networks automatically from the digital elevation model(DEM)data(Li et al.,2013).This module has various functions,such as flow direction definition,flow accumulation calculation,and drainage network generation,and the order of runoff simulation can be analyzed conveniently using a special approach(Xu,2009).

    2.2.Interpolation

    The Thiessen method,gradient plus inverse distance squared method(GIDS),and inverse distance squared method(IDS)were used in the model for spatial interpolation(Nalder and Wein,1998;Teegavarapu and Chandramouli,2005).

    2.3.Evapotranspiration estimation

    Eco-GISMOD integrates eight methods(Allen et al.,1998;Hargreaves and Samani,1985;Priestley and Taylor,1972;Turc,1961)to estimate potential evapotranspiration(PET)and also to enable users to import their own PET data.

    Actual evapotranspiration(AET)can be assumed to be present,depleting water in the root zone of the soil(Beven et al.,1995;Kennen et al.,2008),and is expressed as a simple function of PET and the water content of the surface layer.

    2.4.Runoff simulation

    There are two types of grids in the runoff simulation module:the ordinary grid,which is vertically discretized into three layers:the surface layer,soil layer,and groundwater layer;and the river grid,which looks like a single reservoir for channel routing.In an ordinary grid,each portion of the outflow pours into its corresponding parts when the downstream grid is also an ordinary grid.Otherwise,yielding water will come together as streamflow into the river grid,then move among these river grids and eventually drain out of the watershed(Fig.1).In Fig.1,R is the daily precipitation,E is evapotranspiration,Qr0is the inflow from the upstream neighboring surface layer,Qr1is the inflow from the soil layer,Qr2is the inflow from the groundwater layer,Qr3is the inflow of river grid,Qs1is the saturation excessflow,Qs2is the infiltration excessflow,Qb0is the recharge flow from the soil layer,Qb1is rechargeflow from the groundwater layer,Qc0is the lateral outflow of the surface layer,Qc1is the lateral outflow of the soil layer,Qc2is unconfined groundwater flow,Qc3is confined groundwaterflow,Qx0is the infiltration flow of the surface layer,Qx1is the infiltration flow of the soil layer,Qout0is the outflow of surface layer,Qout1is the outflow of the soil layer,Qout2is the outflow of the groundwater layer,and Qout3is the discharge of river grid.

    Fig.1.Runoff generation structure.

    It is assumed that evapotranspirationfirst takes place in the surface layer.The amount of water used for evapotranspiration is mainly dependent upon the water content of the surface layer.If there is no water in the surface layer,water from the soil layer will be supplied to the surface layer.

    2.5.Eco-simulation

    Plants are crucial in terrestrial hydrological cycles,with their leaves playing key roles during plant development in terms of their functions of photosynthesis and transpiration,but soil moisture,solar radiation,and nutrient elements restrain the plant growth.

    In Eco-GISMOD,the simplified crop growth modules of the environmental policy impact climate(EPIC)model,which was developed by Jimmy Williams from Texas A&M University(Williams et al.,1989),is coupled with a new kind of distributed hydrological model(Li et al.,2015)through two key parameters:leaf area index and soil moisture.It is assumed that the concentration of carbon dioxide in the air is constant,and the plant growth is not restricted by the root system and nutrient substance.The main limiting factors of plant growth are temperature and soil moisture in the model.Suitable temperature and rainfall promote plant growth and development,which increase the leaf area,leading to a larger amount of precipitation interception,leaf evaporation,and transpiration,and a decrease in the soil infiltration.Conversely,the deficiency of soil water restricts plant growth and yield.

    2.5.1.Natural plants

    There are two kinds of natural plants in the model,forest and grass.Their plant roots absorption is different from one another due to the root distribution functions.Forests can fully use the water in the soil with a developed root system,but grass can only uptake water from the surface soil.

    First,the daily cumulative effective heat(heat unit,simplified as N)is computed according to the daily maximum and minimum temperature,and the plant growth reference temperature.

    Then,the possible daily increase of biomass is estimated based on energy-biomass conversion parameters,the day length,latitude and longitude,and effective photosynthesis radiation,which are obtained from the relationship between solar radiation and the leaf area index.

    Third,the impact of soil moisture and temperature is considered,and the greatest impact factor is selected as the plant growth-limiting factor.

    The heat unit formula is as follows:

    where i is number of days;j is the plant type;Tmaxand Tminare the daily maximum temperature and minimum temperature,respectively;and Tbaseis the minimum temperature required for plant growth.

    The heat unit index(Z),which reflects the status of plant development throughout the growing season,can be obtained by comparing N and the total heat for plant needs from growth to maturity(M).The formula is as follows:

    Plant growth is also closely related to solar radiation.The photosyntheticallyactiveradiation(X,MJ·m-2·d-1)intercepted byplantscanbecalculatedfromsolarradiation(S,MJ·m-2·d-1)and leaf area index(I)according to Beer's Law:

    Daily increased biomass(ΔB,kg·d-1)is affected by X,and the energy-biomass conversion parameter(E,kg·MJ-1):

    where D is day length(h)and ΔD is the change in day length(h·d-1),which is computed as

    where L is the latitude(o)and V is the daily solar elevation angle(o):

    The final increased daily biomass(ΔF,kg·d-1)is mainly affected by the daily plant stress index(G)and ΔB:

    2.5.1.1.Leaf area index

    Leaf surface is the most direct and rapid channel for changing the substance and energy between plants and the environment,and is involved in ecosystem water circulation through evapotranspiration and interception(Running and Coughlan,1988;Tian et al.,2003;Wasseige et al.,2003).In Eco-GISMOD,rainfall is divided into two parts.One part is net rainfall that falls into the soil,and the other part is what is intercepted on the leaf surface,and will evaporate quickly.Therefore,leaf area index is a very important indicator for estimating and evaluating the water consumption of vegetation.The leaf area index is estimated for growing and declining periods individually.

    (1)Growing period:

    Plants are considered in the growing period when the heat unit is less than or equal to 1(Zi≤ 1),the leaf area index of which is computed as

    where Iiis the leaf area index of day i,Imaxis the maximum leaf area index,ΔI is the variation of the leaf area index,and ΔU is the variation of the heat unit factor(U),which is calculated as follows:

    The default values of θ1and θ2are set as 6.5 and 10 in this formula,respectively,which can be modified under different situations.

    G ranges between 0 and 1 in the model,which considers the effect of soil moisture and temperature on plant growth,but ignores the impact of nitrogen,phosphorus,and root development.The formula is as follows:

    where Ws,iis the water stress factor of day i,and Ts,iis the temperature factor of day i.

    As one of the key limiting factors in the vegetation,Ts,iis obtained from the following formula:

    where Tg,jis the average temperature required for the growth of plant j(°C),Tb,jis the basic temperature for the growth of plant j(°C),and To,jis the optimum temperature for the growth of plant j(°C).

    Another limiting factor,Ws,i,is determined by the relationship between the available water supply from the soil layer(Wsupply,i)and the water requirement of the plant(Wdemand,i).Wswill be 0 if Wsupply,iis greater than Wdemand,i.Otherwise,Ws,iis computed as follows:

    where h0is the height of the soil layer(m);Agis grid area(m2);Sais the actual water content of soil;k0is the hydraulic conductivity of the soil layer(m3·d-1);Peis the potential evaporation;and j0is a factor that depends on the relationships between Sa,the critical soil water content(Sc),and the wilting point of the soil water content(Sw):

    where v is an empirical coefficient,the value of which is between 0.80 and 0.95.The available water supply comes mainly from the soil layer when the plant type is forest,and from the surface layer when the plant type is grass.(2)Declining period:

    Plants are regarded as in the declining period when the heat unit index is greater than 1(Zi>1),the leaf area index of which is computed as

    where I0and Z0are the corresponding values of leaf area index and heat unit index when plants begin declining,and Cjis the fading rate of plant j.

    2.5.1.2.Dry matter

    The quality of dry matter from growth to maturity is determined by the plant harvest index of day i(Ai):

    (1)Without considering water shortage:

    where Hjis the harvest index of plant j during the whole growth period,and Ukis the heat unit factor of day k.

    (2)Considering water shortage:

    where Yjis the drought sensitivity index of plant j and Wirepresents the water consumption on day i:

    Finally,the harvested dry matter from growth to maturity Ydis obtained as follows:

    where n is the total number of days.

    2.5.2.Crops

    In contrast to natural plants,crops mainly depend on irrigation,and their growth and development have regularity.For instance,spring wheat is sown in the middle of March,and harvested in May or June.The water consumption of summer corn from jointing to heading periods accounts for 65%of the whole growth period.However,summer corn is usually sown in May and harvested in September or October.The water consumption of each growth stage is not very different.

    According to the characteristics of different crop types,the crop water consumption is estimated by combining the crop coefficient method and potential evapotranspiration method.The crop growth period is divided into four stages(seeding,jointing,heading,and maturity),and the water consumption of each stage can be obtained by multiplying the potential evapotranspiration by the crop coefficient.

    2.5.2.1.Leaf area index

    Variation of the crop leaf area is simulated by the photosynthesis module of the crop,which is a kind of terrestrial net primary productivity(NPP)estimation model developed by Huang et al.(2006).The formula is as follows:

    where Jiis the ratio of yellow leaves to green leaves on day i,Kiis the transfer coefficient of photosynthetic products(0-0.48 for spring wheat,and 0-0.54 for corn)on day i,Pn,iis the net photosynthesis on day i,and Slis specific leaf area(m2/g).

    The photosynthetic products are not transferred to the leaf surface during the heading and maturity period.The leaf area index can be computed as follows:

    where δ and ε are key parameters that are determined by crop type and soil moisture,with values ranging from 0 to 8 and 0.5 to 1.0,respectively;and O is the standard fertility index:

    where Atis the effective accumulated temperature from planting to harvesting(°C),B is the initial air temperature of crop growth(°C),Tiis daily mean temperature on day i(°C),and x is the cumulative number of days from sowing.

    2.5.2.2.Actual water consumption

    Water consumption of crop j on day i(Wi,j)is determined by the following formula:

    where Ks,jis the coefficient of crop j in the growing stage s;and efis the irrigation efficiency,the value of which can be as low as 0.4 to 0.5 withflood irrigation,but as high as 0.6 to 0.8 with micro-sprinkle irrigation and drip irrigation.

    Irrigation water is assumed mainly to be taken from rivers.Therefore,the model requires some water intakes(also called control points),which can be used to calculate the total amount of water consumption in the watershed.

    where t is the number of crop types and m is the total number of days of crop growth.

    2.5.2.3.Dry matter

    Crop dry matter can be obtained without considering water shortage(setting Wsas 0),using the same formula as that of natural plants.

    3.Case study

    In order to evaluate the performance of Eco-GISMOD,the upper and middle parts of the Heihe River Basin were selected as the study area,where plant growth is severely restricted by seasonaldroughtand watershortage.The period of 1990-1993 was chosen for the case study.The calibration period was set from 1990 to 1992,and 1993 was selected for validation.Observed precipitation and discharge data were provided by the Cold and Arid Regions Environmental and Engineering Research Institute(CAREER)of the Chinese Academy of Sciences.DEM data with 1 km×1 km horizontal resolution was used and other spatial distribution data with the same resolution,such as land use data,soil data,and geological data,were all obtained from the Environmental and Ecological Science Data Center for West China.

    3.1.Parameter setting

    Requirements for air temperature,solar radiation,and soil moisture vary according to vegetation type and the spatial distribution of each kind of plant.The vegetation was divided intofive types(coniferous forest,broad-leaved forest,shrub,farmland,and grass)from mountainous areas to plains in the Heihe River Basin.The initial values of plant growth parameters were chosen according to the parameters recommended in the EPIC model(Table 1).The values of parameters in the hydrological module were obtained from the results of previous studies(Li et al.,2015).

    3.2.Results and discussion

    3.2.1.Runoff

    The accuracy of runoff simulation plays an important role in plant growth and the development of the eco-hydrological model.Therefore,a river grid near the Yingluoxia Hydrological Station was used to make a comparison with the measured runoff in this study.It can be seen from Fig.2 that Eco-GISMOD reproduces the monthly runoff series with Nash-Sutcliffe efficiency coefficient(NSE)values between 0.88 and 0.93,the corresponding deviation of which ranges from-2.1%to 2.3%.The performances of Eco-GISMOD at theYingluoxia Hydrological Station are acceptable in both the calibration and validation periods.

    Table 1 Initial values of plant growth parameters.

    3.2.2.NPP

    The totalbiomass amountconverted from energy throughout the growth period is calculated by considering the limit factor of plant growth and potential biomass,through subtraction of the self-energy consumption of vegetation.NPP cannot only reflect the production capacity of ecosystems in the natural environment,but is also the key factor in evaluating the quality of the terrestrial ecosystem.

    For the purpose of comparison with other study results,we assume that the daily increased biomass is only used for plant growth,development,and reproduction.Other energy consumption was not considered in this study.Therefore,a comparative analysis of the results of Eco-GISMOD,the data obtained from high-resolution Satellite Pour l’Observation de la Terre(SPOT)remote sensing images using the C-FIX model(Lu et al.,2005),and the data calculated from the TESim model(Peng,2007)was used for model validation in this study(Fig.3).In general,results show that the NPP values of shrub, grass, and farmland (223.2, 246.4, and 394.6 g·m-2·year-1)from Eco-GISMOD are close to those of the C-FIX model(268,244,and 364 g·m-2·year-1),but the NPP values of coniferous forest and broad-leaved forest from Eco-GISMOD(29.2 and 91.8 g·m-2·year-1)are much smaller than those calculated by the C-FIX model(266 and 184 g·m-2·year-1).For the TESim model,in addition to the simulated NPP of broad-leaved forest that is close to the result of the C-FIX model(266.7 g·m-2·year-1),there are signifi-cant differences between the C-FIX model and Eco-GISMOD with the other four vegetation types.

    Fig.2.Comparison between simulated and observed monthly streamflow at Yingluoxia Hydrological Station.

    Fig.3.Comparison of simulated net primary productivity from three models.

    Although the results of these three models are different from one another,Eco-GISMOD shows a strong performance of plant growth,especially for shrub,grass,and farmland.For instance,Hu et al.(1994)and Wang et al.(1988)found that the value of NPP of high mountain meadows from the Qilian Mountains varies from 150 to 240 g·m-2·year-1,which is consistent with the results of Eco-GISMOD.Chen et al.(2008)showed that the value of grass NPP ranges from 128 to 370 g·m-2·year-1,based on statisticsofaboveground biomass obtained using the harvesting method.According to the estimation of farmland biomass in arid areas of western China from 1992 to 1998,some studies have pointed out that the value of farmland NPP is about 390.42 g·m-2·year-1(Zhang et al.,2006),which is in accordance with the simulated results of Eco-GISMOD.

    Fig.4.Spatial distribution of annual average net primary productivity(units:g·m-2·year-1).

    Table 2 Variation of net primary productivity with different simulation conditions.

    Spatial distribution of annual average NPP was analyzed in detail(Fig.4).There was a strong relationship between NPP and topography.A high value of NPP is mainly distributed in the middle part of the agricultural oasis.The lower values are mainly found at the southern part of the Qilian Mountains.Without consideration of water stress,the increase of NPP mainly occurs in the middle part of the Heihe River Basin,where it is dominated by grass and shrub(Table 2).The value of NPP of grassland increases significantly,with a range from 41 to 57 g·m-2·year-1,and the NPP increase for shrub is between 6.9 and 17.4 g·m-2·year-1.The broad-leaved forest and coniferous forest experience less water stress,and the values of NPP of these two types of plants change within 3 g·m-2·year-1.

    3.2.3.Water productivity

    Thewaterproductivitiesofthesefivetypesofvegetationwere compared with those of previous studies,in order to assess the performance of Eco-GISMOD in the simulation of plant growth.Inviewoftheoverallsimulationresults,thewaterproductivityof farmlandishighest,ranging from1.04to1.44kg/m3andwithan average annual value of 1.23 kg/m3.Wu and Hu(2009)also found that the water productivity of crops is between 0.93 and 1.53 kg/m3,based on research in Gansu Province from 1995 to 2007.The water productivity of coniferous forest ranks second,as simulated by Eco-GISMOD,with values changing between 0.19 and 0.41 kg/m3.Through the analysis of a forest plantation in Qinghai Province,Liu et al.(2006)pointed out that the water productivity of Qinghai spruce and larch coniferous forest varies from 0.19 to 0.26 kg/m3,as estimated according to the water efficiency rate(about 44%)and water potential productivity(0.44-0.58 kg/m3).Water productivity of shrub and grassland simulatedbyEco-GISMODisthelowestamongthesevegetation types,with annual mean values of 0.016 and 0.014 kg/m3,respectively.Chen et al.(2006)simulated the water productivity ofartificialgrasslandontheLoessPlateauusingtheEPICmodel.Results show that the water potential productivity of alfalfa and other pasture vegetation is between 0.05 and 0.58 and kg/m3.These facts illustrate that Eco-GISMOD performs well with regard to water productivity simulation,especially for the simulation of crops and coniferous forests.

    4.Conclusions

    In orderto evaluate the eco-hydrologicalprocesses focusing on water exchange between plants and the surrounding natural environment at a catchment scale,Eco-GISMOD,anew kind ofeco-hydrologicalmodelis proposed.It combines the simplified EPIC model and a new kindofdistributedhydrologicalmodelGISMOD.The advantage of this model is that water interaction between plants and soil water can be simply described by a generalized physical process in various situations.Eco-GISMOD can reflect the plant growth status through leaf area index.Dry matter and the spatial and temporal variations of water exchange can also be estimated using this model.

    A case study was carried out in the Heihe River Basin in northwestern China to evaluate the system's performance.The results show that forest and crops generally grow well with enough water supply,but water shortages,especially in summer,inhibit the growth of grass and cause grass degradation.This demonstrates that the simulation system is appropriate for eco-hydrological study in arid and semi-arid areas and will be helpful for the planning and management of water resources in the future.

    Allen,R.G.,Pereira,L.S.,Raes,D.,Smith,M.,1998.Crop evapotranspiration:guidelines for computing crop water requirements-FAO Irrigation and Drainage Paper 56.Food and Agriculture Organization(FAO),Rome.

    Beven,K.J.,Lamb,R.,Quinn,P.F.,Romanowicz,R.,Freer,J.,1995.TOPMODEL.In:Singh,V.P.(Ed.),Computer models of watershed hydrology.Water Resources Publications,Highlands Ranch,pp.627-668.

    Chen,B.,Li,J.,Li,X.F.,2006.Long term simulation of water potential productivity of alfalfa on rain fed highland in south of Loess Plateau.Agric.Res.Arid Areas 24(3),31-35.https://doi.org/10.3321/j.issn:1000-7601.2006.03.006(in Chinese).

    Chen,Z.H.,Ma,Q.Y.,Wang,J.,Qi,Y.,Li,J.,Huang,C.L.,Ma,M.G.,Yang,G.J.,2008.Estimation of Heihe Basin net primary productivity using the CASA model.J.Nat.Resour.23(2),263-273.https://doi.org/10.11849/zrzyxb.2008.02.011(in Chinese).

    Hargreaves,G.H.,Samani,Z.A.,1985.Reference crop evapotranspiration from temperature.Appl.Eng.Agric.1(2),96-99.https://doi.org/10.13031/2013.26773.

    Hu,Z.Z.,Sun,J.X.,Li,Y.,Long,R.J.,Yang,F.L.,1994.The characteristics of biomass and conversion efficiency of solar radiation for principal types of alpine grasslands in Tianzhu,Gansu Province,China.Chin.J.Plant Ecol.18(2),121-131(in Chinese).

    Huang,Y.,Wang,Y.,Zhang,W.,Yu,Y.Q.,Wang,P.,2006.Simulating net primary production of agricultural vegetation in China(I):Model establishment and sensitivity analysis.J.Nat.Resour.21(5),790-801.https://doi.org/10.11849/zrzyxb.2006.05.013(in Chinese).

    Kennen,J.G.,Kauffman,L.J.,Ayers,M.A.,Wolock,D.M.,Colarullo,S.J.,2008.Use of an integratedflow model to estimate ecologically relevant hydrologic characteristics at stream biomonitoring sites.Ecol.Model.211(1-2),57-76.https://doi.org/10.1016/j.ecolmodel.2007.08.014.

    Li,L.,Xu,Z.X.,Li,Y.L.,Liu,W.F.,Zhang,L.Y.,2013.A preprocessing program for a distributed hydrological model:Development and application.J.Hydroinf.15(4),1258-1275.https://doi.org/10.2166/hydro.2013.140.

    Li,L.,Xu,Z.X.,Zuo,D.P.,Zhao,J.,2015.A grid-based integrated surfacegroundwater model(GISMOD).J.Water Clim.Change 7(2),296-320.https://doi.org/10.2166/wcc.2015.006.

    Liu,S.,He,K.N.,Wang,Z.N.,Chang,G.L.,Li,S.R.,2006.Extinction characteristics of artificial forest in removal lands in Datong of Qinghai Province.Sci.Soil Water Conserv.4(3),59-64.https://doi.org/10.3969/j.issn.1672-3007.2006.03.012(in Chinese).

    Lu,L.,Li,X.,Frank,V.,2005.Estimation of net primary productivity of Heihe River Basin using remote sensing.J.Desert Res.25(6),823-830(in Chinese).

    Moret,D.,Braud,I.,Arru′e,J.L.,2007.Water balance simulation of a dryland soil during fallow under conventional and conservation tillage in semiarid Aragon,Northeast Spain.Soil Tillage Res.92(1-2),251-263.https://doi.org/10.1016/j.still.2006.03.012.

    Nalder,I.A.,Wein,R.W.,1998.Spatial interpolation of climate normals:Test of a new method in the Canadian boreal forest.Agric.For.Meteorol.92(4),211-225.https://doi.org/10.1016/S0168-1923(98)00102-6.

    Peng,H.C.,2007.Study on Dynamic Simulation of Ecosystem in Heihe River Basin.Ph.D.Dissertation.Chinese Academy of Sciences,Lanzhou(in Chinese).

    Priestley,C.H.B.,Taylor,R.J.,1972.On the assessment of surface heatflux and evaporation using large scale parameters.Mon.Weather Rev.100(2),81-92.https://doi.org/10.1175/1520-0493(1972)100<0081:OTAOSH>2.3.CO;2.

    Running,S.W.,Coughlan,J.C.,1988.A general model of forest ecosystem processes for regional applications,I:Hydrologic balance,canopy gas exchange and primary production processes.Ecol.Model.42(2),125-154.https://doi.org/10.1016/0304-3800(88)90112-3.

    Teegavarapu,R.S.V.,Chandramouli,V.,2005.Improved weighting methods,deterministic and stochastic data-driven models for estimation of missing precipitation records.J.Hydrol.312(1-4),191-206.https://doi.org/10.1016/j.jhydrol.2005.02.015.

    Tian,Y.H.,Wang,Y.J.,Zhang,Y.,Knyazikhin,Y.,Bogaert,J.,Myneni,R.B.,2003.Radiative transfer based scaling of LAI retrievals from reflectance data of different resolutions.Rem.Sens.Environ.84(1),143-159.https://doi.org/10.1016/S0034-4257(02)00102-5.

    Turc,L.,1961.Estimation of irrigation water requirements,potential evapotranspiration:A simple climatic formula evolved up to date.Ann.Agron.12(1),13-49.

    Vertessy,R.A.,Hatton,T.J.,Benyon,R.G.,Dawes,W.R.,1996.Long-term growth and water balance predictions for a mountain ash(Eucalyptus regnans)for catchment subject to clear-felling and regeneration.Tree Physiol.16(1-2),221-232.https://doi.org/10.1093/treephys/16.1-2.221.

    Wang,Q.J.,Yang,F.T.,Shi,S.H.,1988.A preliminary study on the formation of underground biomass in alpine meadow.In:Proceedings of the International Symposium Alpine Meadow Ecosystem,pp.73-80(in Chinese).

    Wasseige,C.D.,Bastin,D.,Defourny,P.,2003.Seasonal variation of tropical forest LAI based onfield measurements in Central African Republic.Agric.For.Meteorol.119(3-4),181-194.https://doi.org/10.1016/S0168-1923(03)00138-2.

    Williams,J.R.,Jones,C.A.,Kiniry,J.R.,Spanel,D.A.,1989.The EPIC crop growth model.Transactions of the ASABE 32(2),497-511.https://doi.org/10.13031/2013.31032.

    Wu,J.X.,Hu,G.L.,2009.Interannual change on water productivity of crop irrigation in Gaotai County.Developing(8),116-117(in Chinese).

    Xu,Z.X.,2009.Hydrological Models.Science Press,Beijing(in Chinese).

    Xu,Z.X.,Zhao,J.,2016.Development and applications of eco-hydrological models:Past and future.J.Hydraul.Eng.47(3),346-354.https://doi.org/10.13243/j.cnki.slxb.20151160(in Chinese).

    Zhang,J.,Pan,X.L.,Gao,Z.Q.,Shi,Q.D.,Lu¨,G.H.,2006.Satellite estimates and change detection of net primary productivity of oasis-desert based on ecosystem process with remotely sensed forcing in arid western China.Acta Geograph.Sin.61(1),15-25.https://doi.org/10.11821/xb200601002(in Chinese).

    Zhao,J.,Xu,Z.X.,Zuo,D.P.,Wang,X.M.,2015.Temporal variations of reference evapotranspiration and its sensitivity to meteorological factors in Heihe River Basin,China.Water Sci.Eng.8(1),1-8.https://doi.org/10.1016/j.wse.2015.01.004.

    欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲av一区麻豆| 巨乳人妻的诱惑在线观看| 黄色毛片三级朝国网站| 午夜两性在线视频| 18禁裸乳无遮挡免费网站照片 | 亚洲真实伦在线观看| 脱女人内裤的视频| 久99久视频精品免费| 久久亚洲精品不卡| 搡老熟女国产l中国老女人| 一本久久中文字幕| 老汉色∧v一级毛片| 日日摸夜夜添夜夜添小说| 亚洲成人久久爱视频| 亚洲激情在线av| 午夜日韩欧美国产| 青草久久国产| 日韩成人在线观看一区二区三区| 中文字幕精品免费在线观看视频| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品成人综合色| 日韩免费av在线播放| 一区二区三区激情视频| 精品久久久久久久久久免费视频| 婷婷丁香在线五月| 精品欧美国产一区二区三| 99国产精品一区二区蜜桃av| 免费看美女性在线毛片视频| 在线观看免费视频日本深夜| 国产精品 欧美亚洲| aaaaa片日本免费| 亚洲无线在线观看| 亚洲国产精品sss在线观看| 动漫黄色视频在线观看| 国产成人av激情在线播放| 中国美女看黄片| 国产一区二区在线av高清观看| 12—13女人毛片做爰片一| 国内精品久久久久久久电影| 国产精品久久久久久精品电影 | 人人妻,人人澡人人爽秒播| 日日夜夜操网爽| 中文字幕人妻熟女乱码| 免费高清在线观看日韩| 国产av一区二区精品久久| 亚洲全国av大片| 婷婷亚洲欧美| 脱女人内裤的视频| 国产亚洲欧美精品永久| 又大又爽又粗| 19禁男女啪啪无遮挡网站| 99国产极品粉嫩在线观看| 真人做人爱边吃奶动态| 日韩高清综合在线| 久久久水蜜桃国产精品网| 久久精品国产亚洲av香蕉五月| 亚洲国产看品久久| 波多野结衣av一区二区av| 韩国精品一区二区三区| 亚洲中文日韩欧美视频| 久久性视频一级片| 免费高清视频大片| 欧美av亚洲av综合av国产av| 国产精品1区2区在线观看.| 黄色a级毛片大全视频| 午夜老司机福利片| 精品免费久久久久久久清纯| 哪里可以看免费的av片| 啦啦啦免费观看视频1| 男女午夜视频在线观看| 亚洲中文字幕日韩| 白带黄色成豆腐渣| 亚洲午夜理论影院| 亚洲国产欧美网| 午夜免费鲁丝| 久久性视频一级片| 性色av乱码一区二区三区2| 免费看a级黄色片| 国产熟女xx| 日日摸夜夜添夜夜添小说| 欧美日韩一级在线毛片| 免费电影在线观看免费观看| 欧美av亚洲av综合av国产av| 久久婷婷成人综合色麻豆| 国产精品 国内视频| 美女国产高潮福利片在线看| 色尼玛亚洲综合影院| 欧美一级毛片孕妇| 日本免费a在线| 久久久国产欧美日韩av| 国产亚洲av高清不卡| 亚洲av片天天在线观看| 欧美色欧美亚洲另类二区| 亚洲欧美一区二区三区黑人| 亚洲av电影在线进入| 女人高潮潮喷娇喘18禁视频| 色婷婷久久久亚洲欧美| 人人妻人人看人人澡| 国产单亲对白刺激| 一个人观看的视频www高清免费观看 | 日日夜夜操网爽| 久久久国产成人免费| a级毛片a级免费在线| 99在线人妻在线中文字幕| 一级黄色大片毛片| 国产又黄又爽又无遮挡在线| 久久午夜亚洲精品久久| 欧美在线一区亚洲| avwww免费| 动漫黄色视频在线观看| 国产精品乱码一区二三区的特点| 国产99久久九九免费精品| 一级作爱视频免费观看| 久久久水蜜桃国产精品网| 亚洲成人国产一区在线观看| 丝袜在线中文字幕| 看黄色毛片网站| 色av中文字幕| 国产单亲对白刺激| 1024视频免费在线观看| 国产午夜精品久久久久久| 国产欧美日韩一区二区三| 12—13女人毛片做爰片一| 无限看片的www在线观看| 成年女人毛片免费观看观看9| 久久这里只有精品19| 成人午夜高清在线视频 | 亚洲成国产人片在线观看| 大型av网站在线播放| 日本精品一区二区三区蜜桃| 午夜a级毛片| 亚洲国产欧美日韩在线播放| 757午夜福利合集在线观看| 9191精品国产免费久久| 亚洲成人久久爱视频| 亚洲一码二码三码区别大吗| 男女之事视频高清在线观看| 免费电影在线观看免费观看| 欧美又色又爽又黄视频| 亚洲中文字幕日韩| 国产乱人伦免费视频| 国产激情欧美一区二区| 国产成人啪精品午夜网站| 中出人妻视频一区二区| 在线观看免费日韩欧美大片| www.www免费av| 精品免费久久久久久久清纯| 91麻豆精品激情在线观看国产| 一进一出抽搐gif免费好疼| 熟妇人妻久久中文字幕3abv| 成人精品一区二区免费| 老司机午夜福利在线观看视频| 香蕉国产在线看| 高清毛片免费观看视频网站| 精品久久久久久久毛片微露脸| 丰满的人妻完整版| 18禁黄网站禁片免费观看直播| 好男人电影高清在线观看| 美女午夜性视频免费| 男女做爰动态图高潮gif福利片| 美女高潮喷水抽搐中文字幕| 99精品欧美一区二区三区四区| 国产三级黄色录像| 色综合站精品国产| 久久中文字幕一级| 男女之事视频高清在线观看| 亚洲黑人精品在线| 精品欧美国产一区二区三| www国产在线视频色| 无人区码免费观看不卡| 久久久国产成人精品二区| 亚洲天堂国产精品一区在线| 在线观看日韩欧美| 别揉我奶头~嗯~啊~动态视频| 中国美女看黄片| 久久婷婷人人爽人人干人人爱| 韩国av一区二区三区四区| 香蕉久久夜色| 亚洲一区二区三区不卡视频| 精品久久久久久久久久免费视频| 精品第一国产精品| 国产成人av激情在线播放| 日日摸夜夜添夜夜添小说| 亚洲午夜理论影院| 欧美 亚洲 国产 日韩一| 在线十欧美十亚洲十日本专区| 精品不卡国产一区二区三区| 国产成人欧美| 韩国精品一区二区三区| 久久久国产欧美日韩av| 香蕉av资源在线| 老汉色av国产亚洲站长工具| 日韩av在线大香蕉| 麻豆成人午夜福利视频| 国产人伦9x9x在线观看| 国产亚洲精品av在线| 国产亚洲精品一区二区www| 国产精品99久久99久久久不卡| 午夜免费成人在线视频| 亚洲精品久久成人aⅴ小说| 视频区欧美日本亚洲| 男人的好看免费观看在线视频 | www.999成人在线观看| 搡老妇女老女人老熟妇| 国产真实乱freesex| 国产精品久久电影中文字幕| 国产1区2区3区精品| 脱女人内裤的视频| 一级a爱视频在线免费观看| 婷婷精品国产亚洲av在线| 亚洲第一青青草原| 久久中文字幕一级| 给我免费播放毛片高清在线观看| 99久久无色码亚洲精品果冻| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久九九精品二区国产 | 熟妇人妻久久中文字幕3abv| 色老头精品视频在线观看| 99久久久亚洲精品蜜臀av| 久久天堂一区二区三区四区| 久久久久久久精品吃奶| 国产麻豆成人av免费视频| 成人国语在线视频| 两人在一起打扑克的视频| 久久伊人香网站| 最近最新中文字幕大全电影3 | 国产亚洲av高清不卡| 伊人久久大香线蕉亚洲五| 久久人人精品亚洲av| 男女床上黄色一级片免费看| 久久亚洲真实| 欧美色欧美亚洲另类二区| 我的亚洲天堂| 亚洲最大成人中文| 天天躁夜夜躁狠狠躁躁| 亚洲精品美女久久久久99蜜臀| 亚洲五月天丁香| 国产单亲对白刺激| 精品卡一卡二卡四卡免费| 亚洲精品中文字幕在线视频| 欧美激情久久久久久爽电影| 精品一区二区三区视频在线观看免费| 变态另类成人亚洲欧美熟女| 少妇粗大呻吟视频| 免费在线观看影片大全网站| 国产成人系列免费观看| 黄色视频不卡| 欧美激情 高清一区二区三区| 久久精品91无色码中文字幕| av超薄肉色丝袜交足视频| 久久久久国产一级毛片高清牌| 免费看日本二区| 人妻丰满熟妇av一区二区三区| 我的亚洲天堂| 色综合婷婷激情| 国产熟女xx| 少妇熟女aⅴ在线视频| 精品国产超薄肉色丝袜足j| 变态另类成人亚洲欧美熟女| 啦啦啦免费观看视频1| 日韩精品免费视频一区二区三区| 久久精品影院6| 国产黄a三级三级三级人| 满18在线观看网站| 男女视频在线观看网站免费 | 国产精品 国内视频| 母亲3免费完整高清在线观看| 很黄的视频免费| 国产视频内射| 国产亚洲欧美精品永久| 日韩欧美国产在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产精品国产高清国产av| 国产黄片美女视频| 黄片小视频在线播放| 久久精品影院6| 国产精品二区激情视频| 成人亚洲精品一区在线观看| 精品久久久久久久毛片微露脸| 国产亚洲精品久久久久5区| 成年人黄色毛片网站| 久久久国产成人免费| 国产精品一区二区免费欧美| 国产精品电影一区二区三区| 国产精品国产高清国产av| 国产视频一区二区在线看| 伊人久久大香线蕉亚洲五| 搡老岳熟女国产| 久久亚洲精品不卡| 日日干狠狠操夜夜爽| 视频在线观看一区二区三区| 成熟少妇高潮喷水视频| 国产久久久一区二区三区| 精品欧美一区二区三区在线| 少妇的丰满在线观看| 丰满人妻熟妇乱又伦精品不卡| 无限看片的www在线观看| 日本熟妇午夜| 免费无遮挡裸体视频| 侵犯人妻中文字幕一二三四区| 一区二区三区激情视频| 搡老岳熟女国产| 欧美午夜高清在线| 亚洲av第一区精品v没综合| 欧美黑人巨大hd| 国产单亲对白刺激| 中文字幕最新亚洲高清| 欧美日本视频| 老熟妇乱子伦视频在线观看| 亚洲国产欧美日韩在线播放| 美女大奶头视频| 91麻豆精品激情在线观看国产| 亚洲一码二码三码区别大吗| 禁无遮挡网站| 国产在线精品亚洲第一网站| 日日爽夜夜爽网站| 国产亚洲精品久久久久久毛片| 好男人电影高清在线观看| 最好的美女福利视频网| 亚洲精品久久成人aⅴ小说| 亚洲av成人av| 国产精品久久视频播放| 久久香蕉激情| 人妻丰满熟妇av一区二区三区| 久久精品亚洲精品国产色婷小说| 日本五十路高清| 久久天堂一区二区三区四区| 人成视频在线观看免费观看| 欧美人与性动交α欧美精品济南到| 亚洲国产中文字幕在线视频| 亚洲 欧美一区二区三区| 亚洲最大成人中文| 老汉色∧v一级毛片| 日韩 欧美 亚洲 中文字幕| 久久久久久久午夜电影| a级毛片在线看网站| 成熟少妇高潮喷水视频| 亚洲欧美日韩高清在线视频| 亚洲中文字幕一区二区三区有码在线看 | 欧美+亚洲+日韩+国产| 黄色女人牲交| 少妇裸体淫交视频免费看高清 | 美女午夜性视频免费| 满18在线观看网站| 正在播放国产对白刺激| 欧美日韩黄片免| 欧美性猛交黑人性爽| 97超级碰碰碰精品色视频在线观看| 91九色精品人成在线观看| 无遮挡黄片免费观看| 18禁美女被吸乳视频| 一卡2卡三卡四卡精品乱码亚洲| 国产免费男女视频| 亚洲成人精品中文字幕电影| 91成年电影在线观看| 中文字幕精品亚洲无线码一区 | 精品一区二区三区视频在线观看免费| 欧美黑人欧美精品刺激| 国产成人影院久久av| 亚洲欧美日韩高清在线视频| 日韩精品青青久久久久久| 国产黄色小视频在线观看| 国产片内射在线| 国产精品98久久久久久宅男小说| 午夜免费观看网址| 久久中文字幕一级| 色综合欧美亚洲国产小说| 99国产精品一区二区三区| 亚洲欧美日韩无卡精品| 19禁男女啪啪无遮挡网站| 首页视频小说图片口味搜索| 成人永久免费在线观看视频| 在线观看一区二区三区| 免费搜索国产男女视频| 黄色丝袜av网址大全| 国产黄片美女视频| 久久精品夜夜夜夜夜久久蜜豆 | 波多野结衣高清无吗| 天堂动漫精品| 成人国语在线视频| 一本大道久久a久久精品| 精华霜和精华液先用哪个| 日日干狠狠操夜夜爽| 91麻豆av在线| 美女国产高潮福利片在线看| 一区福利在线观看| 午夜福利在线观看吧| 国产aⅴ精品一区二区三区波| 日本黄色视频三级网站网址| 男女床上黄色一级片免费看| 国产高清有码在线观看视频 | 女性生殖器流出的白浆| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕高清在线视频| 亚洲精品一区av在线观看| 狠狠狠狠99中文字幕| 亚洲精品国产一区二区精华液| 一级黄色大片毛片| 1024香蕉在线观看| 久99久视频精品免费| 精品卡一卡二卡四卡免费| 日本 av在线| 中文字幕人妻熟女乱码| 日本熟妇午夜| 国产精品一区二区三区四区久久 | 看免费av毛片| 此物有八面人人有两片| 老鸭窝网址在线观看| 黑人欧美特级aaaaaa片| 欧美日韩精品网址| 男人舔女人下体高潮全视频| 夜夜爽天天搞| 天堂影院成人在线观看| 亚洲成人久久性| 性色av乱码一区二区三区2| 中文在线观看免费www的网站 | 婷婷精品国产亚洲av| 亚洲精品国产精品久久久不卡| 十分钟在线观看高清视频www| 国产精品亚洲一级av第二区| 亚洲无线在线观看| 国产成人啪精品午夜网站| 久久久久久九九精品二区国产 | 色综合站精品国产| 精品午夜福利视频在线观看一区| 国产午夜福利久久久久久| 亚洲性夜色夜夜综合| 欧美成人免费av一区二区三区| 好看av亚洲va欧美ⅴa在| 他把我摸到了高潮在线观看| 中亚洲国语对白在线视频| 亚洲欧洲精品一区二区精品久久久| av免费在线观看网站| av欧美777| 一边摸一边抽搐一进一小说| 伦理电影免费视频| 亚洲国产欧美日韩在线播放| 国产av不卡久久| 午夜成年电影在线免费观看| 色婷婷久久久亚洲欧美| 国产免费av片在线观看野外av| 久久香蕉激情| 可以在线观看毛片的网站| 久久精品人妻少妇| 亚洲av日韩精品久久久久久密| 桃色一区二区三区在线观看| 亚洲国产日韩欧美精品在线观看 | 天天躁夜夜躁狠狠躁躁| 91成年电影在线观看| 精品无人区乱码1区二区| 欧美乱码精品一区二区三区| 精华霜和精华液先用哪个| 99久久无色码亚洲精品果冻| 级片在线观看| 999精品在线视频| 欧美成人性av电影在线观看| 日本在线视频免费播放| 动漫黄色视频在线观看| 黄色视频不卡| xxxwww97欧美| 亚洲av五月六月丁香网| 黄色女人牲交| 国产极品粉嫩免费观看在线| 在线观看舔阴道视频| 日韩欧美一区二区三区在线观看| 黄片大片在线免费观看| 精品国产乱子伦一区二区三区| 久久久久亚洲av毛片大全| 欧美日韩亚洲综合一区二区三区_| 法律面前人人平等表现在哪些方面| 变态另类成人亚洲欧美熟女| 欧美zozozo另类| 精品国产国语对白av| 丰满的人妻完整版| 成年人黄色毛片网站| 国产成人精品久久二区二区免费| 少妇被粗大的猛进出69影院| 青草久久国产| 好看av亚洲va欧美ⅴa在| 久热这里只有精品99| 男人的好看免费观看在线视频 | 伦理电影免费视频| 成人av一区二区三区在线看| 极品教师在线免费播放| 级片在线观看| 男人舔女人的私密视频| 亚洲aⅴ乱码一区二区在线播放 | √禁漫天堂资源中文www| 精品一区二区三区视频在线观看免费| videosex国产| 色综合欧美亚洲国产小说| 亚洲中文字幕日韩| 欧美日韩亚洲综合一区二区三区_| 亚洲avbb在线观看| 欧美又色又爽又黄视频| 一个人观看的视频www高清免费观看 | 亚洲免费av在线视频| 手机成人av网站| 男女床上黄色一级片免费看| 国产又黄又爽又无遮挡在线| 免费看a级黄色片| 中国美女看黄片| 国产伦在线观看视频一区| 久热这里只有精品99| 欧美日韩福利视频一区二区| 国产亚洲精品av在线| 国产一区在线观看成人免费| 国产成人系列免费观看| 一区二区三区精品91| 搡老熟女国产l中国老女人| 丰满的人妻完整版| 女人被狂操c到高潮| 一级黄色大片毛片| 露出奶头的视频| 国产高清videossex| 男人舔奶头视频| 好男人在线观看高清免费视频 | 99精品欧美一区二区三区四区| 免费看日本二区| 一区二区三区激情视频| 最近最新中文字幕大全免费视频| 淫妇啪啪啪对白视频| 免费看a级黄色片| 后天国语完整版免费观看| 9191精品国产免费久久| 天堂动漫精品| 婷婷丁香在线五月| 少妇粗大呻吟视频| 免费看十八禁软件| 免费观看精品视频网站| avwww免费| 老司机午夜十八禁免费视频| 在线观看免费视频日本深夜| 色播亚洲综合网| 丁香六月欧美| 亚洲激情在线av| 91麻豆av在线| 91字幕亚洲| a在线观看视频网站| 一本大道久久a久久精品| 国产欧美日韩精品亚洲av| 看免费av毛片| 国产高清激情床上av| 欧美日韩一级在线毛片| 一级黄色大片毛片| 老鸭窝网址在线观看| 亚洲av片天天在线观看| 国产三级黄色录像| 精品国内亚洲2022精品成人| 亚洲久久久国产精品| 伦理电影免费视频| 久久亚洲精品不卡| 精品久久久久久成人av| 男女下面进入的视频免费午夜 | 欧美日本视频| 国产精品爽爽va在线观看网站 | 亚洲中文字幕日韩| 丝袜人妻中文字幕| 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说| 国产精品亚洲一级av第二区| 亚洲精品久久国产高清桃花| 午夜两性在线视频| 欧美乱妇无乱码| 亚洲一区二区三区色噜噜| 无人区码免费观看不卡| 久久久久亚洲av毛片大全| 老司机午夜十八禁免费视频| 久久久久亚洲av毛片大全| bbb黄色大片| 免费看a级黄色片| 亚洲欧美精品综合久久99| 色在线成人网| 色av中文字幕| 精品乱码久久久久久99久播| 久久 成人 亚洲| 听说在线观看完整版免费高清| 91在线观看av| 中文字幕另类日韩欧美亚洲嫩草| 日韩国内少妇激情av| 亚洲精华国产精华精| 在线观看免费视频日本深夜| 黄片大片在线免费观看| 久久婷婷成人综合色麻豆| 久9热在线精品视频| 悠悠久久av| 黄色片一级片一级黄色片| 色哟哟哟哟哟哟| 精品欧美国产一区二区三| 欧美成人午夜精品| 母亲3免费完整高清在线观看| 国产精品1区2区在线观看.| 日韩大码丰满熟妇| 又紧又爽又黄一区二区| 两性夫妻黄色片| 一进一出抽搐gif免费好疼| 黄色片一级片一级黄色片| 国产一级毛片七仙女欲春2 | 久久精品亚洲精品国产色婷小说| 久久中文字幕一级| 欧美绝顶高潮抽搐喷水| 亚洲av第一区精品v没综合| 欧美国产精品va在线观看不卡| 两个人视频免费观看高清| 老熟妇乱子伦视频在线观看| 亚洲精华国产精华精| 成人特级黄色片久久久久久久| 成人手机av| 国产精品综合久久久久久久免费| 欧美黑人巨大hd| 午夜两性在线视频| 一本一本综合久久| 国内少妇人妻偷人精品xxx网站 | 精品电影一区二区在线|