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

    Impacts of reduced wind speed on physiology and ecosystem carbon flux of a semi-arid steppe ecosystem

    2014-03-31 00:28:36DongYanJinQiongGaoYaLinWangLiXu
    Sciences in Cold and Arid Regions 2014年6期

    DongYan Jin, Qiong Gao, YaLin Wang, Li Xu

    1. Academy of Disaster Reduction and Emergency Management, Ministry of Civil Affairs & Ministry of Education, Beijing Normal University, Beijing 100875, China

    2. College of Natural Sciences, Beijing Normal University, Beijing 100875, China

    Impacts of reduced wind speed on physiology and ecosystem carbon flux of a semi-arid steppe ecosystem

    DongYan Jin1, Qiong Gao2*, YaLin Wang1, Li Xu1

    1. Academy of Disaster Reduction and Emergency Management, Ministry of Civil Affairs & Ministry of Education, Beijing Normal University, Beijing 100875, China

    2. College of Natural Sciences, Beijing Normal University, Beijing 100875, China

    Decreasing wind speed is one aspect of global climate change as well as global warming, and has become a new research orientation in recent decades. The decrease is especially evident in places with frequent perennially high wind speeds. We simulated decreased wind speed by using a steel-sheet wind shield in a temperate grassland in Inner Mongolia to examine the changes in physical environmental variables, as well as their impacts on the photosynthesis of grass leaves and net ecosystem exchange (NEE). We then used models to calculate the variation of boundary layer conductance (BLC) and its impact on leaf photosynthesis, and this allowed us to separate the direct effects of wind speed reduction on leaf photosynthesis (BLC) from the indirect ones (via soil moisture balance). The results showed that reduced wind speed primarily resulted in higher moisture and temperature in soil, and indirectly affected net assimilation and water use efficiency of the prevalent bunch grassStipa krylovii. Moreover, the wind-sheltered plant community had a stronger ability to sequester carbon than did the wind-exposed community during the growing season.

    wind speed; reduction; boundary layer conductance; photosynthesis; net ecosystem exchange; semi-arid steppe

    1 Introduction

    Wind plays an important part in the resource dynamics of semi-arid grasslands (Sterk, 2003; Munsonet al., 2011; Raviet al., 2011; Mooreet al., 2013). Much research has shown that the trend of wind speed is decreasing worldwide, and evidences have been reported from the United States (Pryoret al., 2009; Pryor and Ledolter, 2010), Canada (Tuller, 2004), Australia (McVicaret al., 2008), Italy (Pirazzoli, 2003), and Czechoslovakia (Brazdilet al., 2009). The same trend was also found in China, where annual, spring, and summer mean wind speeds have declined by 0.018, 0.021, and 0.015 m/(s·a), respectively. More wind speed reduction was found in the north, the Tibetan Plateau, and in coastal areas in the east and southeast parts of China (Guoet al., 2011).

    Several causes have been suggested to explain the wind speed reduction. The major hypotheses are changed patterns of atmospheric circulation at high altitudes caused by climate change, increased surface roughness caused by thriving vegetation, and increased urban density and numbers of tall buildings (Iacono, 2009; Joseph, 2010). Overgrown former agricultural land, afforestation, and the change of landscape management practices could account for 25%?60% of the observed stilling in wind speed. Reduced wind speed definitely impacts ecosystem processes, functions, andstructures, especially in cold or arid areas with excessive evapotranspiration induced by extreme wind speed. Reduced wind speed increases the boundary layer thickness of the Earth’s surface and slows down the exchange of energy and water between the atmosphere and biosphere (McVicaret al., 2012). Decreased evapotranspiration and heat exchange rate are conducive to maintaining surface soil moisture and increasing the canopy temperature (Cleugh, 1998, 2002).

    Previous studies on the effects of reduced wind speed were mostly focused on agriculture (Cleugh, 1998; BrandIeet al., 2004) and desert systems (Banget al., 2010). Much less is known about the impacts of reduced wind speed on environments, plant photosynthesis, and carbon sequestration of temperate grasslands. Typical steppe, a particular type of temperate grassland, is widely distributed in the eastern Eurasian steppe zone, and occupies a total area of 4.1×107ha, more than half of which located in China (Haoet al., 2008). Wind-induced excessive surface evapotranspiration intensifies water stress on the physiology and growth of plants, and extreme wind speed also causes soil erosion and land degradation (Raviet al., 2011). Thus, decreasing wind speed is considered to favor ecosystem production and resources/land conservation.

    In this paper, we describe a wind shield experiment conducted in a temperate grassland in Inner Mongolia. Wind speed, soil temperature and moisture, plant photosynthesis, and net ecosystem exchanges were measured on both sides of the windbreak to test the hypothesis that reducing wind speed favors grass physiology and community carbon sequestration, and to explore the mechanisms of the direct and indirect impacts of wind speed reduction on ecosystem processes.

    2 Materials and methods

    2.1 Experiment design

    The study was conducted in Taips Banner (114°51′E ?115°49′E, 41°35′N ?42°10′N), an administrative division occupying most of southern half of the Xilingol League of Inner Mongolia. The site is administered by Beijing Normal University. Due to the control of Mongolian high atmospheric pressure, the predominant wind blows from the northwest at an annual average speed of 3.41 m/s. The highest wind speeds occur in April, reaching 10 m/s and above. According to the latest 35 years of meteorological data, the annual, July, and January mean temperatures are 2.1, 34, and ?18 °C, respectively. The length of the frost-free season is about 100 days. The annual mean precipitation is 387 mm, with about 64.6% occurring from June to August. The soils are classified as chestnut soil, with low soil organic material content. The temperate grassland is dominated byStipa krylovii Roshev(Zhuet al., 2007).

    In the summer of 2010, the Beijing Normal University landscape simulation and analysis laboratory established an experimental plot for wind field manipulation. The plot is 160 m long and 50 m wide, with the width aligned to the northwest, the prevailing wind direction. To reduce the wind speed, a 150-m-long and 3-m-high steel-sheet wind shield was erected in the middle of the plot, with its normal pointing to the prevailing wind direction of 290°49′41.2′ (north equal to 0 or 360 deg, clockwise positive). The wind shield divided that plot into two identical parts, each of which was 150m×25m. The northwest side of the wind shield was exposed to the prevailing wind and was used as the control (EX); the other side was used as the sheltered treatment (SH). Areas within 5 m of the wind shield were left as buffer zones on both sides, and the rest of the 20m×150m at each side of the wind shield was divided into 2 rows and 15 columns, thus forming 30 subplots 10m×10m on each side. One row of subplots facing the major wind on either side was used for clipping defoliation treatments, but those are not described in this paper. This paper addresses only the unclipped rows, leaving the foliage clipping experiment for later analysis.

    2.2 Experimental measurements

    Soil water contents and temperature of 20 randomly selected subplots were sampled every 30 min at a depth of 10 cm with EM50 devices (ECH2O EC-TM, Decagon Devices, Inc., Pullman, WA), so that each treatment had 10 replicates. We measured the leaf photosynthesis rate of the dominant speciesStipa krylovii, a perennial bunch grass, on June 30, July 23, August 31, and September 23, which spanned the whole growth season of 2012. For each of the four months, diurnal stomatal conductance, leaf transpiration, and net assimilation rates, as well as microclimate variables (e.g., leaf temperature, incident photon flux density, and vapor pressure), were measured with six replicates on both sides of wind shield at every 2 h from 08:00 to 18:00, with a portable photosynthesis system (model LI-6400, Li-Cor, Inc., Lincoln, NE).

    The diurnal dynamics of the net ecosystem exchange (NEE) were also measured continuously for 24 hours in the same plots of the leaf photosynthesis measurement from May to September, at half-hour intervals by a fully automatic closed dynamic chamber (CDC) system (model LI-8150, Li-Cor, Inc.).

    2.3 Model descriptions

    The stomatal conductance model we used was a simplified version by Gaoet al.(2005), which calcu-lates stomatal conductance as a function of light intensity, soil water potential, and normalized vapor deficit in air. The leaf photosynthesis models for C3 plants were formulated by Thornley and Johnson (1990) (see Appendix A), which describe a process of carboxylase activated by light, then reacting with carbon dioxide to produce carbohydrates (carboxylation), or reacting with oxygen to produce carbon dioxide (photorespiration) (Gaoet al., 2004). To evaluate the impact of wind speed reduction on grass leaf physiology, we coupled a wind-speed-governed boundary layer conductance model with the stomatal conductance model. This enabled us to translate the effects of wind speed reduction into a calculation of photosynthesis. The details of the model descriptions are given in Appendix B. Model parameters were evaluated from the observed data by means of the nonlinear least-square technique (Table 1).

    Table 1 List of model parameter values evaluated by the nonlinear least-square technique

    3 Results

    3.1 Effects of wind shield on wind velocity

    Wind speed behind the wind shield was reduced significantly, and the closer to the wind shield the greater wind speed reduction (Figure 1). In particular, at 5 m away from the middle of the wind shield, the wind speed of the sheltered side was decreased on average by 58%, 29%, and 48% for April, August, and the entire growth season, respectively. The corresponding average wind speed reductions at 15 m away from the middle of the windbreak were 40%, 2%, and 30%, respectively. These differences were statistically significant except the ones in August, when wind directions diverted most from the prevailing direction.

    3.2 Soil moisture and temperature

    The amount of total precipitation in the growing season (from May to September) was 344 mm. During this period, the maximum and minimum monthly precipitation rates were 172 mm and 16 mm in July and August, respectively (Figure 1). The monthly variation of topsoil moisture was closely related to the amount of precipitation (Table 2).

    Soil moisture in the shelter area was significantly greater than that on the windward side throughout the whole growing season (P<0.001) (Table 2). The mean value on the leeward side, from June to September, was 0.127±0.037 cm3/cm3, 14.4% greater than that on the windward side. Moreover, the maximum difference of soil moisture between sheltered and exposed, 18.6%, occurred in August, and the minimum of 8.8% occurred in September. Mean sheltered and exposed soil temperatures during the growing season of 2012 were 20.33±5.70 °C and 19.31±4.99 °C, respectively, so that the sheltered side was about 1 °C warmer than the exposed side (P<0.001). The higher soil moisture and temperature in the sheltered side implied decreased surface evaporation and sensible heat loss because of the reduced wind speed.

    Figure 1 Wind speed and precipitation on both sides of the windbreak. Averaged daily wind speed measured along the middle of the windbreak; four days were measured in each month

    Table 2 Topsoil temperature and moisture (0-20 cm deep) on both sides of the wind shield

    3.3 Observed leaf photosynthesis

    Figure 2 plots the averages of the observed diurnal net assimilation rate (AN), stomatal conductance (GS), transpiration rate (Tr), and water use efficiency (WUE) ofStipa kryloviiunder the sheltered and exposed treatments from June to September of 2012, respectively. The results show that shelteredANwas 45% greater than exposed, because of warmer and wetter soil. The mean values of sheltered and exposedANwere 7.76±4.18 and 5.35±3.96 μmol/(m2·s), respectively.

    The differences ofGSbetween the leeward and windward were insignificant, except those at 08:00, 14:00, and 18:00. The exposedGSwas about 14% and 20% higher than sheltered in the morning and nightfall, respectively, but was 13% lower at 14:00. The diurnal profiles of the leaf transpiration on both sides of the shield were similar, but the overall average of transpiration in the sheltered side was 3% lower that of the exposed side.

    Figure 2 Observed diurnal dynamics of photosynthesis, stomatal conductance, transpiration, and water use efficiency ofStipa kryloviion both leeward and windward sides

    We calculated intrinsicWUEleafin the form ofANdivided byGS(Waterhouseet al., 2004), which presents substantial intra- and interspecific variation and obtains a better description of plant characteristics than the conventionalWUE. TheWUEleafat the leeward side was obviously greater than that of the windward side during the entire growing season, mostly due to the difference inAN. The average dailyWUEleafin the leeward was 34.5±17.6 μmolCO2/molH2O, 40% higher than that of the windward.

    3.4 Effects of wind speed on leaf boundary layer conductance and gas exchanges

    Predicted leaf boundary layer conductance (BLCleaf) was calculated for the control, as well as at 5 m and 15 m away from the wind shield (leeward) based on observed climate data. For an average surface roughness of 3.7×10?3m (Huet al., 2006), we calculated theBLCleafwith 0.1 m as the characteristic length, using the wind speed on both sides of the shield. The details of the model descriptions are given in Appendix C.

    The calculated meanBLCleafat the sheltered side decreased by 25% and 17% at 5 m and 15 m away from the windbreak, for a given wind speed reduction of 48% and 30%, respectively (Figure 3a). Moreover, the total conductance (Gt), such asGSin series withBLCleaf, varied with wind speed and, more evidently, with moisture conditions in air (vapor pressure deficit,VPD) and soil (soil water potential,Si).

    TheGtis primarily sensitive to wind speed when the latter is smaller than a particular value. In other words,Gtasymptotically approaches a constant value when the wind speed is beyond this range. The range decreases quickly with moisture in soil and air. When the soil water potential varied from ?0.033 (the field capacity) to ?1 MPa, and the vapor pressure deficit increased from 0.5 to 2 kPa, the asymptoticGtdecreased from 1.24 to 0.12 molH2O/(m2·s) at 4 m/s wind speed. However, when wind speed gradually decreased,Gtunder all of the water conditions dropped quickly or slowly. For instance,Gtunder 0.2 kPa vapor pressure deficit and ?0.033 MPa soil water potential decreased when wind speed reduced from 4 m/s to 1 m/s, and sharply deceased when the speed was less than 1 m/s. However,Gtremained constant under 0.5 kPa vapor pressure deficit and ?0.5 MPa soil water potential when the wind speed declined from 4 to 1 m/s, and then decreased quickly when the wind speed was less than 1 m/s. The behavior ofGtwas thus primarily controlled by the stomatal conductance, and the boundary layer conductance came into play when the wind speed was small; the drier the water condition, the smaller the value of wind speed at whichGtsharply decreased.

    We also calculatedANas a function of wind speed, and found that the direct effect of wind speed reduction onAN(i.e., the boundary layer change) was negligible for wind speed greater than 0.5 m/s.

    Figure 3 Calculated sheltered and exposedBLCleafbased on observed wind speed from June to September (a), and the variation ofGtwith wind speed under different vapor pressure deficits (VPD) and soil water potentials (Si) (b). Leaf temperature, incident photon flux density, and pressure were fixed to 20 °C, 1,000 μmol/(m2·s), and 101 kPa, respectively, in theBLCleafcalculation

    3.5 Net ecosystem exchange

    Because of the reduced transpiration and increased carbon assimilation, theNEE(negative for downward flux, or carbon input to the ecosystem) at the leeward side was ?1.05±0.59, ?0.35±2.26, ?1.87±3.57,?0.53±1.03, and ?0.92±3.28 μmol/(m2·s) lower than the windward side for May, June, July, August, and September, respectively. We summarized the negative (carbon input or sequestration) and positive (carbonloss from ecosystem) values ofNEEseparately for the five months and then converted their units into gC/(m2·d) (Figure 4). Paired-t tests (two-tailed) were used to test the mean difference between sheltered and exposedNEE. The results showed that carbon emission of the sheltered side was 0.35±0.40 gC/(m2·d), or 24%, less than the exposed, but the carbon sequestration of the sheltered side was 0.64±0.37 gC/(m2·d), or 39%, greater than the exposed from May to September (Figure 5). The maximum and minimum values of carbon emission were 3.08 and 0.21 gC/(m2·d) in the exposed, and 2.04 and 0.18 gC/(m2·d) with sheltered treatments, respectively, which implied the rates of heterotrophic and autotrophic respiration of the ecosystem in the growth season. Moreover, the maximum and minimum values of carbon sequestration were 1.66 and 0.18 gC/(m2·d) (exposed) and 2.32, 0.48 gC/(m2·d) (sheltered) in August, respectively. On average, the shelteredNEEsequestered 0.98 gC/(m2·d) more than the exposedNEEfrom May to September in 2012.

    It is important to note that the measurements ofNEE, as well as leaf photosynthesis, on both sides of windbreak were conducted in the chamber with the same air velocity, and the different performances of plants between sheltered and exposed were mainly due to the variation of ambient conditions such as soil moisture and temperature. Therefore, we consider the differences inANandNEEbetween the two sides of the windbreak to mainly result from the indirect effects of reduced wind speed on the ecosystem.

    Figure 4 Diurnal dynamics of net ecosystem CO2exchange (NEE) under the sheltered and exposed treatments fromMay to September in 2012. Empty circles indicateNEEin exposed, and solid black circles indicateNEEin sheltered

    Figure 5 Sum of dailyNEE, carbon sequestration, and emission on the both sides of the shield from May to September

    4 Discussion

    4.1 Direct effects of reduced wind speed

    As a result of wind speed reduction and changes in turbulent transfer rates, the microclimate in the sheltered zone was modified (Cleugh, 2002; Brandleet al., 2004). In this study, the sheltered topsoil water content tended to be 15% wetter than in the open, and topsoil temperature tended to be about 1 °C warmer than in the exposed. According to Zhang and Brandle (1996), the combination of increased moisture and temperature in topsoil has a positive effect on plant physiology and growth in semi-arid grasslands, and our findings confirmed this.

    ReducedBLCleafresulting from decreased wind velocity contributes to the decoupling of stomatal exchange with ambient conditions and lessening the relationship between changes in stomatal conductance and the rate of transpiration (Schuepp, 1993). Our results indicated that the magnitude of shelteredBLCleafwithin 5H(whereHis the height of the wind shield) was reduced by 21% when the wind speed decreased by 39%. Also, the calculatedGtvaried with wind speed and water conditions; we found thatGtincreased with wind speed only at wind speeds in the range between 0 and 0.5 m/s in moderately dry conditions. However, the control of wind speed on the total conductance was stronger when moisture stresses from soil and air were low, and the range of the wind speed increased to between 0 and 2.0 m/s when the soil water potential and vapor pressure deficit were?0.5 MPa and 0.5 kPa, respectively, and even increased to 4 m/s when the soil water potential and vapor pressure deficit were ?0.033 MPa and 0.2 kPa, respectively.

    In addition, we found that the direct influences of reduced wind speed onANwere negligibly small unless the wind speed was very close to zero. Previous studies on the impact of wind speed on net photosynthesis rate and transpiration rate showed similar results (Daudetet al., 1999). For example, Kitayaet al.(2003) found that the net photosynthetic rate and the transpiration rate increased sharply as the air current speeds increased from 0.01 to 0.2 m/s, and then slowed down gradually in the range from 0.2 to 1.0 m/s. In contrast to transpiration, the net photosynthetic rate was almost constant at air current speeds ranging from 0.5 to 1.0 m/s. Thus, the influence of wind speed on the plant photosynthesis process could be ignored unless it was below 0.2 m/s.

    4.2 Indirect impact of reduced wind speed

    Reduced wind speed indirectly affects plant physiological activities by changing the ambient microclimate (Brandleet al., 2004). The sheltered meanANofStipa kryloviiin this study was about 45% more than the exposed in 2012 because of the wetter and warmer topsoil in the leeward side. However, there were slight differences in measuredGSandTrbetween sheltered and exposed. One thing to note here is that the diurnal variation of the photosynthesis process was averaged from June to September. The shelteredGSwas significantly greater than the exposedGSin June and July, but was converted by the abnormal rainfall in September. However, greater exposedGSdid not represent greaterAN, especially in the end of the growing season, probably because of the non-stomatal limitation.

    WUEleafis described as the net carbon uptake per amount of water lost from the ecosystem (Bacelaret al., 2012). The intrinsicWUE, defined as net carbon assimilation divided by stomatal conductance in this experiment (Figure 2), showed the advantage of the reduced wind speed. The result of higherWUEleafon the leeward side than the windward side was also consistent with previous studies. For example, Davis and Norman (1988) suggested that under some conditions, sheltered plants made more efficient use of available water. However, our result differs from Monteith (1981), who suggested that water use efficiency in shelter was unlikely to increase unless there was a significant decrease in the saturation vapor pressure deficit.

    Changes of soil moisture and temperature caused by the reduced wind speed also affected ecosystem carbon sequestration and emission. Our experiment results show that carbon sequestration in the leeward plots was 0.64±0.37 gC/(m2·d) more than in the exposed plots. The higher carbon assimilation rate and lower plant respiration rate in the sheltered region was caused by warmer temperature than in the exposed region in daytime but cooler at night (BrandIeet al., 2004). On the other hand, Reyet al. (2012) suggested that the variable pattern and magnitude of the relationship between net ecosystem carbon balance (NECB) and wind speed under different atmospheric stability conditions implies the role of convective turbulence. If a stable atmospheric stratification precludes CO2effluxes, unstable atmospheric conditions, characterized by large sensible heat fluxes, are more efficient than neutral conditions at a given wind speed. Hence, the ecosystem carbon exchange process is not only affected by wind speed but is also determined by atmospheric stability conditions on some level, which needs further study in the future.

    5 Conclusion

    This paper presents the direct and indirect impacts of windbreak-reduced wind speed on soil moistureand temperature, plant photosynthesis, and carbon sequestration in a semi-arid ecosystem during the growing season of 2012. As a result of wind speed reduction and changes in the turbulent transfer rate, topsoil moisture and temperature and boundary layer resistance on the sheltered side were all increased. The lower boundary layer conductance, coupled with the stomatal conductance, brought down leaf transpiration, but this effect of reduced BLC depended on the moisture stresses from air and soil, and was negligible for wind speed greater than about 2.0 m/s. The net assimilation rate ofStipa kryloviiand community carbon sequestration were also indirectly affected by the reduced wind speed: Warmer and wetter soil caused by the reduced evapotranspiration and sensible heat flux enhanced photosynthesis and carbon sequestration. Our results, limited to a single plot within one growth season, imply that the indirect impact of reduced wind speed on plant physiology and ecosystem carbon sequestration via soil moisture and temperature was much stronger than the direct effect with altered boundary layer conductance. In conclusion, reduced wind speed exhibited a positive effect on mitigation of water stress and it enhanced grass growth in a semi-arid steppe at a regional scale.

    Acknowledgments:

    This research was jointly supported by the National Science Foundation of China (Nos. 41171445 and 41321001), the State Key Laboratory of Earth Surface Processes and Resources Ecology (No. 2012-TDZY-31), and the National Program on Key Basic Research Project (No. 2014CB954303).

    Bacelar EL, Moutinho-Pereira JM, Gon?alves BM,et al., 2012. Water use strategies of plants under drought conditions. In: Plant Responses to Drought Stress. Springer, New York, pp. 145-170.

    Bang C, Sabo JL, Faeth SH, 2010. Reduced wind speed improves plant growth in a desert city. PloS One, 5(6): e11061.

    Brandle JR, Hodges L, Zhou XH, 2004. Windbreaks in North American agricultural systems. In: Nair PKR, Rao MR, Buck LE (eds.). New Vistas in Agroforestry. Springer: Dordrecht, The Netherlands, pp. 65-78.

    Brazdil R, Chroma K, Dobrovolny P,et al., 2009. Climate fluctuations in the Czech Republic during the period 1961-2005. International Journal of Climatology, 29(2): 223-242.

    Cleugh H, 1998. Effects of windbreaks on airflow, microclimates and crop yields. Agroforestry Systems, 41(1): 55-84.

    Cleugh H, 2002. Field measurements of windbreak effects on airflow, turbulent exchanges and microclimates. Animal Production Science, 42(6): 665-677.

    Daudet FA, Le Roux X, Sinoquet H,et al., 1999. Wind speed and leaf boundary layer conductance variation within tree crown: Consequences on leaf-to-atmosphere coupling and tree functions. Agricultural and Forest Meteorology, 97(3): 171-185.

    Davis J, Norman J, 1988. Effects of shelter on plant water use. Agriculture, Ecosystems & Environment, 22: 393-402.

    Gao Q, Yu M, Zhang XS,et al., 2005. Modelling seasonal and diurnal dynamics of stomatal conductance of plants in a semi-arid environment. Functional Plant Biology, 32(7): 583.

    Gao Q, Zhang XS, Huang YM,et al., 2004. A comparative analysis of four models of photosynthesis for 11 plant species in the Loess Plateau. Agricultural and Forest Meteorology, 126(3-4): 203-222.

    Guo H, Xu M, Hu Q, 2011. Changes in near-surface wind speed in China: 1969-2005. International Journal of Climatology, 31(3): 349-358.

    Hao Y, Wang Y, Mei X,et al., 2008. CO2, H2O and energy exchange of an Inner Mongolia steppe ecosystem during a dry and wet year. Acta Oecologica, 33(2): 133-143.

    Hu X, Liu LY, Yan P,et al., 2006. The effect of different soil surface on wind erosion in Taipusi County, Inner Mongolia. Research of Soil and Water Conservation, 13(4): 116-119.

    Iacono MJ, 2009. Why is the Wind Speed Decreasing? Miltion, MA: Blue Hill Meteorological Observatory.

    Joseph M, 2010. Why winds are slowing. [http://www.nature.com/ news/2010/101017/full/news.2010.543.html, accessed May 5, 2014]

    Kitaya Y, Tsuruyama J, Shibuya T,et al., 2003. Effects of air current speed on gas exchange in plant leaves and plant canopies. Advances in Space, 30(1): 177-182.

    McVicar TR, Roderick ML, Donohue RJ,et al., 2012. Global review and synthesis of trends in observed terrestrial near-surface wind speeds: Implications for evaporation. Journal of Hydrology, 416: 182-205.

    McVicar TR, Van Niel TG, Li LT,et al., 2008. Wind speed climatology and trends for Australia, 1975-2006: Capturing the stilling phenomenon and comparison with near-surface reanalysis output. Geophysical Research Letters, 35(20): L20403.

    Monteith J, 1981. Coupling of plants to the atmosphere. In: Grace J, Ford ED, Jarvis PG (eds.). Plants and Their Atmospheric Environment. Blackwell Scientific Publications, Oxford, pp. 1-29.

    Moore JR, Manley BR, Park D,et al., 2013. Quantification of wind damage to New Zealands planted forests. Forestry, 86(2): 173-183.

    Munson SM, Belnap J, Okin GS, 2011. Responses of wind erosion to climate-induced vegetation changes on the Colorado Plateau. Proceedings of the National Academy of Sciences of the United States of America, 108(10): 3854-3859.

    Pirazzoli PA, 2003. Recent near-surface wind changes in the central Mediterranean and Adriatic areas. International Journal of Climatology, 23(8): 963-973.

    Pryor SC, Barthelmie RJ, Young DT,et al., 2009. Wind speed trends over the contiguous United States. Journal of Geophysical Research, 114(D14): D14105.

    Pryor SC, Ledolter J, 2010. Addendum to "Wind speed trends over the contiguous United States.'' Journal of Geophysical Research: Atmosphere, 115(D10): D10103.

    Ravi S, D'Odorico P, Breshears DD,et al., 2011. Aeolian processes and the biosphere. Reviews of Geophysics, 49(3): RG3001.

    Rey A, Belelli-Marchesini L, Were A,et al., 2012. Wind as a main driver of the net ecosystem carbon balance of a semi-arid Mediterranean steppe in the southeast of Spain. Global Change Biology, 18: 539-554.

    Schuepp PH, 1993. Leaf boundary layers. New Phytologist, 125(3): 477-507.

    Sterk G, 2003. Causes, consequences and control of wind erosion in Sahelian Africa: A review. Land Degradation & Development, 14(1): 95-108.

    Thornley JH, Johnson IR, 1990. Plant and Crop Modelling. Clarendon Press, Oxford, UK.

    Tuller SE, 2004. Measured wind speed trends on the west coast of Canada. International Journal of Climatology, 24(11): 1359-1374.

    Waterhouse JS, Switsur VR, Barker AC,et al., 2004. Northern European trees show a progressively diminishing response to increasing atmospheric carbon dioxide concentrations. Quaternary Science Reviews, 23(7-8): 803-810.

    Zhang HH, Brandle JR, 1996. Windbreak effect on biomass and grain mass accumulation of corn: A modeling approach. Agronomy Journal, 88(4): 607-613.

    Zhu YJ, Gao Q, Liu JS,et al., 2007. Aggregation of plant functional types based on models of stomatal conductance and photosynthesis. Journal of Plant Ecology (Chinese version), 31(5): 873-882.

    Appendix A: Thornley and Johnson (1990) leaf model for C3 plants

    The model takes the following mathematical form:

    where:

    gscandgsoare stomatal conductance for CO2and O2, respectively;CaandOa, the partial CO2and O2pressure on leaf surface, respectively;gχandgr(μmol/(m2·s·kPa)), the carboxylation and photorespiration conductance, respectively. These parameters were considered dependent on temperature, so that

    where constant parametersα25,gχ25,gr25andRd25were determined by fitting the model to the field data;Tis the leaf temperature (°C);R= 8.314 kJ/(mol·K), the gas constant;Eaj= 37.0 kJ/mol, andEar= 66.4 kJ/mol, are the temperature coefficients.H= 220 kJ/mol, andS= 0.71 kJ/(mol·K). Constant 298 has units of K.

    Appendix B: Gaoet al. (2005) stomatal conductance model for C3 species

    The model by Gaoet al. (2005) define stomatal conductance a function of incident photon flux density, vapor pressure deficit in air, and soil water potential. The model takes the following form:

    where:

    Ip

    0

    ψs(MPa) is the soil water potential;Ip(mmol/(m2·s)) is the incident photon flux density;π0(MPa), the osmotic pressure at zero light intensity (dark);β((kPa·m2·s)/mol), the apparent elastic modulus of the guard cell structure;rz0((kPa·m2·s)/mol), soil-to-leaf hydraulic resistance;kβgis thus approximately the sensitivity of stomatal conductance to vapor pressure deficit (D).λ(dimensionless) is the constant parameter;πp(MPa), the maximum osmotic pressure inducible by incident light;C(kPa), CO2on leaf surface;C0(kPa), the reference CO2.kI(mmol/(m2·s)), the half-saturation light intensity. The parameters to be fitted areπ0,πp,β,rz0andkI.

    Appendix C: Boundary layer conductance model

    The boundary layer conductance model was based on the principle of hydromechanics and the Maxwell-Gilliland function. The expressions are:

    where:

    gb(molH2O/(m2·s)) is the boundary layer conductance;Ha(m) is the thickness of the boundary layer at a specific characteristic lengthL(m);DW(m2/s) is vapor diffusivity,R(m3·kPa)/(K·mol)) is the gas constant;MA,B(g/mol) is the molar mass of water and air;VA,Bis the molecular volume of water and air;μ(Pa·s) is the dynamic viscosity; ρ (kg/m3) is fluid density;n(dimensionless) is a constant fitted by the model; andu(m/s) is the wind velocity.

    : Jin DY, Gao Q, Wang YL,et al., 2014. Impacts of reduced wind speed on physiology and ecosystem carbon flux of a semi-arid steppe ecosystem. Sciences in Cold and Arid Regions, 6(6): 0556-0565.

    10.3724/SP.J.1226.2014.00556.

    Received: February 20, 2014 Accepted: May 26, 2014

    *Correspondence to: Dr. Qiong Gao, Professor of College of Natural Sciences, Beijing Normal University, 19 Xinjiekouwai Avenue, Beijing 100875, China. Tel: +86-10-58806050; Email: gaoq@bnu.edu.cn

    免费黄色在线免费观看| av线在线观看网站| 久久av网站| 欧美xxⅹ黑人| 99热网站在线观看| 人妻人人澡人人爽人人| 亚洲国产毛片av蜜桃av| 黄色一级大片看看| 欧美bdsm另类| 色播在线永久视频| 精品亚洲成a人片在线观看| 熟女av电影| 老鸭窝网址在线观看| 各种免费的搞黄视频| 一区二区三区四区激情视频| 999久久久国产精品视频| 高清视频免费观看一区二区| 亚洲精品,欧美精品| √禁漫天堂资源中文www| 美国免费a级毛片| 亚洲美女搞黄在线观看| 欧美日韩综合久久久久久| 99香蕉大伊视频| 亚洲情色 制服丝袜| 97精品久久久久久久久久精品| 日韩视频在线欧美| 欧美精品一区二区大全| 欧美另类一区| 黄色视频在线播放观看不卡| 免费看不卡的av| 久久女婷五月综合色啪小说| 久久精品国产亚洲av高清一级| 日韩一本色道免费dvd| 中文欧美无线码| 国产麻豆69| 欧美精品国产亚洲| 青春草亚洲视频在线观看| 美女福利国产在线| 天天躁日日躁夜夜躁夜夜| 777久久人妻少妇嫩草av网站| 久久精品久久久久久久性| 国产无遮挡羞羞视频在线观看| 在线观看人妻少妇| 亚洲精品aⅴ在线观看| 日韩中文字幕视频在线看片| 欧美日韩av久久| 欧美97在线视频| 一级毛片黄色毛片免费观看视频| 啦啦啦在线观看免费高清www| 热99国产精品久久久久久7| 丝袜喷水一区| 欧美精品人与动牲交sv欧美| 亚洲情色 制服丝袜| 久久久久久免费高清国产稀缺| 亚洲视频免费观看视频| av免费观看日本| 校园人妻丝袜中文字幕| 久久久久网色| 天堂俺去俺来也www色官网| 国产一级毛片在线| 国产一区二区激情短视频 | 蜜桃国产av成人99| 日韩熟女老妇一区二区性免费视频| 亚洲精品,欧美精品| 一级毛片 在线播放| 亚洲国产精品一区三区| 女人被躁到高潮嗷嗷叫费观| 久久国产精品大桥未久av| 精品一区二区三区四区五区乱码 | 日韩中字成人| 日韩免费高清中文字幕av| av在线老鸭窝| 大话2 男鬼变身卡| 午夜激情久久久久久久| 麻豆精品久久久久久蜜桃| 熟女av电影| 777久久人妻少妇嫩草av网站| 老女人水多毛片| 久久人妻熟女aⅴ| 免费不卡的大黄色大毛片视频在线观看| 最近手机中文字幕大全| 国产淫语在线视频| 青草久久国产| 男男h啪啪无遮挡| 免费黄网站久久成人精品| 亚洲成色77777| 性少妇av在线| 久久久久精品性色| 国产成人一区二区在线| 国产精品秋霞免费鲁丝片| 亚洲国产欧美网| 丁香六月天网| a级毛片在线看网站| 欧美少妇被猛烈插入视频| 亚洲五月色婷婷综合| 18+在线观看网站| 一区二区三区精品91| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 天天操日日干夜夜撸| 夫妻午夜视频| 婷婷成人精品国产| 在线天堂最新版资源| 大码成人一级视频| 日韩中文字幕欧美一区二区 | 日本欧美视频一区| 大香蕉久久成人网| 好男人视频免费观看在线| 国产野战对白在线观看| 91午夜精品亚洲一区二区三区| 两个人看的免费小视频| 国产在线一区二区三区精| 日本免费在线观看一区| 搡老乐熟女国产| 热re99久久精品国产66热6| 精品少妇黑人巨大在线播放| 久久这里只有精品19| 久久国内精品自在自线图片| 亚洲男人天堂网一区| 五月伊人婷婷丁香| 国产成人a∨麻豆精品| 涩涩av久久男人的天堂| 日韩av在线免费看完整版不卡| 成人18禁高潮啪啪吃奶动态图| 亚洲成人一二三区av| 欧美激情极品国产一区二区三区| 桃花免费在线播放| 欧美成人午夜精品| 亚洲av电影在线进入| 一本—道久久a久久精品蜜桃钙片| 亚洲精品日本国产第一区| 午夜免费男女啪啪视频观看| 国产精品秋霞免费鲁丝片| 亚洲欧美清纯卡通| a级片在线免费高清观看视频| 欧美精品一区二区免费开放| 欧美精品av麻豆av| 久久久久网色| www.自偷自拍.com| 尾随美女入室| 中文字幕人妻丝袜制服| 亚洲欧美成人综合另类久久久| 伦理电影免费视频| 欧美国产精品va在线观看不卡| 美女视频免费永久观看网站| 秋霞伦理黄片| 国产有黄有色有爽视频| 少妇被粗大猛烈的视频| 久久午夜综合久久蜜桃| 女性被躁到高潮视频| 男女午夜视频在线观看| 精品少妇黑人巨大在线播放| 成人黄色视频免费在线看| 黑人欧美特级aaaaaa片| av在线app专区| 高清黄色对白视频在线免费看| 在线观看国产h片| 国产乱来视频区| 亚洲精品美女久久av网站| 色吧在线观看| 国产成人aa在线观看| 亚洲在久久综合| 亚洲av成人精品一二三区| 色婷婷久久久亚洲欧美| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 三级国产精品片| 国产在线免费精品| 国产精品欧美亚洲77777| 免费看av在线观看网站| 亚洲欧美一区二区三区黑人 | 国产精品.久久久| 成人国语在线视频| 国产精品久久久久久精品古装| 大陆偷拍与自拍| 国产亚洲午夜精品一区二区久久| 久久精品国产a三级三级三级| 交换朋友夫妻互换小说| 国产成人a∨麻豆精品| 亚洲三级黄色毛片| 国产精品秋霞免费鲁丝片| 亚洲欧美中文字幕日韩二区| 午夜福利在线免费观看网站| 两个人看的免费小视频| 亚洲少妇的诱惑av| 国产精品一国产av| 在线 av 中文字幕| 看十八女毛片水多多多| 亚洲欧洲精品一区二区精品久久久 | 成年女人毛片免费观看观看9 | 中国国产av一级| 熟妇人妻不卡中文字幕| 国产精品国产三级专区第一集| 大片免费播放器 马上看| 精品一品国产午夜福利视频| 少妇人妻久久综合中文| 国产精品99久久99久久久不卡 | 日韩中字成人| 男人爽女人下面视频在线观看| 欧美+日韩+精品| 日日撸夜夜添| 蜜桃国产av成人99| 男女边摸边吃奶| 亚洲成人手机| 免费在线观看完整版高清| 亚洲综合色惰| 国产日韩欧美视频二区| 免费观看在线日韩| 免费观看性生交大片5| 午夜影院在线不卡| 日韩不卡一区二区三区视频在线| 在线观看www视频免费| 日本免费在线观看一区| 91精品伊人久久大香线蕉| 国产av一区二区精品久久| 精品国产乱码久久久久久小说| 美女国产高潮福利片在线看| 最近中文字幕2019免费版| 精品一区在线观看国产| 国产亚洲精品第一综合不卡| 国产精品无大码| 热99国产精品久久久久久7| 亚洲精品美女久久久久99蜜臀 | 五月伊人婷婷丁香| 黄片播放在线免费| 麻豆精品久久久久久蜜桃| 亚洲,欧美精品.| 日韩av不卡免费在线播放| 男女高潮啪啪啪动态图| 久久久久网色| 日韩精品免费视频一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲一级一片aⅴ在线观看| 国产在视频线精品| 丁香六月天网| 国产熟女欧美一区二区| 亚洲av.av天堂| 91精品国产国语对白视频| 伊人久久大香线蕉亚洲五| 国产白丝娇喘喷水9色精品| 久久这里只有精品19| 欧美日韩国产mv在线观看视频| 免费黄频网站在线观看国产| 天堂中文最新版在线下载| 国产视频首页在线观看| 在线观看一区二区三区激情| 美女国产视频在线观看| 另类精品久久| 18在线观看网站| 一本久久精品| 波野结衣二区三区在线| 一区二区三区四区激情视频| 国产成人一区二区在线| 国产视频首页在线观看| h视频一区二区三区| 又大又黄又爽视频免费| 婷婷色av中文字幕| 汤姆久久久久久久影院中文字幕| 国产精品99久久99久久久不卡 | 一级毛片 在线播放| 婷婷色av中文字幕| 国产极品天堂在线| 欧美日韩精品网址| 人体艺术视频欧美日本| 色播在线永久视频| 女人久久www免费人成看片| av天堂久久9| videos熟女内射| 制服诱惑二区| 国产日韩欧美视频二区| 美女大奶头黄色视频| 午夜91福利影院| 伦精品一区二区三区| 国产精品免费大片| 极品少妇高潮喷水抽搐| 亚洲国产av新网站| 久久精品久久久久久久性| 精品福利永久在线观看| 欧美 亚洲 国产 日韩一| 国产精品.久久久| 成人漫画全彩无遮挡| 国产色婷婷99| 成年女人毛片免费观看观看9 | 久久韩国三级中文字幕| 国产在线视频一区二区| 丰满乱子伦码专区| 午夜福利乱码中文字幕| 亚洲av福利一区| 九九爱精品视频在线观看| 久久99精品国语久久久| 一级a爱视频在线免费观看| 国产av码专区亚洲av| 国精品久久久久久国模美| 亚洲精华国产精华液的使用体验| 国产精品欧美亚洲77777| 波多野结衣av一区二区av| 成人国产麻豆网| 久久久久久免费高清国产稀缺| 老鸭窝网址在线观看| 高清欧美精品videossex| av有码第一页| 青青草视频在线视频观看| 赤兔流量卡办理| 亚洲av福利一区| 精品人妻一区二区三区麻豆| 日本vs欧美在线观看视频| 晚上一个人看的免费电影| 亚洲av中文av极速乱| 不卡av一区二区三区| 亚洲经典国产精华液单| 久久国产亚洲av麻豆专区| 亚洲在久久综合| 国产精品一区二区在线不卡| 久久婷婷青草| 亚洲欧洲日产国产| 亚洲精品久久成人aⅴ小说| 女人精品久久久久毛片| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 久久久国产欧美日韩av| 亚洲av.av天堂| 熟妇人妻不卡中文字幕| 搡女人真爽免费视频火全软件| 婷婷色av中文字幕| 一区在线观看完整版| 国产亚洲一区二区精品| 欧美精品一区二区免费开放| 亚洲成人一二三区av| 免费高清在线观看视频在线观看| 一边摸一边做爽爽视频免费| 婷婷色麻豆天堂久久| 久久久久网色| 男的添女的下面高潮视频| 老司机亚洲免费影院| 99精国产麻豆久久婷婷| 午夜老司机福利剧场| 90打野战视频偷拍视频| 亚洲婷婷狠狠爱综合网| 99国产综合亚洲精品| 国产精品久久久久久av不卡| 久久久久国产网址| 少妇猛男粗大的猛烈进出视频| 99国产精品免费福利视频| 久久ye,这里只有精品| 人妻 亚洲 视频| 午夜久久久在线观看| 国产又爽黄色视频| 国产精品免费大片| 另类亚洲欧美激情| 欧美精品av麻豆av| 欧美黄色片欧美黄色片| 久久久久久久精品精品| 亚洲国产精品999| 亚洲欧洲精品一区二区精品久久久 | 18禁观看日本| 精品国产一区二区久久| 99久久人妻综合| 日韩伦理黄色片| 如何舔出高潮| 26uuu在线亚洲综合色| 如何舔出高潮| 亚洲美女视频黄频| 日日啪夜夜爽| 可以免费在线观看a视频的电影网站 | 久久这里有精品视频免费| 国产男人的电影天堂91| 亚洲精品久久成人aⅴ小说| 熟女少妇亚洲综合色aaa.| 亚洲国产看品久久| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 久久久精品94久久精品| 午夜福利,免费看| 久久综合国产亚洲精品| 十八禁网站网址无遮挡| 91精品伊人久久大香线蕉| 女性生殖器流出的白浆| 久久午夜福利片| 成人亚洲欧美一区二区av| 日韩欧美精品免费久久| 肉色欧美久久久久久久蜜桃| 在线 av 中文字幕| 91午夜精品亚洲一区二区三区| 国产精品免费视频内射| 成人手机av| 亚洲国产欧美日韩在线播放| 99香蕉大伊视频| 午夜精品国产一区二区电影| 亚洲欧洲国产日韩| 丝袜喷水一区| 亚洲av福利一区| 日韩精品免费视频一区二区三区| 国产亚洲最大av| 国产精品不卡视频一区二区| 最近中文字幕高清免费大全6| 亚洲国产日韩一区二区| 色播在线永久视频| 久久人妻熟女aⅴ| 久久久久网色| 这个男人来自地球电影免费观看 | 欧美黄色片欧美黄色片| 黄片播放在线免费| 伦理电影大哥的女人| 午夜日韩欧美国产| 国产av精品麻豆| 久久国内精品自在自线图片| 精品人妻偷拍中文字幕| 国产精品国产三级国产专区5o| 中文欧美无线码| 国产精品久久久久久久久免| 国产精品秋霞免费鲁丝片| 精品人妻一区二区三区麻豆| 国产男人的电影天堂91| 欧美97在线视频| 可以免费在线观看a视频的电影网站 | 亚洲精品国产av成人精品| 国产在线免费精品| 久久精品国产综合久久久| 日本免费在线观看一区| 水蜜桃什么品种好| 91国产中文字幕| 国产精品蜜桃在线观看| 亚洲av.av天堂| 久久久久久久亚洲中文字幕| 国产精品 欧美亚洲| 欧美老熟妇乱子伦牲交| 免费观看性生交大片5| videosex国产| 亚洲av免费高清在线观看| 交换朋友夫妻互换小说| 激情视频va一区二区三区| 成人黄色视频免费在线看| 久久人人97超碰香蕉20202| 性少妇av在线| 日韩电影二区| 国产精品女同一区二区软件| 国产欧美亚洲国产| 一区二区三区乱码不卡18| 亚洲国产精品一区二区三区在线| 一级黄片播放器| 制服诱惑二区| 高清在线视频一区二区三区| 久久久国产一区二区| 国产午夜精品一二区理论片| 人人妻人人添人人爽欧美一区卜| 看十八女毛片水多多多| 老司机亚洲免费影院| 1024视频免费在线观看| 久久久久久人人人人人| 99热网站在线观看| 国产淫语在线视频| 国产成人av激情在线播放| 国产一区二区激情短视频 | 亚洲激情五月婷婷啪啪| 飞空精品影院首页| 成人国语在线视频| 国产精品国产av在线观看| 日韩成人av中文字幕在线观看| 国产成人精品在线电影| 国产片特级美女逼逼视频| 日产精品乱码卡一卡2卡三| 晚上一个人看的免费电影| 美女国产视频在线观看| 精品人妻在线不人妻| 青春草亚洲视频在线观看| 国产成人精品福利久久| 亚洲av免费高清在线观看| 亚洲精品日韩在线中文字幕| 国产爽快片一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 一区二区日韩欧美中文字幕| 交换朋友夫妻互换小说| 久久久久国产精品人妻一区二区| 亚洲人成电影观看| 国产深夜福利视频在线观看| 国产男女内射视频| 捣出白浆h1v1| 免费在线观看视频国产中文字幕亚洲 | 我的亚洲天堂| 黑丝袜美女国产一区| 欧美日韩av久久| 午夜av观看不卡| 国产日韩欧美亚洲二区| 91久久精品国产一区二区三区| 成年人免费黄色播放视频| 夜夜骑夜夜射夜夜干| 精品人妻偷拍中文字幕| 欧美成人精品欧美一级黄| 国产成人免费观看mmmm| 18+在线观看网站| 一二三四中文在线观看免费高清| 免费观看av网站的网址| 一区二区三区激情视频| 久久久久久伊人网av| 在线观看免费视频网站a站| 9色porny在线观看| 老司机影院毛片| 欧美日韩精品网址| 99香蕉大伊视频| 亚洲欧美色中文字幕在线| 成人18禁高潮啪啪吃奶动态图| 久久久久久久久久久免费av| 欧美日韩精品网址| 日韩伦理黄色片| 国产色婷婷99| 看免费av毛片| 日本-黄色视频高清免费观看| 三级国产精品片| 精品一区二区免费观看| 另类亚洲欧美激情| 精品卡一卡二卡四卡免费| 欧美精品一区二区大全| 中文字幕制服av| 成年女人毛片免费观看观看9 | 天堂俺去俺来也www色官网| 菩萨蛮人人尽说江南好唐韦庄| 成年女人在线观看亚洲视频| 看免费av毛片| 亚洲一码二码三码区别大吗| 国产精品久久久久久精品古装| 成人国语在线视频| 亚洲,欧美精品.| 精品久久久精品久久久| 久久99精品国语久久久| 国产亚洲av片在线观看秒播厂| 中文字幕最新亚洲高清| 99国产综合亚洲精品| 少妇熟女欧美另类| 亚洲av男天堂| 热99久久久久精品小说推荐| 午夜日韩欧美国产| 丝袜人妻中文字幕| 亚洲少妇的诱惑av| 亚洲av电影在线观看一区二区三区| 久久久久久久久久人人人人人人| 侵犯人妻中文字幕一二三四区| 久久久精品国产亚洲av高清涩受| 少妇被粗大的猛进出69影院| 久久午夜福利片| 这个男人来自地球电影免费观看 | 美女脱内裤让男人舔精品视频| 亚洲欧美日韩另类电影网站| 亚洲av男天堂| 免费在线观看完整版高清| 亚洲综合精品二区| 香蕉国产在线看| 亚洲第一av免费看| 日本爱情动作片www.在线观看| 国产成人精品久久久久久| 婷婷色av中文字幕| 啦啦啦在线观看免费高清www| 久久久久国产网址| 午夜激情av网站| 精品国产一区二区久久| kizo精华| 我要看黄色一级片免费的| 亚洲精品,欧美精品| 欧美av亚洲av综合av国产av | 女人被躁到高潮嗷嗷叫费观| 熟女少妇亚洲综合色aaa.| 五月开心婷婷网| 一级a爱视频在线免费观看| 国产麻豆69| 日韩在线高清观看一区二区三区| 精品久久久精品久久久| 日本av免费视频播放| 国产色婷婷99| 午夜福利视频精品| 亚洲精品自拍成人| 十八禁高潮呻吟视频| 亚洲av电影在线观看一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av综合色区一区| 久久久国产一区二区| 国产免费现黄频在线看| 国产日韩欧美视频二区| 少妇的丰满在线观看| 久久鲁丝午夜福利片| av女优亚洲男人天堂| 一区二区av电影网| 黄片小视频在线播放| 欧美 亚洲 国产 日韩一| 日韩中字成人| 国产午夜精品一二区理论片| 亚洲欧美清纯卡通| 美女主播在线视频| 国产精品香港三级国产av潘金莲 | 久久久久久久久免费视频了| 最近最新中文字幕大全免费视频 | 成人漫画全彩无遮挡| 免费观看在线日韩| 亚洲精品国产av蜜桃| 久久ye,这里只有精品| 黄频高清免费视频| 欧美亚洲日本最大视频资源| 午夜久久久在线观看| 99精国产麻豆久久婷婷| 老汉色∧v一级毛片| 多毛熟女@视频| 亚洲人成77777在线视频| 热re99久久精品国产66热6| 高清av免费在线| 老司机影院毛片| 蜜桃国产av成人99| 国产精品 欧美亚洲| av有码第一页| 三级国产精品片| 大片电影免费在线观看免费| 成人亚洲精品一区在线观看| 麻豆精品久久久久久蜜桃| 又粗又硬又长又爽又黄的视频| 成人亚洲精品一区在线观看| 在线观看www视频免费| 熟女少妇亚洲综合色aaa.| 欧美日韩综合久久久久久|