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

    Optimized nitrogen application methods to improve nitrogen use efficiency and nodule nitrogen fixation in a maize-soybean relay intercropping system

    2018-03-07 11:39:59YONGTaiwenCHENPingDONGQianDUQingYANGFengWANGXiaochunLlUWeiguoYANGWenyu
    Journal of Integrative Agriculture 2018年3期

    YONG Tai-wen, CHEN Ping, DONG Qian, DU Qing, YANG Feng, WANG Xiao-chun, LlU Wei-guo,YANG Wen-yu

    1 College of Agronomy, Sichuan Agricultural University/Key Laboratory of Crop Eco-physiology and Farming System in Southwest,Ministry of Agriculture, Chengdu 611130, P.R.China

    2 Shehong Farm Bureau, Suining 629200, P.R.China

    1.lntroduction

    In China, the continuous increase of food production largely depends on chemical nitrogen (N) fertilizer inputs(Zhu and Chen 2002).However, the overuse of N fertilizer and low N use efficiency (NUE) leads to wasted resources and environmental pollution and is contrary to sustainable agricultural production (Kwong et al.2002; Xing and Zhu 2002; Wang et al.2017).Surplus N fertilization results in substantial soil acidification, with N leaching and emissions as direct triggers of water and air pollution (Guo et al.2010;Zhao et al.2012; Zhang et al.2013), which is contrary to environmental-friendly agricultural production.

    As the main food and economic crops in China, maize and soybean are widely cultivated.However, during maize cropping seasons from 1980 to 2010 in China, the NUE declined from 30.2 to 29.9 kg grain kg-1N; in contrast, NUE improved from 39.4 to 53.2 kg grain kg-1N in the USA (Yu et al.2015).To achieve both food security and agricultural sustainable development, optimal N management practices are urgently needed.Suitable N nutrient management and cropping systems can improve nutrient use efficiency and utilization of solar energy while being environmentally friendly.

    Efficient N management technology has become one of the most urgent requirements for sustainable agriculture in China.Researchers are trying to reduce N input directly by decreasing chemical N fertilization rates (Liu et al.2011;Ruan et al.2013; Yan et al.2013) and combining chemical N with organic fertilizer or new fertilizer types (slow release or controlled release) to meet crop demand (Shoji et al.2001;Fernández-Escobar et al.2004; Meng et al.2005; Noellsch et al.2009).Previous studies have shown that optimal nutrient management strategies can significantly reduce N fertilization rate and increase crop yield, with multiple benefits for agriculture and the environment (Chen et al.2006; Ju et al.2009; Constantin et al.2010; Ruan et al.2013).In wheat-maize rotation cropping system, the total N fertilizer reduced from 430 to 360 kg N ha-1with improved treatments,resulting in a maize yield increase by 7-14%, and a reduction in wheat yield, N2O and NO emissions by 1-2, 7 and 29%,respectively (Liu et al.2011).In addition, reduced N with optimal fertilization practices improve NUE, thereby reducing seasonal cumulative N2O emission (Yan et al.2013).

    Global food demand is increasing with economic development and population growth.Especially, China feeds approximately 20% of the world’s population with 7% of the world’s farm land (Godfray et al.2010; Tilman et al.2011; Larson 2013).In China, the average unit grain yield increases from 1.09 tonne ha-1in 1949 to 6.51 tonne ha-1in 2014 (Zhou 2017).However, large amounts of soybeans are imported annually to meet the grain and soybean oil consumption in China.Relay intercropping and intercropping systems can be high-yielding, land-use efficient, and resource-efficient, and can efficiently control weeds, diseases and pests (Li et al.2006; Corre-Hellou et al.2011; Rodríguez-Navarro et al.2011; Zuo et al.2015; Wang et al.2017; Yang et al.2017).In particular, those systems including soybean crops can increase soybean yield with limited arable land.On one hand, legume-based cropping system can reduce N leaching by 50% compared with conventional systems (Drinkwater et al.1998).Soybeans reduce N input requirements by biological nitrogen fixation,which meets 50-60% of N demand (Salvagiotti et al.2008).On the other hand, fertilization methods have a significant impact on intercropped crops, with better performance when heterogeneous N was supplied at interspecific rather than intraspecific rows (Wu et al.2014).Interspecific facilitation can increase resource use efficiency and land productivity without negative impacts on the environment (Oljaca et al.2000; Li et al.2001, 2003; Hochman et al.2011).Therefore,soybean-based relay intercropping and intercropping systems, with improved fertilization methods can increase soybean yield and decrease environment costs.However,environmental factors, e.g., light, heat, and rainfall, can limit cropping systems.Intercropping systems are employed in areas with two (or three) crops a year, such as the Huang-Huai-Hai region in China (Liu et al.2017), whereas relay intercropping systems are used in areas with one crop a year(or three crops in two years), such as in Southwest China(Yang et al.2017).

    Maize-soybean relay intercropping system can increase NUE, light use efficiency, and land productivity, which has become a major planting pattern in Southwest China (Gao et al.2010; Yan et al.2010; Xiang et al.2012; Yang et al.2014, 2017; Yong et al.2015; Wang et al.2017).Despite these benefits, few studies are focused on the lower N application rate and its application methods to maizesoybean relay intercropping system.The early study showed that crops performed better at lower N application rate (270 kg N ha-1), and the fertilizer placement can increase grain yield (Dong et al.2014).However, the lower N applied rate and fertilization methods on NUE and soybean biological nitrogen fixation are still unclear in the maize-soybean relay intercropping system.Therefore, this study was performed to test the optimized N management strategy for maize-soybean relay intercropping system in Southwest China.It was hypothesized that lower N input with optimized fertilization locations for maize-soybean relay intercropping system would increase NUE, N agronomy efficiency (NAE), and total yield compared with conventional N application rate and location.The objectives of this study were to: (1) investigate NUE, NAE and yield responses of crops to lower N and fertilization locations and to (2)determine ureide, leghemoglobin content, and nitrogenase activity responses of soybean to lower N and fertilization locations.

    2.Materials and methods

    2.1.Experimental site

    The experiment was performed in Renshou County (29°40′-30°16′N(xiāo), 104°00′-104°30′E), Sichuan Province, Southwest China, during the 2012 and 2013 cropping seasons(April-November).The climate of the experimental site is subtropical monsoon humid, with a mean annual temperature of 17.4°C and mean annual rainfall of 1 009.4 mm.The frost-free period lasts approximately 312 days, and the average annual sunshine is 1 196.6 h.The soil is anthrosol with a clay loam texture, and the soil (0-20 cm) nutrient contents are shown in Table 1.

    A semi-compact maize variety (Zea mays L.Chuandan 418)and shade-tolerant soybean variety (Glycine max L.Nandou 12) were chosen as the experimental crops.In the maize-soybean relay intercropping system, a wide-narrow row planting was adopted, with wide rows of 160 cm and narrow rows of 40 cm.The total ratio of maize to soybean rows was 2:2.Maize plants were in the narrow rows with a row spacing of 40 cm, soybean plants in the wide rows with a row spacing of 40 cm.The distance between maize and soybean rows was 60 cm.The plant density was 52 500 plants per ha for maize and 105 000 plants per ha for soybean.The plant spacing was 19 cm for all treatments,with a post-emergence density of one maize plant and two soybean plants per hole (Fig.1).

    The total number of plots was 15, and the area of each plot was 6 m×6 m.To avoid edge effects, border rows were planted along the edges of plots.Maize was sown on 6 April 2012 and 2 April 2013 and harvested on 3 August 2012 and 29 July 2013.Soybean was sown on 16 June 2012 and 14 June 2013, with simultaneous application of maize topdressing and soybean base fertilizer, and harvested on 27 October 2012 and 2 November 2013.

    2.2.Experimental design and treatments

    Five fertilizer treatments were investigated in a singlefactor randomized block design with three replications in the maize-soybean relay intercropping system.Three total nitrogen application rates included no nitrogen application(NN: 0 kg N ha-1), conventional nitrogen application rate employed by the local farmer (CN: 330 kg N ha-1), and lower nitrogen (LN: 270 kg N ha-1) with three topdressing distances(LND), e.g., 15 cm (LND1), 30 cm (LND2) and 45 cm (LND3)from maize rows.The NN and CN treatments were set as the control group (CK).Urea (N content: 46%), CaP2H4O8(P2O5content: 12%) and KCl (K2O content: 60%) were used as N, P and K fertilizer.All maize treatments received 105 kg P2O5ha-1and 112.5 kg K2O ha-1as base fertilizer,and all soybean treatments received 63 kg P2O5ha-1and 52.5 kg K2O ha-1as base fertilizer.

    All fertilizers for crops were closely placed to the plant stems in NN and CN treatments, respectively.In CN treatment, the total N fertilizer (330 kg N ha-1) was applied at 270 kg N ha-1for maize and 60 kg N ha-1for soybean.The N fertilizer for maize was applied at 120 kg N ha-1for base fertilizer and 150 kg N ha-1for topdressing, and N fertilizer for soybean was applied as base fertilizer.In optimized fertilization for LN, all base fertilizers for maize were closely placed to the plant’s stems.The total N fertilizer (270 kg N ha-1) was applied at 120 kg N ha-1for maize base fertilizer,and 150 kg N ha-1for maize topdressing.Topdressing for maize was integrated with soybean base fertilizer, and strips placed at distances of LND1, LND2 and LND3 from maize rows (Fig.2).

    2.3.Determination of crop shoot biomass, grain yield,and N content

    Four representative plants from each plot were separately sampled at five fully developed trifoliate leaf stage (V5),full bloom stage (R2), beginning seed stage (R5), and physiological maturity stage (R8) of soybeans (Fehr et al.1971), respectively.Crops were harvested from 12 m2areas of each plot to calculate grain yield and shoot biomass.Samples were divided into different organs and exposed to 105°C for 30 min to destroy tissues; they were oven-dried at 80°C to a constant weight before weighing.Nitrogen concentration in plant tissues was determined using BUCHI Auto Kjeldahl Unit K-370 (Buchi Labortechnik AG, Switzerland) according to the Kjeldahl method (Hou et al.2006).

    2.4.Determination of symbiotic nitrogen fixation by soybean

    Soybean samples were collected at the V5, R2 and R5 stages (Fehr et al.1971).Shoots were cut at ground level and divided into different organs.Root systems attached to the soil block (20 cm long×40 cm width×30 cm deep) were dug out and gently shaken to remove bulk soil.The entire root system was placed in nylon nets, soaked in water,and carefully washed and separated for counting nodules and lateral roots and for weighing nodules.The roots and nodules of four representative plants from each plot were stored for morphology and physiology determination.Root volume was measured with WinRHIZO System (WinRHIZO Pro LA2400; Regent Instruments Inc., Quebec, Canada);the main root length was manually determined with a ruler.Roots were then exposed to 105°C for 30 min and ovendried at 80°C to a constant weight before weighing.

    Table 1 Properties of the experimental soil in 2012 and 2013

    Fig.1 Maize-soybean relay intercropping in Renshou County,Sichuan Province, Southwest China, in August 2013.

    The acetylene (C2H2) reduction assay (ARA) was used to assess nitrogen fixation.The nodule samples were collected in 10-mL penicillin bottles in the presence of 20% (v/v)acetylene and then incubated for 30 min in a 30°C constant temperature water bath.Next, 0.1 mL of gas sample from each bottle was analyzed for ethylene content by a Shimadzu GC-2010 Plus Gas Chromatograph with Flame Ionization Detector (FID) (Siczek and Lipiec 2011).To measure the content of leghemoglobin, 2 g fresh soybean root nodules were ground into homogenate with approximately 5 volumes of phosphate buffer solution (5°C; pH 6.8;6.8 mol L-1).The samples were centrifuged at 100 r min-1for 15 min at 15°C, the sediment was discarded, and the supernatant was centrifuged at 39 000 r min-1for another 20 min at 5°C.The final supernatant was measured using a spectrophotometer (Mapada Instruments, Shanghai,China) at 540 nm (Riley and Dilworth 1985).Ureide content was analyzed in dry leaves, stems, pods, shells,and grains.The steps were modified according to Young and Conway (1942): 0.4 g tissue samples were added to 8 mL 0.1 mol L-1phosphate buffer (pH 7.0) containing 50%alcohol (v/v), and added to a water bath (80°C) for 5 min.The samples were centrifuged at 5 000 r min-1for 10 min,and clarified supernatant was collected and stored at-20°C.Then, 1.5 mL extracted samples mixed with 0.5 mL of 0.5 mol L-1NaOH were boiled for 6 min.After chilling to room temperature, 0.5 mL of 0.65 mol L-1HCl was added to the mixture and boiled for 6 min.After chilling, 0.5 mL of 0.4 mol L-1phosphate buffer (pH 7.0) and 0.5 mL of 0.5%C6H8N2·HCl were added and samples were kept for 6 min at room temperature.Samples were chilled to 0°C and then 2.5 mL concentrated HCl and 0.5 mL of 2% K3[Fe(CN)6] were added; samples were prechilled to 0°C and kept at room temperature for 15 min.Finally, the ureide concentration was measured by a spectrophotometer at 535 nm.

    Fig.2 Diagrammatic sketch of maize-soybean relay intercropping system and fertilizer application locations in 2012 and 2013.In the maize-soybean relay intercropping system, strip width is 2 m.Maize and soybean were sown in two separate rows with 40 cm width, and 60 cm between maize and soybean rows.Nitrogen application rates were 330 kg ha-1 (CN: 270 kg N ha-1 for maize and 60 kg N ha-1 for soybean); 270 kg N ha-1 (LN, shared by maize and soybean).Fertilization locations under LN were 15 cm(LND1), 30 cm (LND2), and 45 cm (LND3) from maize rows, CN and NN were applied between two plants in the same planting row.Post-emergence densities were 1 maize plant and 2 soybean plants per hole.Maize planting density was 52 500 plants ha-1,and soybean planting density was 105 000 plants ha-1.indicates fertilizer application locations.

    2.5.Calculation

    Plant organ N accumulation (NA, kg ha-1), total N accumulation (TNA, kg ha-1), NUE (%), NAE (kg kg-1) and nitrogen fixation (ARA, mL h-1plant-1) were calculated using the following formulas (Bellaloui et al.2014; Duan et al.2014).

    Total grain yield and shoot N uptake of maize-soybean relay intercropping system were the sum of component crops grain yield and shoot N uptake, respectively.The NUE and NAE of the maize-soybean relay intercropping system were determined by the total N uptake and total yield relative to the amount of N supplied.

    2.6.Statistical analysis

    The study was performed using one-factor design, namely,lower N application with three different topdressing positions.Conventional N application and no N application were set as the control.Therefore, data were analyzed with oneway ANOVA of SPSS (ver.15, SPSS, Chicago, IL, USA)in several steps: treatments were considered the fixed factor, and data were the dependent variable.Significant differences among treatments were determined by the least significant difference (LSD) (P≤0.05).

    3.Results

    3.1.Crop yield and nitrogen uptake

    Yield data of all five treatments showed that different N fertilization treatments had remarkable effects on crop grain yield (P<0.05) (Table 2).Compared with CN, LN increased the system grain yield by 5.7% in 2012 and by 3.3% in 2013.In addition, the system grain yield of LN was increased by 34.8% in 2012 and by 25.6% in 2013 compared with NN.Under different fertilization locations, crops performed better in LND1 and LND2 than that in LND3.The highest grain yield of maize and soybean was in LND1 and LND2,at 23.6 and 16.2% higher than LND3, respectively.In addition, the system grain yields of LND1 and LND2 were higher than LND3.Maize and soybean shoot N uptakes increased by 2.9 and 24.5% in LN compared with CN,respectively.Compared with LND3, shoot N uptake of maize and soybean in LND1 increased by 24.3 and 25.1%, and in LND2 increased by 14.9 and 13.6%, respectively (Table 3).

    3.2.Nitrogen utilization

    NAE and NUE were significantly affected by N fertilization treatments (Table 4).Compared with CN, the NAE and NUE of maize, soybean, and the relay intercropping system increased by 35.2, 98.3 and 51.6%, and 32.2, 160.4 and 72.5% in LN, respectively.The optimization of fertilization location increased the NAE and NUE of maize, soybean and the relay intercropping system by 242.6, 40.1 and 155.8%in LND1, and 77.0, 85.2 and 81.2% in LND2 compared with LND3, respectively.

    3.3.Soybean nodule nitrogen fixation

    The number and dry weight of soybean nodulesNodulenumber and dry weight gradually increased from the V5 to R5 stage (Fehr et al.1971) (Fig.3).Across the growth stages of soybeans, lower N fertilization promoted nodule growth.Compared with CN, the nodule number and dry weight in LND1 increased by 32.0 and 28.5% at the V5 stage, by 41.4 and 42.3% at the R2 stage, and by 47.5 and 35.8% at the R5 stage, respectively.In addition,compared with CN, the nodule number and dry weight in LND2 increased by 56.0 and 44.2% at the V5 stage, by 63.2 and 49.4% at the R2 stage, and by 40.4 and 29.1% at the R5 stage, respectively.Optimization of fertilization location promoted soybean nodule growth, and the nodule number and dry weight of LND1 and LND2 were higher compared with LND3, respectively.

    Table 2 Grain yield under different N application treatments (kg ha-1)

    Table 3 Shoot nitrogen uptake under different N application treatments (kg ha-1)

    Nitrogen fixation and leghemoglobin contentNitrogen fixation was assessed using ARA (Fig.3).With the advancement of growth stages, ARA first increased and subsequently decreased, peaking at the R2 stage.Under various N levels, the ARA of LN was higher than CN and NN,respectively.Fertilization location significantly affected N fixation ability.Compared with LND3, the ARA of LND1 and LND2 was respectively increased by 51.6 and 16.7% at the V5 stage, by 6.2 and 15.4% at the R2 stage, and by 24.0 and 32.9% at the R5 stage.Leghemoglobin content increased with the advancement of growth stages, peaking at the R5 stage (Table 5).Under LN treatments, the leghemoglobin content of LND1 and LND2 was significantly higher than LND3, respectively.Compared with CN, leghemoglobin content in LND1 and LND2 was significantly increased by 101.1 and 55.8% at the V5 stage, and by 73.6 and 91.7%at the R2 stage, and by 63.3 and 32.3% at the R5 stage,respectively.

    Ureide content and N accumulation of soybeanUreides are the dominant long-distance transport forms of N that are fixed in soybean nodules (Collier and Tegeder 2012).At different growth stages of soybean, ureide storage in soybean organs differed (Table 6).During vegetative stages, most ureides are stored in stems, while during reproductive stages, ureide was transferred from stems to pods and stored in grain.At all sampling stages, ureide content was higher in LN compared with CN and NN.Compared with LND3, ureide content of LND1 and LND2 increased by 23.1 and 22.6% in leaves at the V5 stage,by 22.2 and 11.6% in stems at the V5 stage, and by 17.1 and 23.9% in grain at the R8 stage, respectively.Across the growth stages of soybean, N accumulation of leaf was higher than in stems (Table 7).N fertilization influenced the N accumulation of soybean organs, and LN had greater N accumulation than CN and NN.Compared with LND3,N accumulation of LND1 and LND2 increased by 3.9 and 2.3% in leaves at the V5 stage, by 3.1 and 2.8% in stems at the V5 stage, and by 20.7 and 11.9% in pods at the R8 stage, respectively.

    Morphological characteristics of soybean rootsHigh levels of N application suppress soybean root system growth(Table 8).The root volume, main root length, and lateral root number were higher in LN compared with CN.At the V5,R2, and R5 stages of soybean, the root volume of CN was decreased by 43.73, 29.19, and 26.35% compared with LN,respectively.Additionally, at the V5, R2, and R5 stages of soybean, the main root length and lateral root number of LN were increased by 29.29, 9.77 and11.16%, and 24.63, 20.99 and 8.43%, respectively.Optimized fertilizer application methods promoted soybean root system growth.At V5,R2 and R5 stages, the root volume, main root length and lateral root number of LND1 and LND2 were greater than those of LND3, respectively.

    4.Discussion

    4.1.Optimizing N management to improve soybean biological nitrogen fixation

    n stem.3.9.9.9 Sy a a b c±4.4.2.6.3±5±5±2-plicatio 61 54 38 28 ap nal N 13 a b c 20 IS.6.4.1±2-ve±5±1 0.5 d ntio.0.9.2 33 26 19 9.6±N,con a a b b.8.4.8.5±0 ely;C IM±3 E (%)2)±4±2-ctiv 28.5.3.5.4 27 19 18 pe NU),res stem.2.5.3.9 D3a b c c Sy.6.2.8.7±1±3±3±0-62 47 30 28 m (LN 5 c 5.nts a b c.9.7.0),4 D2<0.0 atme 12 20 IS±2-.0.4.1±2±1 1.8 d t P n tre 28 21 14 8.6±m(xù) (LN 0 c D)a plicatio a b c b c IM±2.4.7.0.0±2±4±1-1), 3 ND nce (LS re ap 34.7.8.7.1 25 16 20(L.t d cm iffe ntN re-ding 1.7 a stem 0.6 b iffe 1.0 b System 1.1 a 15 sy nificanig erd 9.7±9.5±5.2±5.6±lu ping sts nd E)u b c n rate,inc bylea cy(NU 13 20 IS 0.9 a 2.9±0.4 a 3.1±0.2 b 2.0±0.2 c 1.2±-atio y intercrop la nces plic an iffere re ien IM 1.0 a ffic use e)2)1N ap erN 0.7 a 0.4 c 0.9 b--soybe t d nifican en 6.5±grain kg-6.9±3.2±4.5±erlow aize nitrog nd wsu m,m n sigea E (kg; syste n m nd a b c±1 E)a NA System.5.1.2.0.5±1±0 0.8 c 5.8±-ro aize n.s to m itro 13 10 4.00 ge an olum cy(NA o n oybe e cthin fficien 12 20 IS 0.2 b 0.2 a 0.3 c 0.2 c-sition po r;NN,n d spe rsntlette mye 2.3±2.8±1.7±1.3±sing re iffe no a l farme ro.3 dres m top lo±1 0.9 b 0.3 d 0.8 c ca; IS,intercrop.D ag IM-SD en.8 10 7.7±2.3±4.7±Nitrog bythe aize ed pe nce fro d m s±, intercrop meanre t1)ista ble 4 en, d ploy tm em ata a Ta Trea D1 LN D2 LN D3 LN CN NN1)L raND te2) IM e dTh

    Fig.3 Effects of different N application treatments on nodulation and biological nitrogen fixation at different growth stages.Nitrogen application rates were 0 (NN),270 kg N ha-1 (LN; shared by maize and soybean), and 330 kg ha-1 (CN; 270 kg N ha-1 for maize and 60 kg N ha-1 for soybean).Fertilization locations under LN were 15 cm (LND1), 30 cm (LND2), and 45 cm (LND3) from maize rows, CN and NN were applied between two plants in the same planting row.V5, five fully developed trifoliate leaf stage; R2, the full bloom stage; R5, the beginning seed stage.Nodules were separately sampled at V5, R2 and R5 stages of soybean.The data are means±SE.Different letters showed significant differences between different treatments under different growth stages by least significant difference (LSD) at P<0.05.

    The present study clearly demonstrated that the maizesoybean relay intercropping system with LN treatments presented an advantage over CN treatment.Moreover,optimization of fertilization methods further contributes to this advantage.The biological nitrogen fixation of soybean can meet almost half of N demand in the lifespan (Salvagiotti et al.2008).However, the grain yield of soybean with full-N is average 11% greater than zero-N, namely, N fertilizer is necessary to increase yield when indigenous N sources are insufficient to meet the growing demand(La Menza et al.2017).Intercropping of cereals and legumes are efficient cropping systems, e.g., promoting nitrogen fixation, enhancing biodiversity,increasing and stabilizing yields, improving resource use and reducing the damage caused by diseases and pests (Awal et al.2006; Hauggaard-Nielsen et al.2008; Gao et al.2010; Ghanbari et al.2010; Wang et al.2017; Yang et al.2017).Previous studies pointed out that compared with sole cereal,intercropped cereal with legume reduces N fertilizer input through biologicalnitrogen fixation, while planting patterns and fertilization treatments significantly affect legume nitrogen fixation (Li et al.2001; Ghosh et al.2006; Salvagiotti et al.2008; Hauggaard-Nielsen et al.2009; Wu et al.2014).

    Table 5 Leghemoglobin content of soybean under different nitrogen application treatments in 2013 (mg plant-1)1)

    In the maize-soybean relay intercropping system, soybean was sown at the wide rows when maize was at the 12th leaf stage.The topdressing for maize was integrated with soybean base fertilizer, and strips-placed at interspecific rows.The optimized fertilizer application methods increased NUE by sharing N fertilizer.Traditionally, fertilizers are hole-placed or strip-placed near plants stems.However, fertilizer placement near to plants may not be advantageous to nutrients efficient utilization in intercropping systems.In this study, LN treatments were advantageous to the soybean root system, nodule growth, and biological nitrogen fixation.Substantially, compared with intraspecific interactions, the interspecific interactions and facilitations are more important, namely, balancing interspecific and intraspecific interactions is the key fator to maximize the nutrients use efficiency of intercropping systems (Zhang and Long 2003).Compared with fertilizers placement at the plant stems, fertilizer placed at interspecific rows can promote the interspecific interactions and facilitation (Wu et al.2014).However, if the fertilizers are placed far away from plants, the availability of nutrients is low before the roots reaching the fertilizers (Nakamura et al.2008).

    After maize harvest, soybean performed better in LND1 and LND2 compared with LND3.It demonstrated the previous finding that starter N fertilizer input can influence soybean growth (Osborne and Riedell 2006).Additionally, the previous study also indicated that high N input has a negative effect on nodule nitrogen fixation (Gan et al.2002).High soil mineral N concentration is adverse to legume biological N fixation and thereby decreasing N transfer in legume-nonlegume mixtures (Thilakarathna et al.2016).Lower N input significantly increased soybean nitrogen fixation and promoted fixed N transfer (Stern 1993; Yong et al.2015).

    Fertilizer-placed strips were far from the soybean rows in LN treatments, which decreased soil N concentration and mitigated the N suppression of soybean nitrogen fixation compared with CN.The variation of biological nitrogen fixation could be explained by reduced N competition in the intercropping system (Li et al.2011).Moreover, compared with NN, the increased soil N concentration in LND1 and LND2 promoted soybean root growth compared with NN, thereby enhancing the potential for N acquisition.Maize nitrogen uptake led to decreased available N in the soybean rhizosphere, which forced soybeans to fix more nitrogen from the air (Li et al.2006).It also meets the soybean starter nitrogen needs during early growth stages and effectively promoted soybean root growth (Osborne and Riedell 2006).Therefore, crops performed better when fertilizers were placed at interspecific rows in the maize-soybean relay intercropping system.Namely, lower N input with 15-30 cm between the fertilization locations and the maize rows promoted competitive nitrogen uptake and met crop N demand in the growing stage.

    n Grain.2.2.3.1.9c d a b b c.6.6.1.1.1±0±0±2 plicatio±1±0 45 43 42 41 35 ap na ntio ell.4.4.2.2.3a ab ve±0 b b c±0.4.2.8.8.1 Sh±0±0 27 26±0 R8 27 26 l N 23 N,con Stem.7.2.1.9.6bc a ab ab c ely;C ctiv±0.4±3.8.1.6.9±2±1±3 pe 32 37 36 34 28),res 0.1 b D3 Leaf 2.4±0.3 a 3.1±2.6±b b b 0.6 a 2.8±0.1 a 0.3 a 2.5±m(xù) (LN),45 c 5..0<0 b c d..0.5.5.8.5a b t P.9.5.7.0.4 D)a d±1±0 55±1 olg-1)1)Po±1 49 45 ha±0 ge 51 D2 51 rvest s 0 c1), 3 nce (LS 8,the(L 3 (μm re Stem.3.4.9.9.5a ab iffe b c c±0 01 R5 t d.5.2.2.0 ta m (LN±0±1±1 cm.4±0 tsin 2 57 57 50 49 59 ge; R 15 ding nifican ig sts en af 4.5±Le ND tm trea 0.3 a 0.4 a 4.6±0.6 a 0.4 a 4.5±4.4±0.3 b lu g s 3.3±in bylea ginnee tion n rate,inc be atio nces plica re d sta ap d a.5.3.4.4.2b c e a p iffe t d en±1 Stem±0 26 e;R±0±0 32.5.5.8.4.6 itrog 28 35 30 5,the plic erN stag nifican R2 ntn erlow ig n s om nd re b lo 0.1 a af iffe 0.6 a ea 8.1±d c 0.4 c Le 8.6±0.5 b 6.9±7.0±1.0 d 5.8±llb erd fu ro n m aize n.nd nsu ge e colum±0 orga in o n th.6.1.2.9.1 gea b c c d; R2,the±0 Stem±0 s to m itro±0±0±0.0.8.3.0.8 le po r;N an N,n rs 26 oybe 21 22 wsu 19 23 afsta te V5 ntlette re af.5.4f s l farme t o a a iffe±0 0.1 b nten Le±0 0.5 b trifolia singsition dres 0.3 c ca.D.3.3 co 10 10 8.4±7.1±8.8±ed m top lo s±Ureide evelop bythe SD nce fro ed an ploy mere ble 6 t2)en1) V5 ND tetm2)L ra e d D1 ista em ata a Ta, d Trea LN D2 D3, five fullyd LN LN CN NN Th

    4.2.Saving fertilizer and improving yield by optimizing N fertilization rates and locations

    ?

    ?

    Abuse and low efficiency of chemical N fertilizer pollute the environment.Recently,new theories and technologies for N reduction that maintain crop yield and improve NUE have been reported (Caliskan et al.2008; Constantin et al.2010; Ruan et al.2013; Gaudin et al.2015; Li et al.2017).There are two ways to reduce chemical N input: breeding N-efficient cultivars and optimizing N management strategies.N management strategies are based in part on the use of efficient cultivars.Therefore,both approaches are essential to developing sustainable agriculture.However, breeding efficient cultivars may take many years,while optimizing nutrient management can be implemented quickly.

    It is well documented that excessive N input results in a reduction of maize yield(Abbasi et al.2012), yield loss may result from soil acidification.The previous study confirmed that the soil pH of maize rows significantly decreased in the no partition condition compared with a partition in a maize-soybean relay intercropping system(Wang et al.2017).In addition, soil pH significantly decreased with the increase in N applied rate (Zhao et al.2014).In this study,lower N application rate increased both crop yield and NUE, which supported previous findings by Ju et al.(2015) and Caliskan et al.(2008).Under CN treatment, maize yield loss may result from soil acidification,and soybean vigorous growth results in decreased grain yield.Under LN treatment,N topdressing for maize, and soybean N,P and K fertilizers were mixed and shared with each other.Although the N input of LN treatments decreased compared with CN,the base fertilizers of soybean that shared by maize could be another reason for the maize yield increase in LN treatments.This confirmed the previous finding that intercropped crops perform better when fertilizer is supplied at interspecific rather than intraspecific rows (Wu et al.2014).

    Zhang et al.(2003) pointed out that adjusting the relative roles of interspecific and intraspecific interaction is the key factor to improved NUE and yield in intercropping systems.A recent study shows that fertilizers placed at interspecific rows can promote interspecific interaction and facilitation,and can increase yield compared with intraspecific placement in an intercropping system (Wu et al.2014).In the present study, under LN treatments, crops performed better in LND1 and LND2 than in LND3.Optimized fertilization locations further promoted root growth, improving the potential to utilize soil resources and thereby increasing resource use efficiency and crop yield.The increased distance between the fertilization locations and maize rows resulted in decreased maize yield, crop NUE, and NAE.However, soybean yield presented a unimodal curve and peaked at LND2.This suggested that the optimized fertilization location could be 15 or 30 cm, or between 15 and 30 cm, because optimized fertilization locations can alleviate interspecific nutrient competition, and can promote interspecific interaction and facilitation.Therefore, the optimized methods can increase fertilizer use efficiency and grain yield.It is well documented that the N competition ability of maize is stronger than that of soybean in maizesoybean relay intercropping system (Wang et al.2017).Interspecific competitive (Li et al.2001; Zhang et al.2011)and complementary N use (Chalk 1997; Xiao et al.2004)in cereal-legume intercropping system, decrease soybean“nitrogen suppression” effects and promote N fixation.In addition, maize root exudates increased faba bean biological nitrogen fixation by promoting nodule formation and N2fixation (Li et al.2016).Barley promoted nodule formation of pea and decreased interspecific N competition(Hauggaard-Nielsen et al.2009).In the present study, N accumulation of soybean organs in LND1 and LND2 was higher than that in other treatments, respectively, suggesting a more favorable balance between fixed atmospheric N and absorbed fertilizer N.

    5.Conclusion

    The results showed that the NAE, NUE, and yield of maize peaked at LND1, which were 242.60, 76.97, and 23.60%respectively higher than LND3.Similarly, the NUE of soybean peaked at LND1, was 85.23% higher than that of LND3.The NAE and yield of soybean peaked at LND2,were increased by 59.85 and 16.20% compared with LND3.Moreover, the nitrogen fixation ability of soybean was greater under LND1 and LND2.On the whole, crop performed better under LND1 and LND2 treatments.The NUE and total grain yield of the maize-soybean relay intercropping system were significantly higher in LND1 and LND2 than that in LND3 and CN, respectively.Those findings suggested that lower N application at 15-30 cm from the fertilizer application location to the maize row was optimal.This study contributes to improved nutrient management for legume-cereal relay intercropping system and helps to develop an efficient, sound agricultural system which is environmentally friendly.Nevertheless, this N application method requires assay using more legume-cereal relay intercropping systems and fertilizer types and rates.

    Acknowledgements

    The research was supported by the National Key Research and Development Program of China (2016YFD0300202),and the National Natural Science Foundation of China(31671625, 31271669).Thanks also go to two anonymous reviewers for their constructive comments.

    Abbasi M K, Tahir M M, Sadiq A, Iqbal M, Zafar M.2012.Yield and nitrogen use efficiency of rainfed maize response to splitting and nitrogen rates in Kashmir, Pakistan.Agronomy Journal, 104, 448-457.

    Awal M A, Koshi H, Ikeda T.2006.Radiation interception and use by maize/peanut intercrop canopy.Agricultural & Forest Meteorology, 139, 74-83.

    Caliskan S, Ozkaya I, Caliskan M E, Arslan M.2008.The effects of nitrogen and iron fertilization on growth, yield and fertilizer use efficiency of soybean in a mediterranean-type soil.Field Crops Research, 108, 126-132.

    Chalk P M.1997.Dynamics of biologically fixed N in legumecereal rotations: A review.Australian Journal of Agricultural Research, 49, 303-316.

    Chen X, Zhang F, R?mheld V, Horlacher D, Schulz R, B?ning-Zilkens M, Wang P, Claupein W.2006.Synchronizing N supply from soil and fertilizer and N demand of winter wheat by an improved Nminmethod.Nutrient Cycling in Agroecosystems, 74, 91-98.

    Collier R, Tegeder M.2012.Soybean ureide transporters play a critical role in nodule development, function and nitrogen export.The Plant Journal, 72, 355-367.

    Constantin J, Mary B, Laurent F, Aubrion G, Fontaine A,Kerveillant P, Beaudoin N.2010.Effects of catch crops, no till and reduced nitrogen fertilization on nitrogen leaching and balance in three long-term experiments.Agriculture Ecosystems & Environment, 135, 268-278.

    Corre-Hellou G, Dibet A, Hauggaard-Nielsen H, Crozat Y,Gooding M, Ambus P, Dahlmann C, von Fragstein P, Pristeri A, Monti M.2011.The competitive ability of pea-barley intercrops against weeds and the interactions with crop productivity and soil N availability.Field Crops Research,122, 264-272.

    Dong Q, Yong T W, Liu X M, Liu W Y, Xu T, Song C, Wang X C,Yang W Y.2014.Effect of nitrogen application methods on crop yield and grain filling characteristics of maize in maizesoybean relay strip intercropping system.Acta Agronomica Sinica, 40, 2028-2039.(in Chinese)

    Drinkwater L E, Wagoner P, Sarrantonio M.1998.Legumebased cropping systems have reduced carbon and nitrogen losses.Nature, 396, 262-265.

    Fehr W R, Caviness C E, Burmood D T, Pennington J S.1971.Stage of development descriptions for soybeans, Glycine Max (L.) Merrill.Crop Science, 11, 929-931.

    Fernández-Escobar R, Benlloch M, Herrera E, García-Novelo J M.2004.Effect of traditional and slow-release N fertilizers on growth of olive nursery plants and N losses by leaching.Scientia Horticulturae, 101, 39-49.

    Gan Y, Stulen I, Posthumus F, Keulen H V, Kuiper P.2002.Effects of N management on growth, N2fixation and yield of soybean.Nutrient Cycling in Agroecosystems, 62, 163-174.

    Gao Y, Duan A, Qiu X, Sun J, Zhang J, Liu H, Wang H.2010.Distribution and use efficiency of photosynthetically active radiation in strip intercropping of maize and soybean.Agronomy Journal, 102, 1149-1157.

    Gaudin A C M, Janovicek K, Deen B, Hooker D C.2015.Wheat improves nitrogen use efficiency of maize and soybeanbased cropping systems.Agriculture, Ecosystems &Environment, 210, 1-10.

    Ghanbari A, Dahmardeh M, Siahsar B A, Ramroudi M.2010.Effect of maize (Zea mays L.)-cowpea (Vigna unguiculata L.)intercropping on light distribution, soil temperature and soil moisture in arid environment.Journal of Food Agriculture& Environment, 8, 102-108.

    Ghosh P K, Manna M C, Bandyopadhyay K K, Ajay, Tripathi A K, Wanjari R H, Hati K M, Misra A K, Acharya C L, Rao A S.2006.Interspecific interaction and nutrient use in soybean/sorghum intercropping system.Agronomy Journal,98, 1097-1108.

    Godfray H C J, Beddington J, Crute I, Haddad L, Lawrence D,Muir J F, Pretty J, Robinson S, Thomas S, Toulmin C.2010.Food security: The challenge of feeding 9 billion people.Science, 327, 812-818.

    Guo J H, Liu X J, Zhang Y, Shen J L, Han W X, Zhang W F,Christie P, Goulding K W T, Vitousek P M, Zhang F S.2010.Significant acidification in major Chinese croplands.Science, 327, 1008-1010.

    Hauggaard-Nielsen H, Gooding M, Ambus P, Correhellou G,Crozat Y, Dahlmann C, Dibet A, Von F P, Pristeri A, Monti M.2009.Pea-barley intercropping for efficient symbiotic N2-fixation, soil N acquisition and use of other nutrients in European organic cropping systems.Field Crops Research,113, 64-71.

    Hauggaard-Nielsen H, Jornsgaard B, Kinane J, Jensen E S.2008.Grain legume-cereal intercropping: The practical application of diversity, competition and facilitation in arable and organic cropping systems.Renewable Agriculture and Food Systems, 23, 3-12.

    Hochman Z, Carberry P S, Robertson M J, Gaydon D S,Bell L W, Mcintosh P C.2011.Prospects for ecological intensification of Australian agriculture.European Journal of Agronomy, 44, 109-123.

    Hou Z, Li P, Li B, Gong J, Wang Y.2006.Effects of fertigation scheme on N uptake and N use efficiency in cotton.Plant and Soil, 290, 115-126.

    Ju C, Buresh R J, Wang Z, Zhang H, Liu L, Yang J, Zhang J.2015.Root and shoot traits for rice varieties with higher grain yield and higher nitrogen use efficiency at lower nitrogen rates application.Field Crops Research, 175, 47-55.

    Ju X T, Xing G X, Chen X P, Zhang S L, Zhang L J, Liu X J, Cui Z L, Yin B, Christie P, Zhu Z L.2009.Reducing environmental risk by improving N management in intensive Chinese agricultural systems.Proceedings of the National Academy of Sciences of the United States of America,106, 3041-3046.

    Kwong K F N K, Bholah A, Volcy L, Pynee K.2002.Nitrogen and phosphorus transport by surface runoff from a silty clay loam soil under sugarcane in the humid tropical environment of Mauritius.Agriculture Ecosystems & Environment, 91,147-157.

    Larson C.2013.Losing arable land, china faces stark choice:Adapt or go hungry.Science, 339, 644.

    Li B, Li Y, Wu H, Zhang F, Li C, Li X, Lambers H, Li L.2016.Root exudates drive interspecific facilitation by enhancing nodulation and N2fixation.Proceedings of the National Academy of Sciences of the United States of America,113, 6496-6501.

    Li H, Hu B, Chu C.2017.Nitrogen use efficiency in crops:Lessons from Arabidopsis and rice.Journal of Experimental Botany, 68, 2477-2488.

    Li L, Sun J, Zhang F, Guo T, Bao X, Smith F A, Smith S E.2006.Root distribution and interactions between intercropped species.Oecologia, 147, 280-290.

    Li L, Sun J H, Zhang F S, Li X L, Rengel Z, Yang S C.2001.Wheat/maize or wheat/soybean strip intercropping.II.Recovery or compensation of maize and soybean after wheat harvesting.Field Crops Research, 71, 123-137.

    Li L, Zhang F, Li X, Christie P, Sun J, Yang S, Tang C.2003.Interspecific facilitation of nutrient uptake by intercropped maize and faba bean.Nutrient Cycling in Agroecosystems,65, 61-71.

    Li Q Z, Sun J H, Wei X J, Christie P, Zhang F S.2011.Overyielding and interspecific interactions mediated by nitrogen fertilization in strip intercropping of maize with faba bean, wheat and barley.Plant and Soil, 339, 147-161.

    Liu C, Wang K, Meng S, Zheng X, Zhou Z, Han S, Chen D, Yang Z.2011.Effects of irrigation, fertilization and crop straw management on nitrous oxide and nitric oxide emissions from a wheat-maize rotation field in northern China.Agriculture, Ecosystems & Environment, 140, 226-233.

    Liu X, Rahman T, Yang F, Song C, Yong T, Liu J, Zhang C,Yang W.2017.PAR interception and utilization in different maize and soybean intercropping patterns.PLoS ONE,12, e0169218.

    Meng L, Ding W, Cai Z.2005.Long-term application of organic manure and nitrogen fertilizer on N2O emissions, soil quality and crop production in a sandy loam soil.Soil Biology &Biochemistry, 37, 2037-2045.

    La Menza N C, Monzon J P, Specht J E, Grassini P.2017.Is soybean yield limited by nitrogen supply? Field Crops Research, 213, 204-212.

    Nakamura R, Kachi N, Suzuki J I.2008.Root growth and plant biomass in Lolium perenne exploring a nutrient-rich patch in soil.Journal of Plant Research, 121, 547-557.

    Noellsch A, Motavalli P, Nelson K A, Kitchen N R.2009.Corn response to conventional and slow-release nitrogen fertilizers across a claypan landscape.Agronomy Journal,101, 607-614.

    Oljaca S, Cvetkovic R, Kovacevic D, Vasic G, Momirovic N.2000.Effect of plant arrangement pattern and irrigation on efficiency of maize (Zea mays) and bean (Phaseolus vulgaris) intercropping system.The Journal of Agricultural Science, 135, 261-270.

    Osborne S, Riedell W.2006.Starter nitrogen fertilizer impact on soybean yield and quality in the Northern Great Plains.Agronomy Journal, 98, 1569-1574.

    Riley I T, Dilworth M J.1985.Cobalt requirement for nodule development and function in Lupinus angustifolius L.New Phytologist, 100, 347-359.

    Rodríguez-Navarro D N, Oliver I M, Contreras M A, Ruiz-Sainz J E.2011.Soybean interactions with soil microbes,agronomical and molecular aspects.Agronomy for Sustainable Development, 31, 173-190.

    Ruan W B, Ren T, Chen Q, Zhu X, Wang J G.2013.Effects of conventional and reduced N inputs on nematode communities and plant yield under intensive vegetable production.Applied Soil Ecology, 66, 48-55.

    Salvagiotti F, Cassman K G, Specht J E, Walters D T, Weiss A, Dobermann A.2008.Nitrogen uptake, fixation and response to fertilizer N in soybeans: A review.Field Crops Research, 108, 1-13.

    Shoji S, Delgado J, Mosier A, Miura Y.2001.Use of controlled release fertilizers and nitrification inhibitors to increase nitrogen use efficiency and to conserve air and water quality.Communications in Soil Science and Plant Analysis, 32,1051-1070.

    Siczek A, Lipiec J.2011.Soybean nodulation and nitrogen fixation in response to soil compaction and surface straw mulching.Soil and Tillage Research, 114, 50-56.

    Stern W R.1993.Nitrogen fixation and transfer in intercrop systems.Field Crops Research, 34, 335-356.

    Thilakarathna M S, Mcelroy M S, Chapagain T, Papadopoulos Y A, Raizada M N.2016.Belowground nitrogen transfer from legumes to non-legumes under managed herbaceous cropping systems.A review.Agronomy for Sustainable Development, 36, 58.

    Tilman D, Balzer C, Hill J, Befort B L.2011.Global food demand and the sustainable intensification of agriculture.Proceedings of the National Academy of Sciences of the United States of America, 108, 20260-20264.

    Wang X C, Deng X Y, Pu T, Song C, Yong T W, Yang F, Sun X, Liu W G, Yan Y H, Du J, Liu J, Su K, Yang W Y.2017.Contribution of interspecific interactions and phosphorus application to increasing soil phosphorus availability in relay intercropping systems.Field Crops Research, 204, 12-22.

    Wu B, Fullen M A, Li J, An T, Fan Z, Zhou F, Zi S, Yang Y,Xue G, Liu Z.2014.Integrated response of intercropped maize and potatoes to heterogeneous nutrients and crop neighbours.Plant and Soil, 374, 185-196.

    Xiang D B, Yong T W, Yang W Y, Wan Y, Gong W Z, Cui L,Lei T.2012.Effect of phosphorus and potassium nutrition on growth and yield of soybean in relay strip intercropping system.Scientific Research and Essays, 7, 342-351.

    Xiao Y, Li L, Zhang F.2004.Effect of root contact on interspecific competition and N transfer between wheat and fababean using direct and indirect15N techniques.Plant and Soil,262, 45-54.

    Xing G X, Zhu Z L.2002.Regional nitrogen budgets for China and its major watersheds.Biogeochemistry, 57–58,405-427.

    Yan G, Zheng X, Cui F, Yao Z, Zhou Z, Deng J, Xu Y.2013.Two-year simultaneous records of N2O and NO fluxes from a farmed cropland in the northern China plain with a reduced nitrogen addition rate by one-third.Agriculture, Ecosystems& Environment, 178, 39-50.

    Yan Y H, Gong W Z, Yang W Y, Yan W, Chen X L, Chen Z Q,Wang L Y.2010.Seed treatment with uniconazole powder improves soybean seedling growth under shading by corn in relay strip intercropping system.Plant Production Science 72, 367-374.

    Yang F, Huang S, Gao R, Liu W, Yong T, Wang X, Wu X,Yang W.2014.Growth of soybean seedlings in relay strip intercropping systems in relation to light quantity and red:far-red ratio.Field Crops Research, 155, 245-253.

    Yang F, Liao D, Wu X, Gao R, Fan Y, Raza M A, Wang X,Yong T, Liu W, Liu J, Du J, Shu K, Yang W.2017.Effect of aboveground and belowground interactions on the intercrop yields in maize-soybean relay intercropping systems.Field Crops Research, 203, 16-23.

    Yong T W, Liu X M, Feng Y, Song C, Wang X C, Liu W G,Su B Y, Li Z, Yang W Y.2015.Characteristics of nitrogen uptake, use and transfer in a wheat-maize-soybean relay intercropping system.Plant Production Science, 18,388-397.

    Young E G, Conway C F.1942.On the estimation of allantion by the rimini-schryver reaction.Journal of Biological Chemistry,142, 839-853.

    Yu P, Li X, White P J, Li C.2015.A large and deep root system underlies high nitrogen-use efficiency in maize production.PLoS ONE, 10, e0126293.

    Zhang F S, Long L.2003.Using competitive and facilitative interactions in intercropping systems enhances crop productivity and nutrient-use efficiency.Plant and Soil,248, 305-312.

    Zhang G, Yang Z, Dong S.2011.Interspecific competitiveness affects the total biomass yield in an alfalfa and corn intercropping system.Field Crops Research, 124, 66-73.

    Zhang W, Dou Z, He P, Ju X T, Powlson D, Chadwick D, Norse D, Lu Y L, Zhang Y, Wu L.2013.New technologies reduce greenhouse gas emissions from nitrogenous fertilizer in China.Proceedings of the National Academy of Sciences of the United States of America, 110, 8375-8380.

    Zhao S, Qiu S, Cao C, Zheng C, Zhou W, He P.2014.Responses of soil properties, microbial community and crop yields to various rates of nitrogen fertilization in a wheatmaize cropping system in north-central China.Agriculture,Ecosystems & Environment, 194, 29-37.

    Zhao X, Zhou Y, Min J, Wang S, Shi W, Xing G.2012.Nitrogen runoff dominates water nitrogen pollution from rice-wheat rotation in the Taihu Lake region of China.Agriculture Ecosystems & Environment, 156, 1-11.

    Zhou L.2017.Grain production and food security in China.In: Zhou L, ed., Reform and Development of Agriculture in China.Springer, Singapore.pp.141-169.

    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.

    Zuo Y, Zhang Z, Liu C, Zhang W.2015.Achieving food security and high production of bioenergy crops through intercropping with efficient resource use in China.Frontiers of Agricultural Science & Engineering, 2, 134-143.

    天堂动漫精品| 男人和女人高潮做爰伦理| 精品一区二区三区视频在线 | 麻豆成人午夜福利视频| 真实男女啪啪啪动态图| 又紧又爽又黄一区二区| 亚洲欧美日韩无卡精品| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| cao死你这个sao货| 国产成人啪精品午夜网站| 日韩国内少妇激情av| 最好的美女福利视频网| 偷拍熟女少妇极品色| 热99re8久久精品国产| 国产在线精品亚洲第一网站| 麻豆久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 99国产精品一区二区蜜桃av| 久久伊人香网站| 婷婷亚洲欧美| 男女做爰动态图高潮gif福利片| 免费观看人在逋| aaaaa片日本免费| 看免费av毛片| 黄色 视频免费看| 丰满人妻一区二区三区视频av | 久久久国产成人免费| 最近最新中文字幕大全电影3| 成年版毛片免费区| 宅男免费午夜| 精品一区二区三区视频在线 | 99久久久亚洲精品蜜臀av| 在线免费观看的www视频| 欧美xxxx黑人xx丫x性爽| 成人三级做爰电影| 黄色日韩在线| 中文字幕人成人乱码亚洲影| a级毛片在线看网站| 美女高潮的动态| 91在线观看av| 美女午夜性视频免费| 51午夜福利影视在线观看| 两人在一起打扑克的视频| 亚洲av片天天在线观看| 亚洲美女黄片视频| 亚洲av第一区精品v没综合| 成人鲁丝片一二三区免费| 国产熟女xx| 国产av在哪里看| 国产精品爽爽va在线观看网站| 亚洲人成电影免费在线| 久久久久免费精品人妻一区二区| 男女午夜视频在线观看| 成在线人永久免费视频| 夜夜夜夜夜久久久久| 亚洲无线观看免费| 男人舔女人下体高潮全视频| 91久久精品国产一区二区成人 | a在线观看视频网站| 真人一进一出gif抽搐免费| 国产高清三级在线| 日韩人妻高清精品专区| 欧美性猛交╳xxx乱大交人| 中文字幕精品亚洲无线码一区| 一区二区三区国产精品乱码| 久久人人精品亚洲av| 免费大片18禁| 午夜成年电影在线免费观看| 中文亚洲av片在线观看爽| 美女高潮喷水抽搐中文字幕| 操出白浆在线播放| 国产精品亚洲美女久久久| 国产野战对白在线观看| 国产成人精品久久二区二区91| 99久国产av精品| 亚洲欧美日韩东京热| 美女午夜性视频免费| 啪啪无遮挡十八禁网站| 国产精品一区二区免费欧美| 一区二区三区高清视频在线| or卡值多少钱| 精品午夜福利视频在线观看一区| 国产精品久久久久久久电影 | 久久久色成人| 黑人操中国人逼视频| 午夜影院日韩av| 91av网站免费观看| 欧美国产日韩亚洲一区| 999精品在线视频| www.999成人在线观看| 别揉我奶头~嗯~啊~动态视频| 午夜精品一区二区三区免费看| 两个人的视频大全免费| 听说在线观看完整版免费高清| x7x7x7水蜜桃| 亚洲午夜精品一区,二区,三区| 久久这里只有精品中国| 两性夫妻黄色片| 久久久久久久久免费视频了| h日本视频在线播放| 美女大奶头视频| 黄色视频,在线免费观看| 国产成人av教育| 丁香六月欧美| 色哟哟哟哟哟哟| 1024香蕉在线观看| 免费搜索国产男女视频| 一本一本综合久久| 两个人看的免费小视频| 国产精品1区2区在线观看.| 久久欧美精品欧美久久欧美| 视频区欧美日本亚洲| 国产精品美女特级片免费视频播放器 | 男插女下体视频免费在线播放| 蜜桃久久精品国产亚洲av| 搡老熟女国产l中国老女人| 久久精品人妻少妇| 亚洲精品美女久久av网站| 搡老妇女老女人老熟妇| 好男人在线观看高清免费视频| 国产精品永久免费网站| 中文字幕av在线有码专区| 亚洲一区二区三区色噜噜| 成人国产一区最新在线观看| 国产精品久久久人人做人人爽| 哪里可以看免费的av片| 俺也久久电影网| 亚洲一区高清亚洲精品| 日韩大尺度精品在线看网址| 亚洲欧美精品综合一区二区三区| 日本免费一区二区三区高清不卡| 99热6这里只有精品| 日韩欧美在线二视频| 欧美日韩国产亚洲二区| 我的老师免费观看完整版| 特大巨黑吊av在线直播| 国产真实乱freesex| 精品久久久久久成人av| 熟妇人妻久久中文字幕3abv| 欧美乱码精品一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 少妇丰满av| 久久久水蜜桃国产精品网| 男女之事视频高清在线观看| 在线观看午夜福利视频| x7x7x7水蜜桃| 精品国内亚洲2022精品成人| 国产视频一区二区在线看| 成人av一区二区三区在线看| 午夜福利欧美成人| 精品免费久久久久久久清纯| 国产视频内射| 国产综合懂色| 全区人妻精品视频| 午夜免费观看网址| а√天堂www在线а√下载| 中出人妻视频一区二区| 久久久精品大字幕| 变态另类成人亚洲欧美熟女| 午夜免费观看网址| cao死你这个sao货| 成人鲁丝片一二三区免费| 精品久久久久久久久久久久久| 看片在线看免费视频| 全区人妻精品视频| 国产一区在线观看成人免费| 午夜福利免费观看在线| 欧美乱妇无乱码| 成年女人毛片免费观看观看9| 精品久久久久久久久久免费视频| 精品福利观看| 国产三级中文精品| 亚洲精品美女久久久久99蜜臀| a级毛片在线看网站| 午夜成年电影在线免费观看| 女同久久另类99精品国产91| 国产精品香港三级国产av潘金莲| 69av精品久久久久久| 国产成人av教育| 伦理电影免费视频| 成人永久免费在线观看视频| 亚洲国产高清在线一区二区三| 美女午夜性视频免费| 淫秽高清视频在线观看| 精品一区二区三区四区五区乱码| 成人一区二区视频在线观看| 一二三四在线观看免费中文在| av天堂在线播放| 深夜精品福利| 1024香蕉在线观看| 老鸭窝网址在线观看| 桃色一区二区三区在线观看| 国产综合懂色| 国产免费男女视频| 国产三级在线视频| 综合色av麻豆| 日本撒尿小便嘘嘘汇集6| 脱女人内裤的视频| 日本黄色视频三级网站网址| 亚洲精品美女久久av网站| 国产亚洲av嫩草精品影院| 麻豆一二三区av精品| 亚洲自拍偷在线| 天天躁日日操中文字幕| 看免费av毛片| 久久久国产欧美日韩av| 欧美成人一区二区免费高清观看 | 九九在线视频观看精品| 久久人人精品亚洲av| 九九热线精品视视频播放| 这个男人来自地球电影免费观看| 婷婷丁香在线五月| 免费看光身美女| 欧美又色又爽又黄视频| 成人特级av手机在线观看| 久久精品91无色码中文字幕| 中文字幕最新亚洲高清| 欧美激情久久久久久爽电影| 久久人人精品亚洲av| 最近视频中文字幕2019在线8| 午夜激情欧美在线| 少妇的逼水好多| 999久久久国产精品视频| 亚洲熟妇熟女久久| 成人av在线播放网站| 国产av麻豆久久久久久久| 人妻夜夜爽99麻豆av| 国产高清视频在线观看网站| 欧美日韩一级在线毛片| 国产伦精品一区二区三区四那| 日本与韩国留学比较| 亚洲国产精品成人综合色| 欧美中文日本在线观看视频| 亚洲中文字幕一区二区三区有码在线看 | 99热这里只有精品一区 | 国产成人影院久久av| 国产伦在线观看视频一区| 国产精品免费一区二区三区在线| 99久国产av精品| 日本免费一区二区三区高清不卡| 午夜福利视频1000在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久精品欧美日韩精品| 国产又黄又爽又无遮挡在线| avwww免费| 亚洲成人久久爱视频| 国产精品99久久99久久久不卡| 一个人免费在线观看的高清视频| 欧美大码av| 美女大奶头视频| 精品熟女少妇八av免费久了| 亚洲国产高清在线一区二区三| 精品免费久久久久久久清纯| 免费在线观看日本一区| 日本五十路高清| 欧美激情在线99| 亚洲第一欧美日韩一区二区三区| 一级毛片精品| 亚洲乱码一区二区免费版| 欧美成人性av电影在线观看| 成年女人毛片免费观看观看9| 搞女人的毛片| 国产精品久久久久久亚洲av鲁大| 欧美日韩中文字幕国产精品一区二区三区| 在线看三级毛片| 搞女人的毛片| 亚洲成人中文字幕在线播放| av视频在线观看入口| 亚洲avbb在线观看| 在线十欧美十亚洲十日本专区| 久久久成人免费电影| www.www免费av| 亚洲欧美日韩高清在线视频| 亚洲电影在线观看av| 99在线人妻在线中文字幕| 性色avwww在线观看| 夜夜躁狠狠躁天天躁| 男人舔女人的私密视频| 天天添夜夜摸| aaaaa片日本免费| 不卡av一区二区三区| 亚洲精品在线观看二区| 观看美女的网站| 国产综合懂色| 日韩免费av在线播放| 制服丝袜大香蕉在线| 久久久久精品国产欧美久久久| 亚洲乱码一区二区免费版| 国产淫片久久久久久久久 | www.精华液| 村上凉子中文字幕在线| 黄色片一级片一级黄色片| 日韩中文字幕欧美一区二区| 性色avwww在线观看| 色老头精品视频在线观看| 最近视频中文字幕2019在线8| 免费在线观看日本一区| 国产av不卡久久| 国产精品永久免费网站| 久久99热这里只有精品18| 亚洲一区高清亚洲精品| 色在线成人网| 欧美3d第一页| 精品久久久久久,| 狂野欧美白嫩少妇大欣赏| 亚洲一区二区三区色噜噜| 一边摸一边抽搐一进一小说| 激情在线观看视频在线高清| 国产成人精品久久二区二区91| 成人永久免费在线观看视频| 我要搜黄色片| 无限看片的www在线观看| 国产精品日韩av在线免费观看| 99热这里只有是精品50| 中文字幕人妻丝袜一区二区| 在线播放国产精品三级| 欧美午夜高清在线| 中文在线观看免费www的网站| avwww免费| 欧美成人免费av一区二区三区| 亚洲美女黄片视频| 又爽又黄无遮挡网站| 亚洲av成人不卡在线观看播放网| 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 欧美乱妇无乱码| 亚洲成人免费电影在线观看| 99riav亚洲国产免费| 曰老女人黄片| 亚洲七黄色美女视频| 亚洲专区字幕在线| 精品国产乱码久久久久久男人| 麻豆一二三区av精品| 久久久水蜜桃国产精品网| 岛国在线观看网站| 国产精品精品国产色婷婷| 一区二区三区国产精品乱码| 女生性感内裤真人,穿戴方法视频| 级片在线观看| 久久久久久久久久黄片| 日本一本二区三区精品| 天堂动漫精品| 成人欧美大片| 亚洲精品美女久久久久99蜜臀| 国产高清三级在线| 国产精品一区二区三区四区久久| 精品一区二区三区av网在线观看| 日本精品一区二区三区蜜桃| 脱女人内裤的视频| 欧美一区二区国产精品久久精品| 亚洲av成人一区二区三| 国产精品女同一区二区软件 | 久久久久久久精品吃奶| 亚洲色图av天堂| 欧美激情在线99| 久久久久国产精品人妻aⅴ院| 91av网站免费观看| www日本在线高清视频| 久久久久久大精品| 香蕉久久夜色| 亚洲av第一区精品v没综合| 国产免费男女视频| 亚洲精品中文字幕一二三四区| 特级一级黄色大片| 久久久久久久精品吃奶| 欧美日韩精品网址| 精品久久久久久久末码| 久久香蕉精品热| 俺也久久电影网| 亚洲国产高清在线一区二区三| 男人的好看免费观看在线视频| 小说图片视频综合网站| 毛片女人毛片| 精品电影一区二区在线| 国产视频内射| 欧美在线黄色| 99久久成人亚洲精品观看| 最新中文字幕久久久久 | 神马国产精品三级电影在线观看| 最新在线观看一区二区三区| 巨乳人妻的诱惑在线观看| 欧美一级毛片孕妇| 黄色女人牲交| 国产又黄又爽又无遮挡在线| 天堂√8在线中文| 一级毛片女人18水好多| a级毛片a级免费在线| 91在线观看av| 最近最新中文字幕大全免费视频| 又黄又爽又免费观看的视频| 亚洲五月天丁香| 黄色成人免费大全| 亚洲五月婷婷丁香| ponron亚洲| 国产精品美女特级片免费视频播放器 | 99久久无色码亚洲精品果冻| 久久国产精品影院| 久久性视频一级片| 午夜视频精品福利| 国产日本99.免费观看| 色在线成人网| 法律面前人人平等表现在哪些方面| 婷婷六月久久综合丁香| 久久久久久久久免费视频了| 成年女人永久免费观看视频| 韩国av一区二区三区四区| 观看美女的网站| 欧美激情久久久久久爽电影| 丝袜人妻中文字幕| 国产高潮美女av| 午夜激情欧美在线| 欧美日韩黄片免| 性色avwww在线观看| 欧美成人免费av一区二区三区| 天天躁日日操中文字幕| 午夜久久久久精精品| 精品福利观看| 少妇裸体淫交视频免费看高清| 成人精品一区二区免费| 免费无遮挡裸体视频| 亚洲国产欧美人成| 日本在线视频免费播放| 此物有八面人人有两片| 亚洲狠狠婷婷综合久久图片| 97超级碰碰碰精品色视频在线观看| 在线十欧美十亚洲十日本专区| 久久久水蜜桃国产精品网| 伦理电影免费视频| 伊人久久大香线蕉亚洲五| 亚洲成av人片免费观看| 久久精品亚洲精品国产色婷小说| 日本五十路高清| 国产欧美日韩一区二区精品| www日本在线高清视频| 国产免费av片在线观看野外av| 午夜激情欧美在线| 麻豆成人av在线观看| 宅男免费午夜| 99久久精品国产亚洲精品| 国产精品一区二区三区四区免费观看 | 日本黄色片子视频| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 又黄又粗又硬又大视频| 国产精品一区二区三区四区久久| 欧美中文综合在线视频| 亚洲美女黄片视频| 日本一二三区视频观看| 午夜日韩欧美国产| 国产精品香港三级国产av潘金莲| 精品一区二区三区视频在线 | 精品国产超薄肉色丝袜足j| 夜夜爽天天搞| 动漫黄色视频在线观看| 美女 人体艺术 gogo| 亚洲美女黄片视频| av黄色大香蕉| 国产精品九九99| 国产精品99久久久久久久久| 色精品久久人妻99蜜桃| 国产高潮美女av| 两人在一起打扑克的视频| 国产高清激情床上av| 国产aⅴ精品一区二区三区波| 十八禁网站免费在线| 99久久精品一区二区三区| 少妇裸体淫交视频免费看高清| 国产精品98久久久久久宅男小说| 麻豆成人午夜福利视频| 91av网一区二区| 国产人伦9x9x在线观看| 亚洲精品美女久久av网站| 香蕉久久夜色| 久久久久国内视频| 97碰自拍视频| 国产成人啪精品午夜网站| 十八禁网站免费在线| 国产熟女xx| 国产精品久久视频播放| 亚洲专区中文字幕在线| 人妻久久中文字幕网| 国产91精品成人一区二区三区| 国产欧美日韩精品一区二区| 国产高清有码在线观看视频| 欧美日韩黄片免| 国产视频内射| 亚洲激情在线av| 精品久久久久久久人妻蜜臀av| 日本免费a在线| 亚洲人成网站在线播放欧美日韩| 非洲黑人性xxxx精品又粗又长| 色噜噜av男人的天堂激情| 91老司机精品| 俄罗斯特黄特色一大片| 黄色成人免费大全| 一进一出抽搐动态| 午夜福利欧美成人| 精品国内亚洲2022精品成人| а√天堂www在线а√下载| 欧美三级亚洲精品| xxxwww97欧美| 99re在线观看精品视频| 亚洲人成伊人成综合网2020| 午夜免费成人在线视频| 很黄的视频免费| 日本免费一区二区三区高清不卡| 亚洲成人久久爱视频| 又粗又爽又猛毛片免费看| 18禁黄网站禁片午夜丰满| 男人舔女人的私密视频| 在线观看免费午夜福利视频| 国产精品综合久久久久久久免费| 一本一本综合久久| 日本黄色视频三级网站网址| 男人的好看免费观看在线视频| av福利片在线观看| 一个人观看的视频www高清免费观看 | 亚洲av成人一区二区三| 久久久久久国产a免费观看| 国产成人啪精品午夜网站| 巨乳人妻的诱惑在线观看| 国产97色在线日韩免费| 久久久久亚洲av毛片大全| x7x7x7水蜜桃| 999久久久国产精品视频| 午夜福利视频1000在线观看| 一区二区三区激情视频| 桃红色精品国产亚洲av| 日本成人三级电影网站| 国产精品久久久久久人妻精品电影| 亚洲最大成人中文| 亚洲欧美日韩卡通动漫| 亚洲专区字幕在线| 在线观看美女被高潮喷水网站 | 国产精品女同一区二区软件 | 亚洲一区高清亚洲精品| 亚洲精品色激情综合| 香蕉久久夜色| 给我免费播放毛片高清在线观看| 日本五十路高清| 国产aⅴ精品一区二区三区波| 久久中文字幕人妻熟女| 99视频精品全部免费 在线 | 香蕉av资源在线| 99久久99久久久精品蜜桃| 丝袜人妻中文字幕| 无遮挡黄片免费观看| 真实男女啪啪啪动态图| 男女之事视频高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 99热这里只有精品一区 | 亚洲av成人精品一区久久| 日本免费a在线| 麻豆久久精品国产亚洲av| 精品日产1卡2卡| АⅤ资源中文在线天堂| 岛国在线观看网站| 国产美女午夜福利| 在线永久观看黄色视频| 国产亚洲精品av在线| 俺也久久电影网| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看日本一区| 老鸭窝网址在线观看| 亚洲精品在线观看二区| 男女那种视频在线观看| 日日干狠狠操夜夜爽| 亚洲成人久久性| 国产视频内射| 国产成人系列免费观看| 日韩欧美在线二视频| 黄片小视频在线播放| 国产精品野战在线观看| 日本免费a在线| 99视频精品全部免费 在线 | 999精品在线视频| 三级毛片av免费| 日本黄大片高清| 免费高清视频大片| 亚洲中文日韩欧美视频| 69av精品久久久久久| 国产极品精品免费视频能看的| 久久久久性生活片| 少妇丰满av| 怎么达到女性高潮| 观看美女的网站| 亚洲人成网站在线播放欧美日韩| 麻豆成人午夜福利视频| 无限看片的www在线观看| 他把我摸到了高潮在线观看| www.精华液| 亚洲精品美女久久av网站| 小蜜桃在线观看免费完整版高清| 国内毛片毛片毛片毛片毛片| 亚洲va日本ⅴa欧美va伊人久久| 欧美激情在线99| 久久香蕉精品热| 俺也久久电影网| 神马国产精品三级电影在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美黄色片欧美黄色片| 国产亚洲精品av在线| 啪啪无遮挡十八禁网站| 欧美三级亚洲精品| 在线观看免费午夜福利视频| 欧美在线黄色| 精品一区二区三区av网在线观看| 一个人免费在线观看的高清视频| 老司机在亚洲福利影院| 麻豆国产av国片精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清videossex| 亚洲中文av在线| 亚洲av五月六月丁香网|