• <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.

    久久中文看片网| 婷婷六月久久综合丁香| 亚洲欧美精品综合一区二区三区| 亚洲在线观看片| 久久久国产成人精品二区| 99在线视频只有这里精品首页| 美女cb高潮喷水在线观看 | 精品日产1卡2卡| 国产97色在线日韩免费| 国产成人一区二区三区免费视频网站| 亚洲avbb在线观看| 怎么达到女性高潮| 男女那种视频在线观看| 午夜久久久久精精品| 国产亚洲精品久久久com| 偷拍熟女少妇极品色| www.999成人在线观看| 欧美午夜高清在线| 国产野战对白在线观看| 久久久国产欧美日韩av| 国产欧美日韩一区二区三| av视频在线观看入口| aaaaa片日本免费| 亚洲欧美日韩高清在线视频| 午夜视频精品福利| 久久中文看片网| 岛国在线免费视频观看| 国产一区二区在线av高清观看| 在线免费观看的www视频| 免费高清视频大片| 久久久久久久精品吃奶| 成熟少妇高潮喷水视频| 国产精品野战在线观看| 亚洲av成人不卡在线观看播放网| 国产精品女同一区二区软件 | 久久久久亚洲av毛片大全| 国产黄a三级三级三级人| 黄片大片在线免费观看| 怎么达到女性高潮| 97碰自拍视频| 我要搜黄色片| 欧美极品一区二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 美女午夜性视频免费| 亚洲九九香蕉| 国产激情久久老熟女| 亚洲最大成人中文| 久久这里只有精品中国| АⅤ资源中文在线天堂| 最新中文字幕久久久久 | 欧美性猛交╳xxx乱大交人| 俄罗斯特黄特色一大片| xxxwww97欧美| 国产亚洲精品综合一区在线观看| 少妇人妻一区二区三区视频| 母亲3免费完整高清在线观看| 老司机福利观看| 免费看光身美女| 国产三级在线视频| 亚洲七黄色美女视频| 性欧美人与动物交配| 国产亚洲精品一区二区www| 神马国产精品三级电影在线观看| 国内毛片毛片毛片毛片毛片| 首页视频小说图片口味搜索| 曰老女人黄片| 在线永久观看黄色视频| 1024手机看黄色片| 色综合婷婷激情| 此物有八面人人有两片| 欧美日韩中文字幕国产精品一区二区三区| 国产日本99.免费观看| 9191精品国产免费久久| 男人舔女人下体高潮全视频| 亚洲aⅴ乱码一区二区在线播放| 欧美一级毛片孕妇| www日本在线高清视频| 18禁黄网站禁片免费观看直播| 91在线精品国自产拍蜜月 | 岛国在线观看网站| 久久精品国产清高在天天线| 国产亚洲精品久久久com| 亚洲欧美日韩高清在线视频| 人人妻人人澡欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 亚洲男人的天堂狠狠| 欧美日韩精品网址| АⅤ资源中文在线天堂| aaaaa片日本免费| 欧美日本视频| 欧美av亚洲av综合av国产av| 两个人视频免费观看高清| 欧美中文综合在线视频| 一级黄色大片毛片| 99在线人妻在线中文字幕| 在线观看日韩欧美| 国产熟女xx| 天天躁狠狠躁夜夜躁狠狠躁| 欧美xxxx黑人xx丫x性爽| 国产aⅴ精品一区二区三区波| 成人特级黄色片久久久久久久| 一二三四社区在线视频社区8| 熟女人妻精品中文字幕| 黑人操中国人逼视频| 亚洲自拍偷在线| 精品免费久久久久久久清纯| 国产高清videossex| 黑人巨大精品欧美一区二区mp4| 国产精品,欧美在线| 久9热在线精品视频| 亚洲中文日韩欧美视频| 精品国产三级普通话版| 亚洲性夜色夜夜综合| а√天堂www在线а√下载| 精品久久久久久久末码| a在线观看视频网站| 久久这里只有精品19| 制服人妻中文乱码| 国产黄片美女视频| 精品久久久久久,| 色综合亚洲欧美另类图片| 一二三四社区在线视频社区8| 国产精品女同一区二区软件 | 久久中文字幕一级| 国产成人影院久久av| 一卡2卡三卡四卡精品乱码亚洲| 看黄色毛片网站| 日韩中文字幕欧美一区二区| 国产真人三级小视频在线观看| 舔av片在线| 青草久久国产| 免费在线观看亚洲国产| 中文在线观看免费www的网站| 国产成+人综合+亚洲专区| 曰老女人黄片| 嫩草影视91久久| 亚洲av美国av| 成在线人永久免费视频| 成人av在线播放网站| 国产aⅴ精品一区二区三区波| 夜夜看夜夜爽夜夜摸| 亚洲国产高清在线一区二区三| 色噜噜av男人的天堂激情| 免费看美女性在线毛片视频| 一级作爱视频免费观看| 亚洲欧美日韩高清在线视频| 国产一级毛片七仙女欲春2| 国产91精品成人一区二区三区| 91在线精品国自产拍蜜月 | 久久这里只有精品中国| 狠狠狠狠99中文字幕| 一级作爱视频免费观看| 一级毛片精品| 不卡一级毛片| 一二三四社区在线视频社区8| 成年人黄色毛片网站| 国产成人啪精品午夜网站| 国产精品亚洲美女久久久| 久久久国产欧美日韩av| 极品教师在线免费播放| 美女高潮喷水抽搐中文字幕| 久久久久国内视频| 亚洲成人久久爱视频| 日本在线视频免费播放| 曰老女人黄片| 少妇人妻一区二区三区视频| 老鸭窝网址在线观看| 熟女人妻精品中文字幕| 国产av麻豆久久久久久久| 国产午夜福利久久久久久| 人人妻人人澡欧美一区二区| 好男人电影高清在线观看| 亚洲中文av在线| 中文字幕精品亚洲无线码一区| 99久久99久久久精品蜜桃| 97碰自拍视频| 99国产综合亚洲精品| 99久久综合精品五月天人人| 在线播放国产精品三级| 不卡av一区二区三区| 国产精品亚洲一级av第二区| 三级国产精品欧美在线观看 | 亚洲色图 男人天堂 中文字幕| 精品久久久久久成人av| 国产探花在线观看一区二区| 国产精品影院久久| 亚洲av电影在线进入| 国产成人av教育| 欧美最黄视频在线播放免费| 亚洲 欧美 日韩 在线 免费| 午夜精品久久久久久毛片777| 国产激情久久老熟女| 变态另类成人亚洲欧美熟女| 国产爱豆传媒在线观看| 香蕉av资源在线| 母亲3免费完整高清在线观看| 欧美极品一区二区三区四区| 老司机深夜福利视频在线观看| 日本黄大片高清| 国产精品一区二区三区四区久久| 给我免费播放毛片高清在线观看| 法律面前人人平等表现在哪些方面| 亚洲国产高清在线一区二区三| 波多野结衣高清作品| 国模一区二区三区四区视频 | 91麻豆av在线| 一级毛片精品| 国产精品1区2区在线观看.| 搡老岳熟女国产| 欧美zozozo另类| 国产成人啪精品午夜网站| 成人特级黄色片久久久久久久| 亚洲精品美女久久久久99蜜臀| 这个男人来自地球电影免费观看| 少妇丰满av| 最近最新免费中文字幕在线| 首页视频小说图片口味搜索| 亚洲午夜精品一区,二区,三区| 男人舔女人下体高潮全视频| 天天添夜夜摸| 婷婷丁香在线五月| 久久久久久久久中文| 国产av在哪里看| 成年女人永久免费观看视频| 两人在一起打扑克的视频| 国产精品野战在线观看| 日本黄大片高清| 国产私拍福利视频在线观看| 久久久久国内视频| 天堂影院成人在线观看| 超碰成人久久| 男女下面进入的视频免费午夜| 露出奶头的视频| 99riav亚洲国产免费| 香蕉国产在线看| 久久热在线av| 亚洲欧美日韩高清在线视频| 黑人巨大精品欧美一区二区mp4| 97人妻精品一区二区三区麻豆| 国产高清视频在线观看网站| 这个男人来自地球电影免费观看| 亚洲第一电影网av| av在线蜜桃| 制服丝袜大香蕉在线| 亚洲av电影不卡..在线观看| 女同久久另类99精品国产91| 国产欧美日韩精品亚洲av| 国产成人av教育| 91av网一区二区| 欧美午夜高清在线| 亚洲黑人精品在线| 老司机福利观看| 国产精品久久久久久久电影 | 亚洲激情在线av| 黄频高清免费视频| 午夜精品久久久久久毛片777| 99热精品在线国产| 神马国产精品三级电影在线观看| 国产精品精品国产色婷婷| 日本撒尿小便嘘嘘汇集6| 超碰成人久久| 久久天堂一区二区三区四区| 一本久久中文字幕| 国产 一区 欧美 日韩| 人妻夜夜爽99麻豆av| 九九热线精品视视频播放| 欧美黄色淫秽网站| 小蜜桃在线观看免费完整版高清| 国产真人三级小视频在线观看| 欧美成狂野欧美在线观看| 国产高清三级在线| 久久久久久人人人人人| 偷拍熟女少妇极品色| 18禁裸乳无遮挡免费网站照片| 级片在线观看| 在线观看舔阴道视频| 一二三四在线观看免费中文在| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 亚洲自偷自拍图片 自拍| av欧美777| 男插女下体视频免费在线播放| 91在线观看av| 搡老岳熟女国产| 俺也久久电影网| 99热这里只有是精品50| 每晚都被弄得嗷嗷叫到高潮| 九九热线精品视视频播放| 久久久久性生活片| 免费看光身美女| 国产精品久久久久久精品电影| 国产成人aa在线观看| 亚洲av五月六月丁香网| 欧美黄色片欧美黄色片| 国产成人av教育| 亚洲中文日韩欧美视频| 精品久久久久久,| 桃红色精品国产亚洲av| 少妇的逼水好多| 99国产极品粉嫩在线观看| 国产精品日韩av在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 两个人的视频大全免费| 最近最新中文字幕大全免费视频| 成人三级做爰电影| 日本黄大片高清| 免费大片18禁| 99久国产av精品| 日韩中文字幕欧美一区二区| 日韩三级视频一区二区三区| 18禁黄网站禁片午夜丰满| 99国产精品99久久久久| 亚洲成人免费电影在线观看| 老汉色av国产亚洲站长工具| tocl精华| 脱女人内裤的视频| av国产免费在线观看| 亚洲五月婷婷丁香| 亚洲av五月六月丁香网| 夜夜爽天天搞| av福利片在线观看| 亚洲黑人精品在线| 亚洲欧美精品综合一区二区三区| a级毛片在线看网站| 欧美绝顶高潮抽搐喷水| 国产精品精品国产色婷婷| 一级毛片高清免费大全| 亚洲中文av在线| 一级毛片女人18水好多| 国产精品一区二区三区四区久久| 欧美一级a爱片免费观看看| 18禁国产床啪视频网站| 国产aⅴ精品一区二区三区波| 久久婷婷人人爽人人干人人爱| 色综合欧美亚洲国产小说| 国产男靠女视频免费网站| 在线观看免费午夜福利视频| 精品电影一区二区在线| 国产真实乱freesex| 999精品在线视频| 国产亚洲精品久久久com| 精品不卡国产一区二区三区| 黄片大片在线免费观看| 国产一区二区三区视频了| 午夜福利视频1000在线观看| 亚洲一区二区三区色噜噜| 禁无遮挡网站| 我的老师免费观看完整版| 欧美成狂野欧美在线观看| 夜夜看夜夜爽夜夜摸| 18禁裸乳无遮挡免费网站照片| 在线观看66精品国产| 嫩草影院入口| 在线观看66精品国产| 色精品久久人妻99蜜桃| 久久香蕉精品热| 18禁裸乳无遮挡免费网站照片| 小说图片视频综合网站| 国产亚洲av嫩草精品影院| 国产精品久久久久久久电影 | 亚洲精品乱码久久久v下载方式 | 国产成人影院久久av| 在线a可以看的网站| 夜夜爽天天搞| 国产视频一区二区在线看| 中国美女看黄片| 欧美日韩综合久久久久久 | 国产免费男女视频| 岛国在线免费视频观看| 国语自产精品视频在线第100页| cao死你这个sao货| 在线观看66精品国产| 国模一区二区三区四区视频 | 亚洲人与动物交配视频| 亚洲真实伦在线观看| av片东京热男人的天堂| 精华霜和精华液先用哪个| 中出人妻视频一区二区| 一a级毛片在线观看| 亚洲男人的天堂狠狠| 久久这里只有精品中国| 观看免费一级毛片| 99在线人妻在线中文字幕| 日日夜夜操网爽| 久久久水蜜桃国产精品网| 一级毛片精品| 极品教师在线免费播放| 岛国在线观看网站| 中国美女看黄片| 全区人妻精品视频| 精品一区二区三区视频在线 | 老汉色av国产亚洲站长工具| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 国产熟女xx| or卡值多少钱| 黑人欧美特级aaaaaa片| 一区福利在线观看| 成年版毛片免费区| 露出奶头的视频| 国产成人啪精品午夜网站| 亚洲成人精品中文字幕电影| 男人的好看免费观看在线视频| 熟妇人妻久久中文字幕3abv| 男人舔女人下体高潮全视频| 免费看十八禁软件| 窝窝影院91人妻| 久久久成人免费电影| 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| 少妇丰满av| x7x7x7水蜜桃| 一本综合久久免费| 久久精品综合一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 亚洲乱码一区二区免费版| 亚洲国产色片| h日本视频在线播放| 我要搜黄色片| www.熟女人妻精品国产| 欧美在线黄色| 成人鲁丝片一二三区免费| 日本黄大片高清| 久久亚洲真实| 欧美精品啪啪一区二区三区| 淫秽高清视频在线观看| 日本三级黄在线观看| 美女大奶头视频| 国产精品亚洲av一区麻豆| 18禁黄网站禁片免费观看直播| www.精华液| 国产91精品成人一区二区三区| 午夜精品一区二区三区免费看| 亚洲精品国产精品久久久不卡| 黑人操中国人逼视频| 亚洲精品中文字幕一二三四区| 国产乱人伦免费视频| 久久久久九九精品影院| 好男人电影高清在线观看| 亚洲一区高清亚洲精品| 黄色女人牲交| 99热精品在线国产| 91av网站免费观看| 在线a可以看的网站| 好男人在线观看高清免费视频| 国产三级中文精品| 亚洲av电影不卡..在线观看| 亚洲av美国av| 国产亚洲精品av在线| 狠狠狠狠99中文字幕| 女同久久另类99精品国产91| 国产美女午夜福利| www日本黄色视频网| 国产精品永久免费网站| 国产黄a三级三级三级人| 亚洲国产日韩欧美精品在线观看 | 桃红色精品国产亚洲av| 老汉色∧v一级毛片| 99riav亚洲国产免费| 国语自产精品视频在线第100页| 偷拍熟女少妇极品色| 日韩欧美在线二视频| 久久久国产欧美日韩av| 国产97色在线日韩免费| 波多野结衣高清无吗| 久久久久国内视频| 亚洲熟女毛片儿| 99在线视频只有这里精品首页| 麻豆久久精品国产亚洲av| 一级作爱视频免费观看| 日韩成人在线观看一区二区三区| 黄色丝袜av网址大全| 中文字幕人妻丝袜一区二区| 中文字幕高清在线视频| 日韩av在线大香蕉| 亚洲精品在线美女| 在线免费观看不下载黄p国产 | 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区久久| 搡老熟女国产l中国老女人| 男女视频在线观看网站免费| 成人av一区二区三区在线看| 久久精品国产99精品国产亚洲性色| 久久精品夜夜夜夜夜久久蜜豆| 亚洲在线观看片| 天天添夜夜摸| 变态另类成人亚洲欧美熟女| 亚洲人与动物交配视频| 最新中文字幕久久久久 | 黄色女人牲交| 国产1区2区3区精品| 免费观看人在逋| 国产精品亚洲一级av第二区| x7x7x7水蜜桃| 国产蜜桃级精品一区二区三区| 啦啦啦韩国在线观看视频| 国产精品精品国产色婷婷| 国产黄片美女视频| 亚洲成人久久爱视频| 欧美在线黄色| 看免费av毛片| а√天堂www在线а√下载| 欧美三级亚洲精品| 成年免费大片在线观看| 香蕉av资源在线| 国产成人精品久久二区二区免费| 母亲3免费完整高清在线观看| 亚洲成人久久性| 99久国产av精品| 欧美黑人欧美精品刺激| 国产免费男女视频| 禁无遮挡网站| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久久久久免费视频| 国产精品自产拍在线观看55亚洲| 精品不卡国产一区二区三区| aaaaa片日本免费| 丰满的人妻完整版| 成人精品一区二区免费| 最近最新免费中文字幕在线| 伊人久久大香线蕉亚洲五| 精品久久久久久成人av| 久久久久久国产a免费观看| 亚洲午夜精品一区,二区,三区| av女优亚洲男人天堂 | 真实男女啪啪啪动态图| 久久中文字幕人妻熟女| а√天堂www在线а√下载| 亚洲中文字幕日韩| 99久久精品国产亚洲精品| 久久精品国产亚洲av香蕉五月| 看片在线看免费视频| 国产成人精品久久二区二区91| 国产午夜福利久久久久久| 国产精品亚洲av一区麻豆| 在线免费观看的www视频| 成人永久免费在线观看视频| 亚洲成av人片免费观看| 精华霜和精华液先用哪个| 一级毛片女人18水好多| 两个人视频免费观看高清| 麻豆久久精品国产亚洲av| 91麻豆av在线| 久久香蕉国产精品| 国产精品98久久久久久宅男小说| 美女 人体艺术 gogo| 精品国产美女av久久久久小说| 亚洲国产欧美网| a在线观看视频网站| 欧美日本亚洲视频在线播放| 亚洲天堂国产精品一区在线| 香蕉丝袜av| 亚洲中文字幕一区二区三区有码在线看 | 小蜜桃在线观看免费完整版高清| 色精品久久人妻99蜜桃| 女人高潮潮喷娇喘18禁视频| 国产高清三级在线| а√天堂www在线а√下载| 国产高潮美女av| 蜜桃久久精品国产亚洲av| 亚洲无线在线观看| 亚洲成人免费电影在线观看| 欧美中文日本在线观看视频| 波多野结衣高清作品| 成人国产综合亚洲| 日韩有码中文字幕| 久久精品国产亚洲av香蕉五月| 9191精品国产免费久久| 免费观看的影片在线观看| 无遮挡黄片免费观看| 91老司机精品| 亚洲精品中文字幕一二三四区| 午夜视频精品福利| 亚洲第一电影网av| 国产精品女同一区二区软件 | 日本黄色片子视频| 亚洲精品美女久久久久99蜜臀| 婷婷亚洲欧美| 精品无人区乱码1区二区| 亚洲,欧美精品.| 亚洲最大成人中文| 黄色片一级片一级黄色片| 国产蜜桃级精品一区二区三区| 搡老熟女国产l中国老女人| 男人和女人高潮做爰伦理| 最新中文字幕久久久久 | 波多野结衣高清作品| 老熟妇乱子伦视频在线观看| 很黄的视频免费| 深夜精品福利| 亚洲av五月六月丁香网| 欧美色欧美亚洲另类二区| 五月玫瑰六月丁香| 国产黄片美女视频| 色视频www国产| 国产av不卡久久| 三级男女做爰猛烈吃奶摸视频| 免费在线观看亚洲国产| 五月玫瑰六月丁香| 国产黄片美女视频| 最新在线观看一区二区三区| 一二三四社区在线视频社区8| 一个人看的www免费观看视频| 精品一区二区三区四区五区乱码| 国产成年人精品一区二区| 老司机在亚洲福利影院| 亚洲成a人片在线一区二区| 老汉色av国产亚洲站长工具| 一二三四在线观看免费中文在| 女人被狂操c到高潮| 亚洲国产欧美人成| 亚洲国产高清在线一区二区三| 熟女少妇亚洲综合色aaa.| 欧美性猛交╳xxx乱大交人| 99久久久亚洲精品蜜臀av|