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

    Effect of subsoil tillage depth on nutrient accumulation,root distribution,and grain yield in spring maize

    2014-03-13 05:50:32HonggungCiWeiXiuzhiZhngJieqingPingXiogongYnJinzhoLiuJinghoYunLihunWngJunRen
    The Crop Journal 2014年5期

    Honggung Ci,Wei M,Xiuzhi Zhng,Jieqing Ping,Xiogong Yn,Jinzho Liu,Jingho Yun,Lihun Wng,*,Jun Ren,*

    aInstitute of Agricultural Resource and Environment,Jilin Academy of Agricultural Sciences,Changchun 130333,China

    bInstitute of Crop Science,Chinese Academy of Agricultural Sciences,Beijing 100081,China

    cDepartment of Agronomy,Purdue University,915 W State Street,West Lafayette,IN 47907,USA

    1.Introduction

    Northeast China is the largest spring maize production area in China.In 2011,the spring maize area in northeast China was 12.53 million hectares,accounting for 37.3% of the national spring maize planting area.The overall yield of spring maize in this area was 80.07 million tons,accounting for 41.5% of the national yield [1].Spring maize in the northeastern region is almost completely rain fed and the main factor limiting the yield is lack of soil water [2,3].At present,its farming system of long-term continuous cropping is dominated by small-sized four-wheeled tractors,which are used for stubble removal,soil preparation,sowing,fertilization,cultivation,and other operations [3,4].However,owing to years of overexploitation of the soil as well as improper mechanical manipulation of the soil,the effective topsoil depth has gradually decreased and the plow pan has thickened [5,6].The average effective topsoil depth is only 15.1 cm in northeast China (against an average of 16.5 cm in China) [7],much smaller than that in North America,which reaches 35 cm on average [7].The shallow and compacted topsoil not only restricts the root development of plants,but hinders their absorption of nutrients and water.The poor soil properties also reduce their tolerance to abiotic stress,especially resistance against natural disasters[5,8].

    Subsoil tillage is one of the most effective ways to break up a plow pan in farming management [9–11].It can break up compacted layers,loosening the soil and deepening the topsoil without inverting it,increasing soil permeability [12].As reported by many authors,subsoil tillage can play important roles in efficiently breaking up the plow pan[9,13],promoting water storage in the soil [14,15],adjusting the proportion of solid,liquid,and gas of soil,improving the structure and characteristics of topsoil[15–18],and improving the ecological environment for root development and root activities that enhance the anti-stress capacity of plants[19–21].Varsa et al.[22] concluded that deeper tillage could result in higher corn grain yield,especially in dry seasons.Currently,the time and method for subsoil tillage vary in different regions based on local production operations and included inter tillage [23–25],autumn tillage [26,27],and spring tillage [26–28].In northeastern China,most farmers are accustomed to inter tillage,which is integrated with fertilizer topdressing at the V6 or V7 stage.Wang et al.[5]developed a technique of subsoil inter tillage and water conservation in maize production with high yield and high efficiency.Based on this technique,a local standard was established (DB22/T1237-2011).However,most studies have been performed over only one or two years,with subsoil tillage depths mainly around 30 cm,and the results could not accurately reflect the contribution of the technique to grain yield and nutrient accumulation.

    In the present study,a four-year(2009–2012)experiment in the middle region of northeastern China was conducted to evaluate the impacts of different subsoil tillage depths on soil properties,nutrient (N,P,K) uptake,root morphology,and grain yield of spring maize.The aim was to characterize the effects of subsoil tillage on spring maize and thereby provide useful guidance for soil management coupled with appropriate machinery operations for sustainable crop production in the northeastern region.

    2.Materials and methods

    2.1.Site description

    The experiment was conducted on the research farm of Jilin Academy of Agricultural Sciences,Gongzhuling city,Jilin province (43°29′55″N,124°48′43″E) during 2009–2012.The field had been under continuous conventional maize management for many years before the experiment.During the maize growth season,the meteorological parameters were as follows:average temperature was 19.6 °C,frost-free days were 125–140 days,effective cumulative temperature was 2600–3000 °C,and total sunshine hours were 1220 h.The properties of the black soil in the 0–20 cm plow layers were as follows: organic matter,26.4 mg kg-1; available nitrogen,244 mg kg-1; available phosphorus,35.9 mg kg-1; available potassium,140 mg kg-1; and pH 6.59.The precipitation totals during the maize growing seasons in the years 2009–2012 were 234.2,628.2,320.6,and 519.3 mm,respectively.

    2.2.Design of experiment

    Three tillage treatments were established,consisting of conventional soil management (CK),subsoil tillage to 30 cm depth (treatment T1),and subsoil tillage to 50 cm (treatment T2).The experiment was laid out in a randomized block design with four replicates of each treatment,and each plot was of 140 m2.Conventional soil management was ridge tillage,a long-term continuous maize system,which is dominated by small-sized four-wheeled tractors for soil preparation before sowing.Subsoil tillage was performed with a subsoiling chisel plow in combination with inter tillage in mid-to-late June(V6 stage).Three treatments were applied with basal fertilizer,which comprised 90 kg ha-1N,90 kg ha-1P2O5,and 90 kg ha-1K2O.Pure nitrogen of 135 kg ha-1was added at the 6-expanded-leaves stage (urea with N 46%),phosphate fertilizer as diammonium phosphate (18-46-0),and potassium chloride (K2O 60%).Maize was overseeded on April 25,2009,April 24,2010,April 26,2011,and April 25,2012.At the V3 stage,seedlings were thinned to a density of 60,000 plants ha-1,which is the optimum density for maize hybrids grown in the experimental area.The hybrid was Xianyu 335,which was harvested on September 25,2009,September 24,2010,September 26,2011,and September 24,2012.The experimental area was kept free of weeds,insects and diseases with chemicals based on standard practices.No irrigation was applied.

    2.3.Test items and methods

    Soil samples from the 0–20 cm plow layer were collected before sowing and conventional chemical methods for determining soil nutrient content were used.At the stage of maize physiological maturity,three representative maize plants for each treatment were collected; leaves,stalks,kernels and cobs were divided,dried and crushed;and N,P and K contents for each fraction were determined.Total N content was determined by the micro-Kjeldahl method,total P content was obtained with method of molybdenum–antimony–D-iso-ascorbic-acidcolorimetry (MADAC) and total K content was tested by flame photometry[29].The middle two rows of each plot were harvested at maturity and grain yield was corrected to 14%moisture content.

    A maize root sample was dug with the section sampling method.At the 12-leaf stage (July 4) and early filling stage(August 3),three plants with uniform appearance were selected from each plot for root sampling.Boundaries were separately defined based on the area occupied by each plant,with row distances of 65 cm and plant spaces of 26 cm.Soil and root samples were collected from each 10-cm layer to 80 cm depth.All roots in each soil layer were carefully removed and rinsed with water to remove adhering soil.A 0.05 mm sieve was used to prevent the loss of fine roots during washing.Roots were placed into a zip-locking bag to soak up water and stored at-20 °C.The roots in each layer were scanned with a scanner (Epson V700,Germany) to an image file.The WinRhizoPro5.0 software (Pro2004b,Canada)was used to evaluate root length,surface area,and diameter.The root dry weight of each layer was evaluated after oven drying at 70 °C to constant weight.

    At the 12-leaf and early filling stages,soil samples from the soil layers were collected,and treated with 0.01 mol L-1CaCl2.A TRACCS2000 continuous flow analyzer was used to determine the ammonium and nitrate nitrogen contents of the soil.The Olsen method was used to test readily available phosphorus of the soil and the water content was also measured at the 12-leaf stage [29].A soil hardness tester (Yamanaka type,Japan) was used to measure the soil compaction of the 0–80 cm soil layer at the 12-leaf stage.

    2.4.Data analysis

    Microsoft Excel 2007 software was used for data processing and drawing,and SAS 8.0 statistical software was used for variance analysis and multiple comparisons.

    Table 1-Biomass,grain yield,and yield components under different subsoiling treatments.

    3.Results

    3.1.Dry matter accumulation and yield in maize

    Significant differences in biomass and grain yields were found among the three treatments (Table 1).Under the T1 and T2 treatments,grain yields were increased by 4.2–23.0% with an average of 12.8%and dry biomass was increased by 9.2–24.5%with an average of 14.6%.Based on the yield components,subsoiling was responsible for an increase in grain weight,which,comparing T1 and T2 treatments with the control(CK),were increased by 12.7%and 15.2%,respectively.The number of ears was increased by-0.2–0.7% with an average of 0.4%compared with the control (CK).The kernel number was increased by-0.5–6.3%with an average of 2.7%.There was no significant difference between T1 and T2 treatments.

    Environment(year)had a significant effect on biomass and grain yield and the interaction between year and treatment was also significant (Table 2).There were significant differences in precipitation and rainy period during 2009–2012(Fig.1),which influenced mainly the slight annual differences in yield components.Although rainfall was sufficient in early 2009,the grain weight was reduced by severe drought in later months of that year,resulting in no significant difference between treatments.Heavy precipitation events occurred mainly in late 2010,resulting in lower kernel number and significantly higher grain weight.Under the T1 and T2 treatments,grain weights were increased by 23.7 and 26.7%,respectively,compared to CK treatment.Grain yield and biomass showed a slight difference between treatments owing to increased rainfall in July and August,masking the effect of subsoil tillage.Although precipitation was abundant throughout the development period in 2012,there was severe lodging caused by Typhoon “Bolaven” at the end of August.Grain filling was thereby affected and 100-kernel weight was reduced,in particular under the CK treatment.It was concluded from the results of the four-year experiment that there were no significant differences between different subsoiling depth treatments with respect to dry biomass,yield,or yield components.However,significant differences were observed in 2012,when dry biomass and yield for subsoil tillage to 50 cm were increased by 8.6%and 8.8%respectively,compared with subsoil tillage to 30 cm.

    3.2.Maize nutrient uptake

    As with grain yield and biomass,the year also affected N,P,and K accumulations,and there was significant interaction between year and subsoil tillage treatment (Table 2).Drought inhibited the accumulation of N,P,and K in plants,resulting in lower uptake by plants in 2009.In 2010,the nutrients in soil moved down with heavy rainfall in July and August,leading to reduced N and K absorption by the plant.With respect to nutrient distribution,the increased N and P accumulation under T1 and T2 treatments were dominated by grain(Table 3).Compared to CK,N accumulation in kernels under subsoiling treatments increased by 11.4–29.1% with an average of 16.9%,whereas P accumulation in the grains increased by an average of 10.7%,ranging from 2.0 to 31.9%.Interestingly,there was only a slight difference in K accumulation among the three treatments.Although K accumulations in straw in 2010,2011,and 2012 under subsoil tillage(T1 and T2)were higher than those in CK,there was no significant difference in the grain among the three treatments.N,P,and K accumulations of the maize plant under T1 and T2 treatments were both significantly higher than those under CK treatment in 2010,2011,and 2012 except for the P accumulation in 2012 (P <0.05),which increased by 9.9–22.1%,1.7–20.5%,and 2.1–25.5%,respectively.

    The N,P,and K accumulations under subsoil tillage up to 50 cm increased by 2.7-2.8%,5.0-8.3%,and 1.6-5.2%,respectively,compared to nutrient accumulation under subsoiling to 30 cm,but there were no significant differences between two treatments.With respect to nutrient distribution,the N,P and K contents in the straws under subsoil tillage to 50 cm increased by 4.0%,-1.7%,and-0.7% respectively,compared to those under 30 cm depth;the N,P,and K content in grains under subsoil tillage 50 cm increased by-1.7%,0.2%,and 1.8%respectively,compared to those under 30 cm depth,but no significant differences were detected between two treatments(Table 3).

    Fig.1-Climatological data(active monthly cumulative temperature and precipitation)for the 2009,2010,2011,and 2012 maize growing seasons.

    3.3.Root morphological characteristics around flowering

    The subsoil tillage had no significant effect on root morphology,especially after flowering (Figs.2–5).At the V12 stage,total root length,root surface area,root diameter,and root dry weight in 0–80 cm soil under subsoil tillage treatment increased by 22.9–23.9%,13.9–17.8%,7.4–26.1%,and 1.2–8.8%compared to those under CK treatment.However,at early filling stage,total root length,root surface area,root diameter,and root dry weight in 0–80 cm soil in subsoil treatment were higher than those in CK treatment,with differences of 43.8–49.8%,28.8–36.5%,13.3–21.3%,and 9.1–13.3% compared to those of CK treatment.Between subsoiling depths there were no significant differences in root length,surface area,diameter,or dry weight,but there were significant differences between some soil layers at different depths,especially in deeper soil layers.

    At the 12-leaf stage,the maximum root length was recorded in the 0–10 cm soil layer under CK treatment and was significantly greater than those in subsoil tillage treatments; as deeper soil was sampled,total root length decreased under CK treatment.For example,the root length in the 40–80 cm soil layer accounted for only 9.7%of total root length and was significantly less than those under T1 and T2 treatments (Fig.2).The maximum percentage for the root length reached 19.6% under subsoil tillage to 50 cm,significantly greater than that under subsoiling to 30 cm.Also,at the early filling stage,root length in the 40–80 cm soil layer accounted for 27.3% of the total length under subsoiling to 50 cm.Significant differences were found among the three treatments.

    The distribution of root surface areas in different soil layers was correlated with root length (Fig.3).At the 12-leaf stage,the distribution of root surface areas in different soil layers were as follows: in the CK treatment,66.0% for the 0–20 cm soil layer,21.1%for the 20–40 cm soil layer,and 12.9%for the 40–80 cm soil layer;for the T1 treatment,57.1%for the 0–20 cm soil layer,28.3%for the 20–40 cm soil layer,and 14.6%for the 40–80 cm soil layer;for the T2 treatment,52.0%for the 0–20 cm soil layer,29.1%for the 20–40 cm soil layer and 18.9%for the 40–80 cm soil layer.At the early filling stage,the root surface areas from the 40–80 cm soil layers had increased,in the order T2 >T1 >CK.

    Table 3-Characteristics of nutrient accumulation and distribution in maize under different subsoiling treatments(kg ha-1).

    Fig.2-Root length in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    The trend of proportions of root dry weights in different soil layers was consistent with those for root length and root surface area.But the proportion of root dry weight in the top soil layer (0–20 cm) was higher and the root dry weight in deeper soil layers was lower (Fig.4).At the 12-leaf stage,the percentages of root dry weights in various soil layers were as follows: for CK,72.2% in the 0–20 cm soil layer,17.5% in the 20–40 cm soil layer,and 10.3% in the 40–80 cm soil layer,for subsoiling to 30 cm,66.0% in the 0–20 cm soil layer,20.9% in the 20–40 cm soil layer,and 13.1% in the 40–80 cm soil layer;for subsoiling to 50 cm,60.9%in the 0–20 cm soil layer,22.8%in the 20–40 cm soil layer,and 16.2% in the 40–80 cm soil layer.At the early filling stage,the percentages of root dry weights in the 0–20 cm soil layers under CK,T1 treatment and T2 treatment increased to 82.1%,75.1%,and 74.0%,respectively,but decreased to 8.5%,10.0%,and 12.5%,respectively,in the 40–80 cm soil layer.The percentages of root dry weights also decreased in the 20–40 cm soil layers.

    Based on the comparisons among different treatments,the maximum value for root dry weight was found in the 0–10 cm soil layer under the CK treatment at the 12th leaf and early filling stages,10.6–31.2%greater than those under the T1 and T2 treatments.Significant differences were observed among the three treatments.For the soil layers in the three treatments,the deeper the subsoiled layer,the lower was the root dry weight; however,the root dry weight in CK treatment began to be significantly lower than those under the T1 and T2 treatments in the 30-cm soil layer.No significant differences were found between the root dry weight in the 0–40 cm soil layer under the T1 and T2 treatments,though that under the T1 treatment was slightly higher than that under the T2 treatment.The maximum root dry weight was identified in the 40–80 cm soil layer under the T2 treatment,and was 15.2% and 20.9% higher than those under the T1 treatment at the 12th leaf stage and early filling stages,respectively.There were significant differences between treatments at the early filling stage(Table S1).

    Fig.3-Root surface area in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    Fig.4-Root dry weight in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    Root diameter is an important root morphological parameter and reflects soil influence on the root system.The maximum root diameter under the three treatments was found in the 0–10 cm layer (Fig.5).The root diameter decreased with increasing soil depth.In the top soil layer,the maximum root diameter was found under the CK treatment; in the soil below 20 cm,the maximum value was found under the T2 treatment; at the 12-leaf stage,the variations among root diameters in the 0–80 cm soil layer under the CK,T1 and T2 treatments were 23.7%,13.8%,and 10.0%,respectively.At the early filling stage,the variations were slightly higher,with values of 28.4%,16.9%,and 11.3%for the CK,T1,and T2 treatments.The smallest variation was found under subsoiling to 50 cm,suggesting that subsoiling efficiently breaks up the plow pan,reduces soil resistance to root penetration into deeper soil layers,and promotes root downward growth and uptake of water and nutrients in deeper soil.

    Fig.5-Root diameter in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    3.4.Soil physicochemical properties

    Significant differences in soil compaction in different soil layers across different subsoiling treatments were found(Table 4).Under the CK treatment,lower compaction was found in the 0–10 cm soil layer,but soil compaction significantly decreased in the 10–20 cm soil layer; under the T1 treatment,lower compaction was found in the 0–20 cm soil layer and the soil compaction began to increase significantly below the 30 cm soil layer.Under the T2 treatment,soil compaction gradually increased with soil depth and remained stable to the 40–50 cm soil layer.There was a significant difference in soil compaction between subsoil tillage management and conventional management in the top 0–40 cm soil layer,indicating that subsoiling efficiently breaks up the plow layer and loosens the soil.Thus subsoil tillage management helps reduce soil compaction in deep soil,in turn facilitating plant growth and development.

    After subsoil tillage,the soil was less compact and water content was significantly increased (Fig.6).At the 12-leaf stage,the maximum water content in the 0–40 cm soil layer was found under the T1 treatment,whereas the maximum water content in the 40–80 cm soil layer was found under the T2 treatment and there were significant differences between the CK and T2 treatments.In the 0–80 cm soil layer,the water content of each soil layer under the T1 and T2 treatments was 6.1%higher in average than that under the CK treatment.The difference was increasingly significant with soil depth.At the early filling stage,the advantages of subsoiling were more significant.In the 0–80 cm soil layer,the water contents for the T1 and T2 treatments were both significantly greater than that for the CK treatment.In the 0–80 cm soil layer,the maximum was found under the T2 treatment,and the average water contents under the T1 and T2 treatments were respectively 7.7% and 6.5% greater than that under the CK treatment.

    The total amounts of mineralized N and readily available phosphorus in the 0–80 cm soil layer showed no significant differences across treatments (Fig.7).However,the nutrient distribution in each soil layer differed.Under CK treatment,mineralized N accumulated mostly in the top 0–20 cm soil layer,whereas under the T1 and T2 treatments,soil N mineralization decreased with increasing depth.In the 20–40 cm soil layer,the mineralized N content under subsoiling treatments was markedly higher than that under CK treatment at the 12-leaf stage.In the 40–80 cm soil layer,although the maximum mineralized N content was found under subsoil tillage,no significant differences were found among the three treatments(Fig.7).

    Table 4-Soil compaction in the 0-80 cm soil layer under different subsoiling treatments(kg cm-2).

    Fig.6-Soil water content in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    Significant differences in soil OlsenP under the three treatments were found in the 0–40 cm soil layer,whereas the content was not different across treatments in the lower 40 cm soil layer (Fig.8).At the 12-leaf stage,the maximum value under CK treatment was 33.1 mg kg-1in the top 0–10 cm soil layer.At the early filling stage,the content of OlsenP under CK treatment was substantially decreased,given that roots were distributed mainly in the 20–30 cm top soil layer,which the OlsenP content under CK treatment was markedly higher than those under the subsoil tillage treatments.Up to the 40–80 cm soil layer,readily available phosphorus reached its maximum at the 12-leaf stage under CK treatment,whereas at the early filling stage,no significant differences were found among three treatments.

    Fig.7-Mineralized nitrogen in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    4.Discussion

    A tillage method is an important management strategy in an agricultural production system [30].In recent years,soil management dominated by small-sized four-wheeled tractors has resulted in soil compaction,severe soil water shortage,and soil degradation [3,5].Our four-year study indicated that inter tillage and subsoiling loosen the soil,break up the plow pan caused by multiyear conventional soil management,and enhance root penetration to depth.Subsoil tillage management also reduces soil bulk density[22,28],deepens the active soil layer,and effectively increases soil water storage capacity[15,31].After subsoiling tillage,the proportions of root length and surface area in deeper soil were significantly increased,especially under subsoil tillage to 50 cm (Figs.2–3),owing largely to the increased depth of the subsoil,which promotes root proliferation during the growing season.Two main contributions are root length and root diameter,which result in increased root surface for water and nutrient absorption[32].Dai et al.[33]emphasized that the root distribution under the plow pan may also play a key role in the uptake and utilization of nutrients and water in deep soil,especially after flowering,for the reason that the active layer for nutrient uptake by the root system is then below the 30 cm soil layer[34].At the early filling stage,the uptake capacities for nutrients and water in the soil under the subsoil tillage treatments were greater than that under the CK treatment(Table 3,Fig.6).Subsoil tillage also had positive effect on soil moisture,especially in deep soil,and soil water content was significantly increased below 40 cm,even during a dry season(Fig.S1).Thus,subsoil tillage not only enhances soil water storage capacity but enhances crop uptake of nutrients and water,increasing grain weight [21] and ultimately,grain yield of maize[35,36].

    Fig.8-Olsen-P in the 0-80 cm soil layer under different subsoiling treatments.Left:12-leaf stage;right:early filling stage.

    Table 5-Effect of subsoiling tillage on grain yield,root system,and soil properties in northeastern China.

    The depth of subsoiling is an important cost consideration for farmers.Most of the published papers concerning northeastern China were reviewed and the results suggested no significant difference between 30 and 40 cm subsoiling depths(Table 5).Most studies have been performed over a single year with too-small differences in subsoiling depth to reflect the actual situation.In the present study,no significant differences were observed in N,P,and K accumulations,biomass,yield and components in maize under different subsoil tillage treatments except in 2012.Environment(year)and interaction with subsoiling treatment showed a significant effect on nutrient uptake,plant growth,and grain yield (Table 1).An accurate evaluation of subsoil tillage should be obtained by a long term experiment[15].However,the deeper the subsoiling layer,the more roots developed in deeper soil under the T2 treatment,and root diameter under the T2 treatment was significantly higher than that under the T1 treatment.Our analysis suggests that subsoil tillage as deep as the 50 cm soil layer improves soil physical behavior and reduces soil mechanical resistance to root penetration [22].More importantly,the aboveground stability of plants is enhanced by a well-developed root system.The experiment field experienced Typhoon Bolaven around August 30,2012,and the lodging rates of plants under the CK,T1 and T2 treatments were 14.8%,4.7%,and 0,respectively.Thus lodging resistance and resistance to environmental stress in maize can be markedly improved by deep subsoil tillage,an advantage to be weighed in view of the trend of increasingly frequent natural disasters in the recent years.

    5.Conclusion

    Inter tillage and subsoiling loosened the soil,significantly increased root length,surface area,dry weight,and diameter,and increased the proportion of roots in the 40–80 cm soil layer.The advantages of inter tillage and subsoiling were the delivery of sufficient nutrients for plant growth,facilitation of N,P,and K accumulations in aboveground plant parts,increase in grain weight,and ultimate increase in maize yield.Moreover,subsoiling to increased depths may improve maize root morphology and resistance to environmental stress,especially lodging resistance.

    This study was supported by the National Key Technology R&D Program of China (2012BAD04B02,2013BAD07B02,and 2011BAD16B10),the Special Fund for Agro-Scientific Research in the Public Interest (201103003 and 201303126-4),and the Key Technology R&D Program of Jilin province,China(20126026).

    Supplementary material

    Supplementary material to this article can be found online at http://dx.doi.org/10.1016/j.cj.2014.04.006.

    [1] MOA (Ministry of Agriculture),China Agricultural Yearbook,China Agricultural Press,Beijing (China),2012.

    [2] C.T.Wang,S.K.Li,Assessment of limiting factors and techniques prioritization for maize production in China,Agric.Sci.China 43(2010)1136–1146(in Chinese with English abstract).

    [3] L.P.Zhao,H.B.Wang,H.Q.Liu,Y.L.Wang,S.X.Liu,Y.Wang,Mechanism of fertility degradation of black soil in corn belt of Songliao plain,Acta Pedol.Sin.43(2006) 79–84(in Chinese with English abstract).

    [4] L.P.Zhao,Z.D.Zhang,H.B.Wang,Q.Gao,Y.Wang,S.H.Wang,J.Feng,S.X.Liu,H.Q.Liu,Q.F.Zheng,Fertility evolution characteristics and fertilizing technologies of black soil in the corn belt of the Songliao plain,J.Jilin Agric.Univ.30(2008)511–516 (in Chinese with English abstract).

    [5] L.C.Wang,H.Ma,J.Y.Zheng,Research on rational plough layer construction of spring maize soil in northeast China,J.Maize Sci.16(2008)13–17(in Chinese with English abstract).

    [6] H.B.Wang,L.M.Chen,L.P.Zhao,H.Q.Liu,Y.Wang,Influence of present farming system of corn belt on fertility degradation in Jilin province,Trans.CSAE 25 (2009) 301–305(in Chinese with English abstract).

    [7] S.H.Zhang,S.K.Li,Domestic and Foreign Corn Industrial Technology Development Report,China Agricultural Science and Technology Press,Beijing,2010.

    [8] Y.Z.Wang,W.Y.Li,Z.R.Yin,Progress in relationship both maize root and its yields,J.Jilin Agric.Sci.24(1999) 6–8(in Chinese with English abstract).

    [9] J.He,H.W.Li,H.W.Gao,Subsoiling effect and economic benefit under conservation tillage mode in northern China,Trans.CSAE 22(2006) 62–67(in Chinese with English abstract).

    [10] H.L.Qin,W.S.Gao,Y.C.Ma,L.Ma,Effects of subsoiling on soil moisture under conservation tillage mode,Agric.Sci.China 41(2008) 78–85(in Chinese with English abstract).

    [11] X.J.Yang,C.J.Li,Impacts of mechanical compaction on soil properties,growth of crops,soil-borne organisms and environment,Agric.Sci.China 41 (2008) 2008–2015(in Chinese with English abstract).

    [12] M.Huang,Y.J.Li,J.Z.Wu,M.C.Chen,J.K.Sun,Effects of subsoiling and mulch tillage on soil properties and grain yield of winter wheat,J.Henan Univ.Sci.Technol.27(2006) 74–77(in Chinese with English abstract).

    [13] E.H.Wang,X.W.Chen,Effect of heavy machinery operation on soil three phases and available nutrient in Phaeozem region,J.Soil Water Conserv.21 (2007) 98–102 (in Chinese with English abstract).

    [14] K.L.Ding,M.J.Hann,Influences of subsoiling on physical characters of soil water storage capacity and plant growth,China Rural Water Hydropower 11 (1997) 13–17(in Chinese with English abstract).

    [15] S.D.Evans,M.J.Lindstrom,W.B.Voorhees,J.F.Moncrief,G.A.Nelson,Effect of subsoiling and subsequent tillage on soil bulk density,soil moisture,and corn yield,Soil Tillage Res.38(1996) 35–46.

    [16] L.C.Wang,S.F.Bian,J.Ren,W.R.Liu,Study on technique way to increase unit area yield in the main production zone of spring maize,J.Maize Sci.15 (2007) 133–134 (in Chinese with English abstract).

    [17] D.Sidhu,S.W.Duiker,Soil compaction in conservation tillage: crop impacts,Agron.J.98(2006) 1257–1264.

    [18] N.H.Abu-Hamdeh,Compaction and subsoiling effects on corn growth and soil bulk density,Soil Sci.Soc.Am.J.67(2003) 1213–1219.

    [19] R.Song,C.S.Wu,J.M.Mu,K.Z.Xu,Effect of subsoiling on root growth of maize,J.Jilin Agric.Univ.22(2000) 73–75(in Chinese with English abstract).

    [20] Y.J.E,J.Y.Dai,W.L.Gu,Studies on the relationship between root growth and yield in maize,Acta Agron.Sin.14(1988)149–154 (in Chinese with English abstract).

    [21] H.Qi,M.Liu,W.J.Zhang,Z.P.Zhang,X.F.Li,Z.W.Song,J.L.Yu,Y.N.Wu,Effect of deep loosening mode on soil physical characteristics and maize root distribution,Acta Agric.Boreali-Occid.Sin.27(2012)191–196(in Chinese with English abstract).

    [22] E.C.Varsa,S.K.Chong,J.O.Abolaji,D.A.Farquhar,F.J.Olsen,Effect of deep tillage on soil physical characteristics and corn(Zea mays L.) root growth and production,Soil Tillage Res.43(1997) 219–228.

    [23] J.B.Xiao,Z.X.Sun,J.T.Yang,Y.L.Zhang,J.M.Zhen,Y.Liu,Effect of subsoiling on soil water and crop yield in semi-arid area,Chin.J.Soil Sci.42(2011) 709–713 (in Chinese with English abstract).

    [24] H.J.Chen,S.Y.Gong,J.L.Li,Z.F.Zhang,F.Z.Chen,Y.F.Zhang,Effect of different subsoiling patterns on the yield and soil moisture for early spring maize,Heilongjiang Agric.Sci.10(2011)18–20(in Chinese with English abstract).

    [25] Y.T.Liu,Y.X.Wang,S.Q.Zhang,L.H.Zheng,H.Y.Yang,L.Liu,Y.L.Lian,Effects of different subsoiling patterns on soil moisture and growth of maize,Heilongjiang Agric.Sci.5(2012) 20–24 (in Chinese with English abstract).

    [26] S.F.Bian,H.Ma,F.Xue,S.Yang,G.B.Tan,Research of subsoiling technology in conservation tillage in semi-arid area of Jilin province,J.Maize Sci.8(2000) 67–78(in Chinese with English abstract).

    [27] Q.Q.Meng,J.G.Xie,H.J.Hu,X.C.Wang,X.F.Wang,W.Wu,Effects of subsoiling treatments on maize yield and its component factors,J.Jilin Agric.Sci.25(2000) 25–28(in Chinese with English abstract).

    [28] J.H.Wang,X.J.Gong,Y.Yu,Y.Zhao,Y.B.Jiang,C.R.Qian,Effects of spring subsoiling on soil physical characters and maize yield,Heilongjiang Agric.Sci.12(2011) 16–18(in Chinese with English abstract).

    [29] S.D.Bao,Soil Agricultural Chemical Analysis,China Agriculture Press,Beijing,2000.257–270 (in Chinese).

    [30] Y.Xiong,Q.Xu,X.L.Yao,Z.L.Zhu,Effect of soil tillage systems on soil fertility,Acta Pedol.Sin.17(1980) 101–119(in Chinese with English abstract).

    [31] X.Z.Wang,L.C.Qiu,Influence of subsoiling to moisture of soil and root growth of maize,J.Shenyang Agric.Univ.42(2011)630–633 (in Chinese with English abstract).

    [32] W.G.Liu,L.Shan,X.P.Deng,Effects of soil bulk density on the growth of maize root system under different water condition,Acta Bot.Boreal-Occident.Sin.22(2002) 831–838 (in Chinese with English abstract).

    [33] J.Y.Dai,Y.J.E,W.L.Gu,Studies on the relationship between root growth and yield in maize(Zea mays L.),Acta Agron.Sin.14(1988) 310–314 (in Chinese with English abstract).

    [34] X.D.Zhu,Y.T.Liu,Study on maize root activity and its space distribution in soil,Acta Agric.Nucl.Sin.3(1982) 17–22(in Chinese with English abstract).

    [35] C.H.Li,S.L.Li,Q.Wang,S.P.Hao,J.F.Han,A study on corn root growth and activities at different soil layers with special bulk density,Agric.Sci.China 38(2005)1706–1711(in Chinese with English abstract).

    [36] C.H.Li,P.P.Mei,Q.Wang,S.P.Hao,Influences of soil bulk density in deep soil layers on the uptake and redistribution of nitrogen,phosphorous and potassium in maize,Agric.Sci.China 40(2007) 1371–1378(in Chinese with English abstract).

    [37] W.Zhao,Effects of different subsoiling treatments on maize yield and its component factors,Heilongjiang Agric.Sci.12(2011) 35–37(in Chinese with English abstract).

    [38] J.Huang,A.W.Wang,Y.R.Zhang,Y.Cao,F.Xue,Effects of new cropping system on soil properties of wide and narrow spacing maize rotation planting,stripe deep loosening and leaving high stubble on the ground,Chin.J.Soil Sci.33(2002)168–171 (in Chinese with English abstract).

    [39] J.F.Liang,Q.Z.Qi,X.H.Jia,S.J.Gong,Y.F.Huang,Effects of different tillage managements on soil properties and corn growth,Ecol.Environ.Sci.19(2010) 945–950 (in Chinese with English abstract).

    国产黄频视频在线观看| 精品人妻在线不人妻| 精品99又大又爽又粗少妇毛片| 91成人精品电影| 日韩精品免费视频一区二区三区 | 2022亚洲国产成人精品| 精品人妻一区二区三区麻豆| 99久久人妻综合| 99热这里只有是精品在线观看| 久久午夜福利片| 久久精品国产综合久久久 | 国产免费现黄频在线看| 色网站视频免费| 成人亚洲精品一区在线观看| 99精国产麻豆久久婷婷| 国产成人精品福利久久| 成人免费观看视频高清| 国产白丝娇喘喷水9色精品| 99国产综合亚洲精品| 制服诱惑二区| freevideosex欧美| 国产黄色免费在线视频| 丝袜人妻中文字幕| a级毛片在线看网站| 满18在线观看网站| 精品福利永久在线观看| 国产精品久久久av美女十八| 精品酒店卫生间| 欧美bdsm另类| 永久免费av网站大全| 亚洲av欧美aⅴ国产| 免费高清在线观看日韩| 日本免费在线观看一区| 免费高清在线观看日韩| 国产成人精品无人区| 91午夜精品亚洲一区二区三区| 国产伦理片在线播放av一区| 日日啪夜夜爽| 午夜老司机福利剧场| 欧美国产精品一级二级三级| 国产国语露脸激情在线看| 精品酒店卫生间| 国产午夜精品一二区理论片| 亚洲av在线观看美女高潮| 国产高清国产精品国产三级| 制服诱惑二区| 欧美xxⅹ黑人| 欧美激情国产日韩精品一区| 男女无遮挡免费网站观看| 精品久久久精品久久久| 一区在线观看完整版| 精品国产露脸久久av麻豆| 国产毛片在线视频| 99久久精品国产国产毛片| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 纯流量卡能插随身wifi吗| 亚洲成人手机| 亚洲高清免费不卡视频| 日韩大片免费观看网站| av网站免费在线观看视频| 中文字幕制服av| 国产激情久久老熟女| 成人亚洲精品一区在线观看| 亚洲精华国产精华液的使用体验| 在线亚洲精品国产二区图片欧美| 欧美成人午夜精品| 少妇 在线观看| 咕卡用的链子| 欧美人与性动交α欧美软件 | 免费观看在线日韩| 久久女婷五月综合色啪小说| 久久鲁丝午夜福利片| a级片在线免费高清观看视频| 91成人精品电影| 自线自在国产av| 欧美激情 高清一区二区三区| 欧美精品一区二区大全| 久久久国产一区二区| 成人无遮挡网站| 晚上一个人看的免费电影| 国产精品不卡视频一区二区| 蜜臀久久99精品久久宅男| 啦啦啦在线观看免费高清www| 精品少妇内射三级| 亚洲av在线观看美女高潮| 黄色配什么色好看| 男人操女人黄网站| 综合色丁香网| 丝袜在线中文字幕| 精品少妇黑人巨大在线播放| 免费人妻精品一区二区三区视频| 亚洲情色 制服丝袜| 最近2019中文字幕mv第一页| 97超碰精品成人国产| 亚洲av中文av极速乱| 久久影院123| 亚洲精品成人av观看孕妇| 妹子高潮喷水视频| 亚洲精品久久午夜乱码| 国产黄频视频在线观看| 亚洲内射少妇av| 国产亚洲欧美精品永久| 2018国产大陆天天弄谢| 欧美+日韩+精品| 如何舔出高潮| 午夜福利影视在线免费观看| 校园人妻丝袜中文字幕| 亚洲av福利一区| 人人妻人人添人人爽欧美一区卜| www.熟女人妻精品国产 | 欧美日韩亚洲高清精品| 五月天丁香电影| 亚洲精品国产av成人精品| 丝袜美足系列| 久久热在线av| 久久久精品94久久精品| 中文字幕av电影在线播放| www日本在线高清视频| 国产高清三级在线| av在线观看视频网站免费| 黑人欧美特级aaaaaa片| 乱人伦中国视频| 欧美日韩综合久久久久久| 桃花免费在线播放| 久久人人爽人人片av| 亚洲,欧美精品.| 香蕉国产在线看| 午夜福利视频精品| 精品一区二区三卡| 日本色播在线视频| 国产在视频线精品| 午夜免费男女啪啪视频观看| 五月伊人婷婷丁香| 夜夜骑夜夜射夜夜干| 成人综合一区亚洲| 欧美精品一区二区免费开放| 99九九在线精品视频| 在线观看人妻少妇| 一区二区三区精品91| 久久97久久精品| 久久精品国产亚洲av天美| 99精国产麻豆久久婷婷| 狂野欧美激情性xxxx在线观看| 热re99久久精品国产66热6| 国产一区二区三区综合在线观看 | 国产又爽黄色视频| av播播在线观看一区| 91在线精品国自产拍蜜月| 久久人人97超碰香蕉20202| 亚洲av电影在线进入| 精品久久蜜臀av无| 亚洲精品456在线播放app| 亚洲国产日韩一区二区| 国产精品国产三级国产av玫瑰| 老司机亚洲免费影院| 午夜免费男女啪啪视频观看| 亚洲伊人久久精品综合| 日韩一区二区三区影片| 99精国产麻豆久久婷婷| 中国国产av一级| 亚洲国产欧美日韩在线播放| 永久网站在线| av在线播放精品| 亚洲国产精品国产精品| 1024视频免费在线观看| 日日爽夜夜爽网站| 精品人妻在线不人妻| 免费人妻精品一区二区三区视频| 国产在视频线精品| 亚洲性久久影院| 精品国产乱码久久久久久小说| 18禁在线无遮挡免费观看视频| 少妇精品久久久久久久| 69精品国产乱码久久久| 国产精品嫩草影院av在线观看| 少妇被粗大的猛进出69影院 | 在线免费观看不下载黄p国产| 国产一区二区激情短视频 | 最新的欧美精品一区二区| 日日摸夜夜添夜夜爱| 精品一区二区三卡| 精品一区二区免费观看| 人妻系列 视频| 90打野战视频偷拍视频| 久久久国产一区二区| 国产一级毛片在线| 欧美精品人与动牲交sv欧美| 亚洲婷婷狠狠爱综合网| 尾随美女入室| 国产精品人妻久久久影院| 日韩一区二区三区影片| 亚洲国产精品一区二区三区在线| 一区二区三区四区激情视频| 精品一品国产午夜福利视频| 久久精品久久久久久噜噜老黄| 亚洲av免费高清在线观看| 欧美+日韩+精品| 亚洲精品乱久久久久久| 激情视频va一区二区三区| 欧美最新免费一区二区三区| 青青草视频在线视频观看| 伊人久久国产一区二区| 国产精品国产三级专区第一集| 久久ye,这里只有精品| 亚洲欧洲国产日韩| 91精品伊人久久大香线蕉| av黄色大香蕉| 22中文网久久字幕| 国产片内射在线| 99热6这里只有精品| 国产一区二区三区综合在线观看 | 侵犯人妻中文字幕一二三四区| tube8黄色片| 日韩在线高清观看一区二区三区| 边亲边吃奶的免费视频| 亚洲精品乱久久久久久| 久久午夜综合久久蜜桃| 春色校园在线视频观看| 精品亚洲成a人片在线观看| 午夜影院在线不卡| 男的添女的下面高潮视频| 亚洲av电影在线观看一区二区三区| 久久国内精品自在自线图片| 欧美成人午夜精品| 国产一区有黄有色的免费视频| 亚洲av电影在线进入| 一区二区三区精品91| 欧美国产精品一级二级三级| 在线免费观看不下载黄p国产| 久久国产精品男人的天堂亚洲 | 啦啦啦中文免费视频观看日本| 女性被躁到高潮视频| 狠狠精品人妻久久久久久综合| videosex国产| 最黄视频免费看| 国产在视频线精品| 国产白丝娇喘喷水9色精品| 男女边吃奶边做爰视频| 免费少妇av软件| 哪个播放器可以免费观看大片| 美女国产高潮福利片在线看| 香蕉丝袜av| 免费大片18禁| 两个人免费观看高清视频| 亚洲欧洲国产日韩| 热re99久久精品国产66热6| 男女免费视频国产| 久久精品国产综合久久久 | 欧美日韩国产mv在线观看视频| 18+在线观看网站| 久久久久久伊人网av| a级毛色黄片| 国产成人精品久久久久久| 91在线精品国自产拍蜜月| 两个人免费观看高清视频| 成人国产av品久久久| 久久午夜综合久久蜜桃| h视频一区二区三区| 亚洲av成人精品一二三区| 最近的中文字幕免费完整| 欧美日韩综合久久久久久| 久久久久久久精品精品| 97超碰精品成人国产| xxx大片免费视频| 在线亚洲精品国产二区图片欧美| 亚洲色图 男人天堂 中文字幕 | 中文字幕另类日韩欧美亚洲嫩草| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 女人精品久久久久毛片| 久久久久久久大尺度免费视频| 亚洲激情五月婷婷啪啪| 精品一区二区三区视频在线| 香蕉丝袜av| 啦啦啦视频在线资源免费观看| kizo精华| 日本猛色少妇xxxxx猛交久久| 亚洲美女搞黄在线观看| 国产熟女午夜一区二区三区| 国产无遮挡羞羞视频在线观看| 成人毛片60女人毛片免费| 大片免费播放器 马上看| 99re6热这里在线精品视频| 精品久久蜜臀av无| 亚洲国产毛片av蜜桃av| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线| 亚洲熟女精品中文字幕| 日韩一区二区三区影片| 丝袜人妻中文字幕| 亚洲精品成人av观看孕妇| 满18在线观看网站| 热99国产精品久久久久久7| 成人免费观看视频高清| 91在线精品国自产拍蜜月| a级片在线免费高清观看视频| 一级片'在线观看视频| 精品人妻一区二区三区麻豆| 日韩不卡一区二区三区视频在线| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 久久人人爽人人片av| 永久免费av网站大全| 亚洲欧美清纯卡通| 美女脱内裤让男人舔精品视频| 啦啦啦视频在线资源免费观看| 免费在线观看完整版高清| 黑人猛操日本美女一级片| 日韩视频在线欧美| 日韩中文字幕视频在线看片| 丝袜喷水一区| 桃花免费在线播放| 国产亚洲av片在线观看秒播厂| 欧美日韩精品成人综合77777| 一边摸一边做爽爽视频免费| 超色免费av| 亚洲人与动物交配视频| 1024视频免费在线观看| 亚洲欧美日韩另类电影网站| 久久久亚洲精品成人影院| 欧美 亚洲 国产 日韩一| 久久久久久久久久久免费av| 啦啦啦视频在线资源免费观看| 高清视频免费观看一区二区| 午夜免费鲁丝| 久久鲁丝午夜福利片| 人人澡人人妻人| 国产精品不卡视频一区二区| 欧美性感艳星| 建设人人有责人人尽责人人享有的| 肉色欧美久久久久久久蜜桃| 一级毛片我不卡| 久久女婷五月综合色啪小说| 99九九在线精品视频| 国产白丝娇喘喷水9色精品| 亚洲 欧美一区二区三区| 亚洲,欧美,日韩| 亚洲国产精品一区三区| 欧美激情 高清一区二区三区| 激情五月婷婷亚洲| 少妇被粗大的猛进出69影院 | 日本91视频免费播放| 国产欧美日韩综合在线一区二区| 亚洲国产成人一精品久久久| 亚洲欧美成人综合另类久久久| 亚洲国产成人一精品久久久| 日韩精品有码人妻一区| 九九爱精品视频在线观看| 欧美人与性动交α欧美精品济南到 | 久久99蜜桃精品久久| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 少妇人妻精品综合一区二区| 亚洲av中文av极速乱| 最近的中文字幕免费完整| 22中文网久久字幕| 少妇熟女欧美另类| 免费人妻精品一区二区三区视频| 国产精品久久久久久精品古装| 亚洲美女黄色视频免费看| 韩国av在线不卡| 最新的欧美精品一区二区| 99热6这里只有精品| 午夜日本视频在线| 麻豆乱淫一区二区| 午夜av观看不卡| 999精品在线视频| 国产精品偷伦视频观看了| 久久99热6这里只有精品| freevideosex欧美| 少妇熟女欧美另类| 国产精品一国产av| 久久人人97超碰香蕉20202| 亚洲欧洲国产日韩| 国产亚洲精品第一综合不卡 | 人妻 亚洲 视频| 久久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 在线观看三级黄色| 午夜视频国产福利| 青春草亚洲视频在线观看| 亚洲国产精品999| 人人澡人人妻人| 久久精品国产亚洲av涩爱| 亚洲情色 制服丝袜| 99国产精品免费福利视频| 视频中文字幕在线观看| 伦理电影大哥的女人| 一二三四中文在线观看免费高清| 久久久久久久久久久免费av| 日韩免费高清中文字幕av| 精品国产露脸久久av麻豆| 春色校园在线视频观看| 欧美精品人与动牲交sv欧美| 欧美xxxx性猛交bbbb| 18在线观看网站| 国产69精品久久久久777片| 两性夫妻黄色片 | 伊人亚洲综合成人网| 亚洲欧美一区二区三区国产| 国产极品天堂在线| 成人二区视频| 草草在线视频免费看| 新久久久久国产一级毛片| 一本色道久久久久久精品综合| 午夜影院在线不卡| 男女国产视频网站| 亚洲在久久综合| 久久99热这里只频精品6学生| 免费看光身美女| 亚洲 欧美一区二区三区| videosex国产| 中文字幕精品免费在线观看视频 | 建设人人有责人人尽责人人享有的| 如日韩欧美国产精品一区二区三区| 啦啦啦啦在线视频资源| 国产一区二区在线观看日韩| 国产在线一区二区三区精| 亚洲成国产人片在线观看| 男女无遮挡免费网站观看| 男女高潮啪啪啪动态图| 久久久久精品人妻al黑| 亚洲精品av麻豆狂野| 成人亚洲欧美一区二区av| 国产精品久久久久成人av| 在线看a的网站| 欧美日韩av久久| 欧美亚洲 丝袜 人妻 在线| 性色av一级| 男人操女人黄网站| 十八禁高潮呻吟视频| 在线观看三级黄色| 国产精品不卡视频一区二区| 一二三四中文在线观看免费高清| av.在线天堂| 国产精品99久久99久久久不卡 | 成人毛片60女人毛片免费| 少妇 在线观看| 国产精品无大码| 女的被弄到高潮叫床怎么办| 精品久久久久久电影网| 国产淫语在线视频| 人人澡人人妻人| 日韩熟女老妇一区二区性免费视频| 一本久久精品| 日韩av在线免费看完整版不卡| 国产无遮挡羞羞视频在线观看| 成人国产麻豆网| tube8黄色片| 国国产精品蜜臀av免费| 捣出白浆h1v1| 婷婷色综合www| 天天躁夜夜躁狠狠久久av| 麻豆精品久久久久久蜜桃| 9色porny在线观看| 亚洲精品一二三| xxxhd国产人妻xxx| 色网站视频免费| 亚洲欧美中文字幕日韩二区| 久久99精品国语久久久| 久久人人97超碰香蕉20202| 各种免费的搞黄视频| 国产av一区二区精品久久| 久久综合国产亚洲精品| 亚洲成人av在线免费| 国产精品久久久久成人av| 亚洲一区二区三区欧美精品| 精品99又大又爽又粗少妇毛片| 美女国产视频在线观看| 久久国产精品男人的天堂亚洲 | 亚洲av电影在线进入| www日本在线高清视频| 国产极品天堂在线| 亚洲图色成人| 国产国拍精品亚洲av在线观看| 国产精品人妻久久久久久| 啦啦啦视频在线资源免费观看| 最近中文字幕2019免费版| 我要看黄色一级片免费的| 亚洲精品自拍成人| 在线亚洲精品国产二区图片欧美| 大码成人一级视频| kizo精华| 久久午夜福利片| 国产免费一级a男人的天堂| 欧美丝袜亚洲另类| av免费在线看不卡| 国产av一区二区精品久久| 亚洲 欧美一区二区三区| 狠狠婷婷综合久久久久久88av| 国产成人aa在线观看| 欧美日韩精品成人综合77777| 一个人免费看片子| 成人影院久久| 国产男女超爽视频在线观看| 久久精品夜色国产| 搡女人真爽免费视频火全软件| 亚洲国产看品久久| 青青草视频在线视频观看| 99精国产麻豆久久婷婷| 韩国av在线不卡| 一区二区av电影网| 免费人成在线观看视频色| 校园人妻丝袜中文字幕| 精品熟女少妇av免费看| 久久99热6这里只有精品| www.色视频.com| 赤兔流量卡办理| 美女福利国产在线| 交换朋友夫妻互换小说| 在线亚洲精品国产二区图片欧美| 大片电影免费在线观看免费| 极品少妇高潮喷水抽搐| 欧美少妇被猛烈插入视频| 女性生殖器流出的白浆| 久久久久网色| 制服诱惑二区| 国产又爽黄色视频| 久久精品久久精品一区二区三区| 亚洲av日韩在线播放| 亚洲精品456在线播放app| 日韩三级伦理在线观看| 亚洲在久久综合| 久久国内精品自在自线图片| 香蕉国产在线看| 精品亚洲成国产av| 黄片播放在线免费| 成年美女黄网站色视频大全免费| 婷婷色av中文字幕| 亚洲av福利一区| 欧美+日韩+精品| 亚洲四区av| 波多野结衣一区麻豆| 国产xxxxx性猛交| 色视频在线一区二区三区| 日本-黄色视频高清免费观看| 人妻系列 视频| 一级毛片 在线播放| 中国美白少妇内射xxxbb| 午夜福利影视在线免费观看| 91国产中文字幕| 一个人免费看片子| 黄色怎么调成土黄色| 精品熟女少妇av免费看| 久热久热在线精品观看| 国产一区亚洲一区在线观看| 免费看光身美女| a级毛片黄视频| 国产极品粉嫩免费观看在线| 午夜视频国产福利| 午夜91福利影院| freevideosex欧美| 日本91视频免费播放| 国产成人av激情在线播放| 国产精品久久久久久久电影| 久久久精品区二区三区| 亚洲国产欧美在线一区| 国产成人精品福利久久| 久久久久精品性色| 亚洲国产av新网站| kizo精华| 国产亚洲精品久久久com| 中文欧美无线码| 免费少妇av软件| 日日摸夜夜添夜夜爱| 日韩一本色道免费dvd| 韩国精品一区二区三区 | av有码第一页| 97人妻天天添夜夜摸| 视频区图区小说| 午夜视频国产福利| 精品久久久精品久久久| 欧美日韩视频精品一区| 夜夜爽夜夜爽视频| 国产精品久久久久久av不卡| 蜜臀久久99精品久久宅男| 国产精品麻豆人妻色哟哟久久| 在线观看美女被高潮喷水网站| 在线看a的网站| 内地一区二区视频在线| 一级毛片电影观看| 全区人妻精品视频| 18禁动态无遮挡网站| 男女高潮啪啪啪动态图| 香蕉丝袜av| 欧美日韩亚洲高清精品| 男女国产视频网站| 欧美日韩一区二区视频在线观看视频在线| 蜜桃国产av成人99| 男女边吃奶边做爰视频| 精品国产露脸久久av麻豆| 国产一级毛片在线| 国产精品欧美亚洲77777| 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| 亚洲内射少妇av| 久久久久久久久久人人人人人人| 麻豆精品久久久久久蜜桃| 亚洲精品乱久久久久久| 国产成人精品福利久久| 99热全是精品| 777米奇影视久久| 精品酒店卫生间| 麻豆精品久久久久久蜜桃| 亚洲综合色网址| 国产男女超爽视频在线观看| 国产精品久久久久久精品古装| 男女边吃奶边做爰视频| 一级毛片黄色毛片免费观看视频| 18禁动态无遮挡网站| 国产精品一区二区在线不卡| 最近最新中文字幕大全免费视频 | 大片电影免费在线观看免费| 夫妻性生交免费视频一级片|