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

    Effects of mulch and planting methods on Medicago ruthenica seed yield and soil physical-chemical properties

    2022-08-18 09:32:56WANGYuxiaZHANGJingYUXiaojun
    Journal of Arid Land 2022年8期

    WANG Yuxia, ZHANG Jing, YU Xiaojun

    College of Grassland Science, Gansu Agricultural University, Lanzhou 730070, China

    Abstract: Medicago ruthenica (L.) Trautv., a wild grass species, is commonly grown as a forage crop in arid and semi-arid areas of China. Herein, we explored mulch patterns and planting methods for optimizing M.ruthenica seed production in the loess plateau of the Gansu Province, China from 2017 to 2019. The experiments comprised of six treatments including (1) flat ground without mulch (F0, control); (2) flat ground with a transparent white 0.008 mm thick plastic film mulch (FP); (3) flat ground with 4500 kg/hm2 straw mulch (FS); (4) furrow with 10 cm ridges (R0); (5) furrow with plastic film mulch (RP); and (6)furrow with straw mulch (RS). Results showed that the harvested seed yield of M. ruthenica was the highest under RP treatment, followed by FP and FS treatments. Soil moisture content from mid-May to mid-August in 2017 was the highest under RP and FP treatments, followed by RS and FS treatments. In 2018, soil moisture content was the highest under RS and FS treatments. In 2017 and 2018, soil temperature was the highest under FP and RP treatments, followed by F0 and R0 treatments. Total and available nitrogen, phosphorus, and potassium contents were the highest under RS and FS treatments,followed by RP and FP treatments. Comprehensive analysis result showed that surface mulch improved soil microenvironment and increased seed yield of M. ruthenica. Straw mulch also effectively recycled excess crop straw, thereby encouraging the sustainable development of agriculture in this area. In conclusion, FS treatment was considered the best mode for M. ruthenica seed production in this area.

    Keywords: arid area; plastic film; straw mulch; soil moisture content; soil temperature; soil chemistry

    1 Introduction

    Medicago ruthenica(L.) Trautv., a perennial herbaceous plant of Leguminosae, is a widely grown forage crop in arid and semi-arid areas due to its high potential yield and tolerance to drought and cold (Shu et al., 2018). Its ability to grow in barren environments makes it promising for forage cultivation and reseeding of natural grassland, thereby playing a substantial role in grassland animal husbandry and ecological restoration (Huanget al., 2007; Liet al., 2013; Xiao et al.,2018a). However, matureM. ruthenicaplants exhibit natural seed dispersal, which hinders the large-scale collection of mature seeds for agricultural purposes. To date,M. ruthenicaseed production and collection has received limited scientific attention (Liet al., 2006; Li and Shi,2006). The difficulty in collectingM. ruthenicaseeds prevents the possibility of large-scale propagation of this species; hence, it is important to researchM. ruthenicaseed production and collection.

    In recent years, due to the modification of Chinese agricultural planting structure, forage seed is not only important for the production and development of grassland industry, but also a key factor for forage production. Furthermore, continuous attention to ecological environment protection, animal husbandry development, and urban greening in China has increased the demand for excellent grass varieties and high-quality forage seeds (Shao et al., 2014).Additionally, such demands have revealed the problems associated with low yield and poor benefits of forage seed production in China.

    Water is a major limiting factor forM. ruthenicaseed production in arid and semi-arid areas because of low precipitation and high evaporation. Mulch can majorly maintain soil moisture and improve water utilization in arid and semi-arid areas. At the same time, mulch covers can reduce evaporation and compensate for inconsistent precipitation. In addition, mulch can also improve soil environment and considerably increase seed yield (Zhao et al., 2012; Chenet al., 2015; Yang et al., 2018). At present, water saving measures include film mulch, straw mulch, and the use of ridge and furrow (Han, 2018). Surface mulch with ridge and furrow treatment regulates soil temperature, soil moisture content, and air permeability, thereby supporting microbial communities and effectively promoting material transformation and nutrient decomposition in arid areas (Liu et al., 2014). Surface mulch with plastic film can promote soil organic carbon(SOC) and total nitrogen (TN) mineralization, and can effectively increase available nitrogen(AN), phosphorus (AP), and potassium (AK) contents in soil. However, soil water collected within the shallow layer under plastic mulch is absorbed by the crops from the lower layer when growing, which causes the upward movement of water from deep layer, thereby reducing deep layer water content (Zhou et al., 2012; Kaderet al., 2017). In straw mulch, straw decomposes to produce organic compounds under favorable soil conditions. This, combined with rainwater seepage, leads to a gradual breakdown of mulch, increasing SOC and available nutrient contents.Thus, straw mulch can serve as an important source of soil organic matter (SOM) content (Tang,2019). In a study of straw mulch treatment for wheat and maize, SOC and TN contents at 0-10 cm soil depth increased by 16.9% and 7.7%, respectively, compared with that of control treatment(Shao et al., 2016). However, in previous studies, it has been observed that plants absorb excess soil AN, AP, and AK during growth and development, decreasing available soil nutrients (Hu et al., 2012; Yin et al., 2016).

    Longzhong Area is located in the middle of the Gansu Province, accounting for 16.8% area of the Gansu Province, China. The distribution of annual precipitation is irregular and varies considerably. Most precipitation occurs from July to September. Due to the concentrated rainfall in plant growth season, Longzhong Area is a typical rain-fed agriculture (Liu and Wang, 2018).However, the researches of effects of mulch treatments onM. ruthenicaseed yield and soil nutrient contents are less concerned in the Longzhong Area (Luo et al., 2014). Our aims are to analyze the effects of different mulch treatments onM. ruthenicaseed yield and soil physical-chemical properties in the arid Longzhong Area. The results might provide a theoretical basis for guiding forage cultivation and precision fertilization in the similar areas.

    2 Materials and methods

    2.1 Study site

    The study site was located in the Longzhong Area, Gansu Province, China (35°58′N, 104°62′E)with an altitude of 2000 m a.s.l. The main soil type in this site is yellow soil, with average annual sunshine of 2409 h, average annual precipitation of 427 mm with extremely uneven distribution,and average annual evaporation of 1510 mm. The annual average temperature is 6.3℃, and average annual active accumulated temperatures ≥5.0℃ and ≥10.0℃ are 2782.5℃ and 2239.1℃,respectively. Temperature extremes ranged from -27.1℃ to 34.3℃, and the average annual frost-free period is 141 d (Liu and Wang, 2017). Soil pH is 7.0-8.2. SOM, TN, total phosphorous(TP), and total potassium (TK) contents are 47.81, 1.38, 1.10, and 8.36 g/kg, respectively.Furthermore, AN, AP, and AK contents are 36.32, 7.37, and 159.89 mg/kg, respectively. The average monthly temperature and precipitation from 2017 to 2019 are shown in Figure 1.

    Fig. 1 Monthly temperature and precipitation in the study area from 2017 to 2019

    2.2 Experimental design

    The experiment comprised of six treatments, each with 30 cm row spacing: (1) flat ground without mulch (F0, control); (2) flat ground with a transparent white 0.008 mm thick plastic film mulch (FP); (3) flat ground with 4500 kg/hm2straw mulch (FS); (4) furrow with 10 cm ridge (R0);(5) furrow with plastic film mulch (RP); and (6) furrow with straw mulch (RS; Fig. 2). In the ridge and furrow design, plants were grown in the furrow, not on the ridge. The area of each cell was 2 m×5 m, and the treatment interval was 0.5 m. The treatment design was a random block,and was repeated thrice. We sown the seeds on 2 May, 2017, using the hole seeding method according to seedling emergence appropriate supplementation, or reduction in seedlings, at a depth of 2-3 cm, with 5-10 seeds per hole to ensure a survival rate of 3-5 per hole. The plant spacing was 10 cm, and the row spacing was 30 cm. Before sowing,M. ruthenicaseeds were soaked in 98% concentrated sulfuric acid for 15 min, followed by washing with running water to break the hard outer cover.

    Fig. 2 Schematic diagram of M. ruthenica planting under different mulch treatments. Gray straight line refers to the transparent plastic film, and dotted line refers to the surface straw cover. R0, furrow with 10 cm ridge; RP,furrow with plastic film mulch; RS, furrow with straw mulch; F0, flat ground without mulch; FP, flat ground with a transparent white 0.008 mm thick plastic film mulch; FS, flat ground with 4500 kg/hm2 straw mulch. The abbreviations are the same as the following figures.

    2.3 Measurements

    2.3.1 Performance seed yield Performance seed yield (kg/hm2) is calculated by Equation 1:

    whereYis the performance seed yield (kg/hm2);Pis the number of pods per unit area (m2);Sis the number of seeds per pod; andWis the thousand-seed weight (g) (Xiao et al., 2018b).

    2.3.2 Actual seed yield

    Due to the dispersal of the matured seeds, we selected a 1-m2sample area away from the side line of the treated sample to measure the actual seed yield in the seed maturity stage. This process was repeated thrice per treatment to calculate the average value.

    2.3.3 Soil moisture content, soil bulk density, and soil porosity

    Soil samples were collected once a month on a sunny day during the growing season ofM.ruthenica(soil moisture content was collected from mid-May to mid-September in 2017 and 2018,while soil bulk density and soil porosity soil samples were collected from mid-May to mid-September in 2018 and 2019). Three soil profiles were assessed using the randomly dug soil from each treatment, and five soil samples were collected from the four 10 cm segments (0-10,10-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 70-80 cm) using a 5-cm diameter soil ring knife.Each cylindrical soil sample was loaded into a sample box and transferred to the laboratory. The samples were placed in an electric thermostatic air blowing drying oven (GZX-GF101-3-BS- Ⅱ,Shanghai, China) and dried at 105℃ until a constant weight was obtained for measurements and subsequent calculations.

    Soil moisture content is calculated by Equation 2:

    whereMis the soil moisture content (%);M2is the mass of wet soil in the aluminum box (g);M1is the mass of dry soil in the aluminum box (g); andM0is the weight of an empty aluminum box dried to a constant weight (g) (Wilke, 2005).

    Soil bulk density is calculated by Equation 3:

    whereBis the soil bulk density (g/cm3);W1is the total weight of ring knife and fresh soil (g);W0is the total weight of ring cutter (g);Mis the soil moisture content (%); andVis the volume of soil core (100 cm3).

    Soil porosity is calculated by Equation 4:

    wherePis the soil porosity (kPa);Gis the specific gravity of soil fixed at 2.65 g/cm3.

    2.3.4 Soil temperature

    Soil temperature was determined by selecting one sunny day in the middle of each month during the growing season. A curved tube geothermometer was used to measure temperature every 5 cm at a depth of 5-25 cm. Measurements were recorded every 2 h from 08:00 to 18:00 (LST), and each measurement was repeated thrice to obtain the average value (Zhou et al., 2017).

    2.3.5 Soil compactness

    Soil compactness was measured using a soil compactness meter (SC-900, spectrum, USA). Soil compactness was measured up to a thickness of 30 cm, and measurements were made 20 times for each plot.

    2.3.6 Soil nutrients Soil nutrients change gradually; hence, they were collected only once during mid-October between 2018 and 2019. Soils cores were collected as described in Section 2.3.3 using a 3.5-cm diameter soil auger. After removing the visible debris, we placed samples into zip-top bags and transported to the laboratory. After natural drying, we passed each sample through 0.25- and 1.00-mm soil sieves for soil nutrient analysis (Hou et al.,2018).

    According to the methods described by Bao (2000), we determined SOM content using the potassium dichromate titration. Soil TN content was determined using the Kjeldahl method. Soil TP content was determined using the molybdenum blue method with an ultraviolet-visible spectrophotometer (SP-756P, Shanghai Guangpu Instuments Co., Ltd., China). Soil TK content was determined using the NaOH melt-flame photometry. Soil AN content was determined using the alkali-hydrolyzed diffusion method. Soil AP content was determined using the 0.5 mol/L NaHCO3-molybdenum-antimony resistance colorimetric method. Soil AK content was extracted using the NH4OAC with flame photometry.

    2.4 Statistical analysis

    Microsoft excel 2010 was used for data collation and visualization, and all values were expressed as mean±standard errors (SE). Compared means in SPSS v.22.0 (SPSS Inc., Chicago, IL, USA)was used for conducting a one-way analysis of variance on the experimental data. Generalized lineal model (GLM) was used to analyze the effects of planting years and different mulch treatments onM. ruthenicaseed production and soil properties. Significance was determined atP<0.05 level. The Duncan method was used for multiple comparison analysis.

    3 Results

    3.1 Seed yield of M. ruthenica

    The performance and actual seed yields ofM. ruthenicafrom 2017 to 2019 were the highest under RP treatment, and the lowest under F0 treatment. Seed yield ofM. ruthenicaunder FS treatment was always higher than that of RS treatment. Furthermore, plastic mulch had the highest seed yield, followed by straw mulch (Fig. 3). Seed yield under R0 and F0 treatments were significantly lower than those of other treatments (P<0.05) from 2017 to 2018. There were no significant differences between RP and FP, between RS and FS, and between R0 and F0 treatments for the performance and actual seed yields (P>0.05). The average actual seed yield ofM. ruthenicaduring 2017-2019 was 203-279 kg/hm2(Fig. 3b).

    Under the same treatment, the performance seed yield ofM. ruthenicain 2018 was significantly higher than those in 2017 and 2019 (P<0.05), with no significant difference observed between 2017 and 2019 for the same treatment (Fig. 3a). The actual seed yield of all treatments in 2018 was significantly higher than that in 2017, and there was no significant difference observed between 2018 and 2019 under R0 and F0 treatments (P>0.05; Fig. 3b).

    Fig. 3 Effects of mulch and planting methods of M. ruthenica on performance (a) and actual seed yields (b)from 2017 to 2019. Different uppercase letters within the same treatment indicate significant differences among different years at P<0.05 level; Different lowercase letters within the same year indicate significant differences among different treatments at P<0.05 level.

    3.2 Effect of mulch and planting methods on soil physical properties

    3.2.1 Soil moisture content

    In 2017, from mid-May to mid-August, the order of soil moisture content for each depth under different mulch treatments was R (furrow)>F (flat ground), P (plastic film mulch)>S (straw mulch)>0 (no mulch), and soil moisture content values for 0-10 and 10-20 cm depths in each treatment were slightly higher than those of other soil depths (Fig. 4). However, in mid-September, soil moisture contents of corresponding soil depth under RS and FS treatments were higher than those of RP and FP treatments. Soil moisture contents under R0 and F0 treatments were significantly lower than those of other treatments (P<0.05).

    Fig. 4 Effects of mulch and planting methods of M. ruthenica on soil moisture content in 2017 and 2018.Different lowercase letters within same soil depth indicate significant differences among different treatments at P<0.05 level. (a), May 2017; (b), June 2017; (c), July 2017; (d), August 2017; (e), September 2017; (f), June 2018;(g), July 2018; (h), August 2018.

    In 2018, from mid-June to mid-August, soil moisture content for each soil depth under different mulch treatments was the highest under RS treatment and the lowest under F0 treatment (Fig. 4).Under RP treatment, soil moisture content was always higher than that of FP treatment. Moreover,plastic mulch treatment was more effective than straw mulch. Soil moisture content of each treatment showed a decreasing trend with increasing in soil depths.

    From mid-June to mid-July, soil moisture content under RS and FS treatments were significantly higher than that of F0 treatment at 0-80 cm soil depth, whereas soil moisture content at 40-80 cm depth under RS and FS treatments were significantly higher than those of other treatments (P<0.05).

    3.2.2 Soil temperature

    From mid-May to mid-September in 2017 and 2018, the average daily soil temperature increased rapidly with plastic mulch treatments, followed by control and straw mulch treatments (Fig. 5).Soil temperature for each treatment exhibited a significant daily change at 5-10 cm soil depth, but it did not vary significantly at 10-25 cm depth. Average soil temperature of each treatment presented an inverse trend with soil depths. From 08:00 to 18:00 on a given day, soil temperature for each depth under each treatment first increased and then decreased.

    In mid-May, mid-June, and mid-September of 2017, soil temperature for each treatment reached the highest at 14:00, and then gradually decreased. Whereas, in mid-July, mid-August of 2017, mid-June, mid-July, and mid-August of 2018, soil temperature for each treatment peaked at 16:00.

    From mid-May to mid-September of 2017, soil temperature under RP treatment at 08:00 was the highest, followed by FP treatment (Fig. 5). At 5 cm soil depth, from mid-May to mid-August,soil temperatures under RP and FP treatments were significantly higher than those of other treatments (P<0.05); and soil temperature under FP treatment was significantly higher than that of FP treatment at 10:00, 12:00, 14:00, 16:00, and 18:00. In mid-May, mid-June, and mid-August,soil temperature under RS treatment was significantly lower than those of RP and FP treatments(P<0.05). In mid-September, at 5 cm soil depth, no significant difference was observed in soil temperature between RP and FP treatments (P>0.05).

    3.2.3 Soil bulk density, soil porosity, and soil compactness

    In 2018 and 2019, soil bulk density for each soil depth under different treatments was the lowest for straw mulch and the highest for control (Fig. 6a). Soil bulk density at 0-20 cm depth was the lowest under RP, FP, RS, and FS treatments. Analysis of variance showed that growth year did not significantly affect soil bulk density. Soil bulk density reached an extremely significant level at 0-20 cm depth (P<0.01).

    In 2018 and 2019, soil porosity under FS treatment was the highest and was the lowest under R0 treatment (Fig. 6b). Furthermore, soil porosity was always higher under flat ground than furrow treatment. Soil porosity for each treatment was the highest at 0-10 and 10-20 cm soil depths, then the values of soil porosity decreased with increasing soil depths. Analysis of variance showed that growth year displayed no significant effect on soil porosity. Effects of treatment on soil porosity at 0-10 and 10-20 cm soil depths reached an extremely significant level (P<0.01).

    In 2018 and 2019, straw mulch treatments exhibited the lowest soil compactness, and control treatment had the highest values (Fig. 6c). Soil compactness of all treatments increased and then decreased with soil depths, with the maximum values at 10-20 cm soil depth and a variation range of 110-152 kPa. Analysis of variance showed that growth year did not significantly affect soil compactness. Treatment did significantly affect soil compactness (P<0.01) at 0-10 cm soil depth. Interaction between growth year and treatment had no significant effect on soil bulk density, soil porosity, and soil compactness (P>0.05).

    3.3 Effects of mulch and planting methods on soil chemical properties

    3.3.1 SOM, soil TN, and AN

    Fig. 5 Effects of mulch and planting methods of M. ruthenica on soil temperature in 2017 (a-e) and 2018 (f-h)

    Fig. 6 Effects of mulch and planting methods of M. ruthenica on soil bulk density (a and b), porosity (c and d),and compactness (e and f) in 2018 and 2019

    Fig. 7 Effects of mulch and planting methods of M. ruthenica on soil organic matter (a and b), total nitrogen (c and d), and available nitrogen (e and f) in 2018 and 2019. Different lowercase letters within the same soil depth indicate significant differences among different treatments at P<0.05 level.

    In 2018 and 2019, SOM content under RS treatment was the highest, and was the lowest under FP treatment; SOM under straw mulch was higher than that of no mulch, and was the lowest for plastic film mulch (S>0>P; Fig. 7a and b). Compared with control treatment, straw mulch increased SOM content, and plastic film mulch reduced it. SOM content of each treatment decreased with soil depths. In 2018 and 2019, SOM content at 0-10 cm soil depth under RS and FS treatments was significantly higher than those of RP and FP treatments (P<0.05). No significant difference was observed in SOM content among different treatments at 10-40 cm soil depth (P>0.05). Growth year significantly influenced SOM content (P<0.01). Effect of treatment on SOM content at 0-10 cm soil depth was extremely significant (P<0.01).

    In 2018 and 2019, soil TN content for each soil depth was the highest under straw mulch,followed by plastic film mulch (S>P), and the lowest soil TN content was observed for no mulch treatment (Fig. 7c and d). TN content presented an inverse trend with soil depths. At 0-40 cm soil depth, no significant difference was observed in soil TN content among different treatments in 2018 and 2019 (P>0.05). Analysis of variance showed that growth year, treatment, and their interaction did not significantly affect soil TN content (P>0.05).

    In 2018 and 2019, soil AN content for each treatment presented a decreasing trend with increasing soil depths (Fig. 7e and f). In 2019, soil AN contents under RS and FS treatments were significantly higher than that of F0 treatment (P<0.05), but no significant difference was observed among other treatments. In 2018 and 2019, no significant difference was observed in soil AN content among treatments at 10-40 cm soil depth (P>0.05). Analysis of variance showed that growth year did not significantly affect the variation of soil AN content. Treatment extremely significantly affected soil AN content at 0-10 cm soil depth (P<0.01).

    3.3.2 Soil TP, AP, TK, and AK

    In 2018 and 2019, soil TP content for each soil depth was the highest for straw mulch, followed by plastic film mulch (S>P), soil TP content decreased with increasing in soil depths (Fig. 8a and b). Analysis of variance showed that growth year did not significantly affect soil TP content,whereas treatment significantly affected soil TP content at 0-10 and 10-20 cm soil depths(P<0.05).In 2018 and 2019, soil AP content for each soil depth was the highest under RS treatment and the lowest under F0 treatment (Fig. 8c and d). Soil AP content under straw mulch was higher than that of plastic film mulch. In 2018, soil AP content under RS treatment at 0-10 cm soil depth was significantly higher than that of F0 treatment (P<0.05). In 2019, soil AP contents under RS and FS treatments were significantly higher than those of R0 and F0 treatments (P<0.05). There was no significant difference among treatments at other soil depths (P>0.05). Growth year did not significantly affect soil AP content, whereas treatment extremely significantly affected soil AP content at 0-10 cm soil depth (P<0.01).

    In 2018 and 2019, soil TK content for each soil depth was the highest under straw mulch,followed by plastic film mulch. Soil TK content presented a decreasing trend with increasing soil depths (Fig. 8e and f). Effect of growth year significantly affected soil TK content at 0-10 cm and 10-20 cm soil depths (P<0.01), and treatment significantly affected soil TK content at 0-10 cm soil depth (P<0.01).

    In 2018 and 2019, soil AK contents under RS and FS treatments were significantly higher at 0-10 and 10-20 cm soil depths than those of other treatments (P<0.05; Fig. 8g and h). There were no significant differences in soil AK content at 20-30 and 30-40 cm soil depths (P>0.05). Growth year did not significantly affect soil AK content. Effects of treatment at 0-10 cm and 10-20 cm soil depths were extremely significant (P<0.01). Interaction between growth year and treatment did not significantly affect soil TP, AP, TK, and AK contents (P>0.05).

    4 Discussion

    4.1 Effects of mulch and planting methods on M. ruthenica seed yield

    We found that seed yield ofM. ruthenicaover the three-year experiment was the highest under plastic film mulch and the lowest under no mulch. Since the surface cover treatment and the fact that ridge and furrow planting methods resulted in the accumulation of natural rainfall, the seed yield may have increased by improving the water utilization rate of the plants. The ridge funnels the rainfall into furrow, and in the early stage of crop growth, the water retention effect of plastic film mulch was higher than that of straw mulch (Fig. 4). Studies have shown that increased temperature and water availability are associated with increased crop growth under mulch conditions. Under the effect of increasing temperature and promoting water consumption, mulch cultivation can expand crop growth periods and increase dry weight per plant (Li et al., 2018;Chen et al., 2019). Therefore, seed yield ofM. ruthenicatreated with plastic mulch was higher than that of straw mulch. In the first year of this study (2017), plastic mulch resulted in the highest soil moisture content (Fig. 4) and soil temperature (Fig. 5), which explains the result of the highest seed yield under FP and RP treatments. Additionally,M. ruthenicaseed yield was the highest in the second year (2018).

    In this study, an enormous gap was observed between performance seed yield and actual seed yield. Results showed that potential seed yield of forage crops is extremely high, but actual seed yield is extremely low, with a gap of approximately 12%-20% (Han et al., 1996). Additionally,the splitting rate ofM. ruthenicais 93% when the seeds are fully mature, with a relatively high proportion of splitting pods (Li et al., 2006). During the harvesting process in the present study,seed pod splitting and maturity were inconsistent, resulting in substantial loss and a discrepancy between performance and actual seed yields.

    All mulch treatments displayed desirable effects on soil microenvironment compared with that of control, mainly by increasing soil moisture content and nutrient content. Although plastic film mulch treatment was associated with high seed yield, soil moisture content, and temperature,plastic film is remained in the soil and degraded slowly, which is detrimental to the emergence and normal growth of crops (Yan et al., 2006). Plastic film is potentially detrimental to soil physical and chemical properties, including the white pollution caused by film residue, which has become a source of soil pollution (Qi et al., 2020; Xu et al., 2020). Straw mulch presented similar,if not equal, effects on seed field and soil properties when compared with that of plastic, and in certain cases resulted in higher soil nutrient contents. Therefore, for the sustainable development ofM. ruthenicain the arid Longzhong Area, straw mulch method should be prioritized.

    4.2 Effects of mulch and planting methods on soil physical properties

    4.2.1 Soil moisture content

    In this study, both plastic film and straw mulch effectively increased soil moisture content. In 2017, RP and FP treatments resulted in the highest soil moisture content, followed by RS and FS treatments. The higher soil moisture content for ridge and furrow treatment instead of flat groud method may have been due to the fact that ridge planting effectively collects rainfall and shelters the growth area from wind evaporation, thereby improving water retention compared with that of flat ground method. Soil moisture content and seed yield ofM. ruthenicatreated with plastic film in 2017 were higher than those of other treatments. Plastic film mulch inhibits soil moisture evaporation, creating internal water circulation (Zhou et al., 2009). Thus, the circulation promotes soil water movement and accumulation from the deeper depths to the surface, and improves water retention (Liang et al., 1990). In this study, soil moisture content of straw mulch was higher than that of plastic film mulch in 2018. Straw mulch inhibited soil water evaporation and ensured that water demand ofM. ruthenicais fulfilled in the later growth period, thus promoting seed production. In addition, the roughness of straw covering surface can effectively intercept rainfall,decrease surface runoff, and allow water to effectively infiltrate the soil, thus increasing soil moisture content (Li et al., 2013).

    Due to the high soil evaporation that occurred in the study area, soil moisture content first increased and then decreased during the growth period ofM. ruthenica(Peng et al., 2018).Moreover, cool weather and increased rainfall occurred in mid-September 2017 (Fig. 4), soil moisture content presented a decreasing trend. We found that soil moisture contents under RS and FS treatments were higher than those of RP and FP treatments, and soil moisture content under plastic film mulch decreased with increasing soil depths. Therefore, we can conclude that the effects of plastic film mulch on increasing soil temperature and moisture were better than that of straw mulch, thus significantly improving the growth ofM. ruthenicain 2017. However, rapid growth ofM. ruthenicaconsumes more soil moisture, causing dehydration in the later growth period (Lu et al., 2016). Simultaneously, plastic film mulch encourages water to move upward from deeper soil depths, resulting in the depletion of soil moisture at deeper depths. Therefore,plastic mulch is not conducive to the long-term sustainable development ofM. ruthenica(Zaongo et al., 1997; Liet al., 1999).

    4.2.2 Soil temperature

    Soil temperature directly affects the development of plants and the emergence of seeds (Peng et al., 2018), and can be influenced by land cover materials and planting methods. In this study, the daily soil temperature of plastic film mulch was higher than that of straw mulch, which suggests that plastic film can effectively absorb solar energy and transmit heat into the soil surface. When air temperature drops, the retained heat at the lower soil depths is distributed toward the surface.Plastic mulch traps this heat at the surface, effectively providing thermal insulation (Chen, 2012;Jiang et al., 2016). Thermal insulation effect of plastic film mulch benefited the emergence ofM.ruthenicaseedlings, which was important for the growth and development ofM. ruthenicain the later growth stage. Straw mulch also delays the emission of soil heat into the atmosphere and effectively prevents daily variations in soil temperature, thereby promoting the growth and development of root and stem system, which increases the seed yield ofM. ruthenica(Fig. 3;Tang, 2019).

    4.2.3 Soil bulk density, soil porosity, and soil compactness

    Soil bulk density indirectly reflects the growth of crop roots and SOM accumulation. The lower the soil bulk density, the looser the soil. Soil porosity is also important for the growth of crop roots, because it controls soil moisture, fertility, aeration, and heat (Wang et al., 2015; Li et al.,2019). We found that soil bulk density and soil compactness decreased after mulched by plastic film and straw, whereas soil porosity increased, indicating that surface mulch effectively improved soil structure. Moreover, soil bulk density and soil compactness at 0-20 cm soil depth under all mulch treatments were significantly lower than those under no mulch treatment, and soil bulk density of furrow and ridge treatment was higher than that of flat ground treatment. These results were associated with improved soil moisture content and soil temperature conditions.Under mulch treatment, rainfall directly falling on the ground is blocked, preventing compaction and promoting loose soil structure. However, Wang et al. (2015) found that during heavy rain,water collected in furrow eventually increases soil compactness and soil bulk density, and decreases soil porosity. In this study, soil bulk density under FS treatment was significantly lower than those of R0 and F0 treatments, which is likely attributed to reduced soil porosity caused by the long-term surface mulch. Straw mulch benefits the formation of soil aggregate structure and porosity, and is conducive to promoting proper growth and spatial distribution of plant roots(Bodner et al., 2015; Liaoet al., 2021). Therefore, straw mulch can be used to create an isolated buffer zone for reducing soil bulk density and soil compactness, and increasing soil porosity(Mulumba and Lal, 2008; Xu et al., 2017).

    4.3 Effects of mulch and planting methods on soil chemical properties

    In this study, soil nutrient contents increased under all mulch treatments compared with control,confirming that mulch effectively improved soil fertility and ultimately promoted the crop growth(Donget al., 2017). Moreover, SOM content under straw mulch was higher than that of plastic film mulch. The maintenance of soil moisture and temperature by mulch accelerates the decay and decomposition of straw, releasing organic particles that infiltrate soil and increase SOM at shallow zones (Akhtar et al., 2018). Meanwhile, surface mulch, ridge and furrow treatments improved soil microenvironment.M. ruthenicamight use soil fertility by accelerating the conversion and utilization of soil nutrients. AsM. ruthenicaabsorb nutrients from the soil, SOM content decreased with the increase of mineralization without additional fertilization, but leguminous forage rhizobia nitrogen fixation, in the process of leguminous plant growth can provide the required nitrogen nutrients, thereby reducing the use of nitrogen in the soil. Thus,planting legume forage did not significantly affect TN content in the soil. Meanwhile, nutrient contents decrease over time without additional fertilization, in particular, plants under plastic film mulch consume more nutrients than they do under other treatments due to fast growth in the early growth stage, which reduces soil fertility in the later growth stage (Pu et al., 2006). Therefore,successive years of plastic film mulch could eventually exhaust soil nutrients, which is not conducive to maintaining soil fertility or sustainable agriculture.

    5 Conclusions

    In this study, we found that performance and actual seed yields ofM. ruthenicawere the highest for plastic film mulch, followed by straw mulch. Soil moisture content in the first year of plantingM. ruthenicawas the highest under RP and FP treatments, followed by RS and FS treatments.While, in the second year, soil moisture content under straw mulch was higher than those of other treatments. Soil temperature was the highest under FP and RP treatments, followed by F0 and R0 treatments. Soil nutrient contents were the highest under RS and FS treatments, and generally deceased with soil depths. Meanwhile, mulch had a negligible effect on soil total nutrient contents,but significantly affected available nutrient contents. All treatments improved shallow soil physical-chemical properties as compared with that of control, and straw mulch presented the best effect. Comprehensive data analysis demonstrated that straw mulch was preferred to produceM.ruthenicaseeds in the Longzhong Area and other similar areas with appropriate potassium fertilizer and phosphate fertilizer supplement.

    Acknowledgements

    This study is funded by the earmarked fund for China Agriculture Research System (CARS) (CARS-34).

    美女国产视频在线观看| 能在线免费看毛片的网站| 少妇人妻久久综合中文| 日本免费在线观看一区| 中文资源天堂在线| 国产伦理片在线播放av一区| 观看免费一级毛片| 18禁动态无遮挡网站| 国产精品一区www在线观看| 女人十人毛片免费观看3o分钟| 成人漫画全彩无遮挡| 久久久精品94久久精品| 国产一区二区亚洲精品在线观看| 国产精品人妻久久久影院| 欧美bdsm另类| 亚洲欧洲国产日韩| 99热国产这里只有精品6| 在线观看av片永久免费下载| 日韩欧美 国产精品| 熟女电影av网| 老司机影院毛片| 国产探花极品一区二区| 免费人成在线观看视频色| 80岁老熟妇乱子伦牲交| 亚洲无线观看免费| 欧美高清成人免费视频www| 国产又色又爽无遮挡免| 日本熟妇午夜| 亚洲最大成人中文| 免费大片黄手机在线观看| 女人久久www免费人成看片| 五月开心婷婷网| 久久女婷五月综合色啪小说 | 国产成人freesex在线| 亚洲无线观看免费| 精品久久久久久久人妻蜜臀av| 日韩一区二区视频免费看| 久热这里只有精品99| 一级毛片aaaaaa免费看小| 免费少妇av软件| 国产乱人偷精品视频| 精品亚洲乱码少妇综合久久| 久久精品国产鲁丝片午夜精品| 亚洲美女搞黄在线观看| 干丝袜人妻中文字幕| 久久人人爽人人片av| 高清欧美精品videossex| 少妇猛男粗大的猛烈进出视频 | 男女无遮挡免费网站观看| 成年免费大片在线观看| 婷婷色综合www| 国产亚洲av嫩草精品影院| 能在线免费看毛片的网站| 18禁在线无遮挡免费观看视频| 黄色怎么调成土黄色| 国产精品三级大全| 高清av免费在线| 亚洲精品国产成人久久av| 美女xxoo啪啪120秒动态图| 我要看日韩黄色一级片| 青春草亚洲视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 成人毛片a级毛片在线播放| kizo精华| 成人午夜精彩视频在线观看| 久久久久久久久久成人| 十八禁网站网址无遮挡 | 91在线精品国自产拍蜜月| 91在线精品国自产拍蜜月| 亚洲欧美日韩卡通动漫| 性插视频无遮挡在线免费观看| 能在线免费看毛片的网站| 日本欧美国产在线视频| 老司机影院成人| 久久人人爽人人片av| 亚洲国产成人一精品久久久| 中文字幕免费在线视频6| 国产中年淑女户外野战色| 丝袜脚勾引网站| 亚洲国产精品成人久久小说| 99热6这里只有精品| 国产精品一区www在线观看| 麻豆成人午夜福利视频| 久久久久精品性色| 日韩 亚洲 欧美在线| 国产黄色视频一区二区在线观看| 午夜激情福利司机影院| 最近的中文字幕免费完整| 亚洲欧美一区二区三区国产| 乱系列少妇在线播放| 国产亚洲一区二区精品| 国产精品一二三区在线看| av国产久精品久网站免费入址| 久久久久久久精品精品| 亚洲精品久久午夜乱码| 一二三四中文在线观看免费高清| 亚洲高清免费不卡视频| 亚洲va在线va天堂va国产| 亚洲婷婷狠狠爱综合网| 最后的刺客免费高清国语| 国产精品熟女久久久久浪| 麻豆乱淫一区二区| 免费在线观看成人毛片| 麻豆乱淫一区二区| 蜜桃亚洲精品一区二区三区| 国产乱人偷精品视频| 免费人成在线观看视频色| 高清视频免费观看一区二区| 人妻制服诱惑在线中文字幕| 亚洲色图av天堂| 国产白丝娇喘喷水9色精品| 插逼视频在线观看| 久久99精品国语久久久| 免费不卡的大黄色大毛片视频在线观看| 人妻夜夜爽99麻豆av| 三级国产精品欧美在线观看| 寂寞人妻少妇视频99o| 国产免费一级a男人的天堂| 国产精品国产三级专区第一集| 日本色播在线视频| 一区二区av电影网| 国产精品秋霞免费鲁丝片| 少妇被粗大猛烈的视频| 国产精品人妻久久久影院| 自拍欧美九色日韩亚洲蝌蚪91 | 大又大粗又爽又黄少妇毛片口| 美女被艹到高潮喷水动态| 岛国毛片在线播放| 我的老师免费观看完整版| 在线观看av片永久免费下载| 亚洲自拍偷在线| 亚洲自偷自拍三级| 99九九线精品视频在线观看视频| 国产高清有码在线观看视频| 少妇猛男粗大的猛烈进出视频 | 国产爽快片一区二区三区| 啦啦啦啦在线视频资源| av在线app专区| 国产高清三级在线| 小蜜桃在线观看免费完整版高清| 国产精品人妻久久久久久| 91精品一卡2卡3卡4卡| 久久鲁丝午夜福利片| 国产高清不卡午夜福利| 亚洲av日韩在线播放| 亚洲精品日韩在线中文字幕| 精品一区二区免费观看| 国产伦精品一区二区三区视频9| 久久久国产一区二区| 久久久久久久精品精品| 嫩草影院入口| 99久久精品国产国产毛片| 丝袜脚勾引网站| 欧美日韩亚洲高清精品| 丝瓜视频免费看黄片| 人妻 亚洲 视频| 十八禁网站网址无遮挡 | 国内精品宾馆在线| av线在线观看网站| 亚洲国产成人一精品久久久| 一级av片app| 国语对白做爰xxxⅹ性视频网站| 国产淫片久久久久久久久| 成人亚洲精品av一区二区| 国产色婷婷99| 九草在线视频观看| 国产精品久久久久久久电影| 国产午夜精品久久久久久一区二区三区| 久久女婷五月综合色啪小说 | 亚洲av电影在线观看一区二区三区 | 搡女人真爽免费视频火全软件| 99热网站在线观看| 亚洲精品国产色婷婷电影| 高清欧美精品videossex| 欧美日韩亚洲高清精品| 波多野结衣巨乳人妻| 亚洲电影在线观看av| 亚洲成人中文字幕在线播放| 国产黄a三级三级三级人| 国产亚洲精品久久久com| 亚洲精品456在线播放app| 午夜视频国产福利| 亚洲欧美精品自产自拍| 亚洲自拍偷在线| 亚洲欧洲日产国产| 日本一二三区视频观看| 蜜桃久久精品国产亚洲av| 国产精品国产av在线观看| 麻豆久久精品国产亚洲av| 国产成人91sexporn| 深爱激情五月婷婷| 亚洲va在线va天堂va国产| 国产v大片淫在线免费观看| 精品视频人人做人人爽| 久久女婷五月综合色啪小说 | 免费大片18禁| 亚洲精品国产av成人精品| 亚洲成人中文字幕在线播放| 免费av毛片视频| 一区二区三区四区激情视频| 水蜜桃什么品种好| 精品亚洲乱码少妇综合久久| 国产毛片a区久久久久| 国产精品福利在线免费观看| 七月丁香在线播放| 午夜福利视频1000在线观看| 国产久久久一区二区三区| h日本视频在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国内精品宾馆在线| 中文字幕制服av| 色5月婷婷丁香| 交换朋友夫妻互换小说| 亚洲综合色惰| 啦啦啦中文免费视频观看日本| 亚洲欧美一区二区三区国产| 黄色欧美视频在线观看| 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 久久6这里有精品| 男人狂女人下面高潮的视频| 在线播放无遮挡| 久久99精品国语久久久| 夫妻午夜视频| 国产淫片久久久久久久久| 欧美一级a爱片免费观看看| 夜夜看夜夜爽夜夜摸| 亚洲怡红院男人天堂| 亚洲国产精品成人综合色| 麻豆乱淫一区二区| 联通29元200g的流量卡| 夜夜看夜夜爽夜夜摸| 日韩大片免费观看网站| 天天躁日日操中文字幕| 欧美成人精品欧美一级黄| 久久精品国产亚洲av涩爱| 成人国产麻豆网| 亚洲婷婷狠狠爱综合网| 久久久久久久大尺度免费视频| 欧美老熟妇乱子伦牲交| 91午夜精品亚洲一区二区三区| 爱豆传媒免费全集在线观看| 插阴视频在线观看视频| 直男gayav资源| 国产又色又爽无遮挡免| 午夜福利高清视频| 国产国拍精品亚洲av在线观看| 插阴视频在线观看视频| av线在线观看网站| 免费av不卡在线播放| 国产伦精品一区二区三区视频9| 少妇猛男粗大的猛烈进出视频 | 国产有黄有色有爽视频| 日韩精品有码人妻一区| 少妇猛男粗大的猛烈进出视频 | 国产精品福利在线免费观看| 午夜激情久久久久久久| 中文乱码字字幕精品一区二区三区| 男女那种视频在线观看| 看十八女毛片水多多多| 成人国产麻豆网| 日韩欧美精品免费久久| 国产精品麻豆人妻色哟哟久久| 一级毛片黄色毛片免费观看视频| 精品一区二区免费观看| 亚洲av一区综合| 国产视频内射| 十八禁网站网址无遮挡 | 秋霞在线观看毛片| 一级爰片在线观看| 亚洲欧美精品自产自拍| 五月玫瑰六月丁香| 老司机影院成人| 日韩欧美一区视频在线观看 | 久久精品夜色国产| 午夜免费观看性视频| 国产免费一级a男人的天堂| 高清av免费在线| 久久精品夜色国产| 日日摸夜夜添夜夜添av毛片| 久久99精品国语久久久| 18禁在线播放成人免费| 婷婷色av中文字幕| 日韩成人伦理影院| 精品人妻视频免费看| 国产黄片视频在线免费观看| 女人被狂操c到高潮| 日韩,欧美,国产一区二区三区| 亚洲成人中文字幕在线播放| 综合色丁香网| 亚洲精品国产av蜜桃| 美女脱内裤让男人舔精品视频| 99热网站在线观看| 国产精品嫩草影院av在线观看| 久久久精品免费免费高清| 国产爱豆传媒在线观看| 99九九线精品视频在线观看视频| 99热这里只有是精品在线观看| 久久精品国产亚洲av涩爱| 看非洲黑人一级黄片| 国模一区二区三区四区视频| 直男gayav资源| 亚洲国产高清在线一区二区三| 久久久国产一区二区| 人妻系列 视频| 日本熟妇午夜| 九九久久精品国产亚洲av麻豆| 春色校园在线视频观看| 看黄色毛片网站| 少妇人妻一区二区三区视频| 国产av国产精品国产| 国产乱人偷精品视频| 大话2 男鬼变身卡| 啦啦啦在线观看免费高清www| 免费在线观看成人毛片| 人妻一区二区av| 日韩不卡一区二区三区视频在线| 成年版毛片免费区| 黄片无遮挡物在线观看| 91在线精品国自产拍蜜月| 99久久人妻综合| 午夜福利网站1000一区二区三区| 欧美国产精品一级二级三级 | 少妇人妻 视频| 国内精品宾馆在线| 免费看不卡的av| 在线观看一区二区三区| 亚洲欧美清纯卡通| 99久久精品一区二区三区| 汤姆久久久久久久影院中文字幕| 校园人妻丝袜中文字幕| 精品亚洲乱码少妇综合久久| 纵有疾风起免费观看全集完整版| 深爱激情五月婷婷| 久久精品国产鲁丝片午夜精品| 麻豆久久精品国产亚洲av| 搡老乐熟女国产| 九九爱精品视频在线观看| 欧美三级亚洲精品| 综合色av麻豆| 国产精品蜜桃在线观看| 大陆偷拍与自拍| 亚洲国产欧美人成| 久热久热在线精品观看| 久久久久久久久久人人人人人人| av又黄又爽大尺度在线免费看| 亚洲精品第二区| .国产精品久久| 国产成人精品婷婷| 欧美3d第一页| 在线天堂最新版资源| 噜噜噜噜噜久久久久久91| 免费看日本二区| 国产精品久久久久久久电影| 国产精品爽爽va在线观看网站| 熟妇人妻不卡中文字幕| 国产高潮美女av| 午夜亚洲福利在线播放| 久久久久久久久久久丰满| 成人毛片60女人毛片免费| 久久影院123| 久久久久久久大尺度免费视频| 日本-黄色视频高清免费观看| 久久ye,这里只有精品| 欧美高清性xxxxhd video| 交换朋友夫妻互换小说| 视频中文字幕在线观看| 观看免费一级毛片| 久久99精品国语久久久| 精品久久久久久久末码| 亚洲丝袜综合中文字幕| 26uuu在线亚洲综合色| 精品一区二区三区视频在线| 人妻一区二区av| 99久久中文字幕三级久久日本| 99热6这里只有精品| 亚洲成色77777| 久久ye,这里只有精品| 亚洲美女搞黄在线观看| 日韩一本色道免费dvd| 大片免费播放器 马上看| 亚洲久久久久久中文字幕| 欧美精品一区二区大全| 丝袜美腿在线中文| av专区在线播放| 最新中文字幕久久久久| 乱系列少妇在线播放| 人人妻人人澡人人爽人人夜夜| 一边亲一边摸免费视频| 国产一级毛片在线| 夜夜看夜夜爽夜夜摸| 18禁在线播放成人免费| 最近最新中文字幕大全电影3| 蜜臀久久99精品久久宅男| 一级毛片我不卡| 麻豆成人av视频| 一本色道久久久久久精品综合| 免费在线观看成人毛片| av福利片在线观看| av一本久久久久| 精品人妻一区二区三区麻豆| av在线观看视频网站免费| 亚洲国产最新在线播放| 国产精品人妻久久久久久| 少妇人妻久久综合中文| 少妇熟女欧美另类| 大陆偷拍与自拍| 听说在线观看完整版免费高清| 国产高清三级在线| 哪个播放器可以免费观看大片| 久久国内精品自在自线图片| 亚洲无线观看免费| 亚洲精品日韩av片在线观看| 国产老妇伦熟女老妇高清| 欧美激情久久久久久爽电影| 一边亲一边摸免费视频| 成人亚洲精品av一区二区| 久久热精品热| 一区二区三区免费毛片| 中文资源天堂在线| 日韩三级伦理在线观看| 久热久热在线精品观看| 成人毛片a级毛片在线播放| 国产成人a∨麻豆精品| 久久久精品免费免费高清| www.色视频.com| 欧美3d第一页| 中文天堂在线官网| 好男人视频免费观看在线| 国产伦在线观看视频一区| 一本色道久久久久久精品综合| 婷婷色av中文字幕| 亚洲国产高清在线一区二区三| 日本av手机在线免费观看| 日韩欧美一区视频在线观看 | 老师上课跳d突然被开到最大视频| 2021天堂中文幕一二区在线观| 禁无遮挡网站| 伊人久久精品亚洲午夜| 啦啦啦在线观看免费高清www| 免费黄网站久久成人精品| 在线观看av片永久免费下载| 精品久久久久久久末码| 国产伦在线观看视频一区| 2018国产大陆天天弄谢| 全区人妻精品视频| 我要看日韩黄色一级片| 女人被狂操c到高潮| 精品人妻偷拍中文字幕| 日韩免费高清中文字幕av| 久久久久久久久久成人| 午夜激情久久久久久久| 91精品国产九色| 精品人妻偷拍中文字幕| 色婷婷久久久亚洲欧美| 五月开心婷婷网| 欧美性猛交╳xxx乱大交人| 午夜精品国产一区二区电影 | 久热这里只有精品99| 在线观看一区二区三区激情| 免费人成在线观看视频色| 夜夜爽夜夜爽视频| 在线播放无遮挡| 插逼视频在线观看| 国产亚洲午夜精品一区二区久久 | 永久免费av网站大全| 亚洲aⅴ乱码一区二区在线播放| 日韩在线高清观看一区二区三区| 97超视频在线观看视频| 91精品国产九色| 日本与韩国留学比较| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在| 卡戴珊不雅视频在线播放| 国产成人午夜福利电影在线观看| 亚洲自偷自拍三级| 亚洲不卡免费看| 国产一区二区三区av在线| 欧美xxxx性猛交bbbb| 人妻夜夜爽99麻豆av| 久久国内精品自在自线图片| 人人妻人人澡人人爽人人夜夜| av国产精品久久久久影院| 国产成人一区二区在线| 一级毛片电影观看| 欧美97在线视频| 国产精品av视频在线免费观看| 成人鲁丝片一二三区免费| 婷婷色麻豆天堂久久| 熟女人妻精品中文字幕| 三级国产精品欧美在线观看| av卡一久久| 亚洲av二区三区四区| 天美传媒精品一区二区| 国产毛片在线视频| 精品一区在线观看国产| av专区在线播放| 日韩av免费高清视频| 国产高清三级在线| 精品国产三级普通话版| 一级av片app| 国内少妇人妻偷人精品xxx网站| 毛片女人毛片| 国产伦在线观看视频一区| 国产精品三级大全| 国产色爽女视频免费观看| 少妇被粗大猛烈的视频| 国产精品伦人一区二区| 国产乱人视频| 国产美女午夜福利| 纵有疾风起免费观看全集完整版| 99久久精品国产国产毛片| 全区人妻精品视频| 能在线免费看毛片的网站| 白带黄色成豆腐渣| 亚洲国产精品成人综合色| 亚洲经典国产精华液单| 91在线精品国自产拍蜜月| 亚洲性久久影院| 亚洲精品国产成人久久av| 久久精品熟女亚洲av麻豆精品| 日本wwww免费看| 久久99热6这里只有精品| 中文天堂在线官网| 亚洲精品国产av成人精品| 熟女电影av网| 久久久久九九精品影院| 在线天堂最新版资源| 插阴视频在线观看视频| 国产av码专区亚洲av| 插阴视频在线观看视频| 九草在线视频观看| 免费观看性生交大片5| 免费看光身美女| 日本爱情动作片www.在线观看| 日韩成人伦理影院| 久久综合国产亚洲精品| 亚洲性久久影院| 亚洲精品久久午夜乱码| 日韩国内少妇激情av| 精品少妇久久久久久888优播| 亚洲无线观看免费| 免费大片18禁| 99re6热这里在线精品视频| 国产黄片视频在线免费观看| 亚洲欧美日韩东京热| 女人十人毛片免费观看3o分钟| 女人久久www免费人成看片| 一个人观看的视频www高清免费观看| 国产黄片美女视频| 在线观看免费高清a一片| 成人免费观看视频高清| 精品久久久久久久末码| 色网站视频免费| 国产精品久久久久久av不卡| 免费av不卡在线播放| 精品久久久精品久久久| 成年女人看的毛片在线观看| 国产在视频线精品| 你懂的网址亚洲精品在线观看| 最近手机中文字幕大全| 夜夜看夜夜爽夜夜摸| 大香蕉97超碰在线| 天堂中文最新版在线下载 | 成年av动漫网址| 亚洲美女视频黄频| 亚洲天堂国产精品一区在线| 久久99热这里只有精品18| 真实男女啪啪啪动态图| 亚洲av福利一区| 波多野结衣巨乳人妻| 国产精品精品国产色婷婷| 久久精品国产亚洲av涩爱| 偷拍熟女少妇极品色| 亚洲欧美日韩东京热| 18禁在线无遮挡免费观看视频| 高清午夜精品一区二区三区| 国产老妇伦熟女老妇高清| 日韩人妻高清精品专区| a级毛色黄片| 婷婷色综合大香蕉| 熟妇人妻不卡中文字幕| 久久久久久伊人网av| 午夜精品一区二区三区免费看| 夫妻午夜视频| 久久久久久久久久成人| 免费观看无遮挡的男女| 欧美最新免费一区二区三区| 99久久精品国产国产毛片| 欧美精品一区二区大全| a级毛色黄片| 日韩欧美精品v在线| 亚洲国产成人一精品久久久| 久久精品人妻少妇| 春色校园在线视频观看| 亚洲av日韩在线播放| 大香蕉97超碰在线| 欧美最新免费一区二区三区| 国产av国产精品国产| 欧美日韩综合久久久久久| 在线免费十八禁| 国产精品福利在线免费观看| 成人亚洲精品av一区二区| av免费在线看不卡| 国产淫语在线视频| 美女xxoo啪啪120秒动态图| 成人免费观看视频高清| 婷婷色综合www| 午夜日本视频在线| 99热这里只有精品一区| 免费不卡的大黄色大毛片视频在线观看| 激情 狠狠 欧美| 特级一级黄色大片| 18禁裸乳无遮挡免费网站照片| 交换朋友夫妻互换小说|