Introduction

Wheat-cotton relay intercropping is widely practiced by farmers in the Yellow River cotton producing region of China (Zhang et al. 2007a), because it meets the need of farmers to grow a profitable cash crop as well as to secure food supply. Relay intercropping enables the farmers to grow cotton, in a rotation with winter wheat, when the heat resource does not match the heat requirement of cotton if sown after harvesting winter wheat. The cotton is sown in April, approximately seven weeks before the harvest date of wheat. Strips are left open in the wheat crop at sowing (October/November) to provide space for the cotton plants during their seedling stage (April, May and June). After the wheat harvest in June, cotton plants can exploit the full space, above-ground as well as below-ground. A cotton-wheat relay intercropping system is thus characterized by three main phases: (1) winter wheat (vegetative stage) grown in strips from November till April; (2) intercropping of wheat (reproductive stage) and cotton (seedling stage) from April till June, and (3) sole cotton (vegetative and reproductive stage) from June till October. The two component crops in the system interact directly only during the second phase; however the physiology, ecology and productivity of the relay strip intercropping system are determined by the spatial architecture and temporal dynamics of the leaf canopy and the root systems during the whole growing cycle. There is hardly any information available on the N utilization and requirement of cotton and wheat in relay intercropping systems. This information is needed to develop profitable and sustainable systems.

Willey (1990) suggested that intercrops may well improve the total nutrient capture by taking up nutrients that might be leached in a monocropping system. As a result intercropping could make greater demands on the soil and thus, in the long term, intercropping yield advantages could have to be paid for with higher fertilizer inputs. Still, even though the intercropping system as a whole may consume more nitrogen, each component crop in an intercropping system is likely to take up less nitrogen than in a monocrop situation, due to competition with the other crop. This, for instance, was found by Blaise et al. (2005), who reported that total N uptake in strip intercropping of cotton and pigeon pea in rainfed regions of central India was lower than in sole cotton because of a lower crop productivity.

Efficiency of N-use per each component crop is diagnosed by nitrogen use efficiency (NUE) and is commonly measured by internal (or physiological) efficiency (IE), expressed as kg yield per kg N uptake (Haefele et al. 2003). Several studies show that intercrops use soil nutrients more efficiently than sole crops, because of a higher N-recovery and increased dry matter yields (Hauggaard-Nielsen et al. 2001; Zhang et al. 2004). However, agronomic N-use efficiency may also be decreased by intercropping. For instance, Aggarwal et al. (1992) found that upland rice intercropped with grain legumes took up the same amount of N as a monocrop, but produced significantly less dry matter.

To determine the land use advantage in terms of N yield for intercrops, a modification of the land equivalent ratio (Willey 1979), expressed as a relative N yield total (RNT) (Baumann et al. 2001) or land equivalent ratio (LER) for N yield (NLER) (Szumigalski and Van Acker 2006), can be employed. The RNT is calculated as the sum of the ratios of N-uptakes by component crops in the intercrop to their respective N-uptakes in monoculture. A value of RNT exceeding the LER suggests that intercropping is not nitrogen-efficient.

In many crops, N uptake is related with the accumulation of biomass and growth of leaf area over time (Booij et al. 1996; Lemaire et al. 2007). In case of N shortage, the actual uptake will fall short of N demand, and the concentration of N in the leaves will decrease more rapidly than in a well-fertilized crop. The robust relationships between N uptake, biomass accumulation and leaf area growth thus enable prediction of N uptake in time and provide a basis for the optimization of N application, both in quantity and in timing. Leigh and Johnston (1987) suggested a close association between leaf N concentration, rate of photosynthesis and biomass growth. To assess the balance between N uptake and demand, and to diagnose N deficiency the observed N dilution in the biomass can be compared to the minimal concentration of N in shoots necessary to produce the maximum dry mass (Justes et al. 1994). The trend in N dilution with increasing biomass can be represented by a power relationship (Flenet et al. 2006). The N dilution curve can be used to diagnose N deficiency, to manage N fertilization and for modeling N allocation (Lemaire et al. 2007). N dilution curves have been widely studied for monocrops but not for crops in intercropping systems. Considering the relatively high cost of nitrogen (N) fertilizer, and environmental concerns associated with excessive N application, increasing the N use efficiency of cropping systems is needed urgently.

The N-economy of a crop is the result of many processes occurring in the soil, crop and atmosphere and can be quantified by a wide variety of parameters (Peng and Bouman 2007; Spiertz and Vos 2005). In this study, we determined N-uptake in relation to DM-yield, physiological N use efficiency, and N-dilution in biomass over time. Furthermore, a nitrogen balance sheet analysis was made based on estimates of N-input and -output in the systems. We used different indicators at the crop and system level to analyze the nitrogen economy and use efficiency (Table 1).

Table 1 Indicators for analysis of N use in relay intercropping systems and monocultures

The specific objectives are: (a) to quantify the nitrogen uptake and nitrogen use efficiency of wheat and cotton in relay intercropping systems; (b) to determine the relationships between N uptake and accumulation of biomass and growth of leaf area; (c) to quantify N dynamics of cotton in relation to the effects of intercropping; and (d) to explore the opportunities for a more effective nitrogen management of wheat-cotton intercropping systems.

Materials and methods

Field experiments

Field experiments were conducted in 2001/02, 2002/03 and 2003/04 consecutively on the same field at the Cotton Research Institute of Chinese Academy of Agricultural Sciences (CRI, CAAS), Anyang city, Henan province, China at 36°07’North and 116°22’East.

Soil parameters of the field are: sandy loam, pH 8.0, bulk density 1.36 g cm−3, and organic matter content 13.2 g kg−1. At the start of the experiments, in October 2001, the soil contained 1.02 g kg−1 total N, 0.52 g kg−1 P and 17.3 g kg−1 K. The amount of precipitation in 2002, 2003 and 2004 was 318, 539 and 517 mm, respectively.

For details on crop development we refer to Zhang et al. (2007b)

Experiment 1: N economy of intercropping and monoculture systems

This experiment was conducted in 2001/02, 2002/03 and 2003/04, and comprised six treatments including four different intercropping patterns and monocultures of wheat (Triticum aestivum L.) and cotton (Gossypium hirsutum L.). The four intercropping patterns were all strip intercrops, with strips of cotton and wheat alternating. The number of wheat vs. cotton rows per strip in the four systems was 3:1, 3:2, 4:2 and 6:2. Row distance in wheat was 20 cm and between cotton rows 40 cm, but distance between cotton and wheat was varied (Fig. 1). The experimental systems are characterized by the width of a “minimal combination,” i.e. an adjacent wheat and cotton strip, the row length density (total row length per m2, or simply the number of rows per m, expressed over the whole intercropped area) and homogenized densities, i.e. densities of either crop, expressed per unit total intercrop area (Table 2).

Fig. 1
figure 1

Layout of wheat-cotton intercropping systems and monocultures (reproduced with permission from Zhang et al. 2007a). In the sole wheat crop, the row distance is 20 cm, and in the sole cotton crop, row distance is 80 cm. The sole crops do not have a strip structure. Intercrops, on the contrary, do have a strip structure, with strips of one or more crop rows of wheat and cotton alternating. The 3:1 system, for example, alternates strips of three wheat rows (sown at 20 cm distance from each other) with single rows of cotton. The distance between the cotton and wheat rows in this system is 30 cm. The 3:2 system uses the same design for the wheat strip, but it has a wider path in which two rows of cotton are sown. The distance between the cotton rows is 40 cm, and the distance between wheat and cotton rows is 20 cm. In the 4:2 system, there are 4 wheat rows per strip, and there is 25 cm between the cotton and wheat rows. In the 6:2 system, there are 6 wheat rows per strip, and there is 30 cm between cotton and wheat rows

Table 2 Characteristics of relay intercropping systems and monocultures (Exp. 1)

The six treatments were arranged in four randomized blocks with a plot size of 180 m2. Wheat was sown on 4 November 2001, 2 November 2002 and 3 November 2003. Cotton was sown on 26 April 2002, 25 April 2003 and 25 April 2004. The wheat cultivar was ‘Zhongyu 5’ and cotton cultivars were middle maturing upland Bt cotton ‘Shiyuan 321’ in 2002 and the Verticillium-tolerant variety ‘CCRI45’in 2003 and 2004.

Irrigation water was applied by flooding (350 mm) in 2002 and by drip application in 2003 (342 mm) and 2004 (182 mm). The total amount of nitrogen applied to wheat and cotton ranged from 302 to 412 kg N ha−1 y−1. Information on amount, composition and application dates of fertilizer materials is provided in Table 3. Fertilizers were applied evenly over the whole area of the plots, not discriminating between crop strips and paths.

Table 3 Type, amount and timing of nitrogen applied in Exp. 1

Experiment 2: Response of the intercropping system to nitrogen input

The response of the 3:2 relay intercropping system to three nitrogen doses was investigated in 2002/03 and 2003/04. Cultivars and planting dates were the same as in Expt 1, but the minimum combination width was 20 cm wider. The amount of irrigation water was 177 mm over the whole wheat and cotton growing period. As in Experiment 1, fertilizer doses were applied evenly on strips with and without a crop. Nitrogen was applied at full dose (180 kg ha−1), at half dose, and at 150% dose. Before wheat was sown in 2002, 57 kg N ha−1 was applied as cotton seed cake and 47 kg N ha−1 as compound fertilizer; no base fertilizer was applied in 2003. During the growing season, urea (46% N) was applied on wheat and cotton by drip irrigation according to the designated dose level (50, 100 or 150%).

Measurements

To determine crop yields, plants were harvested at maturity in each plot; the sampling area covered 5 m row length by 2 m width in the 3:1 system, 5 m length by 2.4 m width in the 3:2 system, 3.5 m row length by 3 m width in the 4:2 system and 5 m row length by 2 m width in the 6:2 system.

To assess above-ground dry mass (DM) and nitrogen uptake, 1 m row of each plot was sampled once per 2 weeks. First the number of plants was counted; next a sub-sample was selected for analysis. The subsamples consisted of twenty randomly selected wheat plants, and in the case of cotton of 10 seedlings, or, later on, three plants.

The leaf area of the subsampled cotton plants was determined by measuring length and width of individual leaves, and using the formula (Zhang, unpublished):

leaf area = 0.810 × length × width

The samples were subsequently oven-dried at 65°C to determine dry matter weight.

Nitrogen content of plant and soil samples was determined by the Micro-Kjeldahl method (Ogg 1960). Leaves and stems were weighted and analysed separately in 2002. In 2003 and 2004, analyses were made using the total vegetative aboveground biomass. N-content of reproductive organs (i.e. the spike of wheat and the squares, flowers and bolls of cotton) was determined separately.

Soil samples were taken before wheat sowing in October 2001, 2002 and 2003, taking five samples per plot with a 3 cm diameter auger. Soil samples from 0–30 cm depth and from 30–60 cm depth was analysed separately.

Data analysis

Nitrogen use efficiency

Internal efficiency of nitrogen use (IE) is defined as:

$${\text{IE}}_i = \frac{{Y_i }}{{{\text{N}}_i }}$$
(1)

Where Y i is the yield of crop i (g harvestable product m−2) and N i is nitrogen uptake (g total N-uptake m−2) by crop i.

Relative nitrogen yield

The relative nitrogen uptake is calculated as:

$${\text{RNT}} = {\text{RNW}} + {\text{RNC}} = \frac{{{\text{N}}_{{\text{W,I}}} }}{{{\text{N}}_{{\text{W,S}}} }} + \frac{{{\text{N}}_{{\text{C,I}}} }}{{{\text{N}}_{{\text{C,S}}} }}$$
(2)

In this equation RNW and RNC are the relative values for wheat and cotton N yield, respectively. NW,I, NW,S, NC,I and NC,S represent nitrogen uptake of intercropped and sole wheat, and of intercropped and sole cotton, respectively.

Nitrogen dilution effect

The N dilution curve is described using a power function (Flenet et al. 2006):

$${\text{N}}\% = a_{\text{c}} W^{b_{\text{c}} } $$
(3)

where N% is nitrogen concentration, W is above-ground dry mass, a c is the N percentage at W = 1 g m−2 and b c is a curvature parameter.

Statistical analysis

Data on grain and lint yield, harvest index, nitrogen content, and total nitrogen uptake were analysed by ANOVA in SPSS 11.0, using a randomized block design with cropping system as fixed effect. Least significant differences (LSD) were used to separate treatment means (P < 0.05).

Results

Crop yields and total N-uptake

Crop yields, harvest indices and total nitrogen uptake by the intercropped wheat and cotton, compared to monocultures, are shown in Table 4. The nitrogen yield of wheat in the intercropping systems ranged from 130 to 205 kg ha−1, which was significantly lower (P < 0.01) than the 3-year’s average N yield of sole wheat (253 kg ha−1). Intercropping reduced the nitrogen yield of wheat proportionally to grain yield, reflecting mainly the differences in homogenized plant densities. The reduction in nitrogen yield of wheat by intercropping was greatest in the 3:2 and 4:2 systems, and smallest in the 3:1 system. The response of the 6:2 system varied between years.

Table 4 Homogenized yield, harvest index (HI), nitrogen yield and Internal use efficiency (IE) for three cropping cycles (Exp. 1)

The effect of intercropping on nitrogen uptake of cotton differed between years. In 2002, the nitrogen yield of cotton in intercrops ranged from 120 to 140 kg ha−1 compared to an N yield in the monoculture of 127 kg ha−1, with no significant difference (P > 0.05) between the monocrop and intercrops. However, in 2003 and 2004, cotton N yields in intercrops were significantly lower than in monocrops, except for the 3:2 system. Especially, the total N uptake of cotton in the 6:2 system was extremely low (Table 4), which was caused by an incomplete canopy closure and low biomass and lint yield.

The relative nitrogen yield of cotton (RNC) in intercropping systems was much higher than the relative lint yield (RYC), indicating that the lint yield was much more diminished by intercropping than nitrogen uptake (Fig. 2). RNC values did not differ significantly among the intercropping patterns 3:2, 4:2 and 3:1; but the value was significantly lower for the 6:2 system (P < 0.05). For wheat there were hardly any differences between the relative nitrogen yield of wheat (RNW) and the relative grain yield of wheat (RYW) (Fig. 2). Both parameters, RNW and RYW, were closely associated with row length densities (total row length of one crop species per m2, averaged over the whole field, m m−2). The values were somewhat higher in the 3:1 and 6:2 systems and lower in the 3:2 system.

Fig. 2
figure 2

RNT, relative N yield total (a, c and e) and LER, relative agronomic yield (b, d and f) of wheat and cotton in the intercropping systems, in three years experimentation. SE’s pertain to the contributions of wheat and cotton to the RNT and LER

N-use efficiency of wheat and cotton

The physiological nitrogen use efficiencies (IE) of wheat varied from 24.4 to 30.9 kg grain per kg N uptake for intercrops and from 25.5 to 27.7 kg kg−1 for the monocrop. In 2002 and 2003, the IE of wheat was not significantly affected by intercropping (Table 4), but in 2004, wheat in the 3:2 system showed a significant higher IE-value than in the other systems including the monocrop.

The effect of intercropping on the IE of cotton varied between years. In 2002, it ranged from 4.5 to 6.1 kg lint per kg N uptake in intercropping systems, which was significantly (P<0.02) lower than in sole cotton (9.7 kg kg−1). Due to the large variation, IE did not significantly differ between systems in 2003 and 2004. Averaged over 3 years, the IE values amounted to 6.4, 6.7, 7.5, 8.2 and 9.7 for the 3:1, 3:2, 4:2, 6:2 systems and the monoculture, respectively. Lower IE’s of cotton in intercrops, compared to monoculture, were mainly due to lower DM accumulation and lower (15–30%) harvest indexes (HI).

Crop yields and total N-uptake under reduced input

In Experiment 2 (2002/2003 and 2003/2004), the nitrogen fertilizer input to a 3:2 intercrop was varied between 38 and 150 kg ha−1 during the wheat phase and between 53 and 210 kg ha−1 during the cotton phase in intercropping systems. Grain yields of wheat were somewhat lower under the highest nitrogen application rate (210 kg ha−1) than with the other two N-inputs, due to a lower HI at higher N-input (Table 5). Lint yields of cotton differed strongly between the two seasons, however, did not respond to the reduced N-rates. Despite the substantial difference in N-inputs between treatments, the uptake of nitrogen did not differ significantly (Table 5), indicating ample soil N-reserves. Estimated apparent N recoveries in 2002/2003 and 2003/2004 amounted to 115 and 141 kg ha−1 in the wheat crop, and 80 and 87 kg ha−1 in the cotton crop. Thus matching most of the N requirement of both crops.

Table 5 Homogenized yield, above-ground dry mass (DM) and nitrogen uptake (Nu) under different N application rates (Na) in the 3:2 system (Exp. 2)

Cotton nitrogen dynamics

During the early growth stages, nitrogen uptake was slower in intercropped cotton than in the monocrop (Fig. 3), but uptake caught up in intercrops from 80 days after sowing (DAS). N-uptake leveled off from 110 DAS. In all the 3 years, the rate of N-uptake was most fast and total N-uptake reached the highest level in the 3:2 system and both were lowest in the 6:2 system. In each of the years, cotton in the 3:1 system had a low rate of N-uptake early on, but cotton in this system accumulated N at the fastest rate of all systems at the end of the season. Cotton in the 4:2 system showed similar patterns of N-uptake as cotton in the 3:2 system, but at a somewhat slower rate and lower final level, reflecting the lower density of cotton plants in this system, compared to the 3:2 system (Table 2). Likewise, the low level of nitrogen uptake in the 6:2 system can be explained by the lower canopy density and biomass yield in comparison to other systems (Tables 2 and 4).

Fig. 3
figure 3

Nitrogen uptake by cotton in the intercropping systems and the monoculture in 2002–2004

Major influences of intercropping on the cumulative N-uptake by the cotton seedlings are shown at the end of the intercrop period, i.e. June 14 (49 DAS) in 2002, June 10 (46 DAS) in 2003 and June 17 (53 DAS) in 2004 (Fig. 4). All intercropping systems show large and significant differences in relative N-yield per plant with the monoculture, ranging from 0.18 to 0.48 over 3 years, 2002–2004. There are no significant differences between intercropping systems, except for the significantly lower N-accumulation in the 3:1 system compared to the other intercropping systems in 2003.

Fig. 4
figure 4

Relative nitrogen yield of cotton per plant in the intercropping systems at the time of wheat harvest in 2002–2004

N accumulation in cotton at the end of the wheat growing period is a reflection of the competitive effects of the wheat on the cotton, reflected in the N-content of the above-ground biomass (Fig. 5). The competitive effect is closely related to the width of the space between the wheat strips, which is narrowest in the 3:1 system (60 cm), wider in 3:2 (80 cm) and 4:2 (90 cm), and widest in 6:2 system (100 cm; Fig. 1). Correspondingly, the N-content of cotton seedlings is the least in the 3:1 system, greater in 3:2, still greater in the 4:2 and 6:2 systems, and greatest in the cotton monoculture (Fig. 5). Averaged over 3 years the N content of cotton seedlings in intercrops ranged from about 3.2 to 3.4%, which is significantly lower (P < 0.01) than in monoculture (3.6%). Cotton had a significantly lower N-content in the 3:1 system than in any other intercropping system (P < 0.05). Thus wheat exhibited significant effects on N-uptake of cotton seedlings during intercropping.

Fig. 5
figure 5

Nitrogen content of above-ground cotton biomass in the intercropping systems and the monoculture at the time of the wheat harvest in 2002–2004

Nitrogen dilution curves for wheat and cotton

Nitrogen dilution curves of wheat and cotton are shown in Fig. 6a and b. Wheat showed an ongoing decline in leaf nitrogen content with increase in total dry weight. The dilution was less rapid in the monocrop of wheat than in intercrops, possibly resulting from a lesser dry weight increase as a result of a stronger competition for light compared to wheat grown in intercrop strips, where there is sideways incidence of light. Leaf nitrogen content in cotton decreased rapidly while the crop progressed from the seedling stage to an above-ground dry mass of 200 g m−2, but it remained relatively constant from then on. The level at which N content leveled off was higher in the cotton than in the wheat, and the dilution was somewhat less in the cotton plants of the 3:2 system than of the other systems.

Fig. 6
figure 6

Relationship between above-ground dry mass and nitrogen content for wheat (a) and cotton (b) in the intercrops and the monocultures in 2002

Trends in specific leaf nitrogen content

The seasonal dynamics in crop nitrogen content of intercropped and monocrop cotton were reflected in the specific leaf nitrogen content (SLN); three phases are distinguished (Fig. 7). The first phase of cotton seedling growth; because of retarded growth, SLN at 50 DAS (time of wheat harvest) was significantly lower in intercrops than in sole cotton. During the second phase, from 50 to 75 DAS, the SLN of the intercrops recovered, exceeding the SLN of the monoculture. During the third phase, from DAS 75 to cotton harvest, the SLN of intercrops stayed higher than in the monocrop, because more nitrogen was accumulated in the vegetative organs of intercropped cotton and crop senescence was delayed. The delayed senescence of leaves yielded some increase of dry mass, but did not compensate for the delay in growth during the early stages. These findings show that intercropped cotton accumulates more N per unit of biomass than sole cotton.

Fig. 7
figure 7

Specific leaf nitrogen content (SLN) of cotton before anthesis in the intercropping systems and the monoculture in 2002

Discussion and conclusions

Nitrogen use at the crop level

N yield of wheat ranged from 130 to 205 kg ha−1 in intercrops, which is 62–84% of the monoculture. The total N uptake of wheat ranged from 203 to 288 kg ha−1 in the monoculture which corresponds to the reported 257 kg ha−1 N uptake of wheat at about 8 t ha−1 grain yield, averaged over five regions in China (Liu et al. 2006).

The values for total N uptake of sole cotton, ranging from 110 to 127 kg ha−1 in 2002, 2003 and 2004, were within the range (67–403 kg ha−1) reported by Rochester (2007), but lower than found by Ishaq et al. (2001) and Sainju et al. (2006). Total N uptake did not differ from the monoculture in any intercropping system in 2002; however, in 2003 similar N uptake as in monoculture was only found in intercropping systems 3:2 and 4:2 and a lower uptake in the 3:1 and 6:2 systems. In intercropping system 6:2 the crop failed to close the canopy and therefore light interception and crop growth were reduced; as a consequence the N yield in 2003 and 2004 was significantly less (48%) than in the monoculture.

The different responses between seasons and within intercropping systems are strongly associated with the degree of competition during the cotton seedling stage and with canopy development after the wheat harvest when the full space can be used by the cotton crop. The nitrogen uptake of intercropped cotton was only 15–45% of the sole cotton at the time of the wheat harvest. After the wheat harvest, the nitrogen uptake in the intercropping systems recovered. These findings suggest that during the reproductive phase, the delayed fruiting and boll formation resulted in a weaker sink, and therefore more assimilates were retained in the vegetative parts. As a consequence the leaves could stay green for a longer period of time.

The physiological nitrogen use efficiency (IE) of wheat was not much affected by intercropping, but the N and grain yield per unit homogenized land area were strongly associated with the plant densities. Internal efficiencies (IE) of wheat in intercrops and the monocrop were similar except in the 3:2 system that showed higher values; IE of wheat ranged from 26.5 to 29.0 kg kg−1 on average over 3 years. IE was higher in 2004 due to a higher harvest index (HI). A higher HI in intercrops may be explained from a border row effect as a result of border plants receiving more light. Liu et al. (2006) reported that IE of wheat ranged from 19.8 to 66.4 kg kg−1 in China based on estimates by the QUEFTS model. Our results were within the range of 25.5–30.5 kg kg−1.

Based on lint yields IE of cotton ranged from 6.4 to 8.2 kg lint per kg N uptake in intercrops in 3 years, which was significantly lower (P < 0.05) than in monoculture (9.7 kg. kg−1) except in the 6:2 system. Compared to the monoculture, the lower IE of intercropped cotton was due to a lower HI (18%) as a result of the delay in growth and development during the early growth stages (Zhang et al. 2007b). The variation in IE of intercropped cotton under various environmental condition (years) indicates that genotype × environment interactions may play an important role. Thus, a more detailed analysis of the N dynamics of various genotypes of cotton under contrasting environmental conditions is needed.

Nitrogen dynamics at the plant level

Plant N uptake is co-regulated by soil N supply and shoot growth, as concluded by Lemaire et al. (2007) based on studies under different environments (temperate and subtropical) with various crop species (C3 and C4). The relationship between N uptake and biomass accumulation reflects the feed-back regulation of N absorption capacity of roots by shoot growth itself under non-limiting N supply. During the intercropping period, the nitrogen content of cotton seedlings was significantly lower than in the monoculture. We conclude that intercropping decreased N uptake of cotton seedlings strongly during the intercropping period. Cotton plants compensate with extra N uptake at later growth stages for a delay in development as a consequence of shading in the seedling phase.

The coefficient b c of the N dilution curve in sole wheat was 0.48, which was consistent with the study of Flenet et al. (2006). The N demand per unit biomass of cotton in intercrops was greater than in the monoculture. The specific leaf nitrogen content (SLN) of cotton in intercrops ranged from 2.0 to 2.3 g m−2. The relationship between SLN and radiation use efficiency (RUE) for cotton is consistent with other species (Milroy and Bange 2003). In cotton a linear increase of RUE was found for leaf N contents ranging from 2 to 5% (Sadras 1996). Milroy and Bange (2003) reported also that leaf assimilate rates were about 3 times higher at SLN-values within a range from 2.78 to 4.32 g m−2 compared to 1.45 g m−2. This corresponds with the findings of Reddall et al. (2004); they reported that leaf assimilate rates rapidly increased within a range from 0.9 to 2.5 g m−2. Thus the lower SLN values in the 3:1 system are likely to have resulted in a reduced light use efficiency. The recovery of the nitrogen content in cotton biomass after the intercropping period, is in good agreement with compensatory effects documented for cotton after loss of reproductive organs (Sadras 1995). Our results confirm the compensatory effect for cotton with a delay in development in wheat-cotton intercropping systems.

Relative N yields and N balance at the system level

We found that the relative nitrogen yield total (RNT) in the intercropping systems ranged from 1.4 to 1.7; only the RNT of the 6:2 system was significantly lower. Compared to the land equivalent ratios (LER), which ranges on average from 1.3 to 1.4 in intercropping systems (Zhang et al. 2007a), the RNT is 8–21% higher. It was concluded that cotton wheat relay intercropping systems utilize N less efficiently than the monocrops at the system level.

In our experiment the N-doses applied in the intercropping systems (range: 302–412 kg ha−1) and in the monocultures (for wheat: 124–211 kg ha−1 and for cotton 121–187 kg ha−1) were less than in farmers’ practice. The conventional N-fertilizer doses in farmers’ practice in the Yellow River region amounted to 375 kg ha−1 in wheat and 360 kg ha−1 in cotton (Zhen et al. 2006). The soil indigenous N estimated in this study ranged from 115 to 141 kg ha−1 for wheat, which corresponds quite well with the 135 kg ha−1 reported for the same region (Liu et al. 2006). The N-excess (Table 6) derived from the N balance sheet – calculated as N input (applied plus soil indigenous N) minus N removal by wheat grain and seed cotton – ranged from 400 to 410 kg ha−1 in the intercropping systems, which was much higher than in sole cotton (147 kg ha−1) and wheat (205 kg ha−1). Thus, more N in the intercropping systems was prone to losses by leaching or other processes than in the sole cotton system. In both systems, excessive N is being applied in practice, resulting in unnecessarily high input costs and environmental pollution.

Table 6 N balance sheet including the N residues and N excess in the intercropping and monocropping systems in 2 years (Exp. 1)

The total N content of the soil in a 60 cm soil profile increased from 0.82–0.89 g kg−1 in Nov. 2001 (the beginning of experiment 1) to 0.91–0.96 g kg−1 in Dec. 2004 at the end of 3-years experimental period. This corresponds with an average increase in total soil N of about 600 kg ha−1 over a period of 3 years. So, a substantial part of the excess N accumulated in the soil. The increase of soil N might be due to the full return of residues and immobilization of fertilizer N. We concluded that the wheat-cotton intercropping system enriched soil N-status more than a sole cotton or wheat crop. Based on these findings, the N management of the intercropping systems as well as of monocrops should be improved by means of proper timing and demand-oriented dosing of N fertilization. Simple and robust crop simulation models and GIS can help to improve the decision support to farmers in developing more profitable and sustainable cotton cropping systems (McKinion et al. 2001).