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

    Influence of short-term experimental warming on heat-water processes of the active layer in a swamp meadow ecosystem of the Qinghai-Tibet Plateau

    2016-10-17 09:06:20GuangShengLiuGenXuWang
    Sciences in Cold and Arid Regions 2016年2期

    GuangSheng Liu, GenXu Wang

    1. College of Environmental Science and Engineering, Xiamen University of Technology, Xiamen, Fujian 361024, China

    2. Key Laboratory of Mountain Surface Processes and Ecological Regulation, Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China

    ?

    Influence of short-term experimental warming on heat-water processes of the active layer in a swamp meadow ecosystem of the Qinghai-Tibet Plateau

    GuangSheng Liu1*, GenXu Wang2

    1. College of Environmental Science and Engineering, Xiamen University of Technology, Xiamen, Fujian 361024, China

    2. Key Laboratory of Mountain Surface Processes and Ecological Regulation, Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China

    ABSTRACT

    Climate change is now evident in the Qinghai-Tibet Plateau (QTP), with impacts on the alpine ecosystem, particularly on water and heat balance between the active layer and the atmosphere. Thus, we document the basic characteristics of changes in the water and heat dynamics in response to experimental warming in a typical alpine swamp meadow ecosystem. Data sets under open top chambers (OTC) and the control manipulations were collected over a complete year. The results show that annual (2008) air temperatures of OTC-1 and OTC-2 were 6.7 °C and 3.5 °C warmer than the control. Rising temperature promotes plant growth and development. The freeze-thaw and isothermal days of OTCs appeared more frequently than the control, owing to comparably higher water and better vegetation conditions. OTCs soil moisture decreased with the decrease of soil depth; however, there was an obviously middle dry aquifer of the control, which is familiar in QTP. Moreover, experimental warming led to an increase in topsoil water content due to poorly drained swamp meadow ecosystem with higher organic matter content and thicker root horizons. The results of this study will have some contributions to alpine cold ecosystem water-heat process and water cycle under climate change.

    experimental warming; open top chamber; Qinghai-Tibet Plateau; soil moisture; soil temperature; alpine swamp meadow ecosystem

    1 Introduction

    IPCC (2013) projected that mean global temperature would increase 1.8-4.0 °C for the 21st century. The rate of temperature increase in high latitude and altitude regions would be even larger (Thomas et al.,2004). The Qinghai-Tibet Plateau (QTP) houses the greatest portion of high altitude permafrost in the world (Zhou et al., 2000) and is extremely sensitive to global warming. The mean annual temperature increased at approximately 0.21 °C/decade over the past 50 years at QTP (Liu and Chen, 2000), which exceeds the Northern Hemisphere and the same latitudinal zone during the same period (IPCC, 2013).

    Climate warming significantly influences permafrost, active layer conditions and seasonal frozen soils and, as a result, alpine swamp meadow ecosystem stability (Osterkamp and Lachenbruch, 1990; Jorgenson et al., 2001; Nelson et al., 2003; Walker et al., 2003). The influence of a changing climate on permafrost, activelayer thermal regime, and vegetation, may activate,through complex mechanisms, both negative and positive feedbacks related to global warming (Walker et al.,2003; Chapin et al., 2005). The QTP's ecosystem has experienced a loss of vegetation, glacial ablation, permafrost degradation and the disappearance of wetlands and lakes during the last 30 years in response to a warmer climate (Wang et al., 2007).

    Soil heat-water processes reflect the spatiotemporal distribution of soil temperature (heat) and moisture (dry vs. wet), thereby constituting an extremely important component in the energy and water cycles of the permafrost soil-atmosphere system (Li et al., 2002; Hu et al., 2009). These cycles affect the soil water and heat transfer processes and water hosting condition, and would also be key to alpine ecosystem stability and could have an important feedback effect on global climate (Woo and Winter, 1993; Carey and Woo, 1999;Zhang et al., 2003). Transfer in energy and water balance, accompanied with the response of alpine frost ecosystems to global climate change, would be likely reflected by changes in water content and the associated heat of frozen soil. Meanwhile, some monitoring of the soil heat-water processes have been performed in the QTP (Zhao et al., 2000; Wu et al., 2003; Wang et al.,2007) in recent years, and the studies illustrate the response of freezing-thawing processes to a cold region ecological environment. Wang et al. (2007) showed that moisture in the active layer in a permafrost region had an important influence on water cycling and the growth of meadow vegetation. Investigations in sub-arctic Canada and Alaska showed that climate changes led to an increase in topsoil water content in the permafrost regions (Jorgenson et al., 2001; Turetsky et al., 2002).

    Therefore, understanding the laws governing variations in soil moisture and associated temperature of the alpine swamp meadow ecosystem is important in understanding regional variations in water and heat cycling arising from global changes. For this reason,an open top chamber (OTC) was adopted to simulate global warming. The influence of short-term experimental warming on heat-water processes of the active layer in an alpine swamp meadow ecosystem of the Qinghai-Tibet Plateau was analyzed through data sets under different manipulations over a complete year. The purpose of this study is to analyze the response of the hydrologic and thermal regime to plot-scale experimental warming at a typical alpine swamp meadow ecosystem underlain by permafrost in QTP.

    2 Materials and methodology

    2.1Study area

    The Fenghuoshan permafrost region on the Qinghai-Tibet Plateau (34°40'N-34°48'N, 92°50'E-93°03'E), located in the upper watershed of the Zuomo Xikong River, a tributary of the Yangtze River, was chosen as the system of study. The region covers a total area of 127.63 km2at elevations ranging from 4,680 to 5,360 m, belonging to an arid climate area and without glacial or multiyear snow-pack. The climate is cold and dry. The mean annual (1973-2005) air temperature, extreme maximum air temperature, extreme minimum air temperature are -5.2, 23.2 and -37.7 °C,respectively. The mean annual precipitation is 290.9 mm, 80% of which falls during the growing season(from May to September). The mean annual evaporation is 1,316.9 mm. The mean annual ground temperature is -1.5 °C to -4.0 °C (Wang et al., 2008), the permafrost thickness ranges from 50 m to 120 m and the active layer thickness is in the range of 0.8 m to 2.5 m. Additionally, these areas belong to an alpine meadow system which are dominated by plant species such as Kobresia tibetica Maximowicz, Stipa aliena Keng and Festuca spp.. The main soil type is alpine meadow soil (Table 1), and the parent material is mostly quaternary sediment, which includes the weathering material of metamorphic rock and intermediate intrusive rock (Wang et al., 2008).

    Table 1 Physical and chemical characteristic of selected alpine swamp meadow soils of the Qinghai-Tibet Plateau

    2.2Methodology

    2.2.1Layout of experimental warming equipment The OTC is an experimental warming structure first used by the International Tundra Experiment(ITEX) to study the influence of global warming on high altitude and latitude region ecosystems (Hollister and Webber, 2000).

    The experimental plots have almost identical vegetation coverage (90%), vegetation type and terrain. Vegetation and soil conditions were classified prior to the onset of the manipulation (Table 1). The results show similar characteristics between all plots where,no significant correlation (p >0.4) was found between the three plots for any physical and chemical characteristics. The vegetation conditions of three plots were all under non-degraded status, and were similar to each other.

    Two OTCs were made of 8 mm thick PMMA fiber material, with a top and basal area of 1.5m×1.5m and a height of 1.7 m. The base was embedded in the soil at a depth of 0.3 m. In order to gain a higher rate of temperature increase, one of the OTCs (OTC-1) included a partial closed top manipulation, through a 0.6 m diameter hole on the top. This type of OTC was used in the study of Havstrom et al. (1993) in subarctic (Abisko)and high-arctic sites (Svalbard). Additionally, the other OTC (OTC-2) represents a traditional open top chamber that has a complete open top manipulation. The control plot, kept the original state with no change. The experiment aimed to include two levels warming, with a mean annual air temperature of 3-6 °C (OTC-1) and 1-3 °C (OTC-2) warmer than the control.

    2.2.2Layout of the observational experiment

    Each manipulation housed a pair of 1.2 m deep soil water and temperature observation wells, five soil moisture and temperature sensors were located at depths of 0.05, 0.20, 0.40, 0.65, and 1.20 m, respectively. Soil moisture was determined by frequency domain reflectometry (FDR), using a calibrated soil moisture sensor equipped with a Theta-probe (Hollend Eijkelammp Co. Giesbeek). Soil moisture was derived from changes in the soil's dielectric constant that was converted to a millivolt signal with an accuracy of ±2%. It should be noted that the soil volume moisture measured by FDR is the unfrozen soil moisture. Soil temperature was monitored using a thermal resistance sensor sensitive to temperature changes in the range of -40 °C to 50 °C, with an overall system precision of ±0.02 °C.

    The thermal hygrometer was placed at 80-cm height above the alpine swamp meadow to measure atmospheric temperature and humidity. All of the probes were connected to a Campbell CR1000 data logger, developed by Campbell Scientific Ltd.. The frequency of data sample was 30 min, and the data set was collected over a one-year period from November 1, 2007 to December 31, 2008.

    The aboveground biomass under different manipulations was sampled during the growing season by clipping the vegetation to ground level in 0.1m×0.1m subplots. Samples were dried for 72 hours at 55 °C and then weighed.

    2.2.3Data analysis

    The mean unfrozen moisture storage was calculated by:

    where, θiand Thiare the soil moisture content and the represent thickness of soil at depths of 0.05, 0.20, 0.40,0.65, and 1.20 m, respectively. For example, the represent thickness at 0.05 m depth is 0.1 m (0-0.1 m).

    In order to describe the thermal regime under different manipulations, some parameters were computed(Guglielmin et al., 2008):

    1) The number of isothermal days, during which all the hourly measuring temperature ranges only between ±0.5 °C, indicating that the energy in the soil is consumed to change the water phase both during the freezing and thawing process, when latent heat exchange produces the "zero curtain".

    2) The number of frozen days, during which all the hourly measuring temperatures are negative and at least one is colder than -0.5 °C, indicating that the soil profile is frozen.

    3) The number of thawed days, during which all the hourly measuring temperatures are positive and at least one is warmer than -0.5 °C, indicating that the soil profile is in a liquid state.

    4) The number of freeze-thaw days, during which there are both negative and positive temperatures and at least one value is greater than ±0.5 °C, indicating that freeze-thaw cycles occur.

    3 Results

    3.1Experimental warming effects of OTCs

    Mean annual air temperatures (2008) under OTC-1 and OTC-2 manipulations were 6.7 °C and 3.5 °C higher than the control, respectively (Figure 1). This almost overlaps with the 1.8-4.0 °C increasing temperature range predicted by IPCC (2013) using Global/General circulation model. Air temperatures observed inside the OTCs have typically been 0.5-3.0 °C warmer than outside (Marion et al., 1997), which is similar with OTCs.

    At the same time, the maximum monthly air temperature of OTC-1 and OTC-2 were 9.2 °C and 4.3 °C warmer than the control, respectively. However, the minimum monthly air temperature of OTC-1 and OTC-2 were just 4.7 and 1.7 °C warmer than the control, respectively (Figure 1). These results are different from other results in that the increasing rate of summer is larger than that of winter (Marion et al., 1997). These results are mainly due to variation in sunshine hours,which is unique to the climate of QTP and differs from other regions previously studied (Figure 1).

    There was significant difference (p <0.001) in the above ground biomass between OTCs and the control(Table 2). The aboveground biomass of OTC-1 and OTC-2 were 1.57 and 2.11 times of the control, respectively, in the plant-growing period. OTCs improve the microclimate environment of the alpine swamp meadow ecosystem and fulfill the energy need of the vegetation, in favor of growth and development of the vegetation, thus causing an increase in aboveground biomass. Moreover, protection from wind may play a more important role than temperature in providing a favorable environment for plant growth within OTCs(Marion et al., 1997). This is consistent with results from previous studies (Li et al., 2004) which showed that the biomass increases in short-term experimental warming. Besides, lower biomass of OTC-1 than OTC-2 is due to drought treatment to OTC-1 and extreme experimental warming.

    Figure 1 Relationship between increasing monthly air temperature and sunshine hours under different manipulations

    Table 2 Underground biomass (g/m2) variation of alpine swamp meadow ecosystem in plant growing period (2007-2008) under different manipulations

    3.2Division of heat-water processes

    As presented in Figure 2, the coupled changes between soil moisture and temperature occurring over the year in the active layer can be divided into four periods based on freeze-thaw states of the permafrost soil (Osterkamp and Romanovsky, 1997; Wang et al.,2007; Hu et al., 2009):

    1) The freezing process period (A) -- from mid-October to early December, when the active layer is in the process of freezing, which shows a decrease of the total unfrozen moisture storage from 350-370 mm to 70-90 mm;

    2) The entirely frozen period (B) -- from mid-December to mid-May, when the active layer is in a fully frozen state and the total unfrozen moisture storage remains at a constant level in 70-90 mm;

    3) The thawing process period (C) -- from mid-May to early July, when the active layer is in the process of thawing, and the total unfrozen moisture storage increases from 70-90 mm to the vertex value;

    4) The fully thawed period (D) -- from early July to mid-October, when the entire soil layer (0-1.2 m) is in a thawed state, and the total unfrozen moisture storage fluctuates from 300 mm to 400 mm.

    Mean annual ground temperatures (2008) of OTC-1, OTC-2 and the control were -0.05, -1.64 and -1.71 °C, respectively (Figure 2a). Mean annual soil temperatures (2008) of OTC-1 were 3.49, 2.30, 1.46,1.34 and 0.45 °C warmer than the control for the 0.05,0.20, 0.40, 0.65 and 1.20 m soil depths, respectively. However, the differences between OTC-2 and the control were -0.07, 0.24, 0.22, 0.20 and 0.41 °C, respectively, which were inconsistent with OTC-1. Clearly, there was a significant difference between the mean ground temperatures of OTC-1 and the control manipulation; however, there is little difference between OTC-2 and the control manipulations.

    There were some differences between the total unfrozen moisture storage (0-1.20 m) under the OTC-1, OTC-2 and the control manipulations (Figure 2b). The mean annual unfrozen moisture storage(2008) of OTC-1, OTC-2 and the control were 230.7,212.2, and 195.4 mm, respectively. These are mainly due to much longer duration of fully thawed period in OTCs and the observed soil moisture was unfrozen water, which lead to relatively high value of unfrozen moisture storage in OTCs.

    Figure 2 Mean soil temperature (a) and total unfrozen moisture storage (b) of the active layers (0-1.20 m) under OTC-1, OTC-2 and the control manipulations

    3.3Thermal regime

    3.3.1Soil thermal condition

    Results show significant differences in the number of isothermal days, frozen days, thawed days and freeze-thaw days at different depths under different manipulations (Table 3).

    In the ground surface (0.05 cm), the isothermal days mainly appeared in May and October in OTC-2 and the control, however, appeared only in November in OTC-1, showing that "zero curtain" only appeared in the freezing process in OTC-1. It was due to greater air temperature fluctuation in OTC-1, which led to severe freeze-thaw cycle on the ground surface. At a soil depth of 0.2 m, the isothermal days of OTC-1 mainly appeared in November, but they have appeared in May, October and November in OTC-2 and the control. The isothermal days of OTCs were more than the control, due to comparably higher water conditions and vegetation conditions in OTCs (Table 2).

    With an increase of soil depths, frozen days of all manipulations generally increased, whereas thawed days show a converse, decreasing trend. The frozen days of OTC-1 was significantly less than the control,whereas OTC-2 is only slightly less than the control. This indicated that frozen days decrease with temperature increase. In contrast, the thawed days of OTC-1 was significantly more than the control,whereas OTC-2 was only slightly more than the control, indicating that thawed days increase with increasingairtemperature.Moreover,the freeze-thaw days near the surface appeared more frequently in OTC-1 than in OTC-2 and the control. The freeze-thaw days appeared only in April and May in OTC-2 and the control, but also appeared in November in OTC-1 due to greater air temperature fluctuation (Guglielmin et al., 2008).

    3.3.2Soil temperature variation process

    The contour chart clearly shows the duration of freeze-thawing, frozen depth, the extension within a certain period and freezing-thawing front movement(Figure 3). With the periodic cycle of solar radiation in one year, soil heat-water processes achieve their own related cycles.

    As presented in Figure 3, in the freezing process period (A), the dates of active layers beginning to freeze at different depths of OTCs were earlier than the control, but all coincide in late October and early November. The period when the freezing front moved from the surface to 0.55 m depth in OTC-1, OTC-2 and the control were 20, 25-30 and 30 days, respectively,and the duration of the whole soil profile becoming entirely frozen are 40, 67 and 70 days, respectively.This indicates that the durations of freezing process shorten and the freezing front moving rate increased with temperature increase.

    In the entirely frozen period (B), days it took for the temperature to decrease from -2 °C to -6 °C in the shallow layer (0.05 m) in OTC-1, OTC-2 and the control were 43, 25 and 20 days, respectively, indicating that rate of temperature change in OTCs were lower than the control. In contrast, days it took for the temperature to decrease from -1 °C to -5 °C in the deeper layer (1.2 m) in OTC-1, OTC-2 and the control were 36, 29 and 18 days, respectively, indicating that rate of temperature change in OTCs were higher than the control. Meanwhile, the areas encompassed by the 0 °C isothermal were OTC-1 < OTC-2 < Control,indicating that entirely frozen durations decreased with temperature increase.

    In the thawing process period (C), the dates of soil layers (0-1.2 m) beginning to thaw at different depths under different manipulations were disperse, ranged from May to August, and were delayed with the increase of soil depth. With the temperature increase, the beginning of thawing days moved ahead, but the differences in duration of thawing process were minor.

    In the fully thawed period (D), contrary to the entirely frozen period, rate of temperature change in OTCs were higher than the control in shallow layers and trended lower in deeper layers. Meanwhile, the areas encompassed by the 0 °C isotherm of OTC-1 were larger than those of OTC-2 and the control, indicating that duration of fully thawed period extended with temperature increase.

    Table 3 Number of isothermal days, frozen days, thawed days and freeze-thaw days (2008)at different depths under different manipulations

    Figure 3 Thermal regime under the control (a), OTC-2 (b) and OTC-1 (c) manipulations

    3.4Moisture distribution

    As presented in Figure 4, in the fully thawed period(D), there were significant differences between moisture distribution in the active layers of OTC-1, OTC-2 and the control. In OTC-1 and OTC-2, soil moisture decreased with the depth increase, varying from 46% and 43% on the ground surface to 28% and 33% at 1.2 m soil depth, respectively. The decreasing rate of OTC-1 was larger than that OTC-2. However, the soil moisture distribution characteristic of the control was similar to the alpine meadow ecosystems study result(Wang et al., 2008; Hu et al., 2009), which showed a dry layer at 0.65 m depth and two higher moisture aquifers at the ground surface and 1.1 m depth. The existence of a middle dry layer is because the upper layer is markedly influenced by precipitation, which provides a water source for the upper soil; the lower layer is mainly influenced by a watertight layer near the bottom of an active layer where water migrates downwards with thawing soil and concentrates relatively. Hence, the middle layer is a dry belt where summer precipitation cannot infiltrate (Zhao et al., 2000).

    In the freezing process period (A) and the thawing process period (C), relative to soil moisture distribution in the period (D), the change rate of soil moisture in OTCs decreased with soil depth increase. Furthermore, soil moisture showed fluctuating change in the shallow layer due to the impact of air temperature and precipitation.

    In the entirely frozen period (B), the unfrozen soil moisture under different manipulations was all about 10%. The unfrozen soil moisture in OTC-1 increased with soil depth increase, ranging from 5% to 11%.

    Figure 4 Soil moisture regime under the control (a), OTC-2 (b) and OTC-1 (c) manipulations

    4 Discussion

    There is a growing body of research in water and heat cycles in permafrost regions. Some monitoring of the soil heat-water processes have been performed in the QTP in recent years (Zhao et al., 2000; Wu et al.,2002; Wang et al., 2008). These studies illustrate soil water circulation in the active layer increased as vegetation cover decreased. Besides, significant response of the alpine frost ecosystems upon the QTP to global climate change would have an important feedback effect on global climate (Zhang et al., 2003). Therefore,further observations and studies are required to quantify the heat-water transfer processes in typical alpine swamp meadow ecosystems and its feedback that occur between climate change and changes to vegetation.

    4.1Effect of short-term experimental warming on soil temperature

    There was a significant difference between mean ground temperatures of OTC-1 and the control manipulation; however, there is little difference between OTC-2 and the control manipulations. The depressed soil temperature in OTC was also found in the Salix site at Alexandra Fjord (Marion et al., 1997). First, it was because OTC-2 has no inclined sides, which dissipates heat through the top by free convection aided by the forced convective flow of the wind. Moreover,they were likely due to vegetation and ground shading,and relative higher moisture in OTC-2 and these may simply reflect the complex physics of energy balance within chambers. Finally, the vegetation condition of OTC-2 was better than the control (Table 2). As reported by Wang et al. (2008) and Hu et al. (2009),vegetation has a warming effect on the freezing process and entirely frozen period; and has a cooling effect on the thawing process and fully thawed period,which led to higher soil temperature in cold season and lower soil temperature in warm season in OTC-2.

    4.2Effect of short-term experimental warming on soil moisture

    As presented in Figure 5, on the ground surface (5 cm), soil moisture was shown as OTC-1 > OTC-2 >control in the fully thawed period. In addition, at 0.20-0.65 m depth, soil moisture was shown as OTC-2 > OTC-1 > control. Experimental warming led to an increase in topsoil water content. Model projections suggest evapotranspiration may increase due to climate warming, resulting in reduce near-surface soil moisture and greater aridity in water-limited systems around the world (Komescu et al., 1998; IPCC, 2013). This is contrary to the results of this study.

    Studies on hydrological processes in the Mackenzie Global Energy and Water Cycle Experiment(Rouse, 2000) and a North American discontinuous permafrost region wetland (Carey and Woo, 1999;Oechel et al., 2000) show that in discontinuous permafrost regions, climate change has led to an increase in topsoil moisture content of the permafrost region(Jorgenson et al., 2001; Turetsky et al., 2002). However, studies in other regions (Sweden) have reported reduced topsoil moisture content and even drought(Christensen et al., 2004). Zavaleta et al. (2003) found that interactions between warming and the dominant biota in a grassland ecosystem produced the reverse effect. Simulated warming increased spring ground surface moisture by 5%-10% under both ambient and elevated CO2. Bremer et al. (1996) reported that surface soil moisture was consistently higher than other treatments under elevated CO2chambers, which is similar to OTCs, tends to be lowest early in the season and then became higher than that of outside. These trends are consistent with our research results, due to evapotranspiration difference.

    Figure 5 Active layer soil moisture regime in the fully thawed period under different warming manipulations in 2008

    Moreover, experimental warming led to an increase in above ground biomass (Table 2). This means better vegetation conditions, thicker root horizons and more soil organic matter content in OTCs. As presented in Figure 6, the soil moisture in the ground surface increase with the soil organic matter content increase. Due to poorly drained swamp meadow ecosystem with higher organic matter content and thicker root horizons, higher water binding capacity occurs in OTCs. Therefore, experimental warming led to an increase in topsoil water content due to poorly drained swamp meadow ecosystem with higher organic matter content and thicker root horizons.

    At a depth of 1.20 m, the soil moisture appeared to be the control > OTC-2 > OTC-1. Experiments of Fukuda et al. (1980) indicated that the amount of migrated moisture in the freezing front has a close relationship with the freezing rate of soils. The duration of freezing process and the frozen period became shorter with a temperature increase (Figure 3, Table 3), so less water migrates to the deeper soil layer, which led to higher moisture in the control.

    These results show that the soil moisture distribution in the profile changes significantly with increasing air temperature. Short-term experimental warming led to an increase in topsoil water content and a deep layer moisture loss, although the unfrozen soil moisture storage increased.

    Figure 6 The relationships between soil organic matter content and soil moisture on the ground surface

    5 Conclusions

    This study describes the soil thermal regime and moisture distribution in an alpine swamp meadow ecosystem due to short-term experimental warming. Air temperatures (2008) of OTC-1 and OTC-2 were 6.70 and 3.48 °C warmer than the control, respectively. Rising temperature promotes plant growth and development, thus making a significantly higher aboveground biomass in OTCs.

    The dates marking the beginning of active layer freezing were postponed and the dates marking the beginning of thawing occurred earlier, shortening the durations of entirely frozen period and frozen days; the fully thawed period extended and the thawed days increased.

    OTCs soil moisture decreased with the soil depth decrease; however, there was an obviously middle dry aquifer of the control, which is familiar in QTP. Moreover, experimental warming lead to an increase in topsoil water content due to poorly drained swamp meadow ecosystem with higher organic matter content and thicker root horizons. The results of this study will have some contributions to alpine cold ecosystem water-heat process and water cycle under climate change.

    Acknowledgments:

    This study was funded by The Natural Science foundation of China (No. 40925002) and The Natural Science foundation of Fujian Province (No. 2015J05114). The authors thank the anonymous reviewers for their fruitful comments, and help in improving the presentation of this research.

    Bremer DJ, Ham JM, Owensby CE, 1996. Effect of elevated atmospheric carbon dioxide and open-top chambers on transpiration in tallgrass prairie. Journal of Environmental Quality, 25: 691-701. DOI: 10.2134/jeq1996.00472425002500040008x.

    Carey SK, Woo MK, 1999. Hydrology of two slopes in subarctic Yukon, Canada. Hydrological Processes, 13: 2549-2562. DOI: 10.1002/(SICI)1099-1085.

    Chapin III FS, Sturm M, Serreze MC, et al., 2005. Role of land-surface changes in Arctic summer warming. Science, 310: 657-660. DOI: 10.1029/2003GL018680.

    Christensen TR, Johansson T, Kerman HJ, et al., 2004. Thawing subarctic permafrost: effects on vegetation and methane emissions.GeophysicalResearchLetters,31:L04501. DOI: 10.1029/2003GL018680.

    Fukuda M, Orhun A, Luthin JN, 1980. Experimental studies of coupled heat and moisture transfer in soils during freezing. Cold Region Science and Technology, 9(3): 223-232. DOI: 10.1016/0165-30232X(80)90028-2.

    Guglielmin M, Evans CE, Cannone N, 2008. Active layer thermal regime under different vegetation conditions in permafrost areas. A case study at Signy Island (Maritime Antarctica). Geoderma, 144: 74-85. DOI: 10.1016/j.geoderma.2007.10.010.

    Havstrom M, Callaghan TV, Jonasson S, 1993. Differential growth responses of Cassiope tetragona, an arctic dwarf-shrub, to environmental perturbations among three contrasting high and subarctic sites. Oikos, 66: 389-402. DOI: 10.2307/3544933.

    Hollister RD, Webber PJ, 2000. Biotic validation of small open-top chambers in a tundra ecosystem. Global Change Biology, 6: 835-842. DOI: 10.1046/j.1365-2486.2000.00363.x.

    Hu HC, Wang GX, Liu GS, 2009. Influences of alpine ecosystem degradation on soil temperature in the freezing-thawing process on Qinghai-Tibet Plateau. Environmental Geology, 57(6): 1391-1397. DOI: 10.1007/s00254-008-1417-7.

    IPCC, 2013. Climate Change 2013: The Physical Science Basis. Cambridge, United Kingdom and New York: Cambridge University Press.

    Jorgenson MT, Racine CH, Walters JC, 2001. Permafrost degradation and ecological changes associated with a warming in central Alaska. Climatic Change, 48: 551-579. DOI: 10.1023/A:1005667424292.

    Komescu AU, Eikan A, Oz S, 1998. Possible impact of climate change on soil moisture availability in the southwest Anatolia development project region (GAP): An analysis from an agricultural drought perspective. Climatic Change, 40: 519-545. DOI: 10.1023/A:1005349408201.

    Li SX, Nan ZT, Zhao L, 2002. Impact of soil freezing and thawing process on thermal exchange between atmosphere and ground surface. Journal of Glaciology and Geocryology, 24(5): 506-511.

    Li YN, Zhao L, Zhao XQ, et al., 2004. Effects of a 5-year mimic temperature increase to the structure and productivity of Kobresia humilis meadow. Acta Agrestia Sinica, 12: 236-239.

    Liu XD, Chen BD, 2000. Climatic warming in the Tibetan Plateau during recent decades. International Journal of Climatology, 20: 1729-1742. DOI: 10.1002/1097-0088(20001130)20:14<1729: AID-JOC556>3.0.CO;2-Y.

    Marion GM, Henry GHR, Freckman DW, 1997. Open-top designs for manipulating field temperatures in high-latitude ecosystems. Global Change Biology, 3(Supplement 1): 20-32. DOI: 10.1111/j.1365-2486.1997.gcb136.x.

    Nelson FE, Lachenbruch AH, Woo MK, et al., 2003. Permafrost and changing climate. In: Proceedings of the Sixth International Conference on Permafrost, Beijing, China, pp. 987-1005.

    Oechel WC, Vourlitis GL, Hastings SJ, et al., 2000. Acclimation of ecosystem CO2exchange in the Alaskan Arctic in response to decadal climate warming. Nature, 406: 978-980. DOI: 10.1038/35023137.

    Osterkamp TE, Lachenbruch AH, 1990. Thermal regime of permafrost in Alaska and predicted global warming. Journal of ColdRegionsEngineering,4(1):38-42.DOI: 10.1061/(ASCE)0887-381X(1990)4:1(38).

    Osterkamp TE, Romanovsky VE, 1997. Freezing of active layer on the coastal plain of the Alaskan Arctic. Permafrost and Periglacial Processes, 8: 23-44. DOI: 10.1002/(SICI)1099-1530.

    Rouse WR, 2000. Progress in hydrological research in the Mackenzie GEWEX study. Hydrological Processes, 14: 1667-1685. DOI: 10.1002/1099-1085.

    Thomas D, Cameron A, Green GE, et al., 2004. Extinction risk from climate change. Nature, 427: 145-148. DOI: 10.1038/nature02121.

    Turetsky MR, Kelman-Wieder R, Vitt DH, 2002. Boreal peatland C fluxes under varying permafrost regimes. Soil Biology and Biochemistry, 34: 907-912. DOI: 10.1016/S0038-0717(02)00022-6.

    Walker DA, Jia GJ, Epstein HE, et al., 2003. Vegetation-soil-thaw-depth relationships along a low arctic bioclimate gradient, Alaska: synthesis of information from the ATLAS studies. Permafrost and Periglacial Processes, 14: 103-123. DOI: 10.1002/ppp.452.

    Wang GX, Li YS, Hu HC, et al., 2008. Synergistic effect of vegetation and air temperature changes on soil water content in alpine frost meadow soil in the permafrost region of Qinghai-Tibet. Hydrological Processes, 22: 3310-3320. DOI: 10.1002/hyp.

    Wang GX, Li YS, Wang YB, et al., 2007. Typical alpine wetland system changes on the Qinghai-Tibet Plateau in recent 40 years. Acta Geographica Sinica, 62(5): 481-491.

    Woo MK, Winter TC, 1993. The role of permafrost and seasonal frost in the hydrology of northern wetlands in North America. Journal of Hydrology, 141: 5-31. DOI: 10.1016/0022-1694(93)90043-9.

    Wu QB, Shen YP, Shi B, 2003. Relationship between frozen soil together with its water heat process and ecological environment in the Tibetan Plateau. Journal of Glaciology and Geocryology,25(3): 250-255.

    Zavaleta ES, Thomas BD, Chiariello NR, et al., 2003. Plants reverse warming effect on ecosystem water balance. PNAS, 100: 9892-9893. DOI: 10.1073/pnas.1732012100.

    Zhang YS, Ohata T, Kadota T, 2003. Land surface hydrological processes in the permafrost region of the eastern Tibetan Plateau.JournalofHydrology,283:41-56.DOI: 10.1016/S0022-1694(03)00240-3.

    Zhao L, Cheng GD, Li SX, et al., 2000. Thawing and freezing processes of active layer in Wudaoliang Region of Tibetan Plateau. Chinese Science Bulletin, 45(23): 2181-2186. DOI: 10.1007/BF02886326.

    Zhou YW, Guo DX, Qiu GQ, et al., 2000. Geocryology in China. Beijing: Science Press, pp. 65-80.

    Liu GS, Wang GX, 2016. Influence of short-term experimental warming on heat-water processes of the active layer in a swamp meadow ecosystem of the Qinghai-Tibet Plateau. Sciences in Cold and Arid Regions, 8(2): 0125-0134.

    10.3724/SP.J.1226.2016.00125.

    *Correspondence to: Ph.D., GuangSheng Liu, College of Environmental Science and Engineering, Xiamen University of Technology, Xiamen, Fujian 361024, China. E-mail: liugs@xmut.edu.cn

    July 27, 2015Accepted: November 19, 2015

    日韩精品有码人妻一区| 日韩欧美国产在线观看| 黄色一级大片看看| 五月伊人婷婷丁香| 在线观看美女被高潮喷水网站| 国产成人aa在线观看| 天堂影院成人在线观看| 中文亚洲av片在线观看爽| 黄片wwwwww| 丰满的人妻完整版| 日韩制服骚丝袜av| 国产高清激情床上av| 国产精品久久久久久精品电影| 日韩制服骚丝袜av| 哪里可以看免费的av片| 久久久精品欧美日韩精品| 两个人视频免费观看高清| 一级毛片我不卡| 麻豆国产av国片精品| 欧美丝袜亚洲另类| 嘟嘟电影网在线观看| 中文字幕人妻熟人妻熟丝袜美| 免费人成在线观看视频色| 欧美+日韩+精品| 欧美日韩综合久久久久久| 国产精品无大码| 我要看日韩黄色一级片| 亚洲人与动物交配视频| 可以在线观看毛片的网站| 成年女人看的毛片在线观看| 哪里可以看免费的av片| 国产亚洲精品久久久com| eeuss影院久久| av在线老鸭窝| 九九爱精品视频在线观看| 大又大粗又爽又黄少妇毛片口| 中文字幕av在线有码专区| 九草在线视频观看| 亚洲第一区二区三区不卡| 2021天堂中文幕一二区在线观| 国产蜜桃级精品一区二区三区| 久久久久久久午夜电影| 女人十人毛片免费观看3o分钟| 日本一二三区视频观看| a级毛色黄片| av黄色大香蕉| 日本-黄色视频高清免费观看| 99热网站在线观看| 久久久国产成人免费| 亚洲在线观看片| 国产一级毛片七仙女欲春2| 亚洲自偷自拍三级| 99在线视频只有这里精品首页| 日韩欧美国产在线观看| 国产午夜福利久久久久久| 丰满人妻一区二区三区视频av| 久久午夜福利片| 99久久人妻综合| 白带黄色成豆腐渣| 亚洲欧美成人综合另类久久久 | 观看美女的网站| 欧美3d第一页| 麻豆成人av视频| 国产精品福利在线免费观看| 日本av手机在线免费观看| 成年免费大片在线观看| 亚洲色图av天堂| 日韩精品有码人妻一区| 亚洲精品影视一区二区三区av| 高清毛片免费看| 人人妻人人澡人人爽人人夜夜 | 在线a可以看的网站| 久久99热这里只有精品18| 国模一区二区三区四区视频| 啦啦啦韩国在线观看视频| 欧美一区二区精品小视频在线| 97在线视频观看| 午夜激情欧美在线| 久久久久国产网址| 99热只有精品国产| 精品无人区乱码1区二区| 亚洲成av人片在线播放无| 99国产极品粉嫩在线观看| 欧美性猛交╳xxx乱大交人| 亚洲一区高清亚洲精品| av福利片在线观看| 国产精品国产三级国产av玫瑰| 禁无遮挡网站| 青春草视频在线免费观看| 在线观看av片永久免费下载| 91在线精品国自产拍蜜月| 久久草成人影院| 九草在线视频观看| 我的女老师完整版在线观看| 51国产日韩欧美| 亚洲国产精品久久男人天堂| av天堂中文字幕网| 麻豆国产97在线/欧美| 最近2019中文字幕mv第一页| 看免费成人av毛片| 久久精品国产亚洲av香蕉五月| 日韩高清综合在线| 国产探花在线观看一区二区| 日韩人妻高清精品专区| 久久综合国产亚洲精品| 欧美极品一区二区三区四区| 一级黄色大片毛片| 精品久久久久久久久av| a级毛色黄片| 成人高潮视频无遮挡免费网站| 一区福利在线观看| 男人的好看免费观看在线视频| 精品国内亚洲2022精品成人| 深夜a级毛片| 99国产精品一区二区蜜桃av| 别揉我奶头 嗯啊视频| 国产熟女欧美一区二区| 欧美性感艳星| 老女人水多毛片| 麻豆久久精品国产亚洲av| 97在线视频观看| 夜夜夜夜夜久久久久| www.av在线官网国产| 丰满的人妻完整版| 日韩一区二区三区影片| 国产精品1区2区在线观看.| 精品国产三级普通话版| 国产精品一区www在线观看| 成人一区二区视频在线观看| av国产免费在线观看| 亚洲精品影视一区二区三区av| 久久久久久久亚洲中文字幕| 国产三级在线视频| 一级毛片我不卡| 国产精品电影一区二区三区| 亚洲成人精品中文字幕电影| 国产黄色小视频在线观看| 精品一区二区三区视频在线| 97人妻精品一区二区三区麻豆| 国产一区二区在线观看日韩| 国产精品久久久久久精品电影| or卡值多少钱| 麻豆成人午夜福利视频| 级片在线观看| 岛国毛片在线播放| av在线观看视频网站免费| 床上黄色一级片| 高清毛片免费观看视频网站| 成人亚洲精品av一区二区| 中国美白少妇内射xxxbb| 精品熟女少妇av免费看| 99久久精品热视频| 亚洲欧美日韩卡通动漫| 观看免费一级毛片| 91aial.com中文字幕在线观看| 亚洲五月天丁香| 十八禁国产超污无遮挡网站| 国产三级中文精品| 毛片一级片免费看久久久久| 亚洲人成网站在线观看播放| 久久精品夜夜夜夜夜久久蜜豆| 久久久色成人| 国内精品一区二区在线观看| 精品久久久久久久久av| 别揉我奶头 嗯啊视频| 国产女主播在线喷水免费视频网站 | 麻豆国产av国片精品| 欧美日韩国产亚洲二区| 久久精品国产亚洲av香蕉五月| 国产淫片久久久久久久久| 亚洲丝袜综合中文字幕| 午夜福利高清视频| а√天堂www在线а√下载| 成人二区视频| 日本熟妇午夜| 国产国拍精品亚洲av在线观看| 婷婷色av中文字幕| 国产免费一级a男人的天堂| 听说在线观看完整版免费高清| 国产91av在线免费观看| 亚洲精品日韩在线中文字幕 | 99视频精品全部免费 在线| 久久久久久国产a免费观看| 如何舔出高潮| 国产女主播在线喷水免费视频网站 | 亚洲乱码一区二区免费版| 性插视频无遮挡在线免费观看| 成人欧美大片| 国产私拍福利视频在线观看| 亚洲国产色片| 真实男女啪啪啪动态图| 国产精品电影一区二区三区| 免费看av在线观看网站| 99久久精品一区二区三区| av.在线天堂| 亚洲美女视频黄频| 男插女下体视频免费在线播放| 中文在线观看免费www的网站| 成人高潮视频无遮挡免费网站| eeuss影院久久| 亚洲国产精品成人综合色| 国产av一区在线观看免费| 色综合色国产| 精品一区二区免费观看| 午夜福利在线观看吧| 中出人妻视频一区二区| 欧美日韩国产亚洲二区| 欧美日本亚洲视频在线播放| 听说在线观看完整版免费高清| 99热这里只有是精品在线观看| 免费av观看视频| 免费看a级黄色片| 亚洲av一区综合| 中出人妻视频一区二区| 国内精品一区二区在线观看| 26uuu在线亚洲综合色| 欧美zozozo另类| 中文字幕制服av| 伦理电影大哥的女人| 国产精品国产三级国产av玫瑰| 少妇熟女欧美另类| 日本黄色视频三级网站网址| 在线观看美女被高潮喷水网站| 国产成人a区在线观看| 少妇人妻一区二区三区视频| 精品人妻熟女av久视频| 国产一级毛片七仙女欲春2| 狠狠狠狠99中文字幕| 欧美最黄视频在线播放免费| 久久久久久久久大av| 桃色一区二区三区在线观看| kizo精华| 91aial.com中文字幕在线观看| 男插女下体视频免费在线播放| 熟妇人妻久久中文字幕3abv| 亚洲精品粉嫩美女一区| 国产 一区 欧美 日韩| 国产一区二区三区av在线 | 人体艺术视频欧美日本| 69av精品久久久久久| 可以在线观看毛片的网站| 草草在线视频免费看| 免费在线观看成人毛片| 日本-黄色视频高清免费观看| 狠狠狠狠99中文字幕| 看片在线看免费视频| 男女做爰动态图高潮gif福利片| 一级黄片播放器| 全区人妻精品视频| 亚洲精品影视一区二区三区av| 亚洲欧美成人精品一区二区| 国产高清视频在线观看网站| 精品久久久久久久久久久久久| 国产三级在线视频| 99热6这里只有精品| 久久精品91蜜桃| 岛国在线免费视频观看| 成人永久免费在线观看视频| 久久久久久久午夜电影| 亚洲在线观看片| 精品国产三级普通话版| 人妻久久中文字幕网| 亚洲精品自拍成人| 亚洲欧美中文字幕日韩二区| 一个人看的www免费观看视频| 国产一区二区在线av高清观看| 又爽又黄a免费视频| 日韩欧美精品v在线| 成人性生交大片免费视频hd| 不卡一级毛片| 亚洲va在线va天堂va国产| 国产成人freesex在线| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 黄片wwwwww| 久久久久性生活片| 欧美+亚洲+日韩+国产| 毛片女人毛片| 精品熟女少妇av免费看| 亚洲va在线va天堂va国产| 中文字幕av在线有码专区| 91aial.com中文字幕在线观看| 久久久久久九九精品二区国产| 2021天堂中文幕一二区在线观| 综合色丁香网| 久久久国产成人精品二区| 国产毛片a区久久久久| 国产精品,欧美在线| 久久午夜亚洲精品久久| 赤兔流量卡办理| 国产一区二区在线观看日韩| 久久久久九九精品影院| 成人一区二区视频在线观看| 91aial.com中文字幕在线观看| 亚洲五月天丁香| 啦啦啦韩国在线观看视频| 最好的美女福利视频网| 夜夜爽天天搞| 欧美最黄视频在线播放免费| 久久99蜜桃精品久久| 国产成人a∨麻豆精品| 免费观看精品视频网站| 日本与韩国留学比较| 六月丁香七月| 日韩一区二区视频免费看| 五月伊人婷婷丁香| 国产片特级美女逼逼视频| 小说图片视频综合网站| 亚洲精品久久国产高清桃花| 黄色日韩在线| 免费搜索国产男女视频| 99久久九九国产精品国产免费| 欧美一区二区亚洲| 麻豆成人av视频| 国产精品爽爽va在线观看网站| 国内精品久久久久精免费| 亚洲av男天堂| 乱码一卡2卡4卡精品| 亚洲美女视频黄频| 亚洲第一电影网av| www.色视频.com| 少妇人妻精品综合一区二区 | 高清毛片免费看| 亚洲国产色片| 日韩在线高清观看一区二区三区| 国产高清激情床上av| 国产老妇伦熟女老妇高清| 欧美激情久久久久久爽电影| 亚洲av免费在线观看| 国产美女午夜福利| 高清午夜精品一区二区三区 | 免费av不卡在线播放| 一级二级三级毛片免费看| 老师上课跳d突然被开到最大视频| 国内精品美女久久久久久| 丰满人妻一区二区三区视频av| 成人性生交大片免费视频hd| 日日摸夜夜添夜夜添av毛片| 18禁黄网站禁片免费观看直播| 国产成人影院久久av| 波多野结衣巨乳人妻| 看非洲黑人一级黄片| 在线观看午夜福利视频| 国产高清三级在线| 村上凉子中文字幕在线| 超碰av人人做人人爽久久| 99热这里只有精品一区| 亚洲经典国产精华液单| 精品久久久久久久久av| 老司机福利观看| 不卡一级毛片| 亚洲av第一区精品v没综合| 久久久精品大字幕| 亚洲成人中文字幕在线播放| 久久午夜亚洲精品久久| 久久鲁丝午夜福利片| 亚洲av成人av| 99久国产av精品| 悠悠久久av| 亚洲国产精品成人综合色| 免费一级毛片在线播放高清视频| eeuss影院久久| 成人av在线播放网站| 老熟妇乱子伦视频在线观看| 嫩草影院入口| 99视频精品全部免费 在线| 日韩欧美精品免费久久| 永久网站在线| 好男人在线观看高清免费视频| 日本成人三级电影网站| 春色校园在线视频观看| 亚洲成人久久性| 国产精品久久久久久精品电影小说 | 两个人的视频大全免费| 三级毛片av免费| 亚洲精品日韩在线中文字幕 | 久久这里只有精品中国| 观看免费一级毛片| 久久精品影院6| 国产精品电影一区二区三区| 亚洲最大成人av| 中国美白少妇内射xxxbb| 18禁在线播放成人免费| 九草在线视频观看| 能在线免费看毛片的网站| 熟女电影av网| 99久久精品国产国产毛片| 寂寞人妻少妇视频99o| 久久久久网色| 亚洲精品日韩在线中文字幕 | 99久久精品热视频| 丰满人妻一区二区三区视频av| 亚洲美女搞黄在线观看| 国产精品国产三级国产av玫瑰| 亚洲av成人精品一区久久| 亚洲熟妇中文字幕五十中出| 精品久久久久久久人妻蜜臀av| 国产精品爽爽va在线观看网站| 久久亚洲国产成人精品v| 免费大片18禁| 99久久九九国产精品国产免费| 六月丁香七月| 中文亚洲av片在线观看爽| 亚洲最大成人手机在线| 精品人妻偷拍中文字幕| 麻豆av噜噜一区二区三区| 国产真实伦视频高清在线观看| 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看| 乱码一卡2卡4卡精品| av又黄又爽大尺度在线免费看 | 日本熟妇午夜| 成人无遮挡网站| 亚洲国产高清在线一区二区三| 1024手机看黄色片| 久久久午夜欧美精品| 国产单亲对白刺激| 男女做爰动态图高潮gif福利片| 亚洲天堂国产精品一区在线| 亚洲精品粉嫩美女一区| 欧美人与善性xxx| 日韩av在线大香蕉| 亚洲中文字幕日韩| 国产麻豆成人av免费视频| 蜜臀久久99精品久久宅男| 国产一级毛片在线| 久久久久久伊人网av| 国产日韩欧美在线精品| 噜噜噜噜噜久久久久久91| 国产探花极品一区二区| 色播亚洲综合网| 国产老妇女一区| 老司机影院成人| 亚洲精品日韩在线中文字幕 | 免费不卡的大黄色大毛片视频在线观看 | 在现免费观看毛片| 国产精品久久久久久亚洲av鲁大| 欧美成人免费av一区二区三区| 日韩欧美一区二区三区在线观看| 欧美最黄视频在线播放免费| 久久综合国产亚洲精品| 在线观看66精品国产| 1024手机看黄色片| 国产精品福利在线免费观看| 直男gayav资源| 亚州av有码| 真实男女啪啪啪动态图| 99热网站在线观看| 亚洲久久久久久中文字幕| 一进一出抽搐gif免费好疼| 欧美一区二区精品小视频在线| 日韩欧美精品v在线| 高清毛片免费观看视频网站| 久久精品人妻少妇| 如何舔出高潮| 在线免费观看的www视频| 成人欧美大片| 亚洲av第一区精品v没综合| 在线a可以看的网站| 好男人在线观看高清免费视频| 欧美又色又爽又黄视频| 精品99又大又爽又粗少妇毛片| 国产私拍福利视频在线观看| 国产在线男女| 99国产精品一区二区蜜桃av| 久久久久国产网址| 日本爱情动作片www.在线观看| 精品少妇黑人巨大在线播放 | 久久精品国产99精品国产亚洲性色| 日韩国内少妇激情av| 波多野结衣高清作品| 一进一出抽搐动态| 九九久久精品国产亚洲av麻豆| 久久久久久久午夜电影| 男女那种视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 最近手机中文字幕大全| 搞女人的毛片| 国产一区二区在线av高清观看| 又爽又黄a免费视频| 亚洲欧洲国产日韩| 久久久国产成人精品二区| 欧美又色又爽又黄视频| av女优亚洲男人天堂| 欧美区成人在线视频| 久久午夜福利片| 国产高清不卡午夜福利| 久久韩国三级中文字幕| 亚洲中文字幕一区二区三区有码在线看| 在线免费观看不下载黄p国产| 看非洲黑人一级黄片| 99热这里只有是精品50| 一级毛片久久久久久久久女| 熟女电影av网| 人妻少妇偷人精品九色| 少妇被粗大猛烈的视频| 精品久久久久久久久久免费视频| 免费大片18禁| av.在线天堂| a级毛色黄片| 丰满的人妻完整版| 日韩欧美在线乱码| 中国美女看黄片| 亚洲精品乱码久久久久久按摩| 校园春色视频在线观看| 天堂√8在线中文| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看的影片在线观看| 国产亚洲5aaaaa淫片| 99久久人妻综合| 欧美一区二区国产精品久久精品| 人妻制服诱惑在线中文字幕| 在现免费观看毛片| av免费观看日本| 99热这里只有是精品50| 搡老妇女老女人老熟妇| 身体一侧抽搐| 老司机福利观看| 国产国拍精品亚洲av在线观看| 搞女人的毛片| 在现免费观看毛片| www日本黄色视频网| 成人美女网站在线观看视频| 一夜夜www| 日本黄大片高清| 中出人妻视频一区二区| 天堂√8在线中文| 少妇熟女aⅴ在线视频| 国产老妇伦熟女老妇高清| 精华霜和精华液先用哪个| 啦啦啦观看免费观看视频高清| 亚洲av第一区精品v没综合| 三级国产精品欧美在线观看| 91精品国产九色| 天堂√8在线中文| 黄色一级大片看看| 可以在线观看的亚洲视频| 色噜噜av男人的天堂激情| 啦啦啦韩国在线观看视频| 69av精品久久久久久| 精品人妻偷拍中文字幕| 久久久久久久久久久丰满| 综合色av麻豆| 全区人妻精品视频| 男人狂女人下面高潮的视频| 久久精品影院6| 少妇高潮的动态图| 亚洲人与动物交配视频| 国产成人91sexporn| 小蜜桃在线观看免费完整版高清| 国产高清有码在线观看视频| 久久久久久九九精品二区国产| 欧美高清性xxxxhd video| 97在线视频观看| 在线a可以看的网站| 夜夜夜夜夜久久久久| 国产成年人精品一区二区| 人妻制服诱惑在线中文字幕| 精品人妻熟女av久视频| 日韩欧美一区二区三区在线观看| 日日摸夜夜添夜夜添av毛片| 99久久精品国产国产毛片| 免费观看a级毛片全部| 欧美潮喷喷水| videossex国产| 欧美日本亚洲视频在线播放| 精品久久久噜噜| av免费在线看不卡| 日本熟妇午夜| 国产美女午夜福利| 高清日韩中文字幕在线| 欧美人与善性xxx| av女优亚洲男人天堂| 草草在线视频免费看| 日韩成人伦理影院| 国产精品嫩草影院av在线观看| 少妇丰满av| 久久久欧美国产精品| 美女cb高潮喷水在线观看| 精华霜和精华液先用哪个| 日韩一区二区视频免费看| 中国美女看黄片| 深夜a级毛片| 99久久精品一区二区三区| 美女内射精品一级片tv| 国产三级在线视频| 联通29元200g的流量卡| 精品久久久久久久末码| 欧美日本亚洲视频在线播放| 联通29元200g的流量卡| 午夜福利在线观看吧| 亚洲欧美日韩无卡精品| 国产欧美日韩精品一区二区| 午夜老司机福利剧场| 免费观看a级毛片全部| 精品久久久久久久人妻蜜臀av| 一本久久中文字幕| 少妇被粗大猛烈的视频| 18+在线观看网站| 久久午夜亚洲精品久久| 婷婷六月久久综合丁香| 色5月婷婷丁香| 国产人妻一区二区三区在| 日日撸夜夜添| 国产精品乱码一区二三区的特点| 亚洲av免费高清在线观看| 国产高清激情床上av| 欧美+日韩+精品| 最近2019中文字幕mv第一页| 99热只有精品国产| 国产白丝娇喘喷水9色精品| a级毛色黄片| 亚洲精品粉嫩美女一区| 国产亚洲精品久久久久久毛片| 午夜精品国产一区二区电影 |