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

    The effects of aeration and irrigation regimes on soil CO2 and N2O emissions in a greenhouse tomato production system

    2018-02-05 07:10:56CHENHuiHOUHuijingWANGXiaoyunZHUYanQaisarSaddiqueWANGYunfeiCAIHuan
    Journal of Integrative Agriculture 2018年2期

    CHEN Hui, HOU Hui-jing, WANG Xiao-yun, ZHU Yan, Qaisar Saddique, WANG Yun-fei, CAI Huan

    jie1

    1 College of Water Resources and Architectural Engineering, Northwest A&F University/Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas, Ministry of Education, Yangling 712100, P.R.China

    2 School of Hydraulic, Energy and Power Engineering, Yangzhou University, Yangzhou 25127, P.R.China

    1. Introduction

    Carbon dioxide (CO2) and nitrous oxide (N2O) are significant greenhouse gases (GHGs) that contribute to stratospheric ozone layer depletion and global warming, with their concentrations increasing by 40 and 20%, respectively,since global industrialization (IPCC 2013). Worldwide,anthropogenic sources of CO2and N2O are dominated by agriculture and have increased by approximately 17% from 1990 to 2005 (Robertson and Grace 2004; Forsteret al.2007). China accounts for 45% of the world’s total vegetable production, occupying 11.6% of the country’s cultivated land(FAOSTAT 2009). Unlike the conventional staple grain crops(i.e., rice, wheat, and corn), vegetable cropping systems are subjected to high nitrogen (N) fertilizer inputs, frequent irrigation, high temperatures, and multiple harvests within a year. For instance, greenhouse vegetable cropping systems in Beijing suburbs utilize N fertilization rates of approximately 1 000 kg N ha–1crop–1(Chenet al.2004), over 2 800 kg N ha–1yr–1in Huiming, Shandong Province (Juet al.2006) and on average 600 kg N ha–1crop–1in the suburban area of Xi’an City, Shaanxi Province (Zhouet al.2006). These values are much greater than the average national N fertilizer rate of approximately 180 kg N ha–1yr–1or 120 kg N ha–1crop–1(Zhu and Chen 2002) that has been the cause of serious pollution (Sunet al.2017). Furthermore, the constant irrigation of these vegetable crops plays an increasingly important role in Chinese agriculture (Sunet al.2013).Irrigation influences the soil organic matter decomposition,microbial biomass and activity, root biomass, soil aeration and nitrogen turnover, which influence soil CO2and N2O production and emission (Rethet al.2005; Huanget al.2007; Scheeret al.2008, 2013). Excessive N fertilization along with abundant soil moisture undoubtedly generates high CO2and N2O emissions measured from vegetable cropping systems. Thus, continuous field measurements of CO2and N2O emissions from vegetable cropping systems are necessary to fully understand crop-specific CO2and N2O emissions (Flynnet al.2005).

    Although CO2and N2O emissions from Chinese grain crops have been extensively documented, research on CO2and N2O emissions from vegetable systems in China are sparse and lead to imprecise values for national CO2and N2O emissions inventory from Chinese croplands(Heet al.2009). Previous studies on GHGs in vegetablefields have extensively analyzed fertilizer application(Riyaet al.2012; Denget al.2013; Liet al.2015a).However, the effect of water management coupled with GHGs in vegetable fields remains uncertain, especially when considering aerated irrigation (AI) practices. As early as the 1994s, AI application has and continues to be employed in root zone soils to overcome hypoxia.The benefits of this technique have been verified in crops including pineapple (Chenet al.2011), soybean and cotton(Bhattarai and Midmore 2009), cucumber (Ehretet al.2010; Niuet al.2013), tomato (Bhattaraiet al.2006; Houet al.2016) and rice (Zhuet al.2012). The majority of AI studies have concentrated on the economic benefits (e.g.,crop production and water use efficiency) and overall fruit quality (Bhattarai and Midmore 2009; Ehretet al.2010;Abuarabet al.2013; Niuet al.2013). Very few studies on AI-derived soil properties (i.e., soil aeration, NO3–and electrolytic conductivity (EC)) have been performed (Niuet al.2012; Ben-Noah and Friedman 2016), especially for soil CO2and N2O emissions. Soil respiration caused by AI has been studied in wheat, cotton and pineapple (Chenet al.2011) but the effect of AI on vegetable-planted soils is insufficiently documented. Furthermore, variations of N2O emissions under AI have been demonstrated in the laboratory (Hwang and Hanaki 2000; Castro-Barroset al.2015) while field experiments are extremely limited.During our previous study, only gas emissions and soil water contents were measured after AI treatment during the autumn-winter season, and the experiment did not capture emissions at the tomato seedling stage and part of the blooming and fruit setting stage (Houet al.2016).Unfortunately, little information was gathered regarding factors influenced by AI and the relationship between gas fluxes and environmental variables. The effects of soil aeration and irrigation regimes on GHGs were not recorded. Thus, a study on the combined effects of soil aeration and irrigation regimes on CO2and N2O emissions from vegetable fields is important to perform a comprehensive assessment and to develop greenhouse gas emission reduction measures.

    In this study, the closed chamber and gas chromatography technique was employed to measure CO2and N2O emissions from the soil of a greenhouse tomato production system in Northwest China. The soil temperature, soil water content, soil organic carbon, soil nitrate content and electrolytic conductivity were monitored concurrently. The objectives of this study were to (1) evaluate soil CO2and N2O emissions under different aeration and irrigation regimes and to (2) discern the primary factors influencing soil CO2and N2O fluxes.

    2. Materials and methods

    2.1. Experimental site

    A field experiment was conducted from 4 April 2015 to 17 January 2016 in a solar greenhouse at the Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of the Ministry of Education, Northwest A&F University, Yangling, Shaanxi Province of China (34°20′N,108°04′E). The greenhouse was 36 m×10.3 m×4 m in dimension. The site was in a semi-humid and liable-drought climate zone with an annual average sunshine duration of 2 163.8 h and frost-free period of 210 d. The experimental site had a Lou soil type in a brown soil zone. The soil texture was silt clay loam and the groundwater depth was at least 100 m below the surface. Physical properties of the top 20 cm soil layer are shown in Table 1.

    Table 1 Physical properties of soil in testing site

    2.2. Field management

    Experimental designTwo irrigation regimes were carried out in the experiment, full irrigation (1.0W) and deficit irrigation (0.6W). Here, W represented the irrigation volume of sufficient water supply which was calculated according to Houet al.(2016). Each irrigation regime contained aeration and non-aeration. Hence, four treatments with three replications were performed with a completely randomized design: (1) aerated deficit irrigation (AI1), (2) non-aerated deficit irrigation (CK1), (3) aerated full irrigation (AI2) and(4) non-aerated full irrigation (CK2).

    Each plot (size: 4.0 m×0.8 m) represented one replication.Eleven plants were planted in each plot with a spacing of 35 cm. During the experimental period, all plots were covered with plastic films. Subsurface drip irrigation was used with drip irrigation tape buried at a depth of 15 cm and a dripper spacing of 35 cm. The Mazzei air injector model 287(a Venturi, Mazzei Injector Company, America) was installed for aeration at the head of each irrigation pipe with an inlet and outlet pressure of 0.1 and 0.02 Mpa, respectively, and set to inject 17% air by volume of water (Zhuet al.2016).Vegetable cropHumus pots of tomato seedling (variety:Feiyue) were adopted in the greenhouse for two consecutive growing seasons. In the spring-summer season, tomatoes were transplanted on 4 April 2015 and harvested on 29 July 2015. In the autumn-winter season, transplantation occurred on 30 August 2015 and the harvest occurred on 17 January 2016.

    Fertilization and irrigationThroughout the tomato growing stage, only base fertilizer was applied in the spring-summer season, which was comprised of organic fertilizer (N-P2O5-K2O≥10%, organic matter≥45%) and compound fertilizer(total nutrients≥45%, including N, P2O5and K2O each at 15%). They were applied on 22 January 2015 at a rate of 5 937.5 and 3 437.5 kg ha–1, respectively. During the autumnwinter season, they were applied on 5 August 2015 at a rate of 3 437.5 and 2 187.5 kg ha–1, respectively. However,the topdressing (spraysven potassium sulfate, K2O≥52%,SO3≥46%) in the autumn-winter season was applied at a rate of 236.8 kg ha–1on 27 October, 2015 and on 22 October,2015 for deficit irrigation and full irrigation, respectively.

    An E601 evaporation pan was placed in the greenhouse to determine the irrigation volume as performed by Houet al.(2016). Irrigation occurred at 8:00 a.m. Throughout the tomato crop growth, with the exception of transplanting irrigation, the irrigation treatments were applied 27 and 35 times at an interval of 3–5 days in the spring-summer season and autumn-winter season, respectively. The total irrigation volume in the spring-summer season for deficit irrigation and full irrigation was 180.1 and 300.1 mm,respectively. In the autumn-winter season, the value for deficit irrigation and full irrigation was 116.2 and 194.6 mm,respectively.

    2.3. Sampling and measurements

    The CO2and N2O fluxes were simultaneously measured in three replications for each treatment by using the static closed chamber method (Houet al.2016). The static chamber with 6-mm thick PVC material covered an area of 25 cm×25 cm, with a height of 25 cm. The bases of the chambers, also made of PVC, were installed in the middle of each plot on the day of transplanting and remained there until tomato harvest. A 3-cm-deep groove on the top edge of the bottom layer and on the base of the chamber was to be filled with water to seal the rim of the chamber.All chambers were wrapped with layers of sponge and aluminum foil to minimize air temperature changes within the chamber that could be caused by sunlight during gas sampling. The inside of each chamber was equipped with an electric fan near the top for mixing air. Gas samples forflux measurements were collected between 10 to 11 a.m. at an interval of approximately 6 and 9 days during the early and late growing stages, respectively. Chamber air samples were collected using a 50-mL syringe following 0, 10, 20 and 30 min after the chambers were placed on pre-fixed plastic frames. A 30-mL air sample was drawn each time with a syringe. Gas samples in the syringes were analyzed within a few hours. Sample sets were discarded unless they yielded anR2linear regression value greater than 0.85.

    Gas samples were analyzed for CO2and N2O concentrations using a gas chromatograph (7890A GC System, Agilent Technologies, America) equipped with aflame ionization detector (FID) and an electron capture detector (ECD) for measuring CO2and N2O, respectively.Fluxes of CO2and N2O from the greenhouse tomato production soils were calculated based on the equation given by Houet al.(2016).

    The soil/air temperatures were measured when air samples were collected. The soil temperature was measured at a depth of 10 cm near the chamber bases using a geothermometer (RM-004). The temperature inside the chamber was measured using mercury thermometers(WNG-01) and the air temperature outside the chamber was recorded using a mercury thermometer (WNG-01) placed 1.5 m above the ground.

    Soil sample cores from 0 to 10 cm were taken between two plants at the head, middle and end of each plot through a diameter gauge at an interval of 8, 18, 10 and 13 days on average for soil water-filled pore space (WFPS), soil organic carbon (SOC), nitrate (NO3–) content and EC,respectively (Chenet al.2016). Part of the fresh soil was used to measure the soil water contentviaoven-drying at 105°C for 12 h and then converted to WFPS using methods by Dinget al.(2007). To determine NO3–, fresh soil samples were mixed thoroughly and sieved through a 5-mm mesh.Subsequently, each fresh soil sample (corresponding to 5 g dry soil) was extracted by 50 mL of 2 mol L–1KCl for 30 min. The NO3–contents were analyzed using an ultraviolet spectrophotometer (EV300PC, Thermo Fisher, England)following the methods of Zhanget al.(2015). Another dried soil sample was passed through a 1-mm sieve to measure SOC using the oil bath heated potassium dichromate oxidation volumetric method (Zheng 2013) and the EC was determined in a 1:5 of soil:H2O (g mL–1) suspension using an EC meter (FE30K Plus, METTLER TOLEDO, China).

    All mature tomato fruits in each plot were manually harvested and individual fruit yields were determined at each picking.

    2.4. Data analyses

    All statistical analyses, including correlations, were computed using SPSS Statistics 22.0. Figures were plotted in Sigmaplot 12.5. A Pearson correlation test was used to further explain the relationships between gas fluxes (CO2and N2O) and environmental variables (i.e., soil temperature,WFPS, SOC, NO3–, and EC).

    3. Results and discussion

    3.1. Tomato yield

    The soil aeration and irrigation regimes influenced tomato yields, with spring-summer season yields ranging from 36.83 to 55.10 t ha–1and autumn-winter season yields ranging from 21.77 to 33.17 t ha–1(Table 2). Aeration under each irrigation regime over the two seasons increased tomato yields compared to the non-aeration treatments and the difference was significant under the full irrigation treatment (P<0.05).This was attributed to the oxygation benefit to enhance availability of dissolved oxygen to the root system (Bhattarai and Midmore 2009), enhancing crop yields for different plant species (Bhattarai and Midmore 2009; Abuarabet al.2013; Niuet al.2013). Moreover, full irrigation during the two seasons enhanced crop yields significantly compared to the deficit irrigation treatments (P<0.05, Table 2). Similarfindings have also been reported in other studies (Patanè and Cosentino 2010; Zhanget al.2017). These results were determined to be an effect of water stress on plant cells’ multiplication and expansion throughout all growth stages (Zhanget al.2017), in turn influencing the crop yield.

    Variance analysis on aeration and irrigation treatments showed that both had a significant effect on total tomato yields (P=0.000 for irrigation,P=0.021 for aeration), while their interaction effect on total crop yield was not significant(P=0.881, Table 2).

    3.2. CO2 and N2O emissions from greenhouse tomato production soils

    Seasonal dynamics and cumulative CO2 emissionsSoil CO2fluxes over all treatments showed fluctuated patterns(Fig. 1-A) that varied from 8.52 to 518.76 mg m–2h–1and from 91.80 to 495.75 mg m–2h–1during the spring-summer season and autumn-winter season, respectively. The maximum total cumulative CO2emission from the soil was 14 836.07 kg ha–1for the AI2 treatment, resulting in a 4.2,29.7 and 20.4% increase compared to AI1, CK1 and CK2 treatments, respectively (Table 2).

    In contrast to non-aerated soils, soil aeration under different irrigation regimes increased the cumulative CO2emissions for each of the two seasons (Table 2) and the difference was found to be significant in the autumn-winter season (P<0.05). Other studies have seen this effect where soil respiration has been shown to increase under aerated irrigation (Chenet al.2011; Houet al.2016).However, throughout typical growing seasons, only a few measurements are conducted and thus some flux peaks might be missed. Such few measurements fail to estimate the variations of greenhouse gases accurately and could severely underestimate cumulative gas emissions (Liuet al.2010), providing inaccuracies in the CO2emissions national inventory from Chinese vegetable cropping systems.Hence, the static closed chamber and gas chromatography technique were used to systematically and continuously study the effect of aerated irrigation on greenhouse gas emissions from greenhouse tomato production soils.Previous studies had found that crop root and soil microbial respiration were the main sources of soil CO2emissions.The higher soil CO2emissions under aerated irrigation were potentially resulted from the effect of aeration on increased oxygen (Guadieet al.2014), soil microbial/enzyme activity(Liet al.2015b) and root respiration (Bhattaraiet al.2008).In addition, a slight higher soil temperature under aerated irrigation may also be responsible for the higher CO2emissions compared to non-aerated irrigation (Fig. 2).

    ?

    Fig. 1 Variation of CO2 fluxes (A) and N2O fluxes (B) for different treatments from greenhouse tomato production soils over two seasons. AI1, aerated deficit irrigation; CK1, non-aerated deficit irrigation; AI2, aerated full irrigation; CK2, nonaerated full irrigation. Solid arrows and dotted arrows denote time of topdressing for full irrigation and deficit irrigation, respectively. Bars indicate standard errors of three replications.

    In our study, the cumulative CO2emissions under the full irrigation regime during the two seasons were slightly higher than those under deficit irrigation(P>0.05, Table 2), which was consistent with other studies (Scheeret al.2013;Zornozaet al.2016). Several explanations may be given for this conclusion.First, deficit irrigation has been shown to reduce soil organic matter breakdown and enhanced SOC accumulation (Zornozaet al.2016). Second, the application of higher water quantities has been shown to increase soil enzyme activity (Zhang and Wang 2006; Wanet al.2008), serving as an indicator of microbial activity and soil biochemical intensity (Mersi and Schinner 1991).

    Variance analysis on aeration and irrigation showed that aeration had a significant effect on total cumulative CO2emissions (P=0.021), while the effect of irrigation on total cumulative CO2emissions was not significant (P=0.447,Table 2). Meanwhile, the aeration and irrigation did not have a significant interaction effect on total cumulative CO2emissions (P=0.881).Seasonal dynamics and cumulative N2O emissions The soil N2O fluxes showed clear seasonal patterns (Fig. 1-B). In the spring-summer season,higher N2O fluxes were observed before 30 April, followed by relatively low N2O flux levels. In the autumn-winter season, soil N2O fluxes slowly decreased before 15 October. Peaks that occurred around 1 November 2015 were strongly ascribed to the topdressing. Additionally, soil N2O fluxes changed minimally after 10 November. Similar patterns were observed by previous researchers who concluded that fertilization coupled with wetted soil resulted in extremely high N2O emissions (Heet al.2009; Liet al.2015a; Zhanget al.2015). This short peak emission and dissipation phenomenon can be strongly attributed to the fact that soil N2O fluxes have been shown to be affected by the life cycle of plants with more and more soil available N taken up by plants as plants grew up gradually (Xiaet al.2013).

    Soil N2O fluxes over all treatments ranged from 2.90 to 311.77 μg m–2h–1and from 5.47 to 136.59 μg m–2h–1over the two seasons, respectively (Fig. 1-B). The highest total cumulative N2O emission was 1.78 kg ha–1measured in the AI2 treatment, a 61.8, 95.6 and 11.3% increase compared to AI1, CK1 and CK2, respectively (Table 2).

    In general, N2O is primarily produced during soil nitrification and denitrification processes (Paulet al.1993;Skiba and Smith 1993), which are highly dependent not only on the water regime, soil aeration and soil temperature but also on fertilizer inputs (Weieret al.1993). The soil oxygen content variation caused by aerated irrigation must inevitably affect the condition of nitrification and denitrification(Hwang and Hanaki 2000; Liikanen and Martikainen 2003;Bhattaraiet al.2006) and has been seen to influence soil N2O emissions (Houet al.2016). In the present study,we found that there was no significant difference between aeration and non-aeration in the cumulative N2O emissions during the spring-summer season (P>0.05, Table 2). This might be because the soil moisture throughout the growing period was relatively low most of the time. The increased soil porosity and the reduced soil moisture due to aeration were bound to reduce soil N2O emissions. Moreover, N2O emissions in the autumn-winter season increased with increasing oxygen (P<0.05, Table 2), which was similar to previous researches (Liikanen and Martikainen 2003;Houet al.2016). Compared to the control treatment,higher WFPS under aerated irrigation (57.6%vs. 57.0%on average, Fig. 2) during this period was beneficial to N2O production. However, the effects of the aeration regimes on N2O emissions predominantly focused on laboratory experiments (Hwang and Hanaki 2000; Castro-Barroset al.2015) and field experiments regarding the effects of different modes of aeration and irrigation combinations on soil N2O emissions have not been reported. In this study,the Mazzei air injector was used for soil aeration and the lack of knowledge about aerated irrigation technology practices could attribute to the differences observed.

    The full irrigation treatment over the two seasons showed significantly increased cumulative N2O emissions relative to the deficit irrigation treatment (P<0.05, Table 2),which was similar to the results reported by Scheeret al.(2013). The full irrigation treatment could have prompted the denitrification processes causing hypoxic conditions that in turn promoted elevated emissions of N2O from the soil. However, the significance of the treatment effects on N2O emissions varied and could be a result of different experimental conditions. Scheeret al.(2013) showed that under all treatments, the highest N2O emissions occurred following heavy rainfall, overriding the effect of irrigation treatments. The present study was conducted in a controlled greenhouse setting and was not skewed by rainfall events.

    Variance analysis on aeration and irrigation showed that irrigation had a significant effect on total cumulative N2O emissions (P=0.000), while the effect of aeration on total cumulative N2O emissions was not significant (P=0.102,Table 2). Meanwhile, the aeration and irrigation did not have a significant interaction effect on total cumulative N2O emissions (P=0.931).

    Fig. 2 Seasonal changes of soil/air temperature (A) and water-filled pore space (WFPS, B) under different treatments in a greenhouse tomato production system. AI1, aerated deficit irrigation; CK1, non-aerated deficit irrigation; AI2, aerated full irrigation; CK2, non-aerated full irrigation; Ta, air temperature.Bars indicate standard errors of three replications.

    3.3. Environmental variables influencing CO2 and N2O emissions

    Soil temperature and WFPSThe soil temperature showed clear seasonal dynamics, increasing as the air temperatures increased (Fig. 2-A). Aeration and full irrigation increased the average soil temperature slightly, but not significantly compared to non-aeration and deficit irrigation over the two seasons (P>0.05).

    The WFPS was slightly influenced by aeration and irrigation regimes (P>0.05, Fig. 2-B). Aeration has been shown to affect WFPS through: (1) changing the soil structure due to the shrinking and movement of soil particles(Ben-Noah and Friedman 2016), inducing evaporation; (2)encouraging soil moisture uptake by the root, a process of water consumption; and (3) affecting animal, microbial and root respiration (Bhattaraiet al.2008), producing water.

    Soil organic carbon, nitrate and ECSimilar patterns of SOC during each season were observed over all treatments(Fig. 3-A). Aeration and irrigation regimes did not have a significant influence on SOC (P>0.05). SOC fluctuated from 7.43 to 9.56 g kg–1(average=8.27 g kg–1) for the AI1 treatment, from 7.31 to 8.80 g kg–1(average=8.10 g kg–1) for the CK1 treatment, from 7.75 to 8.90 g kg–1(average=8.32 g kg–1) for the AI2 treatment and from 7.51 to 8.57 g kg–1(average=8.17 g kg–1) for the CK2 treatment during the two seasons.

    No significant differences in NO3–content were found between the aeration and irrigation regimes (P>0.05,Fig. 3-B). The mean value of the NO3–contents during the two seasons for AI1, CK1, AI2 and CK2 were 118.53, 122.34 112.64 and 113.13 mg kg–1, respectively.

    EC varied slightly throughout the growing tomato period(Fig. 4). The treatment effects between aeration and irrigation regimes on EC during the two seasons were insignificant (P>0.05). In the 0–10 cm layer, the value of EC during the two seasons was, on average, 1.90, 1.87, 1.98 and 1.92 mS cm–1for AI1, CK1, AI2 and CK2, respectively.

    3.4. Correlation of gas fluxes with environmental variables

    Correlation of CO2 with environmental variablesThe Pearson correlation test was used to further explain the relationships between the gas fluxes (CO2and N2O) and environmental variables (i.e., soil temperature, WFPS, SOC,NO3–and EC) (Table 3). The CO2flux was significantly and positively correlated with soil temperature under different irrigation modes (Table 3, Fig. 5). This relationship between CO2emissions and soil temperature has been demonstrated to be a result of increased temperatures causing increased microorganism activities and organic matter decomposition(Tonget al.2015). Similar results were obtained at Three Gorges Reservoir area (Iqbalet al.2010), in California’s Central Valley tomato fields (Kallenbachet al.2010), in a Chinese mountainous area (Liet al.2008), as well as in the humid temperate grasslands of Japan (Wanget al.2009).Furthermore, a significant positive correlation between the CO2flux and WFPS under aerated irrigation was observed(Table 3), calling for further study to investigate the correlations between the CO2flux and other environmental variables.

    Fig. 3 Seasonal changes of soil organic carbon (A) and soil NO3– content (B) under different treatments in a greenhouse tomato production system. AI1, aerated deficit irrigation; CK1, non-aerated deficit irrigation; AI2, aerated full irrigation; CK2, non-aerated full irrigation. Bars indicate standard errors of three replications.

    Correlation of N2O with environmental variablesThe dependence of soil N2O fluxes on WFPS and NO3–showed a positive exponential correlation under all treatments(Table 3, Fig. 6). This trend was in alignment with numerous studies. For example, in an intensively managed vegetable cropping system, Zhanget al.(2015) found that N2O flux was significantly affected by WFPS and the NO3–content.Weslienet al.(2012) reported that N2O fluxes were positively correlated with WFPS (R=0.50) in a carrot cropping system.These findings suggested that nitrogen fertilizer application and moisture are the main factors influencing soil N2O emissions (Xuet al.2004; Gaoet al.2014). Moreover, the highest N2O flux in this study was observed at a WFPS of 58.5 and 69.3% under aerated irrigation and non-aerated irrigation, respectively (Fig. 6-A). This was in accordance with the previous findings reporting that N2O production has a maximum at WFPS of roughly 60% (Linn and Doran 1984;Schmidtet al.2000; Gaoet al.2014). Once, WFPS reaches more than 60%, the availability of O2and CO2substrates for nitrification declines due to severely restricted diffusion rates (Davidson and Schimel 1995). However, these values were higher than previous results of approximately 50%WFPS (Dinget al.2007; Houet al.2016). This could be attributed to fertilization and irrigation field management during transplantation. Abundant substrates based on base fertilizer and large amounts of water irrigated into the soil during transplantation were conducive to gas emissions.

    Fig. 4 Seasonal changes of electrolytic conductivity (EC)for different treatments. AI1, aerated deficit irrigation; CK1,non-aerated deficit irrigation; AI2, aerated full irrigation; CK2,non-aerated full irrigation. Bars indicate standard errors of three replications.

    Fig. 5 Correlations between soil CO2 fluxes and soil temperature under aerated irrigation and non-aerated irrigation among the two seasons. T, soil temperature. * indicates significance at the 0.05 level.

    Table 3 Pearson correlation coefficients between gas fluxes and environmental variables under different irrigation modes1)

    The N2O fluxes from the soil as a function of the soil temperature at a depth of 10 cm was normally distributed,with 33.6% of the total N2O emitted within 16–23°C under all treatments (Fig. 7). Similar results were obtained by Schmidtet al.(2000) and Gaoet al.(2014). Furthermore, N2O fluxes peaked under all treatments when the soil temperature was approximately 18°C (Fig. 7). An exponential positive correlation was observed between the soil N2O fluxes and soil temperature when the soil temperature was below 18°C,while a linear negative correlation was observed when the soil temperature exceeded 18°C (P<0.01).

    Fig. 6 Correlations between soil N2O fluxes and WFPS (A) and between soil N2O fluxes and NO3– (B) under aerated irrigation and non-aerated irrigation among the two seasons. AI, aerated irrigation; CK, non-aerated irrigation. ** indicates significance at the 0.01 level.

    Fig. 7 Correlations between N2O fluxes and soil temperature under aerated irrigation (A) and non-aerated irrigation (B) among the two seasons. T, soil temperature. ** indicates significance at the 0.01 level.

    4. Conclusion

    Our results showed that aeration under full irrigation over the two seasons significantly increased tomato yields by 21.9%on average. Soil CO2and N2O emissions were significantly increased by aeration during the autumn-winter season. Full irrigation over the two seasons increased tomato yields and soil N2O emissions significantly. Moreover, soil temperature was the primary factor influencing CO2emissions, while soil temperature, moisture and NO3–were the primary factors influencing N2O emissions. Our results provide a theoretical foundation and scientific basis for evaluating aerated irrigation techniques on farmland soil ecology. Irrigation coupled with particular soil aeration practices may offer a desirable balance between crop production yields and greenhouse gas mitigation in greenhouse vegetable fields.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (51309192), the National Key Research and Development Program of China (2016YFC0400201)and the Fundamental Research Funds for the Central Universities, China (Z109021510).

    Abuarab M, Mostafa E, Ibrahim M. 2013. Effect of air injection under subsurface drip irrigation on yield and water use efficiency of corn in a sandy clay loam soil.Journal of Advanced Research, 4, 493–499.

    Ben-Noah I, Friedman S P. 2016. Aeration of clayey soils by injecting air through subsurface drippers: Lysimetric andfield experiments.Agricultural Water Management, 176,222–233.

    Bhattarai S P, Midmore D J. 2009. Oxygation enhances growth,gas exchange and salt tolerance of vegetable soybean and cotton in a saline vertisol.Journal of Integrative Plant Biology, 51, 675–688.

    Bhattarai S P, Midmore D J, Pendergast L. 2008. Yield, wateruse efficiencies and root distribution of soybean, chickpea and pumpkin under different subsurface drip irrigation depths and oxygation treatments in vertisols.Irrigation Science, 26, 439–450.

    Bhattarai S P, Pendergast L, Midmore D J. 2006. Root aeration improves yield and water use efficiency of tomato in heavy clay and saline soils.Scientia Horticulturae, 108, 278–288.

    Castro-Barros C M, Daelman M R, Mampaey K E, van Loosdrecht M C, Volcke E I. 2015. Effect of aeration regime on N2O emission from partial nitritation-anammox in a full-scale granular sludge reactor.Water Research,68, 793–803.

    Chen H, Hou H J, Cai H J, Zhu Y, Wang C. 2016. Effects of aerated irrigation on CO2emissions from soils of tomatofields.Scientia Agricultura Sinica, 49, 3380–3390. (in Chinese)

    Chen Q, Zhang X, Zhang H, Christie P, Li X, Horlacher D, Liebig H P. 2004. Evaluation of current fertilizer practice and soil fertility in vegetable production in the Beijing region.Nutrient Cycling in Agroecosystems, 69, 51–58.

    Chen X, Dhungel J, Bhattarai S P, Torabi M, Pendergast L,Midmore D J. 2011. Impact of oxygation on soil respiration,yield and water use efficiency of three crop species.Journal of Plant Ecology, 4, 236–248.

    Davidson E A, Schimel J P. 1995. Microbial processes of production and consumption of nitric oxide, nitrous oxide,and methane. In: Matson P A, Harriss R C, eds.,Biogenic Trace Gases:Measuring Emissions from Soil and Water.Blackwell Science, London. pp. 327–357.

    Deng J, Zhou Z, Zheng X, Li C. 2013. Modeling impacts of fertilization alternatives on nitrous oxide and nitric oxide emissions from conventional vegetable fields in southeastern China.Atmospheric Environment, 81,642–650.

    Ding W, Cai Y, Cai Z, Yagi K, Zheng X. 2007. Nitrous oxide emissions from an intensively cultivated maize-wheat rotation soil in the North China Plain.Science of the Total Environment, 373, 501–511.

    Ehret D L, Edwards D, Helmer T, Lin W, Jones G, Dorais M, Papadopoulos A P. 2010. Effects of oxygen-enriched nutrient solution on greenhouse cucumber and pepper production.Scientia Horticulturae, 125, 602–607.

    Flynn H C, Smith J, Smith K A, Wright J, Smith P, Massheder J. 2005. Climate- and crop-responsive emission factors significantly alter estimates of current and future nitrous oxide emissions from fertilizer use.Global Change Biology,11, 1522–1536.

    FAOSTAT (Food and Agriculture Organization Statistical Data).2009. Crops. [2016-06-13]. http://www.fao.org/faostat/en/#data/QC

    Forster P, Ramaswamy V, Artaxo P, Berntsen T, Betts R, Fahey D W, Haywood J, Lean J, Lowe D C, Myhre G, Nganga J,Prinn R, Raga G, Schulz M, Van Dorland R. 2007. Changes in Atmospheric Constituents and in Radiative Forcing.In:Climate Change 2007: The Physical Science Basis.Contribution of Working Group I to the 4th Assessment Report of the Intergovernmental Panel on Climate Change.Cambridge University Press, Cambridge, United Kingdom.pp. 180-185.

    Gao B, Ju X, Su F, Meng Q, Oenema O, Christie P, Chen X,Zhang F. 2014. Nitrous oxide and methane emissions from optimized and alternative cereal cropping systems on the North China Plain: a two-year field study.Science of the Total Environment, 472, 112–124.

    Guadie A, Xia S, Zhang Z, Zeleke J, Guo W, Ngo H H,Hermanowicz S W. 2014. Effect of intermittent aeration cycle on nutrient removal and microbial community in afluidized bed reactor-membrane bioreactor combo system.Bioresource Technology, 156, 195–205.

    He F, Jiang R, Chen Q, Zhang F, Su F. 2009. Nitrous oxide emissions from an intensively managed greenhouse vegetable cropping system in northern China.Environmental Pollution, 157, 1666–1672.

    Hou H, Chen H, Cai H, Yang F, Li D, Wang F. 2016. CO2and N2O emissions from Lou soils of greenhouse tomato fields under aerated Irrigation.Atmospheric Environment, 132,69–76.

    Huang S, Pant H K, Lu J. 2007. Effects of water regimes on nitrous oxide emission from soils.Ecological Engineering,131, 9–15.

    Hwang S, Hanaki K. 2000. Effects of oxygen concentration and moisture content of refuse on nitrification, denitrification and nitrous oxide production.Bioresource Technology,71, 159–165.

    Intergovernmental Panel on Climate Change (IPCC). 2013.Climate Change 2013:The Physical Science Basis.Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change.Cambridge University Press, Cambridge, United Kingdom and New York, USA.

    Iqbal J, Hu R, Feng M, Lin S, Malghani S, Ali I M. 2010. Microbial biomass, and dissolved organic carbon and nitrogen strongly affect soil respiration in different land uses: A case study at Three Gorges Reservoir Area, South China.Agriculture Ecosystems & Environment, 137, 294–307.

    Ju X T, Kou C L, Zhang F S, Christie P. 2006. Nitrogen balance and groundwater nitrate contamination: Comparison among three intensive cropping systems on the North China Plain.Environmental Pollution, 143, 117–125.

    Kallenbach C M, Rolston D E, Horwath W R. 2010. Cover cropping affects soil N2O and CO2emissions differently depending on type of irrigation.Agriculture Ecosystems &Environment, 137, 251–260.

    Li B, Fan C H, Zhang H, Chen Z Z, Sun L Y, Xiong Z Q. 2015a.Combined effects of nitrogen fertilization and biochar on the net global warming potential, greenhouse gas intensity and net ecosystem economic budget in intensive vegetable agriculture in southeastern China.Atmospheric Environment, 100, 10–19.

    Li H J, Yan J X, Yue X F, Wang M B. 2008. Significance of soil temperature and moisture for soil respiration in a Chinese mountain area.Agricultural and Forest Meteorology, 148,490–503.

    Li Y, Niu W Q, Zhang M Z, Xue L, Wang J W. 2015b. Effects of aeration on rhizosphere soil enzyme activities and soil microbes for muskmelon in plastic greenhouse.Transactions of the Chinese Society for Agricultural Machinery, 46, 121–129. (in Chinese)

    Liikanen A, Martikainen P J. 2003. Effect of ammonium and oxygen on methane and nitrous oxide fluxes across sediment-water interface in a eutrophic lake.Chemosphere,52, 1287–1293.

    Linn D M, Doran J W. 1984. Effect of water-filled pore space on carbon dioxide and nitrous oxide production in tilled and nontilled soils.Soil Science Society of America, 48,1267–1272.

    Liu S, Qin Y, Zou J, Liu Q. 2010. Effects of water regime during rice-growing season on annual direct N2O emission in a paddy rice-winter wheat rotation system in southeast China.Science of the Total Environment, 408, 906–913.

    Mersi W V, Schinner F. 1991. An improved and accurate method for determining the dehydrogenase activity of soils with iodonitrotetrazolium chloride.Biology and Fertility of Soils,11, 216–220.

    Niu W, Guo Q, Zhou X, Helmers M J. 2012. Effect of aeration and soil water redistribution on the air permeability under subsurface drip irrigation.Soil Science Society of America,76, 815–820.

    Niu W Q, Fan W T, Persaud N, Zhou X B. 2013. Effect of postirrigation aeration on growth and quality of greenhouse cucumber.Pedosphere, 23, 790–798.

    Patanè C, Cosentino S L. 2010. Effects of soil water deficit on yield and quality of processing tomato under a Mediterranean climate.Agricultural Water Management,97, 131–138.

    Paul J W, Beauchamp E G, Zhang X. 1993. Nitrous and nitric oxide emissions during nitrification and denitrification from manure-amended soil in the laboratory.Canadian Journal of Soil Science, 73, 539–553.

    Reth S, G?ckede M, Falge E. 2005. CO2efflux from agricultural soils in Eastern Germany - comparison of a closed chamber system with eddy covariance measurements.Theoretical and Applied Climatology, 80, 105–120.

    Riya S, Ju M, Sheng Z, Shi W M, Hosomi M. 2012. Short-term responses of nitrous oxide emissions and concentration profiles to fertilization and irrigation in greenhouse vegetable cultivation.Pedosphere, 22, 764–775.

    Robertson G P, Grace P R. 2004. Greenhouse gas fluxes in tropical and temperate agriculture: The need for a fullcost accounting of global warming potentials. In:Tropical Agriculture in Transiction-opportunities for Mitigating Greenhouse Gas Emissions? Springer, Netherlands. pp.51–63.

    Scheer C, Grace P R, Rowlings D W, Payero J. 2013. Soil N2O and CO2emissions from cotton in Australia under varying irrigation management.Nutrient Cycling in Agroecosystems,95, 43–56.

    Scheer C, Wassmann R, Kienzler K, Ibragimov N, Eschanov R. 2008. Nitrous oxide emissions from fertilized, irrigated cotton (Gossypium hirsutumL.) in the Aral Sea Basin,Uzbekistan: Influence of nitrogen applications and irrigation practices.Soil Biology and Biochemistry, 40, 290–301.

    Schmidt U, Th?ni H, Kaupenjohann M. 2000. Using a boundary line approach to analyze N2O flux data from agricultural soils.Nutrient Cycling in Agroecosystems, 57, 119–129.

    Skiba U, Smith K A. 1993. Nitrification and denitrification as sources of nitric oxide and nitrous oxide in a sandy loam soil.Soil Biology and Biochemistry, 25, 1527–1536.

    Sun S K, Wu P T, Wang Y B, Zhao X N, Liu J, Zhang X H.2013. The impacts of interannual climate variability and agricultural inputs on water footprint of crop production in an irrigation district of China.Science of the Total Environment,444, 498–507.

    Sun S K, Zhang C F, Li X L, Zhou T W, Wang Y B, Wu P T,Cai H J. 2017. Sensitivity of crop water productivity to the variation of agricultural and climatic factors: A study of Hetao irrigation district, China.Journal of Cleaner Production,142, 2562–2569.

    Tong H, Yin K, Giannis A, Ge L, Wang J Y. 2015. Influence of temperature on carbon and nitrogen dynamics during in situ aeration of aged waste in simulated landfill bioreactors.Bioresource Technology, 192, 149–156.

    Wan Z M, Song C C, Guo Y D, Wang L, Huang J Y. 2008.Effects of water gradients on soil enzyme activity and active organic carbon composition underCarex lasiocarpamarsh.Acta Ecologica Sinica, 28, 5980–5986.

    Wang W, Feng J, Oikawa T. 2009. Contribution of root and microbial respiration to soil CO2efflux and their environmental controls in a humid temperate grassland of Japan.Pedosphere, 19, 31–39.

    Weier K L, Doran J W, Walters D T. 1993. Denitrification and the dinitrogen/nitrous oxide ratio as affected by soil water,available carbon, and nitrate.Soil Science Society of America, 57, 66–72.

    Weslien P, Rütting T, Kasimir-Klemedtsson ?, Klemedtsson L.2012. Carrot cropping on organic soil is a hotspot for nitrous oxide emissions.Nutrient Cycling in Agroecosystems, 94,249–253.

    Xia Z W, Xu H, Chen G X, Dong D, Bai E, Luo L G. 2013. Soil N2O production and the δ15N-N2O value: Their relationship with nitrifying/denitrifying bacteria and archaea during a growing season of soybean in northeast China.European Journal of Soil Biology, 58, 73–80.

    Xu Y C, Shen Q R, Li M.L, Dittert K, Sattelmacher B. 2004.Effect of soil water status and mulching on N2O and CH4emission from lowland rice field in China.Biology and Fertility of Soils, 39, 215–217.

    Zhang H M, Xiong Y W, Huang G H, Xu X, Huang Q G. 2017.Effects of water stress on processing tomatoes yield, quality and water use efficiency with plastic mulched drip irrigation in sandy soil of the Hetao Irrigation District.Agricultural Water Management, 179, 205–214.

    Zhang M, Fan C H, Li Q L, Li B, Zhu Y Y, Xiong Z Q. 2015.A 2-yr field assessment of the effects of chemical and biological nitrification inhibitors on nitrous oxide emissions and nitrogen use efficiency in an intensively managed vegetable cropping system.Agriculture,Ecosystems and Environment, 201, 43–50.

    Zhang Y L, Wang Y S. 2006. Soil enzyme activities with greenhouse subsurface irrigation.Pedosphere, 16,512–518.

    Zheng B. 2013.Soil Analysis Technology Guide. China Agriculture Press, Beijing. (in Chinese)

    Zhou J B, Zhai B L, Chen Z J, Xu A M, Feng W H. 2006.Fertilizers application and nutrient accumulations in tomatogrown soils under greenhouse condition in the suburban of Xi’an City.Chinese Journal of Soil Science, 37, 2287–2290.(in Chinese)

    Zhu L F, Yu S M, Jin Q Y. 2012. Effects of aerated irrigation on leaf senescence at late growth stage and grain yield of rice.Rice Science, 19, 44–48.

    Zhu Y, Cai H J, Hou H J, Song L B. 2016. Effects of aerated irrigation on root-zone environment and yield of tomato.Journal of Northwest A&F University(Natural Science Edition), 44, 157–162. (in Chinese)

    Zhu Z L, Chen D L. 2002. Nitrogen fertilizer use in China- Contributions to food production, impacts on the environment and best management strategies.Nutrient Cycling in Agroecosystems, 63, 117–127.

    Zornoza R, Rosales R M, Acosta J A, de la Rosa J M, Arcenegui V, Faz á, Pérez-Pastor A. 2016. Efficient irrigation management can contribute to reduce soil CO2emissions in agriculture.Geoderma, 263, 70–77.

    久久人妻av系列| 亚洲黑人精品在线| av在线天堂中文字幕| a级毛片在线看网站| 母亲3免费完整高清在线观看| 91久久精品国产一区二区成人 | 一个人看的www免费观看视频| 久久久久精品国产欧美久久久| 两人在一起打扑克的视频| 精品无人区乱码1区二区| 久久久精品大字幕| 变态另类成人亚洲欧美熟女| 欧美av亚洲av综合av国产av| 亚洲av第一区精品v没综合| 精品一区二区三区四区五区乱码| 欧美日韩国产亚洲二区| 好男人电影高清在线观看| 88av欧美| 别揉我奶头~嗯~啊~动态视频| 日本在线视频免费播放| 国产成+人综合+亚洲专区| 色哟哟哟哟哟哟| 久久天堂一区二区三区四区| 天天躁日日操中文字幕| 成年女人毛片免费观看观看9| 亚洲午夜理论影院| 日韩免费av在线播放| 国产成年人精品一区二区| 欧美成人一区二区免费高清观看 | 亚洲18禁久久av| 国产成人影院久久av| 精品日产1卡2卡| 成人一区二区视频在线观看| 嫩草影视91久久| 波多野结衣高清作品| 久久久久九九精品影院| 99热这里只有精品一区 | 日本精品一区二区三区蜜桃| 俺也久久电影网| 国产久久久一区二区三区| 不卡av一区二区三区| 亚洲真实伦在线观看| 久久久久久久午夜电影| 国产真实乱freesex| 亚洲人成网站高清观看| 亚洲人成网站在线播放欧美日韩| 亚洲,欧美精品.| 丝袜人妻中文字幕| 精品久久久久久久久久久久久| 成人国产综合亚洲| 一级a爱片免费观看的视频| 叶爱在线成人免费视频播放| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 99国产极品粉嫩在线观看| 国产精品爽爽va在线观看网站| 操出白浆在线播放| 免费一级毛片在线播放高清视频| av中文乱码字幕在线| 久久久水蜜桃国产精品网| 精品欧美国产一区二区三| 中文字幕最新亚洲高清| 变态另类成人亚洲欧美熟女| 99精品欧美一区二区三区四区| a在线观看视频网站| 国产成人啪精品午夜网站| 亚洲午夜精品一区,二区,三区| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一欧美日韩一区二区三区| 黄色日韩在线| 亚洲精品美女久久久久99蜜臀| 久久午夜亚洲精品久久| 99热这里只有是精品50| 亚洲一区二区三区色噜噜| 男女床上黄色一级片免费看| 欧美色欧美亚洲另类二区| 亚洲成av人片在线播放无| 精品人妻1区二区| 听说在线观看完整版免费高清| aaaaa片日本免费| 亚洲成av人片在线播放无| 亚洲美女黄片视频| 欧美日韩乱码在线| 欧美性猛交黑人性爽| bbb黄色大片| 美女被艹到高潮喷水动态| 我的老师免费观看完整版| 免费高清视频大片| 99精品久久久久人妻精品| 免费在线观看日本一区| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 夜夜夜夜夜久久久久| 久久精品国产综合久久久| 久久久久亚洲av毛片大全| 国内毛片毛片毛片毛片毛片| 成人三级做爰电影| www日本在线高清视频| 男人舔奶头视频| 亚洲真实伦在线观看| 在线观看午夜福利视频| 男人和女人高潮做爰伦理| 三级毛片av免费| aaaaa片日本免费| 欧洲精品卡2卡3卡4卡5卡区| 黄色丝袜av网址大全| 在线永久观看黄色视频| 国产男靠女视频免费网站| tocl精华| 99国产综合亚洲精品| 九九热线精品视视频播放| 制服人妻中文乱码| 欧美中文日本在线观看视频| 后天国语完整版免费观看| 在线观看日韩欧美| 两性夫妻黄色片| 热99re8久久精品国产| 国产成人av教育| 精品一区二区三区四区五区乱码| 国产单亲对白刺激| 国产高清videossex| 最新美女视频免费是黄的| 国产激情偷乱视频一区二区| 欧美日本视频| 香蕉久久夜色| 老司机在亚洲福利影院| 三级国产精品欧美在线观看 | 免费看十八禁软件| 午夜激情欧美在线| 亚洲五月天丁香| 国产高清激情床上av| 深夜精品福利| 午夜福利在线观看吧| 在线免费观看不下载黄p国产 | 免费在线观看影片大全网站| 国产伦人伦偷精品视频| 成人18禁在线播放| 全区人妻精品视频| 女生性感内裤真人,穿戴方法视频| 欧美黄色淫秽网站| 白带黄色成豆腐渣| 99热只有精品国产| 午夜激情欧美在线| 亚洲av五月六月丁香网| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 国产午夜精品论理片| 亚洲片人在线观看| 久久久久九九精品影院| 人人妻人人看人人澡| 美女免费视频网站| 成人午夜高清在线视频| 久久久精品大字幕| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区精品| 国产免费男女视频| 色综合亚洲欧美另类图片| 午夜精品久久久久久毛片777| 中文资源天堂在线| 国产精品av视频在线免费观看| 毛片女人毛片| 9191精品国产免费久久| 男人舔女人下体高潮全视频| 日韩有码中文字幕| 国产午夜精品论理片| 在线永久观看黄色视频| 国产美女午夜福利| 一个人看视频在线观看www免费 | av女优亚洲男人天堂 | 桃红色精品国产亚洲av| 9191精品国产免费久久| 国产蜜桃级精品一区二区三区| 男人和女人高潮做爰伦理| 欧美日韩国产亚洲二区| 免费电影在线观看免费观看| 免费看日本二区| 精品国产亚洲在线| 精品国产乱码久久久久久男人| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 免费看美女性在线毛片视频| 99久久成人亚洲精品观看| ponron亚洲| 搡老妇女老女人老熟妇| 夜夜躁狠狠躁天天躁| 久久人人精品亚洲av| 亚洲七黄色美女视频| 久久久久久久久免费视频了| 五月伊人婷婷丁香| 精品欧美国产一区二区三| 欧美乱色亚洲激情| 88av欧美| 日本在线视频免费播放| 男人舔奶头视频| 久久精品91无色码中文字幕| 波多野结衣巨乳人妻| 国产一区二区在线观看日韩 | 夜夜夜夜夜久久久久| 欧美日韩综合久久久久久 | 九色国产91popny在线| 国产精品影院久久| www.www免费av| 亚洲成av人片在线播放无| avwww免费| 成人三级做爰电影| 久久国产精品影院| 中文字幕人妻丝袜一区二区| 欧美3d第一页| 每晚都被弄得嗷嗷叫到高潮| 欧美激情在线99| 国产精品野战在线观看| 亚洲国产看品久久| av女优亚洲男人天堂 | 久久伊人香网站| 色哟哟哟哟哟哟| 观看美女的网站| 级片在线观看| 美女大奶头视频| 俺也久久电影网| 欧美日本视频| 亚洲,欧美精品.| 免费看日本二区| 99久久久亚洲精品蜜臀av| 香蕉丝袜av| 天堂影院成人在线观看| 日本五十路高清| 免费搜索国产男女视频| 亚洲av成人一区二区三| 热99在线观看视频| 国产日本99.免费观看| 黑人巨大精品欧美一区二区mp4| 精品久久久久久久毛片微露脸| 国产精品永久免费网站| 亚洲精品在线观看二区| 性色avwww在线观看| 欧美乱码精品一区二区三区| 国产av不卡久久| 精品日产1卡2卡| 精品不卡国产一区二区三区| 1024香蕉在线观看| 国产精品一区二区三区四区久久| av在线天堂中文字幕| 亚洲成a人片在线一区二区| 淫妇啪啪啪对白视频| 国产黄片美女视频| av国产免费在线观看| 极品教师在线免费播放| 国产淫片久久久久久久久 | av中文乱码字幕在线| 国产av不卡久久| 嫩草影院入口| 狂野欧美激情性xxxx| svipshipincom国产片| 国产乱人伦免费视频| 久久久国产欧美日韩av| 日本精品一区二区三区蜜桃| 日韩高清综合在线| 女人被狂操c到高潮| 丁香六月欧美| 97人妻精品一区二区三区麻豆| 丝袜人妻中文字幕| 国产精品亚洲美女久久久| 国产av不卡久久| 欧美中文综合在线视频| av天堂在线播放| 这个男人来自地球电影免费观看| 亚洲国产精品合色在线| 一二三四社区在线视频社区8| 午夜福利成人在线免费观看| 免费看a级黄色片| 亚洲成av人片免费观看| 超碰成人久久| 国产三级黄色录像| 在线播放国产精品三级| 日韩精品青青久久久久久| 亚洲成av人片免费观看| 欧美日韩瑟瑟在线播放| 老汉色av国产亚洲站长工具| 精品久久久久久成人av| 女同久久另类99精品国产91| 久久久成人免费电影| 精品一区二区三区视频在线观看免费| 好男人电影高清在线观看| 黄色视频,在线免费观看| 成人鲁丝片一二三区免费| 国产熟女xx| av天堂中文字幕网| 午夜福利欧美成人| 久久午夜亚洲精品久久| 国产 一区 欧美 日韩| 亚洲自偷自拍图片 自拍| 97人妻精品一区二区三区麻豆| 日日干狠狠操夜夜爽| 亚洲avbb在线观看| 又紧又爽又黄一区二区| АⅤ资源中文在线天堂| 午夜激情福利司机影院| 成人国产综合亚洲| 国产欧美日韩一区二区三| 夜夜爽天天搞| 亚洲精品乱码久久久v下载方式 | 伦理电影免费视频| 成年版毛片免费区| 国产野战对白在线观看| 亚洲无线在线观看| 免费看十八禁软件| 国产精品久久电影中文字幕| 黄色日韩在线| 国产伦精品一区二区三区视频9 | 97人妻精品一区二区三区麻豆| 最近视频中文字幕2019在线8| 亚洲色图 男人天堂 中文字幕| 亚洲av成人精品一区久久| 两个人的视频大全免费| www.999成人在线观看| 狂野欧美白嫩少妇大欣赏| 国产人伦9x9x在线观看| 日本黄色视频三级网站网址| 舔av片在线| 久久久久精品国产欧美久久久| 村上凉子中文字幕在线| 在线免费观看不下载黄p国产 | 欧美丝袜亚洲另类 | 欧美中文日本在线观看视频| 午夜影院日韩av| 日本一二三区视频观看| 免费无遮挡裸体视频| 精品国产乱码久久久久久男人| 精品久久久久久久人妻蜜臀av| 宅男免费午夜| xxx96com| 国产伦精品一区二区三区视频9 | 国产成人精品久久二区二区免费| 19禁男女啪啪无遮挡网站| 法律面前人人平等表现在哪些方面| 亚洲真实伦在线观看| 99视频精品全部免费 在线 | 全区人妻精品视频| 中国美女看黄片| 色老头精品视频在线观看| 舔av片在线| 天堂av国产一区二区熟女人妻| 免费电影在线观看免费观看| 黑人欧美特级aaaaaa片| 国产精品久久久久久精品电影| 日本三级黄在线观看| 精华霜和精华液先用哪个| 美女高潮的动态| 亚洲美女黄片视频| 两人在一起打扑克的视频| 嫩草影视91久久| 狂野欧美白嫩少妇大欣赏| 黄色 视频免费看| 毛片女人毛片| 久久久久久人人人人人| 99在线视频只有这里精品首页| 精品国产亚洲在线| 男插女下体视频免费在线播放| 99视频精品全部免费 在线 | 欧美日本亚洲视频在线播放| 首页视频小说图片口味搜索| 男人舔奶头视频| 黄片大片在线免费观看| 变态另类成人亚洲欧美熟女| 久久久久国内视频| 国语自产精品视频在线第100页| 男女之事视频高清在线观看| 少妇熟女aⅴ在线视频| 88av欧美| 天堂网av新在线| 亚洲精品美女久久久久99蜜臀| 久久热在线av| 久久这里只有精品中国| 狂野欧美激情性xxxx| 神马国产精品三级电影在线观看| 日本撒尿小便嘘嘘汇集6| 999久久久精品免费观看国产| 国产精品久久视频播放| 成人高潮视频无遮挡免费网站| 两个人视频免费观看高清| 日韩国内少妇激情av| 久久中文字幕人妻熟女| 国产1区2区3区精品| 久久精品影院6| 美女被艹到高潮喷水动态| 欧美不卡视频在线免费观看| 色噜噜av男人的天堂激情| 色av中文字幕| 最新中文字幕久久久久 | 人人妻,人人澡人人爽秒播| 国产成人精品久久二区二区91| 成人高潮视频无遮挡免费网站| 久久香蕉国产精品| 成人精品一区二区免费| 色综合站精品国产| 国产激情偷乱视频一区二区| 偷拍熟女少妇极品色| 少妇的逼水好多| 老汉色av国产亚洲站长工具| 麻豆一二三区av精品| 九色成人免费人妻av| 一级毛片女人18水好多| 久久婷婷人人爽人人干人人爱| 天堂影院成人在线观看| 久久精品国产清高在天天线| 国产精品永久免费网站| 免费电影在线观看免费观看| 在线播放国产精品三级| 国产主播在线观看一区二区| 色尼玛亚洲综合影院| 美女黄网站色视频| 国产高清三级在线| 亚洲,欧美精品.| 中文在线观看免费www的网站| 99精品欧美一区二区三区四区| 夜夜躁狠狠躁天天躁| 亚洲熟妇熟女久久| 国产亚洲精品av在线| 91av网站免费观看| 精品福利观看| 麻豆av在线久日| 亚洲成人免费电影在线观看| 亚洲成人久久爱视频| 国产乱人视频| 亚洲欧美激情综合另类| 巨乳人妻的诱惑在线观看| 国产成人欧美在线观看| 亚洲国产看品久久| 99久久成人亚洲精品观看| 欧美黑人欧美精品刺激| 脱女人内裤的视频| 国产一区二区在线av高清观看| 真人做人爱边吃奶动态| 亚洲九九香蕉| 熟妇人妻久久中文字幕3abv| 一个人观看的视频www高清免费观看 | 手机成人av网站| 18美女黄网站色大片免费观看| 亚洲电影在线观看av| 国产淫片久久久久久久久 | 男女床上黄色一级片免费看| 日本撒尿小便嘘嘘汇集6| 悠悠久久av| 国产精品久久视频播放| 亚洲专区字幕在线| 婷婷六月久久综合丁香| 好男人在线观看高清免费视频| 久久精品国产综合久久久| 久久精品夜夜夜夜夜久久蜜豆| 给我免费播放毛片高清在线观看| 国产成+人综合+亚洲专区| 香蕉av资源在线| 天天躁日日操中文字幕| 在线国产一区二区在线| www日本黄色视频网| 一本一本综合久久| 天天一区二区日本电影三级| 久久香蕉国产精品| 国产精品 欧美亚洲| 国产乱人视频| 十八禁人妻一区二区| 欧美日本视频| 五月玫瑰六月丁香| 熟女少妇亚洲综合色aaa.| 麻豆久久精品国产亚洲av| 成在线人永久免费视频| 中文字幕久久专区| 中文字幕人妻丝袜一区二区| 国产视频一区二区在线看| 一个人免费在线观看电影 | 男人和女人高潮做爰伦理| 国产高清视频在线播放一区| 欧美日韩福利视频一区二区| 久久国产乱子伦精品免费另类| 精品久久久久久成人av| 他把我摸到了高潮在线观看| 午夜精品一区二区三区免费看| 最近视频中文字幕2019在线8| 熟妇人妻久久中文字幕3abv| 国模一区二区三区四区视频 | 欧美成人免费av一区二区三区| 国产主播在线观看一区二区| 可以在线观看毛片的网站| 露出奶头的视频| 国产亚洲精品一区二区www| 国产精品久久久久久人妻精品电影| 动漫黄色视频在线观看| 亚洲成人久久爱视频| 人人妻人人澡欧美一区二区| 国产精品 欧美亚洲| 日韩三级视频一区二区三区| 日韩欧美在线乱码| 亚洲性夜色夜夜综合| 手机成人av网站| 麻豆国产av国片精品| 好男人电影高清在线观看| 99久国产av精品| 亚洲精品色激情综合| 少妇的逼水好多| av视频在线观看入口| 黄色日韩在线| 中国美女看黄片| 久久久久久九九精品二区国产| 日本撒尿小便嘘嘘汇集6| 男女做爰动态图高潮gif福利片| 中文字幕最新亚洲高清| 国产精品九九99| 色老头精品视频在线观看| 国产成人影院久久av| 老熟妇仑乱视频hdxx| av片东京热男人的天堂| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添小说| 黄频高清免费视频| 亚洲电影在线观看av| 无人区码免费观看不卡| avwww免费| 亚洲 欧美一区二区三区| 99国产综合亚洲精品| 成年版毛片免费区| 黄色片一级片一级黄色片| 日韩 欧美 亚洲 中文字幕| 国产又黄又爽又无遮挡在线| 久久中文字幕人妻熟女| 免费av不卡在线播放| 亚洲成人久久爱视频| 男女下面进入的视频免费午夜| av黄色大香蕉| 免费电影在线观看免费观看| 9191精品国产免费久久| 免费看十八禁软件| 日韩人妻高清精品专区| 色综合欧美亚洲国产小说| 深夜精品福利| 国产亚洲精品av在线| 国产精品久久视频播放| 国产av在哪里看| 亚洲国产精品久久男人天堂| 久久久久久国产a免费观看| h日本视频在线播放| 久9热在线精品视频| 国产精品久久久av美女十八| 精品一区二区三区av网在线观看| 成在线人永久免费视频| 熟女电影av网| 色视频www国产| 男女做爰动态图高潮gif福利片| 国产精品野战在线观看| 日本a在线网址| 国内精品美女久久久久久| 女人被狂操c到高潮| 亚洲,欧美精品.| 91在线观看av| 国产精品久久视频播放| 国产男靠女视频免费网站| 成人18禁在线播放| 一个人免费在线观看的高清视频| 又紧又爽又黄一区二区| 国产成人一区二区三区免费视频网站| 久久这里只有精品中国| 黑人操中国人逼视频| 琪琪午夜伦伦电影理论片6080| 波多野结衣巨乳人妻| 亚洲av成人av| 淫妇啪啪啪对白视频| 亚洲激情在线av| 婷婷丁香在线五月| 一区二区三区激情视频| 夜夜爽天天搞| 国产激情偷乱视频一区二区| 综合色av麻豆| 国产伦在线观看视频一区| 亚洲精品中文字幕一二三四区| 在线十欧美十亚洲十日本专区| 亚洲欧美一区二区三区黑人| 一个人看视频在线观看www免费 | 啪啪无遮挡十八禁网站| 麻豆成人av在线观看| 高清毛片免费观看视频网站| av福利片在线观看| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久久久毛片| 丝袜人妻中文字幕| 全区人妻精品视频| 一边摸一边抽搐一进一小说| 变态另类成人亚洲欧美熟女| 欧美成狂野欧美在线观看| 成人国产综合亚洲| 毛片女人毛片| 色在线成人网| 99热这里只有是精品50| 色播亚洲综合网| 日韩欧美国产一区二区入口| 国产aⅴ精品一区二区三区波| 91av网站免费观看| 国产午夜福利久久久久久| 99久久综合精品五月天人人| 久久久久久久午夜电影| 免费无遮挡裸体视频| 两人在一起打扑克的视频| 男女床上黄色一级片免费看| 日韩精品中文字幕看吧| aaaaa片日本免费| 99精品欧美一区二区三区四区| 欧美不卡视频在线免费观看| 99热6这里只有精品| 丰满人妻一区二区三区视频av | 成人av一区二区三区在线看| 法律面前人人平等表现在哪些方面| 欧美日韩中文字幕国产精品一区二区三区| 在线观看66精品国产| 国产v大片淫在线免费观看| 亚洲av成人av| 国产精品98久久久久久宅男小说|