Leveraging cover crop functional traits and nitrogen synchronization for yield-optimized waxy corn production systems

Context: Prolonged monoculture of waxy corn (Zea mays L. var. ceratina Kulesh) exacerbates soil nutrient depletion and compromises soil structural integrity, concomitant with underutilization of photosynthetically active radiation (PAR) resources.

Objective: Implementing cover cropping post-harvest of waxy corn can utilize residual environmental resources for soil quality improvement. Nevertheless, the agronomic consequences of this practice on canopy architecture optimization and resource allocation dynamics in subsequent growing seasons require systematic elucidation.

Methods: The regionally adapted cultivar Wannuo 2000 (W67×W68), predominantly cultivated in northeastern China, was employed for canopy characterization. Three experimental treatments were established, including waxy corn continuous monoculture control (CK), shamrock (Trifolium pratense L.) rotation cropping after waxy corn harvest (ZT) and rapeseed (Brassica napus L.) rotation cropping after waxy corn harvest (ZB). Each treatment incorporated five nitrogen (N) application gradients (0, 60, 120, 180, 240 kg N ha⁻¹) arranged in randomized complete block design (RCBD) with triplicate plots.

Results: Cover crops integration significantly enhanced waxy corn productivity. Mean yields for ZT and ZB systems during 2022–2023 demonstrated 20.74% (8.88 ± 2.50 Mg ha⁻¹) and 22.26% (8.99 ± 3.12 Mg ha⁻¹) increases respectively compared to CK. Remarkably, under 25% N reduction scenarios, ZT and ZB achieved 15.25% (44.58 ± 6.28%) and 20.67% (46.68 ± 7.15%) improvements in nitrogen use efficiency (NUE) relative to conventional practice. The path analysis revealed synergistic interactions between cover crop incorporation and N management mediated through canopy structural optimization. Specifically, enhanced leaf area index (4.56 ± 0.69 m² m⁻²) and elevated pre-silking canopy N content (132.61 ± 26.33 g N plant⁻¹) collectively drove post-silking biomass accumulation (134.88 ± 26.85 g plant⁻¹) and N remobilization efficiency.

Conclusion: The integrated cover crop-nitrogen reduction system enhanced both yield and NUE relative to monoculture benchmarks, demonstrating dual benefits in environmental conservation and agricultural productivity enhancement.

Implication: This study establishes a theoretical framework and provides empirical evidence for the sustainable intensification of waxy corn production systems.

1 Introduction

In developing regions where food security remains precarious, reconciling agricultural productivity with ecological sustainability constitutes a pressing scientific challenge (Sanchez et al., 2020; Najafi Vafa et al., 2022). As a dual-purpose crop serving both dietary and economic functions, waxy corn (Zea mays L. var. ceratina Kulesh) occupies critical positions in food systems, though its yield intensification typically demands escalated agrochemical inputs and labor investments (Feng et al., 2024). The concomitant environmental externalities of nutrient oversupply, particularly nitrogen (N) leaching and phosphorous runoff, progressively compromise soil health parameters and aquatic ecosystem integrity, ultimately diminishing the economic viability and agroecosystem resilience of waxy corn cultivation (Motasim et al., 2024). This dichotomy necessitates innovative cultivation strategies that reconcile yield enhancement with ecological preservation in waxy corn production systems.

Proper N management exerts dual effects on pedospheric systems, demonstrating the ability to enhance soil nutrient availability and stable carbon pool formation. However, excessive application accelerates SOC mineralization, reducing carbon sequestration potential while compromising yield stability (Chen et al., 2009). The environmental ramifications extend beyond terrestrial ecosystems, as surplus N undergoes biogeochemical transformations into nitrous oxide (N₂O), a potent greenhouse gas contributing to climate change (Winnick, 2021; Cui et al., 2024). Sustainable agricultural intensification requires three fundamental approaches: 1) Regional optimization of N application thresholds through precision agricultural techniques, 2) Development of crop cultivars with enhanced N assimilation efficiency, and 3) Implementation of integrated nutrient management systems (Wang et al., 2020; Duan et al., 2024a). Recent investigations further highlight the critical role of legacy N from optimized fertilization regimes in sustaining subsequent crop N demand, thereby improving systemic nitrogen use efficiency (Zhao et al., 2023).

The Soil Science Society of America formally defines cover crops as non-commodity vegetation strategically implemented to maintain continuous soil coverage during agricultural fallow periods. These biological systems provide multidimensional agroecological benefits through canopy closure, organic matter deposition, and rhizosphere engineering, with leguminous species exhibiting particular efficacy in atmospheric N assimilation via symbiotic rhizobia associations (Hubbard et al., 2013). Beyond N enrichment, cover crop systems confer critical ecosystem services including erosion mitigation through aggregate stabilization, hydrological regulation via improved infiltration capacity, pathogen suppression through allelochemical production, and carbon sequestration in stabilized organic fractions (Seitz et al., 2024). The integrated effects of enhanced microbial biodiversity and optimized soil physicochemical properties collectively contribute to agricultural system resilience. Current evidence substantiates that incorporating cover cropping aligns with global sustainability frameworks by addressing soil conservation, climate change mitigation, and yield stability imperatives (Purwanto and Alam, 2019).

Solar radiation constitutes a primary biophysical driver governing plant developmental process, with direct regulatory effects on photosynthetic efficiency, canopy water-use dynamics, and biomass partitioning mechanisms (Durand et al., 2021; Wang et al., 2022). In maize cultivation systems, optimized canopy architecture development enhances photosynthetically active radiation (PAR) interception, subsequently driving carbohydrate assimilation and yield formation (Soto-Gómez and Pérez-Rodríguez, 2022). The practice of maize-legume intercropping demonstrates particular efficacy in this context, where symbiotic interactions stimulate maize population establishment while maintaining radiation capture efficiency. However, suboptimal light environments may induce photomorphogenic adaptations affecting leaf angle distribution and stomatal patterning, necessitating strategic canopy management to balance light harvesting with metabolic expenditure (Yang et al., 2021). Contemporary agronomic research consequently positions canopy structure optimization as a cornerstone strategy for yield intensification, particularly through manipulation of leaf area index (LAI) and vertical foliage arrangement to maximize grain yield under fluctuating radiation regimes (Croce et al., 2024).

Waxy corn production systems confront dual agronomic constraints: suboptimal yield potential and inefficient N utilization, issues that concurrently compromise economic viability and exacerbate environmental externalities (Lu et al., 2024; Elrys et al., 2023). Addressing this sustainability paradox requires innovative cultivation strategies that synergistically enhance productivity while optimizing N economy—a critical nexus for ensuring both food security and ecological stewardship. The post-harvest phase in northern China’s waxy corn systems presents untapped potential, with residual photosynthetic active radiation and accumulated thermal units providing agronomic windows for cover crop integration. This investigation systematically elucidates the mechanistic relationships between cover cropping practices and subsequent crop performance through a comprehensive analysis of N remobilization dynamics and photo assimilate partitioning patterns, ultimately establishing an evidence base for sustainable intensification pathways.

2 Materials and methods

2.1 Experimental location description

The two-year field experiment (2022-2023) was conducted at Gongzhuling Experimental Station (43°52’N, 124°81’E; 206 m asl) within the Northeast China Corn Belt, characterized by a temperate continental monsoon climate with 2,744.32 annual sunshine hours and 4,262.71 MJ m⁻² solar radiation. Daily air temperature and precipitation during the experimental years are presented in Figure 1. Significant interannual precipitation variation was observed, with growing-season totals measuring 598.20 mm in 2022 and 339.80 mm (−43.20%) in 2023, inducing moderate drought conditions. Before the trial initiation, the soil contained 25.22 g kg⁻¹ of organic matter, 78.67 mg kg⁻¹ of alkali N, 197.67 mg kg⁻¹ of available K, and 22.18 mg kg⁻¹ of available P. Comparative soil analysis (0–20 cm depth) between monoculture (CK) and cover crop-integrated systems (ZT&ZB) revealed progressive soil quality enhancement: organic matter content increased by 5.00% (2021) and 9.76% (2022), accompanied by respective 5.88%/18.06%, 8.20%/19.75%, and 16.54%/28.00% improvements in alkaline hydrolysable nitrogen, available potassium, and phosphorus concentrations (Table 1). These findings demonstrate the temporal accumulation effect of cover crops on soil fertility parameters.

Figure 1

Figure 1. The mean daily maximum temperature (Tmax), mean daily minimum temperature (Tmin), mean average temperature (Tavg), precipitation (Pre) during waxy corn growing season in 2022 and 2023.

Table 1

Table 1. Chemical properties of test site soil at 0–20 cm depth. IR, the average increase ratio of ZT and ZB relative to CK.

2.2 Experimental design and agronomic management

This study is a long-term positioned experiment initiated in 2021, therefore, the location of the plots for each treatment has remained unchanged between different years. In the preceding years, the field had been under continuous grain maize cultivation with full-season fertilizer applications of 210 kg N ha⁻¹, 90 kg P ha⁻¹, and 100 kg K ha⁻¹. The study employed a split-plot randomized complete block design with triplicate replications, incorporating three agronomic variables: planting systems (CK: post-harvest fallow; ZT: Trifolium pratense L. rotation cropping; ZB: Brassica napus L. rotation cropping) and N application gradients (0, 60, 120, 180, 240 kg N ha⁻¹). Fifteen treatment combinations were systematically arranged in 38.4 m² experimental plots (5 m × 7.8 m) with 1.5 m isolation buffers, utilizing the regionally adapted waxy corn cultivar ‘Wannuo 2000’ (W67×W68) predominant in Northeast China’s agricultural production systems.

The agronomic protocol implemented a rotational cover crop integration system, cover crops had been planted after the waxy corn harvest in the previous year. Waxy corn was cultivated under wide-narrow row configuration (0.9 m/0.4 m row spacing) during 5 May - 19 August 2022 and 27 April - 17 August 2023 growing windows, as planted at a density of 70,000 plants ha^–1. The cover crops were sown on August 20, 2022, and August 18, 2023, with shamrock planted at a seeding rate of 110 kg ha⁻¹ and rapeseed (cultivar Longyou 7) at a planting density of 330,000 plants ha⁻¹. Before planting waxy corn, a rotary tiller was used to incorporate the cover crops into the soil to a depth of 15 cm, and then plant the waxy corn above the incorporated cover crop. After harvesting the waxy corn, cover crops were planted in wide rows between the waxy corn stubble to avoid root interference (Figure 2). Nutrient management comprised basal application (V3 stage) of 90 kg P₂O₅ ha⁻¹, 90 kg K₂O ha⁻¹, and 50% total N allocation, with the remaining N top-dressed at the V12 stage through precision banding. Hydrological inputs relied exclusively on precipitation events, complemented by integrated pest management strategies to maintain phytosanitary conditions.

Figure 2

Figure 2. Schematic diagram of planting methods for corn and cover crops.

2.3 Sampling and measurements

Yield quantification was conducted at the fresh food period (R3) through destructive sampling of central plot areas (3 m × 3 row section). A representative sampling protocol was implemented, selecting five ears per plot for multivariate analysis: 1000-kernel mass (g), kernel count per ear, ear density (ears ha⁻¹), and grain yield (Mg ha⁻¹ at 14% moisture content).

Plant biomass dynamics were assessed through phased destructive sampling at critical phenological stages: silking (R1) and fresh food period (R3). Triplicate plant specimens per plot were subjected to desiccation (85°C forced-air oven, 48-hour duration) until mass stabilization (± 0.01 g precision). The homogenized tissue samples (≤ 0.50 mm particle size) underwent N quantification via micro-Kjeldahl digestion (CN61 M/KDY-9820, Beijing, China) following AOAC 978.02 protocols. Derived parameters included:

$$\text{Post-silking N uptake amount}=\text{R}3\text{ N content}–\text{R}1\text{ N content}$$

$$\text{Post-silking N uptake percentage in total N}=\text{Post-silking N uptake amount}/\text{R}3\text{ N content}$$

$$\text{N remobilization}=\text{R}1\text{ vegetative organ N content}–\text{R}3\text{ vegetative organ N content}$$

$$\text{Nitrogen use efficiency }(\text{NUE})=(\text{R}3\text{ N content in N fertilized cropland}–\text{R}3\text{ N content in N unfertilized cropland})/\text{Total N input to cropland}$$

$$\text{Post}-\text{silking dry matter }(\text{DM})\text{ weight}=\text{R}3\text{ DM}–\text{R}1\text{ DM}$$

$$\text{Harvest index }(\text{HI})=\text{Grain dry weight}/\text{Aboveground biomass}$$

Canopy radiation dynamics were quantified during the critical silking phase (R1 stage) in both experimental years using a SunScan canopy analysis system (Delta-T Devices Ltd., Cambridge, UK). High-resolution photosynthetically active radiation (PAR) measurements were conducted between 10:00-12:00 solar time under clear-sky conditions (cloud cover< 20%), capturing incident radiation above the canopy and transmitted PAR (TPAR) at four critical canopy strata: ground level, third leaf below ear node, ear node position, and third leaf above ear node. Intercepted PAR (IPAR) was derived through radiative transfer modeling of these vertical attenuation profiles.

Concurrent leaf area characterization employed a non-destructive sampling protocol, measuring length-width dimensions (with 0.75 correction factor) on five representative plants per plot. Leaf area index (LAI) was calculated as the ratio of total foliar surface area (m² plant⁻¹) to projected ground occupancy (m² plant⁻¹). The canopy light extinction coefficient (K) was determined through Beer-Lambert law application:

$$K=-\frac{1}{\text{LAI}}\text{ln}(\frac{{\text{PAR}}_{\text{transmiteed}}}{{\text{PAR}}_{\text{incident}}})$$

Where PAR_transmiteed is the IPAR of waxy corn canopy, and PAR_incident is the TPAR.

2.4 Statistical analysis

Statistical analyses were conducted by using Origin 2022 (OriginLab, Northampton, Massachusetts, USA). The data sources were analyzed through a three-way ANOVA followed by the LSD test to compare the differences among the groups. Pearson correlation coefficients were calculated using SPSS to identify interrelationships among the measured parameters, and path correlation analyses were conducted to better understand causal relationships.

3 Results

3.1 Grain yield and yield components

Agronomic interventions significantly influenced waxy corn productivity through both main and interactive effects (p< 0.01), as detailed in Table 2. The brassica-based rotation cropping system (ZB) under N₁₈₀ fertilization demonstrated peak yields of 12.18 Mg ha⁻¹ (2022) and 11.96 Mg ha⁻¹ (2023), surpassing conventional monoculture (CK) by 22.26% across both seasons. Yield component analysis revealed systemic improvements in ZT and ZB systems: 16.02% higher ear density, 2.71% elevated kernel count per ear, and 14.44% increased 1000-kernel mass relative to CK.

Table 2

Table 2. Yield and yield components of different treatments under different N conditions in 2022 and 2023.

N response curves at CK showed yield increases slowly with N escalation (34.72 kg yield gain per kg N in 2022; 31.73 kg in 2023). Maximum productivity in ZT and ZB occurred at N₁₈₀, with marginal yield gains of 49.12, 53.83 kg per kg N (2022) and 40.22, 61.68 kg (2023) respectively. Quadratic plateau modeling identified optimal N thresholds of 189.25 kg ha⁻¹ for CK, 155.75 kg for ZT, and 172.80 kg for ZB (Figure 3), demonstrating cover crops’ nitrogen-sparing capacity.

Figure 3

Figure 3. The impact of N fertilizer application amount on grain yield and nitrogen use efficiency (NUE) under different cover crop treatments. **Significant at p ≤ 0.01.

3.2 N accumulation, post-silking N uptake, and N remobilization

N dynamics analysis revealed significant treatment effects (p< 0.01) on pre-silking N allocation, post-silking assimilation efficiency, and system-level N utilization (Table 3). Both N application and cover crop integration demonstrated interactive effects across measured parameters. Maximum nitrogen use efficiency (NUE) was achieved at N₂₄₀ for monoculture (CK: 41.24%), while cover crop systems (ZT: 49.89%; ZB: 54.28%) peaked at reduced N levels (N₁₈₀), representing 15.25% and 20.67% efficiency improvements over CK respectively.

Table 3

Table 3. N uptake and remobilization under different treatments in 2022 and 2023.

Cover crop treatments enhanced N acquisition capacity throughout the growth cycle, with ZT and ZB systems exhibiting 30.06-40.18% greater plant N content at silking (R1) and 23.41-34.49% elevation at fresh food period (R3) compared to CK. Post-silking N assimilation proved particularly responsive to cover crop integration, showing 54.55-61.99% increases in absolute uptake (kg N ha⁻¹), and 5.21-26.06% increases proportional contribution to total N under N60-N180 regimes. Notably, this advantage diminished at excessive N levels (N₂₄₀), where soil mineral N saturation likely negated cover crops’ rhizosphere modification benefits.

Mechanistically, Correlation analysis identified post-silking N assimilation (correlation coefficient 0.76) and R3 N status (0.67) as primary NUE determinants (Table 4). The N amplification effect (ΔNUE/N input) for cover crop systems exceeded CK by 6.13-10.40%, demonstrating superior N conversion efficiency through improved soil-plant N synchrony.

Table 4

Table 4. Correlation analysis between N content accumulation or N uptake percentage on total N and yield or yield composition in waxy corn.

3.3 Dry matter accumulation

Interannual biomass partitioning patterns demonstrated significant three-way interactions among growing seasons, N regimes, and cover crop systems (p< 0.01, Table 5). The 43.20% precipitation deficit in 2023 intensified drought stress, amplifying cover crops’ biomass enhancement effects compared to the optimal 2022 conditions. Across both experimental years, cover crop integration systematically improved dry matter (DM) accumulation across growth phases and harvest index (HI), exhibiting N-dependent efficacy gradients. Notably, ZT outperformed ZB under low-moderate N inputs (N₀-N₁₂₀), while this relationship reversed under high N regimes (N₁₈₀-N₂₄₀), suggesting N-mediated interspecific competition dynamics.

Table 5

Table 5. Dry matter production (DM) and harvest index (HI) under different cover crop treatments and various N fertilizers from 2022 and 2023, respectively.

Quantitative analysis revealed, in 2022, the average pre-silking DM of ZT and ZB were 14.19% and 14.97% higher than CK; the average post-silking DM were 28.22% and 37.60% higher than CK; meanwhile, the HI increased 22.33% and 10.68% than CK, respectively. In 2023, the average pre-silking DM of ZT and ZB were 27.48% and 22.47% higher than CK; the average post-silking DM were 50.05% and 64.08% higher than CK; meanwhile, the HI increased 24.26% and 21.32% than CK. Maximum whole-plant biomass (280.08 g plant⁻¹ in 2022; 252.34 g plant⁻¹ in 2023) occurred in ZB under N₁₈₀, mirroring yield optimization patterns. N responsiveness followed quadratic trends, peaking at N₁₈₀ before declining, indicative of metabolic sink limitations under excessive fertilization.

Biometric correlations established strong positive associations between cumulative DM and yield parameters (r = 0.57-0.87, p< 0.01), whereas HI showed no significant relationship (r = 0.11, p > 0.05), confirming biomass accumulation rather than partitioning efficiency as the primary yield determinant in these systems (Table 6). The differential nitrogen × cover crop × environment interactions highlight the necessity for adaptive management strategies balancing N inputs with cover crop functional types under variable precipitation regimes.

Table 6

Table 6. Correlation analysis between dry matter (DM) accumulation or harvest index (HI) and yield or yield composition in waxy corn. .

3.4 Silking and maturity leaf area index and extinction coefficient

Canopy architectural analysis revealed systematic improvements in leaf area index (LAI) under cover crop integration, with ZT and ZB systems exhibiting 17.81 and 17.76% superior LAI enhancement compared to monoculture (CK) across critical growth stages (Figure 4). In the optimal 2022 season, maximum LAI values of 5.58 (ZT) and 5.62 (ZB) at silking under N180 fertilization significantly outperformed CK by 18.54% and 15.66% respectively, demonstrating N-mediated canopy expansion. This advantage persisted through maturity, with ZT and ZB maintaining 17.08 and 19.89% greater LAI, indicative of prolonged photosynthetic capacity.

Figure 4

Figure 4. LAIs of waxy corn under different N fertilizer and cover crop treatments at silking and maturity in 2022 and 2023. *Significant at p ≤ 0.05; **Significant at p ≤ 0.01; ***Significant at p ≤ 0.001; ns, non‐significant (LSD).

The drought-affected 2023 season modified these relationships - while low (N₀) and excessive (N₂₄₀) N inputs negated LAI differences between systems, moderate fertilization (N₆₀-N₁₈₀) sustained ZT/ZB’s canopy superiority (15.59-15.89% over CK). Mechanistically, LAI dynamics directly regulated light extinction coefficients (k) through Beer-Lambert law interactions, where k increased logarithmically with LAI escalation. This light attenuation pattern drove biomass accumulation, evidenced by strong k correlations with pre-silking (R²=0.61) and post-silking (R²=0.46) dry matter production (Figure 5). The synergistic LAI-k relationships underscore cover crops’ dual role in optimizing both light interception efficiency and N utilization precision, particularly under 180 kg N ha⁻¹ fertilization where canopy architecture and radiation use efficiency reached equilibrium.

Figure 5

Figure 5. The correlation between extinction coefficient (k) and pre-silking DM (A), post-silking DM (B), pre-silking N content (C), and post-silking N content (D), respectively.

3.5 Path analysis

Factors affecting the yield and NUE, and the relationship between these causes and effects, are shown in Figure 6. The path analysis revealed that leaf area index (LAI) and extinction coefficient (k) exerted significant direct effects on pre-silking dry matter (DM), with standardized path coefficients (β) of 0.92 and 0.39, respectively. Pre-silking DM subsequently demonstrated substantial direct effects on three yield components: kernel number (β = 0.86), post-silking DM (β = 0.96), and 1000-kernel weight (β = 0.57). These yield components directly influenced final grain yield through distinct pathways, with kernel number showing the strongest association (β = 0.46), followed by post-silking DM (β = 0.42) and 1000-kernel weight (β = 0.19). Regarding nitrogen dynamics, LAI and k significantly affected pre-silking N content through direct pathways (β = 0.63 and 0.37, respectively). Pre-silking N content subsequently influenced three critical nitrogen-related parameters: post-silking DM (β = 0.29), post-silking N content (β = 0.64), and N remobilization (β = 0.61). The NUE was primarily determined by post-silking nitrogen content (β = 0.55) and nitrogen remobilization (β = 0.38), demonstrating their differential contributions to nitrogen utilization.

Figure 6

Figure 6. Path analysis of the relationship among yield, kernel number, 1000-kernel weight, pre-silking dry matter (DM), post-silking DM, leaf area index and extinction coefficient (A), and the relationship among nitrogen use efficiency, post-silking DM, pre-silking N content, post-silking N content, N remobilization, leaf area index and extinction coefficient (B). Arrows point from causes and effects, and solid lines represent direct effects. The values along of the arrows are r of the correlation with ** being p ≤ 0.01.

The results showed further evidence that the interplay of cover crop and N significantly influenced yield and NUE by canopy structure. Canopy architectural improvements mediated by cover crop integration fundamentally enhanced radiation capture efficiency and N utilization dynamics across growth stages. The synergistic interaction between cover crop systems and N management amplified leaf area development, with cover crop treatments consistently outperforming monoculture in both optimal and drought-affected seasons. Enhanced leaf area indices prolonged photosynthetic activity duration while optimizing light distribution patterns within the canopy, driving coordinated improvements in pre-silking resource accumulation and post-silking translocation efficiency.

4 Discussion

Contemporary agroecological research substantiates cover crops as multifunctional biological regulators within agricultural ecosystems, demonstrating tripartite benefits through soil nutrient modulation, erosion mitigation, and systemic soil health enhancement (Sharma et al., 2017; Tribouillois et al., 2022; Goettl et al., 2024). These living mulches operate through three synergistic mechanisms: 1) temporal niche complementarity in nutrient acquisition, 2) rhizosphere engineering via root exudate-mediated microbial consortia, and 3) physical armor against erosional forces through canopy-litter composite layers. Particularly in maize production systems, cover crops exhibit dual-phase optimization—enhancing vegetative growth through microclimate moderation while concurrently implementing nutrient sequestration mechanisms that reduce fertilizer dependency (Sainju et al., 2005; Muñoz et al., 2014). The experimental outcomes substantiate the resilience-enhancing capacity of optimized cover crop systems under climatic variability. Across contrasting growing seasons, cover crop integration with precision N management demonstrated sustained yield superiority over conventional monoculture, particularly under drought-stressed conditions (Table 2). While interannual precipitation fluctuations induced significant reductions in biomass accumulation and N assimilation, the ZB system maintained robust productivity—exhibiting 34.00-39.00% yield advantage under optimal N application (N₁₈₀) compared to control plots. This drought-mitigation efficacy underscores the system’s ability to buffer environmental stressors through enhanced soil moisture retention and N cycling efficiency.

4.1 Effects of cover crops and N application on NUE in waxy corn

The differential performance observed between normal and stress years highlights a critical synergy between cover crop functional traits and N availability to stabilize yield outcomes under hydrological uncertainty. The observed yield enhancement mechanisms align with established principles of ecological intensification, where cover crops functionally compensate for soil N limitation through species-specific strategies. Under N-deficient regimes (N₀-N₁₂₀), the leguminous cover crop (Trifolium pratense L.) likely demonstrated superior efficacy through symbiotic N fixation, effectively bypassing soil N constraints and elevating plant-available N pools during critical growth phases. This biological N provisioning likely translated directly into enhanced kernel initiation and grain filling efficiency. Conversely, brassica cover crops (Brassica napus L.) exhibited maximum performance synergies under optimized N inputs (N₁₈₀-N₂₄₀), where their extensive biomass production likely stimulated microbial-mediated N mineralization and improved edaphic conditions through enhanced aggregate stability. The N-dependent functional complementarity between cover crop types underscores the importance of matching species selection with fertilization regimes to maximize system-level productivity. These findings corroborate global evidence of cover crops’ capacity to mitigate nutrient limitations while amplifying yield components through both direct and indirect pathways (Gabriel and Quemada, 2011; Radicetti et al., 2016). Mechanistic analysis revealed a N-dependent yield optimization threshold governed by Liebig’s law of the minimum, where yield and NUE exhibited quadratic response patterns to fertilization intensity (Figure 3). This relationship delineates an economic-environmental equilibrium threshold—the precise N input level that simultaneously maximizes productivity and minimizes agrochemical surplus. Beyond this critical threshold (N₁₈₀ in cover crop systems), N saturation likely induced metabolic constraints through luxury consumption, triggering assimilates partitioning imbalances that preferentially favored vegetative over reproductive development. The establishment of such N-mediated yield ceilings underscores the imperative for precision N management protocols that harmonize crop physiological demands with ecological carrying capacities, achieving dual optimization of productivity and sustainability in waxy corn production systems.

Strategic N management through agroecological intensification has emerged as a pivotal approach for reconciling agricultural productivity with environmental stewardship (Abdo et al., 2024; Zou et al., 2024). Cover crops function as dynamic biofilters within this paradigm, mediating soil N cycles through three synergistic mechanisms: rhizosphere immobilization of residual N, microbial reactivation of legacy nutrients, and physical-chemical stabilization of soil organic matter. By establishing continuous biological activity during fallow periods, these systems synchronize N mineralization pulses with subsequent crop demand while suppressing nitrification-denitrification cascades that drive greenhouse gas emissions (Baptistella et al., 2020; Hirsh et al., 2021). The resultant enhancement of N residence time and plant availability fundamentally reconfigures traditional input-output relationships, enabling yield stability under reduced fertilization regimes. This ecological engineering strategy exemplifies the transition from linear nutrient management to circular nutrient economies in modern agroecosystems. This study demonstrates that cover crop integration significantly improves NUE across fertilization regimes, with synergistic nitrogen × cover crop interactions driving non-linear response patterns (Table 3). The temporal alignment between enhanced N accumulation during silking—a critical phase for grain differentiation and yield formation (Jacobs and Pearson, 1992)—and cover crop-mediated N provisioning highlights a fundamental physiological linkage. This N assimilation advantage during early reproductive development establishes a cascade effect, optimizing subsequent nutrient partitioning efficiency while minimizing environmental losses through improved N demand-supply synchrony. The diminished NUE enhancement at extreme N levels (N₆₀-N₂₄₀) underscores the importance of balanced fertilization strategies that harmonize agroecological intensification with crop developmental phenology. During critical reproductive phases, the developing grain sink exerts a hierarchical pull-on N resource (Zhang et al., 2016), amplifying translocation efficiency from vegetative tissues—a process modulated by cover crop functional traits. Under moderate N regimes (N₆₀-N₁₈₀), leguminous and brassica systems demonstrated superior N channeling capacity compared to monoculture, reflecting improved source-sink coordination through enhanced rhizosphere N availability. This synergy diminished under excessive N inputs (N₂₄₀), where soil N saturation likely shifted N cycling toward organic matter stabilization pathways rather than immediate plant accessibility. The N flux redirection suggests a threshold-dependent transition in cover crop functionality: from direct N provisioning under limited conditions to long-term soil N reservoir building in high-input scenarios. The drought-affected 2023 season revealed a paradoxical N assimilation pattern: despite marked reductions in vegetative-stage (R1-R3) N content across treatments, post-silking N acquisition and allocation efficiency exceeded 2022 levels, particularly in cover crop systems. This counterintuitive phenomenon stems from drought-induced rhizosphere modification, where early-season water deficits triggered adaptive root architectural remodeling—characterized by enhanced antioxidant enzyme activity and lignification—that improved hydraulic conductivity and nutrient foraging capacity (Wu et al., 2024).

4.2 Effects of cover crops and N application on dry matter accumulation in waxy corn

Drought-adaptive dry matter partitioning mechanisms in waxy corn were fundamentally reconfigured through cover crop integration, demonstrating a compensatory strategy that prioritizes reproductive sink strength under hydrological constraints. While vegetative biomass accumulation (R1-R3) decreased during the drought-affected season, the system exhibited enhanced assimilate remobilization efficiency and harvest index elevation—a paradoxical response indicating metabolic prioritization of grain filling over vegetative persistence (Table 5). This drought-induced source-sink realignment was amplified in cover crop systems, where ZT and ZB treatments increased post-silking photo assimilate allocation rates compared to monoculture, suggesting improved phloem loading capacity and sucrose synthase activity during critical grain formation phases citation. These findings substantiate cover crops’ role as biostimulants of plant resource optimization, effectively decoupling yield formation from environmental stochasticity through improved assimilate partitioning precision. N management analysis revealed distinct dry matter accumulation patterns between production systems (Zhu et al., 2021). The monoculture control (CK) exhibited a linear increase in biomass accumulation with escalating N inputs, while cover crop systems demonstrated quadratic responses, peaking at N₁₈₀ fertilization (Table 5). This N optimization threshold corresponded with maximum yield realization, indicating synchronized source-sink coordination. Notably, supra-optimal N application (N₂₄₀) induced metabolic sink limitations in cover crop systems, reducing translocation efficiency from vegetative tissues to developing grains. The results substantiate that cover crop integration necessitates N input calibration, with 180 kg N ha⁻¹ establishing the equilibrium point between ecological intensification and agronomic productivity in waxy corn systems (Figure 3).

4.3 Effects of cover crops on canopy structure and light interception capacity in waxy corn

Leaf Area Index (LAI) serves as a critical biophysical determinant of canopy light interception efficiency and photosynthetic optimization in cereal cropping systems (Liu et al., 2022; Duan et al., 2024b). This study demonstrates that cover crop integration fundamentally reconfigures canopy architectural dynamics in subsequent waxy corn, with functional divergence between leguminous (Trifolium pratense L.) and brassica (Brassica napus L.) systems across N gradients (Figure 4). The leguminous cover crop exhibited superior LAI enhancement under N-limited regimes (N₀-N₁₂₀), leveraging symbiotic N fixation to sustain canopy expansion. Conversely, brassica dominance emerged under high N availability (N₁₈₀-N₂₄₀), where its extensive biomass production synergized with soil N mineralization to optimize light distribution patterns (Zhang et al., 2020). These N-mediated canopy modifications establish a self-reinforcing cycle—LAI augmentation facilitates radiation-use efficiency improvement, which in turn sustains photosynthetic capacity during reproductive transitions (Rens et al., 2018). The observed species × N interaction underscores the importance of matching cover crop functional types with fertilization strategies to maximize light capture efficiency while maintaining source-sink balance throughout crop developmental phases (Figure 5). Under N-limited conditions, legume-derived N provisioning sustains canopy expansion, while brassica dominance under high-N regimes promotes multi-layered light distribution through biomass-mediated architectural modifications. Path analysis confirms LAI’s pivotal role as a biophysical lever, simultaneously amplifying both carbon assimilation and N utilization precision through optimized source-sink dynamics (Figure 6). The findings establish canopy architectural modulation as a central strategy for decoupling productivity from environmental constraints in sustainable maize production systems.

This study establishes a novel ecological intensification framework for waxy corn production through functional integration of cover crop-mediated canopy optimization and N synchronization strategies. The demonstrated hierarchical mechanism—spanning from rhizosphere N provisioning to photosynthetic efficiency amplification—reveals how strategic species × fertilization pairings reconfigure source-sink coordination, achieving unprecedented yield-NUE synergies. Crucially, the systems approach resolves the historical dichotomy between input reduction and productivity enhancement by engineering self-sustaining light-nutrient coupling effects. These conceptual advances provide a blueprint for sustainable intensification, where managed biodiversity replaces chemical inputs through ecological leverage points in canopy architecture and nutrient cycling. The findings pioneer a paradigm shift from conventional yield-centric management to whole-system resource optimization, positioning cover crop integration as a keystone strategy for decoupling cereal production from environmental externalities while maintaining economic viability.

5 Conclusion

Cover crops (shamrock and rapeseed) integration significantly enhanced waxy corn productivity, and enables sustainable intensification through light-nutrient synergy and resource optimization. This approach reduces N demand while enhancing yield stability via improved canopy-mediated radiation capture and drought-adaptive N remobilization. The demonstrated coordination of ecological processes—soil nutrient stewardship, hydrological buffering, and carbon-N coupling—provides a scalable template for climate-resilient agriculture. By transforming post-harvest periods into productive ecological phases, the system redefines conventional monoculture paradigms, achieving simultaneous productivity gains and environmental stewardship through managed biodiversity.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

MS: Data curation, Formal Analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. LZ: Data curation, Formal Analysis, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. JZ: Investigation, Writing – review & editing. ZL: Investigation, Writing – review & editing. CP: Investigation, Writing – review & editing. ZJ: Investigation, Writing – review & editing. YL: Conceptualization, Formal Analysis, Resources, Supervision, Validation, Writing – review & editing. YW: Conceptualization, Formal Analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. Jilin Province Key Technology R&D Program (20220302004NC)National Maize Industry System Post Scientist Project (CARS-02-19).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abdo, A. I., Sun, D. L., Yang, K., Li, Y. Z., Shi, Z. J., Abd Allah, W. E., et al. (2024). Carbon footprint of synthetic nitrogen under staple crops: A first cradle-to-grave analysis. Global Change Biol. 30, 1365–2486. doi: 10.1111/gcb.17277

PubMed Abstract | Crossref Full Text | Google Scholar

Baptistella, J. L. C. B., Andrade, S. A. L., Favarin, J. L., and Mazzafera, P. (2020). Urochloa in tropical agroecosystems. Front. Sustain. Food Syst. 4. doi: 10.3389/fsufs.2020.00119

Crossref Full Text | Google Scholar

Chen, H. Q., Hou, R. X., Gong, Y. S., Li, H. W., Fan, M. S., and Kuzyakov, Y. K. (2009). Effects of 11 years of conservation tillage on soil organic matter fractions in wheat monoculture in Loess Plateau of China. Soil Tillage Res. 106, 85–94. doi: 10.1016/j.still.2009.09.009

Crossref Full Text | Google Scholar

Croce, R., Elizabete, C. S., Cho, Y. B., Ermakova, M., Harbinson, J., Lawson, T., et al. (2024). Perspectives on improving photosynthesis to increase crop yield. Plant Cell. 36, 3944–3973. doi: 10.1093/plcell/koae132

PubMed Abstract | Crossref Full Text | Google Scholar

Cui, X. Q., Bo, Y., Adalibieke, W., Winiwarter, W., Zhang, X., and Davidson, E. A. (2024). The global potential for mitigating nitrous oxide emissions from croplands. One Earth 7, 401–420. doi: 10.1016/j.oneear.2024.01.005

Crossref Full Text | Google Scholar

Duan, J. K., Liu, H. B., Zhang, X. M., Ren, C. C., Wang, C., Cheng, L. X., et al. (2024a). Agricultural management practices in China enhance nitrogen sustainability and benefit human health. Nat. Food 5, 378–389. doi: 10.1038/s43016-024-00953-8

PubMed Abstract | Crossref Full Text | Google Scholar

Duan, M., Zhang, X. T., Wei, Z., Chen, X., and Zhang, B. Z. (2024b). Effect of maize canopy structure on light interception and radiation use efficiency at different canopy layers. Agronomy 14, 1–20. doi: 10.3390/agronomy14071511

Crossref Full Text | Google Scholar

Durand, M., Murchie, E. H., Lindfors, A. V., Urban, O., Aphalo, P. J., and Robson, T. M. (2021). Diffuse solar radiation and canopy photosynthesis in a changing environment. Agric. For. Meteorol. 311, 1–13. doi: 10.1016/j.agrformet.2021.108684

Crossref Full Text | Google Scholar

Elrys, A. S., Wang, J., Meng, L., Zhu, Q. L., El-Sawy, M. M., Chen, Z. X., et al. (2023). Integrative knowledge-based nitrogen management practices can provide positive effects on ecosystem nitrogen retention. Nat. Food. 4, 1075–1089. doi: 10.1038/s43016-023-00888-6

PubMed Abstract | Crossref Full Text | Google Scholar

Feng, W. J., Xue, W. W., Zhao, Z. Q., Shi, Z. K., Wang, W. J., Bai, Y., et al. (2024). Nitrogen fertilizer application rate affects the dynamic metabolism of nitrogen and carbohydrates in kernels of waxy maize. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1416397

PubMed Abstract | Crossref Full Text | Google Scholar

Gabriel, J. L. and Quemada, M. (2011). Replacing bare fallow with cover crops in a maize cropping system: Yield, N uptake and fertiliser fate. Eur. J. Agron. 34, 133–143. doi: 10.1016/j.eja.2010.11.006

Crossref Full Text | Google Scholar

Goettl, B., Andersen, B., DeSutter, T., Franzen, D., and Wick, A. (2024). Interseeded cover crops in wide-row corn: An opportunity for northern cropping systems. Crop Forage Turfgrass Manag. 10, 1–7. doi: 10.1002/cft2.20268

Crossref Full Text | Google Scholar

Hirsh, S. M., Duiker, S. W., Graybill, J., Nichols, K., and Weil, R. R. (2021). Scavenging and recycling deep soil nitrogen using cover crops on mid-Atlantic, USA farms. Agric. Ecosyst. Environ. 309, 1–15. doi: 10.1016/j.agee.2020.107274

Crossref Full Text | Google Scholar

Hubbard, R. K., Strickland, T. C., and Phatak, S. (2013). Effects of cover crop systems on soil physical properties and carbon/nitrogen relationships in the coastal plain of southeastern USA. Soil Tillage Res. 126, 276–283. doi: 10.1016/j.still.2012.07.009

Crossref Full Text | Google Scholar

Jacobs, B. C. and Pearson, C. J. (1992). Pre-flowering growth and development of the inflorescences of maize. J. Exp. Botany 43, 565–569. doi: 10.1093/jxb/43.4.557

Crossref Full Text | Google Scholar

Liu, G. Z., Yang, Y. S., Liu, W. M., Guo, X. X., Xie, R. Z., Ming, B., et al. (2022). Optimized canopy structure improves maize grain yield and resource use efficiency. Food Energy Secur. 11, 1–11. doi: 10.1002/fes3.375

Crossref Full Text | Google Scholar

Lu, J., Lankhost, J. A., Stomph, T. J., Schneider, H. M., Chen, Y. L., Mi, G. H., et al. (2024). Root plasticity improves maize nitrogen use when nitrogen is limiting: an analysis using 3D plant modelling. J. Exp. Botany 75, 5989–6005. doi: 10.1093/jxb/erae298

PubMed Abstract | Crossref Full Text | Google Scholar

Motasim, A. M., Samsuri, A. W., Nabayi, A., Akter, A., Haque, M. A., Sukor, A. S.A, et al. (2024). Urea application in soil: processes, losses, and alternatives—a review. Discov. Agric. 2, 1–25. doi: 10.1007/s44279-024-00060-z

Crossref Full Text | Google Scholar

Muñoz, J. D., Steibel, J. P., Snapp, S. L., and Kravchenko, A. N. (2014). Cover crop effect on corn growth and yield as influenced by topography. Agric. Ecosyst. Environ. 189, 229–239. doi: 10.1016/j.agee.2014.03.045

Crossref Full Text | Google Scholar

Najafi Vafa, Z., Sohrabi, Y., Mirzaghaderi, G., and Heidari, G. (2022). Soil microorganisms and seaweed application with supplementary irrigation improved physiological traits and yield of two dryland wheat cultivars. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.855090

PubMed Abstract | Crossref Full Text | Google Scholar

Purwanto, B. H. and Alam, S. (2019). Impact of intensive agricultural management on carbon and nitrogen dynamics in the humid tropics. Soil Sci. Plant Nutr. 66, 50–59. doi: 10.1080/00380768.2019.1705182

Crossref Full Text | Google Scholar

Radicetti, E., Mancinelli, R., Moscetti, R., and Campiglia, E. (2016). Management of winter cover crop residues under different tillage conditions affects nitrogen utilization efficiency and yield of eggplant (Solanum melanogena L.) in Mediterranean environment. Soil Tillage Res. 155, 329–338. doi: 10.1016/j.still.2015.09.004

Crossref Full Text | Google Scholar

Rens, L. R., Zotarelli, L., Rowland, D. L., and Morgan, K. T. (2018). Optimizing nitrogen fertilizer rates and time of application for potatoes under seepage irrigation. Field Crops Res. 215, 49–58. doi: 10.1016/j.fcr.2017.10.004

Crossref Full Text | Google Scholar

Sainju, U. M., Whitehead, W. F., and Singh, B. P. (2005). Carbon accumulation in cotton, sorghum, and underlying soil as influenced by tillage, cover crops, and nitrogen fertilization. Plant Soil 273, 219–234. doi: 10.1007/s11104-004-7611-9

Crossref Full Text | Google Scholar

Sanchez, L., Ermolenkov, A., Biswas, S., Septiningsih, E. M., and Kurouski, D. (2020). Raman spectroscopy enables non-invasive and confirmatory diagnostics of salinity stresses, nitrogen, phosphorus, and potassium deficiencies in rice. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.573321

PubMed Abstract | Crossref Full Text | Google Scholar

Seitz, V. A., McGivern, B. B., Borton, M. A., Chaparro, J. M., Schipanski, M. E., Prenni, J. E., et al. (2024). Cover crop root exudates impact soil microbiome functional trajectories in agricultural soils. Res. Square 12, 1–38. doi: 10.21203/rs.3.rs-3956430/v1

PubMed Abstract | Crossref Full Text | Google Scholar

Sharma, P., Laor, Y., Raviv, M., Medina, S., Saadi, I., Krasnovsky, A., et al. (2017). Green manure as part of organic management cycle: Effects on changes in organic matter characteristics across the soil profile. Geoderma 305, 197–207. doi: 10.1016/j.geoderma.2017.06.003

Crossref Full Text | Google Scholar

Soto-Gómez, D. and Pérez-Rodríguez, P. (2022). Sustainable agriculture through perennial grains: Wheat, rice, maize, and other species. A review. Agric. Ecosyst. Environ. 325, 1–14. doi: 10.1016/j.agee.2021.107747

Crossref Full Text | Google Scholar

Tribouillois, H., Constantin, J., Casal, L., Villerd, J., and Therond, O. (2022). Introducing and expanding cover crops at the watershed scale: Impact on water flows. Agric. Ecosyst. Environ. 337, 1–10. doi: 10.1016/j.agee.2022.108050

Crossref Full Text | Google Scholar

Wang, B., Yue, X., Zhou, H., and Zhu, J. (2022). Impact of diffuse radiation on evapotranspiration and its coupling to carbon fluxes at global FLUXNET sites. Agric. For. Meteorol. 322, 1–9. doi: 10.1016/j.agrformet.2022.109006

Crossref Full Text | Google Scholar

Wang, C., Zang, H. D., Liu, J. G., Shi, X. Y., Li, S., Chen, F., et al. (2020). Optimum nitrogen rate to maintain sustainable potato production and improve nitrogen use efficiency at a regional scale in China. A meta-analysis. Agron. Sustain. Dev. 40, 1–14. doi: 10.1007/s13593-020-00640-5

Crossref Full Text | Google Scholar

Winnick, M. J. (2021). Stream transport and substrate controls on nitrous oxide yields from hyporheic zone denitrification. AGU Adv. 2, 1–17. doi: 10.1029/2021AV000517

Crossref Full Text | Google Scholar

Wu, B., Cui, Z. J., Zechariah, E., Guo, L. Z., Gao, Y. H., Yan, B., et al. (2024). Post-anthesis dry matter and nitrogen accumulation, partitioning, and translocation in maize under different nitrate–ammonium ratios in Northwestern China. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1257882

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Y. S., Guo, X. X., Liu, H. F., Liu, G. Z., Liu, W. M., Ming, B., et al. (2021). The effect of solar radiation change on the maize yield gap from the perspectives of dry matter accumulation and distribution. J. Integr. Agric. 20, 482–493. doi: 10.1016/S2095-3119(20)63581-X

Crossref Full Text | Google Scholar

Zhang, D. S., Sun, Z. X., Feng, L. S., Bai, W., Yang, N., Zhang, Z., et al. (2020). Maize plant density affects yield, growth and source-sink relationship of crops in maize/peanut intercropping. Field Crops Res. 257, 1–10. doi: 10.1016/j.fcr.2020.107926

Crossref Full Text | Google Scholar

Zhang, J. S., Yang, L., Guo, S., Chen, Q. W., Chen, F. J., Yuan, L. X., et al. (2016). Use of the stable nitrogen isotope to reveal the source-sink regulation of nitrogen uptake and remobilization during grain filling phase in maize. PloS One 11, 1–16. doi: 10.1371/journal.pone.0162201

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, X., Cai, S. Y., Yang, B. G., Zhao, H. C., Zeng, K., Fan, P. F., et al. (2023). Soil nitrogen dynamics drive regional variation in nitrogen use efficiency in rice: A multi-scale study. Eur. J. Soil Sci. 74, 1–16. doi: 10.1111/ejss.13352

Crossref Full Text | Google Scholar

Zhu, Y. N., Qi, B. F., Hao, Y. W., Liu, H. C., Sun, G. W., Chen, R. Y., et al. (2021). Appropriate NH₄⁺/NO₃^– ratio triggers plant growth and nutrient uptake of flowering chinese cabbage by optimizing the pH value of nutrient solution. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.656144

PubMed Abstract | Crossref Full Text | Google Scholar

Zou, H. Q., Li, D. J., Ren, K. Y., Liu, L. S., Zhang, W. J., Duan, Y. H., et al. (2024). Response of maize yield and nitrogen recovery efficiency to nitrogen fertilizer application in field with various soil fertility. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1349180

PubMed Abstract | Crossref Full Text | Google Scholar