Next Article in Journal
Arbuscular Mycorrhizal Fungi and Associated Microbiota as Plant Biostimulants: Research Strategies for the Selection of the Best Performing Inocula
Next Article in Special Issue
Identification of Wheat Cultivars for Low Nitrogen Tolerance Using Multivariable Screening Approaches
Previous Article in Journal
Effects of Biochar Application in a Sorghum Crop under Greenhouse Conditions: Growth Parameters and Physicochemical Fertility
Previous Article in Special Issue
Nitrogen Application Improved Photosynthetic Productivity, Chlorophyll Fluorescence, Yield and Yield Components of Two Oat Genotypes under Saline Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In-Field Observation of Root Growth and Nitrogen Uptake Efficiency of Winter Oilseed Rape

by
Julien Louvieaux
1,2,*,
Antoine Leclercq
2,
Loïc Haelterman
1 and
Christian Hermans
1
1
Crop Production and Biostimulation Laboratory, Interfacultary School of Bioengineers, Université libre de Bruxelles, Campus Plaine CP 245, Bd du Triomphe, 1050 Brussels, Belgium
2
Laboratory of Applied Plant Ecophysiology, Haute Ecole Provinciale de Hainaut Condorcet, Centre pour l’Agronomie et l’Agro-industrie de la Province de Hainaut, 11 rue Paul Pastur, 7800 Ath, Belgium
*
Author to whom correspondence should be addressed.
Submission received: 3 December 2019 / Revised: 27 December 2019 / Accepted: 8 January 2020 / Published: 10 January 2020

Abstract

:
Field trials were conducted with two nitrogen applications (0 or 240 kg N ha−1) and three modern cultivars of winter oilseed rape (Brassica napus L.) previously selected from a root morphology screen at a young developmental stage. The purpose is to examine the relationship between root morphology and Nitrogen Uptake Efficiency (NUpE) and to test the predictiveness of some canopy optical indices for seed quality and yield. A tube-rhizotron system was used to incorporate below-ground root growth information. Practically, clear tubes of one meter in length were installed in soil at an angle of 45°. The root development was followed with a camera at key growth stages in autumn (leaf development) and spring (stem elongation and flowering). Autumn was a critical time window to observe the root development, and exploration in deeper horizons (36–48 cm) was faster without any fertilization treatment. Analysis of the rhizotron images was challenging and it was not possible to clearly discriminate between cultivars. Canopy reflectance and leaf optical indices were measured with proximal sensors. The Normalized Difference Vegetation Index (NDVI) was a positive indicator of biomass and seed yield while the Nitrogen Balance Index (NBI) was a positive indicator of above-ground biomass N concentration at flowering and seed N concentration at harvest.

1. Introduction

By the year 2050, a societal challenge will be to almost double the food production from existing land areas to feed more than an estimated nine billion people [1,2]. Innovation has therefore to make a step change to achieve environmentally sustainable intensification of agricultural systems. Crop production heavily relies on chemical fertilizers, but a considerable fraction of these gets lost with harmful consequences to the environment. Nitrogen (N) is the macronutrient required in the largest amount by crops. In agricultural ecosystems, N is continually depleted by such processes as exportation of N-containing crop residues, ammonia volatilization, microbial denitrification, erosion or N leaching. Therefore, the soil N reserve must be replenished every so often to sustain crop yield but N fertilizer in excess is detrimental to the human health (methemoglobinemia) and to the environment (eutrophication, greenhouse gas emission) [3,4]. One way to reduce N fertilizer input is to breed crops with greater Nitrogen Use Efficiency (NUE) [5,6]. NUE has two main components: the Nitrogen Uptake Efficiency (NUpE), which describes the plant capacity to acquire N from the soil, and the Nitrogen Utilization Efficiency (NUtE), the capacity to utilize N to produce harvestable organs [6].
Rapeseed (Brassica napus L.) is the second most important oilseed crop in the world [7]. Winter oilseed rape is regarded as a cash crop that diversifies cereal-dominated rotations and that reduces N leaching during autumn and winter [8]. While most of the efforts have been on proteolysis and N remobilization to the seeds [9,10], the root system is still considered as a “black box.” Nonetheless, the root organ is central for N acquisition [11] and its morphology could be optimized [12]. The influence of nitrate on biomass allocation and root morphology is well documented in the model and parented plant Arabidopsis thaliana. During homogeneous low nitrate conditions, root growth is favored overshoot resulting into increased root to shoot biomass ratio [13], and the elongation of main and lateral roots is triggered [14,15,16]. A rationale shared by several authors is that a branched root system that explores a larger soil volume prevents N leaching [17,18,19]. There are indications that at low N supply, the seed yield of rapeseed cultivars is most closely correlated with N uptake (Nupt) [20] and root growth following stem extension [21]. Other reports indicate that N- efficient cultivars are characterized by an important root density during the vegetative growth stage [22]. Therefore, optimizing the root morphology may be a valuable strategy to improve NUpE. Furthermore, genetic differences in NUpE underpin variation in NUE under limiting N conditions [20,23]. This suggests that NUpE may be a valuable target for the creation of N-efficient genotypes.
In this study, a tube-rhizotron system was employed to characterize the spatiotemporal dynamics of the oilseed rape root system. Our aim is to observe the crop root growth at different key developmental stages and to evaluate such experimental setting for discriminating between N fertilization treatments and genotypes. We also tested whether one cultivar with a branched and deep root system would have a greater NUpE. Complementarily, proximal sensors based on canopy optical properties (foliar reflectance and pigment fluorescence) were evaluated.

2. Material and Methods

2.1. Plant Material

The three varieties of winter oilseed rape (Brassica napus L.) are restored hybrids registered in the French catalogue of plant species and varieties with the following names: DK Exception (DKE), DK Expansion (DKA), and Dualis (DUA). Seeds were obtained from the Terres Inovia network for post-registration evaluation of winter oilseed rape cultivars in Northern France. The plant development stages were determined using BBCH-identification keys of oilseed rape [24,25].

2.2. Field Site and Experimental Design

Field trials were conducted from September 2017 to July 2018 at the CARAH experimental farm in Ath, Belgium (N 50°36′52.038″; E 3°45′54.025″), on a Luvisol soil (>2 m depth) with favorable natural drainage. Preceding crop was winter wheat. At sowing, N mineral content in 0–30 cm soil horizon was 78.6 kg ha−1. The experimental design consisted of split-plot randomized blocks with four replicates. On 28 August 2017, microplots (1.5 m × 12 m) were sown at a density of 60 seeds m−2. Two nitrogen fertilization treatments were considered: control/unfertilized plots (N−) and N fertilized plots with 240 kg ha−1 (N+). In the latter case, two doses of ammonium nitrate were applied: 60 kg N ha−1 at BBCH 00 and 180 kg N ha−1 at BBCH 31–32. The purpose of applying fertilizer at sowing was to observe the impact on root development during crop establishment in autumn. Temperature data were obtained from a meteorological station located less than 600 m from the field. Accumulated average daily temperature corresponding to growing degree days was calculated assuming a base temperature of 4.5 °C (GDD4.5) [26].

2.3. Root Scanning and Data Analysis

Root system expansion in the microplots was monitored with a minirhizotron system (CI-600, CID Bio-Science, Camas, USA). Immediately after sowing, transparent plexiglass tubes (105 cm length, 7 cm external diameter) were inserted with a manual auger of the same diameter (Eijkelkamp, Giesbeek, The Netherlands) at an angle of 45° [27]. A total of 18 tubes were installed, with 6 tubes per cultivars and 3 tubes per N treatment. The above-ground part of the tubes (15–20 cm) was covered with dark plastic to block out light. The minirhizotrons reached a depth of approximatively 0.5 m, corresponding to the upper soil layer in which fertilizers are most likely intercepted by roots. The captured images of the roots around the tube outer surface were divided in four zones (A, B, C, D) corresponding to soil horizons with an increasing depth of 12 cm, from 0 cm (ground level) to 48 cm depth (given the angle of the tube at 45° off the vertical) (Figure 1).
Captures of the soil profiles were redrawn with a digital pen and images were then quantified using Optimas v6.1 (Media Cybernetics). Lines were identified with the “automatically find lines” feature. The lengths of individual lines were obtained through the measurement explorer dialog box and summed for each image. The root length surface density (RLSD) was calculated as the length of visible roots per surface (mm cm−2), following to previous studies [28,29].

2.4. Determining Aerial Biomass, Nitrogen Concentration, Nitrogen Uptake, and Nitrogen Uptake Efficiency

Above-ground biomasses covering a plot surface of 0.75 m2 were collected on 27 November 2017 (BBCH 19, 647 GDD4.5) and 25 April 2018 (BBCH 65, 1026 GDD4.5). Samples were dried for 48 h at 70 °C. Nitrogen concentration in shoot tissues (N%shoot) was determined with a Vario MAX Cube (Elementar, Langenselbold, Germany). Nitrogen uptake (Nupt) was calculated as the product of biomass and N%shoot. Nitrogen Uptake Efficiency (NUpE) was calculated as follows [30]: NUpE = (NuptN+ − NuptN−)/N, where NuptN+ is N uptake after fertilization treatment and NuptN− without fertilization treatment and ΔN is the difference of applied nitrogen amount between N+ and N− treatments.

2.5. Canopy Optical Properties

The Normalized Difference Vegetation Index (NDVI) was measured on 30 October 2017 (BBCH 19) with the GreenSeeker® handheld crop sensor (Trimble, Sunnyvale, USA). The sensor was positioned at 60 cm from the top of the canopy to measure the reflectance at 660 nm (red) and 770 nm (near infrared) [31]. Ten measurements were made per microplot. The NDVI is commonly used as an indicator of aerial biomass production and nutrient status of crops [32,33,34]. In canola, an optimal time window to measure NDVI for grain yield prediction is between 50 and 70 days after sowing (DAS), while later readings are poorly predictive because they are influenced by stem elongation and flowering [35,36]. In the present study, NDVI measurements were made at 63 DAS (BBCH 19).
The chlorophyll index (CHL), flavonol index (FLAV), and the Nitrogen Balance Index (NBI) were measured on 23 April 2018 (BBCH 65) with the Dualex® Scientific + (Force-A, Orsay, France), a leafclip sensor. Within each microplot, measurements were conducted on the mature leaves of five plants.

2.6. Yield and Harvest Quality Traits

The central parts of each plot, corresponding to a surface of 10.5 m2, were harvested with a combine harvester (Wintersteiger Delta, Ried im Innkreis, Austria). Subsamples of seeds were analyzed for humidity, OilConc, ProtConc, and GLS (Table 1) by near infra-red spectroscopy (XDS NIR analyzer, Foss, Hilleroed, Denmark). A seed counter was used to determine TSW (Numigral, Chopin Technologies, Villeneuve-la-Garenne, France).

2.7. Statistical Analysis

The split-plot experimental design and the data statistical analyses were done with Statbox v7.0 (FBC Software). Assumptions for analysis of variance (ANOVA) were verified with D’Agostino-Spearman test. During ANOVA, the following variances were calculated: cultivar, fertilization, interaction (cultivar × fertilization), blocks (controlled factor), sub-blocks (controlled factor), and Residuals. ANOVA and Pearson’s correlations were performed at a significance level of α = 0.05. Newman-Keuls post-hoc tests were employed to classify the different results into homogeneous groups with different letters corresponding to homogeneous groups.

3. Results

3.1. Selecting Three Cultivars with Contrasting Root Morphologies

Three winter oilseed rape cultivars with dissimilar root morphologies were identified during a previous laboratory screening [37]. Twelve seedlings of DK Exception (DKE), DK Expansion (DKA), and Dualis (DUA) grew for one week in hydroponics. Cultivars were ranked in order of increasing root biomass production and lateral root length (Figure S1, Table S1). Eventually, these three genotypes were cultivated in field with two N fertilization treatments (N− or N+).

3.2. Observing Root Development Below Ground with a Camera System

Clear tubes were installed between plant rows to observe root development. The roots on the outer surface of the rhizotron (Figure 2) were imaged with a rotating scanner at key developmental stages over the seasons: crop establishment (BBCH 15, 18, and 19) during autumn, as well as stem elongation (BBCH 55) and flowering (BBCH 65) during spring.
The root system progressively explored soil layers and growth was the fastest at the end of autumn. After winter, roots were thinner and less visible. A first comparison was made between the two N fertilization treatments based on average values of the three genotypes. At BBCH 18, root length surface density (RLSD) in the deep soil horizon (36–48 cm) was significantly (P < 0.05) more important in N− compared to N+ conditions (Figure 3b). This clearly indicated a stimulation of root growth in the absence of N fertilization. A second comparison was done between cultivars in response to N treatment. There was an important heterogeneity between rhizotron images, but no obvious difference in rooting depth and placement was observed between cultivars across time points and N treatments (Figure 3f–j).

3.3. Above-Ground Biomass Production and Nitrogen Uptake Efficiency

Biomass and nitrogen concentration in plant tissues were measured following N treatments to calculate NUpE. At BBCH 19, the above-ground dry biomass was 56.0% more important in N+ compared to N− plots and the N concentration was 29.7% greater (Table 2). The total N uptake was approximately double between treatments. Therefore, the NUpE value was close to one, and this indicated the entire N fertilizer amount (60 kg N ha−1) applied in N+ at sowing was absorbed by the culture by autumn end. The NUpE differed between genotypes but not significantly (P > 0.05). At BBCH65, the above-ground dry biomass only increased by +5.9% in N+ compared to N− plots and the N content increased by +27.8%. The calculated NUpE was 0.24, suggesting that the second N fertilizer application was not totally taken up at this developmental stage. Under N+ treatment, the cultivars DKA and DUA showed significant (P < 0.05) greater Nupt compared to DKE (+21.0% and +27.4% respectively), those two cultivars having also the most branched root system in the laboratory screen (Figure S1).

3.4. Complementary Canopy Measurements

Pre-winter development of aerial biomass drives oilseed rape yield formation [38]. The Normalized Difference Vegetation Index (NDVI) is commonly used as an indicator of aerial biomass production and nutrient status in crops [32,33,34]. We then employed a proximal optical sensor (Greeseeker) to non-invasively probe leaf reflectance to assess plant vitality. Currently, this technique is also used in precision agriculture with drones and models to adjust nitrogen fertilizer levels [39]. In canola, an optimal time window to measure NDVI for grain yield prediction is between 50 to 70 days after sowing (DAS), while later readings are poorly predictive because they are influenced by stem elongation and flowering [35,36]. In this study, NDVI was measured at 63 DAS (BBCH 19). Another optical sensor (Dualex® leafclip) was used at flowering (BBCH 65), when nutrient absorption is maximum and leaf reflectance unsuitable. The Nitrogen Balance Index (NBI) is measured by leaf pigment fluorescence and is an indicator of plant N status [40,41,42].
At BBCH 19, the NDVI increased by +4.5% in N+ compared to N− plots. The cultivar DUA had significant (P < 0.05) greater values than DKE and DKA at both N treatments (Table 3). At BBCH 65, the NBI increased by +15.8% in N+ compared to N− plots. That parameter was greatly responding to N fertilization, but no difference was observed between cultivars.

3.5. Yield and Harvest Quality

Upon harvest, yield and seed quality traits were measured (Table 1). SY increased by 28.4% in N+ compared to N− plots (4.64 ± 0.39 vs. 3.61 ± 0.51 t ha−1, respectively) and the NConc increased by 4.7% (3.51 ± 0.06 vs. 3.36 ± 0.11 %) (Figure 4, Table 2). NConc and OilConc were negatively correlated, confirming previous findings [43,44]. The GLS and TSW were essentially cultivar dependent.
NUpE was estimated at harvest, considering the total nitrogen absorbed by harvestable organs (SNU). No marked difference was observed, and values were lower than the one observed at flowering (Table 2), reflecting a poor ability to remobilize nitrogen to harvest organs.

3.6. Predictiveness of Optical Parameters

In order to evaluate the predictiveness of the biosensing methods, indices measured at the canopy level during vegetative growth and flowering stages were correlated with yield and harvest quality parameters (Table S2). At the end of autumn (BBCH 19), NDVI and above-ground dry biomass were positively correlated (r = 0.83), and both negatively with RLSD (r = −0.60; −0.53 respectively). At full flowering (BBCH 65), NBI was only positively correlated (r = 0.67) with N concentration in above-ground plant biomass (N%shoot).
Furthermore, some canopy measurements recorded during vegetative growth in autumn and at flowering were significantly correlated to harvest parameters. The NDVI readings in autumn were correlated to SY (r = 0.75), while NBI at flowering to NConc in seeds and SNU (r = 0.59 and r = 0.47, respectively).

4. Discussion

Improving NUE is essential to ensure the environmental and economic sustainability of oilseed rape production. While the aerial organs are extensively considered in breeding programs, equivalent examinations in roots are receiving much less attention because of the practical challenge of working with the below-ground organs [45]. In situ observation of the root system is not trivial and conventional methods (e.g., soil coring) can be tedious [46]. Here, we employed a minirhizotron system to observe dynamic aspects of root growth with a camera through clear tubes inserted in soil (Figure 1). That non-destructive method allows following root development in soil horizons during the culture cycle with a fair imaging resolution. Creating an intimate contact between transparent tube and soil improves the image quality [27], notwithstanding soil texture and moisture can hamper the further image analysis process. In this study, image processing with the dedicated manufacturer software was not satisfactory enough because of the root overlapping and water condensation outside the tubes, as also reported in other studies [47]. To overcome that problem, root items were redrawn by hand, but it was time consuming. Current tools for processing rhizotron images are often restricted to a limited number of variables and show limitation with highly branched root systems [48]. Improvement has still to be done to take full advantage of the high throughput scanning method. Before wintertime, we observed the full deployment of the root system (Figure 2 and Figure 3). The end of autumn season was a critical time window to observe root development. At BBCH 18, N fertilization treatment reduced RLSD in the deep soil horizon (Figure 3) and stimulated above-ground biomass production (Table 1). During spring, roots were less clearly visible, resulting in lower RLSD values (Figure 2 and Figure 3). This could be attributed to less important resource allocation to the roots, to frost damage during winter, and/or to deeper soil exploration, out of sight of the imaging system (>0.5 m). Indeed, the root system of oilseed rape can reach more than 2.3 m depth, with roots up to 1.5 m before winter [30].
Winter oilseed rape is seen as a poor N− efficient crop that has a low seed production per N unit applied, mainly because of a low N remobilization from vegetative to reproductive organs [39,49]. Nonetheless, the plant can absorb great N amounts during autumn and winter, making it a valuable catch crop to reduce N leaching in soil [50]. Assessing NUpE at the beginning of the growth cycle is crucial, since up to one-third of the total N uptake occurs before winter [39]. The NUpE values were greater in autumn than in spring (Table 2). This can be explained by a fraction of N absorbed before winter being lost with littering of frozen leaves, but still susceptible of being mineralized and reabsorbed during spring [49,51]. An estimation of that lost is between 20 and 100 kg N ha−1 [51,52]. During stem elongation and onwards, occurs the sequential fall of senescent leaves with large N content (2–2.5% of dry weight) [49]. The pod filling is an important phase during which plants are still actively taking up N [39,49]. The sink capacity of the pods is considered as a bottleneck for N remobilization and this could be enhanced through yield improvement [53]. At harvest, we estimated nearly one-fifth of applied N was recovered in the seeds (Table 2), that proportion being lower than in previous reports [49,52,54]. Large N input (240 kg ha−1) and drought stress likely contributed to this low value. Indeed, the year 2018 weather conditions were extreme, and two-thirds less precipitation than annual average during flowering and pod filling, severely depressed yield. Furthermore, the negative impact of drought stress on N recovery to seeds is well documented [52,55].
We appraised two rapid and non-destructive proximal remote sensing methods based on optical properties of plants. The NDVI at the canopy level during autumn was a good estimate of above-ground biomass production in autumn and seed yield at harvest, while the NBI at the leaf level during flowering time, of N concentration in plants tissues and seeds at harvest. Previous studies showed that NUE in oilseed rape was strongly correlated to leaf N content at flowering, especially under low N supply [56]. Therefore, NBI measurements could be convenient for improving rapid screening procedures. The NDVI and NBI delivered complementary information but they were not correlated to each other.
We based our research program on the premise that NUpE could be improved by favoring a more branched root system that explores a larger soil volume to limit N leaching [37]. In a previous report, we challenged a panel of winter oilseed rape cultivars at a young developmental stage (one-week-old seedlings) with divergent nitrate supplies and identified contrasting root morphologies [37]. Some roots traits of young seedlings (e.g., length of the primary root length) were positively correlated to NConc in subsequent field trials. Such a screening strategy in laboratory conditions may speed up for breeders the delivery of genotypes with ideal morphological features. Positive correlation between root traits in laboratory conditions and field performance was also observed in other studies [57]. For example, the primary root length measured in hydroponics was positively correlated to the field yield, and the lateral root density to the concentration of mineral elements in leaves. In this paper, in-field differences between cultivars were not marked regardless of the N fertilization treatments. Nonetheless, at BBCH 65 and following N+ treatment, the DKA and DUA cultivars had greater Nupt than DKE, the latter one showed less lateral root outgrowth in laboratory conditions (Figure S1, Table S1).

5. Conclusions

We have carried a one-year pilot experiment with three hybrid cultivars reflecting the market trends. Hybrids are high-yielding lines with enhanced NUE [58]. In perspective, more genotypes with contrasting root morphologies could be screened to further strengthen the corroboration between root morphological features in laboratory and soil conditions.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4395/10/1/105/s1, Figure S1: Root morphologies of three oilseed rape cultivars cultivated in hydroponics, Figure S2: Global ANOVA of measured traits at different periods in field experiment, Table S1: Root system parameters of three oilseed rape cultivars cultivated in hydroponics, Table S2: Pearson’s correlation coefficients between parameters measured at different time points.

Author Contributions

Conceptualization, J.L. and C.H.; formal analysis, J.L. and A.L.; funding acquisition, J.L. and C.H.; investigation, J.L. and A.L.; methodology, J.L. and L.H.; writing—original draft, J.L.; writing—review and editing, J.L., L.H., and C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MIS and PDR T.0116.19 grants from Fonds de la Recherche Scientifique (F.R.S.-FNRS).

Acknowledgments

C.H. is an F.R.S.-FNRS research associate. The authors thankfully acknowledge the Centre pour l’Agronomie et l’Agro-industrie de la Province de Hainaut (CARAH) for the field trial and analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pradhan, P.; Fischer, G.; van Velthuizen, H.; Reusser, D.E.; Kropp, J.P. Closing yield gaps: How sustainable can we be? PLoS ONE 2015, 10, e0129487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Berners-Lee, M.; Kennelly, C.; Watson, R.; Hewitt, C.N. Current global food production is sufficient to meet human nutritional needs in 2050 provided there is radical societal adaptation. Elem. Sci. Anth. 2018, 6, 52. [Google Scholar] [CrossRef]
  3. Wick, K.; Heumesser, C.; Schmid, E. Groundwater nitrate contamination: Factors and indicators. J. Environ. Manag. 2012, 111, 178–186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. World Health Organization. Chapter 8: Chemical Aspects. In Guidelines for Drinking-Water Quality: Fourth Edition Incorporating the First Addendum; World Health Organization: Geneva, Switzerland, 2017; pp. 155–201. [Google Scholar]
  5. Kant, S.; Bi, Y.M.; Rothstein, S.J. Understanding plant response to nitrogen limitation for the improvement of crop nitrogen use efficiency. J. Exp. Bot. 2011, 62, 1499–1509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Han, M.; Okamoto, M.; Beatty, P.H.; Rothstein, S.J.; Good, A.G. The genetics of nitrogen use efficiency in crop plants. Annu. Rev. Genet. 2015, 49, 269–289. [Google Scholar] [CrossRef]
  7. Stahl, A.; Pfeifer, M.; Frisch, M.; Wittkop, B.; Snowdon, R.J. Recent genetic gains in nitrogen use efficiency in oilseed rape. Front. Plant Sci. 2017, 8, 963. [Google Scholar] [CrossRef] [Green Version]
  8. Laine, P.; Ourry, A.; Macduff, J.; Boucaud, J.; Salette, J. Kinetic-parameters of nitrate uptake by different catch crop species: Effects of low temperatures or previous nitrate starvation. Physiol. Plant. 1993, 88, 85–92. [Google Scholar] [CrossRef]
  9. McAllister, C.H.; Beatty, P.H.; Good, A.G. Engineering nitrogen use efficient crop plants: The current status. Plant Biotechnol. J. 2012, 10, 1011–1025. [Google Scholar] [CrossRef]
  10. Havé, M.; Marmagne, A.; Chardon, F.; Masclaux-Daubresse, C. Nitrogen remobilization during leaf senescence: Lessons from Arabidopsis crops. J. Exp. Bot. 2017, 68, 2513–2529. [Google Scholar]
  11. Hodge, A.; Berta, G.; Doussan, C.; Merchan, F.; Crespi, M. Plant root growth, architecture and function. Plant Soil. 2009, 321, 153–187. [Google Scholar] [CrossRef]
  12. Garnett, T.; Conn, V.; Kaiser, B.N. Root based approaches to improving nitrogen use efficiency in plants. Plant Cell. Environ. 2009, 32, 1272–1283. [Google Scholar] [CrossRef] [PubMed]
  13. De Pessemier, J.; Chardon, F.; Juraniec, M.; Delaplace, P.; Hermans, C. Natural variation the root morphological response to nitrate supply in Arabidopsis thaliana. Mech. Dev. 2013, 130, 45–53. [Google Scholar] [CrossRef]
  14. López-Bucio, J.; Cruz-Ramírez, A.; Herrera-Estrella, L. The role of nutrient availability in regulating root architecture. Curr. Opin. Plant Biol. 2003, 6, 280–287. [Google Scholar] [CrossRef]
  15. Gruber, B.D.; Giehl, R.F.; Friedel, S.; von Wirén, N. Plasticity of the Arabidopsis root system under nutrient deficiencies. Plant Physiol. 2013, 163, 161–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Qin, L.; Walk, T.C.; Han, P.; Chen, L.; Zhang, S.; Li, Y.; Hu, X.; Xie, L.; Yang, Y.; Liu, J.; et al. Adaptation of roots to nitrogen deficiency revealed by 3D quantification and proteomic analysis. Plant Physiol. 2019, 179, 329–347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lynch, J.P. Steep, cheap and deep: An ideotype to optimize water and N acquisition by maize root systems. Ann. Bot. 2013, 112, 347–357. [Google Scholar] [CrossRef] [Green Version]
  18. Li, X.; Zeng, R.; Liao, H. Improving crop nutrient efficiency through root architecture modifications. J. Integr. Plant. Biol. 2016, 58, 193–202. [Google Scholar] [CrossRef] [Green Version]
  19. Thorup-Kristensen, K.; Kirkegaard, J. Root system-based limits to agriculture productivity and efficiency: The farming systems context. Ann. Bot. 2016, 118, 573–592. [Google Scholar] [CrossRef]
  20. Berry, P.M.; Spink, J.; Foulkes, M.J.; White, P.J. The physiological basis of genotypic differences in nitrogen use efficiency in oilseed rape (Brassica napus L.). Field Crops Res. 2010, 119, 365–373. [Google Scholar] [CrossRef]
  21. Kamh, M.; Wiesler, F.; Ulas, A.; Horst, W. Root growth and N-uptake activity of oilseed rape (Brassica napus L.) cultivars differing in nitrogen efficiency. J. Plant Nutr. Soil Sci. 2005, 168, 130–137. [Google Scholar] [CrossRef]
  22. Ulas, A.; Schulte auf’m Erley, G.; Kamh, M.; Wiesler, F.; Horst, W.J. Root-growth characteristics contributing to genotypic variation in nitrogen efficiency of oilseed rape. J. Plant Nutr. Soil. Sci. 2012, 175, 489–498. [Google Scholar] [CrossRef]
  23. Behrens, T.; Ulas, A.; Wiesler, F.; Horst, W.J. Agronomic traits contributing to nitrogen efficiency of winter oilseed rape cultivars. Field Crops Res. 2011, 124, 114–123. [Google Scholar]
  24. Weber, E.; Bleiholder, H. Erläuterungen zu den BBCH-Dezimal-Codes für die Entwicklungsstadien von Mais, Raps, Faba-Bohne, Sonnenblume und Erbse—Mit Abbildungen. Gesunde Pflanz. 1990, 42, 308–321. [Google Scholar]
  25. Lancashire, P.D.; Bleiholder, H.; Van Den Boom, T.; Langelüddeke, P.; Strauss, R.; Weber, E.; Witzenberger, A. A uniform decimal code for growth stages of crops and weeds. Ann. Appl. Biol. 1991, 119, 561–601. [Google Scholar] [CrossRef]
  26. Gabrielle, B.; Denoroy, P.; Gosse, G.; Justes, E.; Andersen, M.N. A model of leaf area development and senescence for winter oilseed rape. Field Crops Res. 1998, 57, 209–222. [Google Scholar] [CrossRef]
  27. Smit, A.L.; George, E.; Groenwold, J. Root observation and measurements at (transparent) interfaces with soil. In Root Methods: A Handbook; Smit, A.L., Bengough, A.G., Engels, C., van Noordwijk, M., Pellerin, S., van de Geijn, S.C., Eds.; Springer: Berlin, Germany, 2000; pp. 235–271. [Google Scholar]
  28. Muñoz-Romero, V.; López-Bellido, L.; López-Bellido, R.J. The effects of tillage system on chickpea root growth. Field Crops Res. 2012, 128, 76–81. [Google Scholar] [CrossRef]
  29. Postic, F.; Beauchêne, K.; Gouache, D.; Doussan, C. Scanner-Based minirhizotrons help to highlight relations between deep roots and yield in various wheat cultivars under combined water and nitrogen deficit conditions. Agronomy 2019, 9, 297. [Google Scholar] [CrossRef] [Green Version]
  30. Dresbøll, D.B.; Rasmussen, I.S.; Thorup-Kristensen, K. The significance of litter loss and root growth on nitrogen efficiency in normal and semi-dwarf winter oilseed rape genotypes. Field Crops Res. 2016, 186, 166–178. [Google Scholar] [CrossRef]
  31. Tremblay, N.; Wang, Z.; Ma, B.-L.; Belec, C.; Vigneault, P. A comparison of crop data measured by two commercial sensors for variable-rate nitrogen application. Precis. Agric. 2009, 10, 145–161. [Google Scholar] [CrossRef]
  32. Serrano, L.; Filella, I.; Peñuelas, J. Remote sensing of biomass and yield of winter wheat under different nitrogen supplies. Crop Sci. 2000, 40, 723–731. [Google Scholar] [CrossRef] [Green Version]
  33. Moges, S.M.; Raun, W.R.; Mullen, R.W.; Freeman, K.W.; Johnson, G.V.; Solie, J.B. Evaluation of green, red and near infrared bands for predicting winter wheat biomass, nitrogen uptake, and final grain yield. J. Plant Nutr. 2004, 27, 1431–1441. [Google Scholar] [CrossRef]
  34. Magney, T.S.; Eitel, J.U.H.; Vierling, L.A. Mapping wheat nitrogen uptake from RapidEye vegetation indices. Precis. Agric. 2017, 18, 429–451. [Google Scholar] [CrossRef]
  35. Holzapfel, C.B.; Lafond, G.P.; Brandt, S.A.; Bullock, P.R.; Irvine, R.B.; Morrison, J.M.; May, W.E.; James, D.C. Estimating canola (Brassica napus L.) yield potential using an active optical sensor. Can. J. Plant Sci. 2009, 89, 1149–1160. [Google Scholar] [CrossRef]
  36. Cowley, R.B.; Luckett, D.J.; Moroni, J.S.; Diffey, S.M. Use of remote sensing to determine the relationship of early vigour to grain yield in canola (Brassica napus L.) germplasm. Crop. Pasture Sci. 2014, 65, 1288–1299. [Google Scholar] [CrossRef]
  37. Louvieaux, J.; De Gernier, H.; Hermans, C. Exploiting genetic variability of root morphology as a lever to improve nitrogen use efficiency in oilseed rape. In Engineering Nitrogen Utilization in Crop Plants; Shrawat, A., Zayed, A., Lightfoot, D.A., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 185–206. [Google Scholar]
  38. Diepenbrock, W. Yield analysis of winter oilseed rape (Brassica napus L.): A review. Field Crops Res. 2000, 67, 35–49. [Google Scholar] [CrossRef]
  39. Bouchet, A.; Laperche, A.; Bissuel-Belaygue, C.; Snowdon, R.; Nesi, N.; Stahl, A. Nitrogen use efficiency in rapeseed. A review. Agron. Sustain. Dev. 2016, 36, 1–20. [Google Scholar] [CrossRef]
  40. Cartelat, A.; Cerovic, Z.G.; Goulas, Y.; Meyer, S.; Lelarge, C.; Prioul, J.-L.; Barbottin, A.; Jeuffroy, M.-H.; Gate, P.; Agati, G.; et al. Optically assessed contents of leaf polyphenolics and chlorophyll as indicators of nitrogen deficiency in wheat (Triticum aestivum L.). Field Crops Res. 2005, 91, 35–49. [Google Scholar] [CrossRef]
  41. Cerovic, Z.G.; Ghozlen, N.B.; Milhade, C.; Obert, M.; Debuisson, S.; Le Moigne, M. Nondestructive diagnostic test for nitrogen nutrition of grapevine (Vitis vinifera L.) based on Dualex leaf-clip measurements in the field. J. Agric. Food Chem. 2015, 63, 3669–3680. [Google Scholar] [CrossRef]
  42. Padilla, F.M.; Gallardo, M.; Peña-Fleitas, M.T.; de Souza, R.; Thompson, R.B. Proximal optical sensors for nitrogen management of vegetable crops: A review. Sensors 2018, 18, 2083. [Google Scholar] [CrossRef] [Green Version]
  43. Grami, B.; Stefansson, B.R.; Baker, R.J. Genetics of protein and oil content in summer rape-heritability, number of effective factors, and correlations. Can. J. Plant Sci. 1977, 57, 937–943. [Google Scholar] [CrossRef]
  44. Jolivet, P.; Deruyffelaere, C.; Boulard, C.; Quinsac, A.; Savoire, R.; Nesi, N.; Chardot, T. Deciphering the structural organization of the oil bodies in the Brassica napus seed as a mean to improve the oil extraction yield. Ind. Crop Prod. 2013, 44, 549–557. [Google Scholar] [CrossRef]
  45. Araus, J.L.; Cairns, J.E. Field high-throughput phenotyping: The new crop breeding frontier. Trends. Plant Sci. 2014, 19, 52–61. [Google Scholar] [CrossRef] [PubMed]
  46. Böhm, W. Methods of studying root systems. In Ecological Studies; Billings, W.D., Golley, F., Lange, O.L., Olson, J.S., Eds.; Springer: Berlin, Germany, 1979; Volume 33. [Google Scholar]
  47. Ohashi, A.Y.P.; de Matos Pires, R.C.; de Oliveira Silva, A.L.; Silva dos Santos, L.N.; Matsura, E.E. Minirhizotron as an in-situ tool for assessing sugarcane root system growth and distribution. Agri. Res. Tech. Open. Access J. 2019, 22, 556182. [Google Scholar]
  48. Lobet, G.; Draye, X. Novel scanning procedure enabling the vectorization of entire rhizotron-grown root systems. Plant Methods 2013, 9, 1–10. [Google Scholar] [CrossRef] [Green Version]
  49. Malagoli, P.; Lainé, P.; Rossato, L.; Ourry, A. Dynamics of nitrogen uptake and mobilization in field-grown winter oilseed rape (Brassica napus) from stem extension to harvest: I. Global N flows between vegetative and reproductive tissues in relation to leaf fall and their residual N. Ann. Bot. 2005, 95, 853–861. [Google Scholar] [CrossRef] [Green Version]
  50. Rossato, L.; Lainé, P.; Ourry, A. Nitrogen storage and remobilization in Brassica napus L. during the growth cycle: Nitrogen fluxes within the plant and changes in soluble protein patterns. J. Exp. Bot. 2001, 361, 1655–1663. [Google Scholar] [CrossRef] [Green Version]
  51. Dejoux, J.-F.; Recous, S.; Meynard, J.-M.; Trinsoutrot, I.; Leterme, P. The fate of nitrogen from winter-frozen rapeseed leaves: Mineralization, fluxes to the environment and uptake by rapeseed crop in spring. Plant Soil 2000, 218, 257–272. [Google Scholar] [CrossRef]
  52. Schjoerring, J.K.; Bock, J.G.H.; Gammelvind, L.; Jensen, C.R.; Mogensen, V.O. Nitrogen incorporation and remobilization in different shoo components of field-grown winter oilseed rape (Brassica napus L.) as affected by rate of nitrogen application and irrigation. Plant Soil 1995, 177, 255–264. [Google Scholar] [CrossRef]
  53. Tilsner, J.; Kassner, N.; Struck, C.; Lohaus, G. Amino acid contents and transport in oilseed rape (Brassica napus L.) under different nitrogen conditions. Planta 2005, 221, 328–338. [Google Scholar] [CrossRef]
  54. Villar, N.; Aranguren, M.; Castellón, A.; Besga, G.; Aizpurua, A. Soil nitrogen dynamics during an oilseed rape (Brassica napus L.) growing cycle in a humid Mediterranean climate. Sci. Rep. 2019, 9, 13864. [Google Scholar] [CrossRef] [Green Version]
  55. Smith, C.J.; Wright, G.C.; Woodroofe, M.R. The effect of irrigation and nitrogen fertilizer on rapeseed (Brassica napus) production in south-eastern Australia. Irrig. Sci. 1988, 9, 15–25. [Google Scholar] [CrossRef]
  56. Stahl, A.; Friedt, W.; Wittkop, B.; Snowdon, R. Complementary diversity for nitrogen uptake and utilization efficiency reveals broad potential for increased sustainability of oilseed rape production. Plant Soil 2016, 400, 245–262. [Google Scholar] [CrossRef]
  57. Thomas, C.L.; Graham, N.S.; Hayden, R.; Meacham, M.C.; Neugebauer, K.; Nightingale, M.; Dupuy, L.X.; Hammond, J.P.; White, J.P.; Broadley, M.R. High-throughput phenotyping (HTP) identifies seedling root traits linked to variation in seed yield and nutrient capture in field-grown oilseed rape (Brassica napus L.). Ann. Bot. 2016, 118, 655–665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Kessel, B.; Schierholt, A.; Becker, H.C. Nitrogen use efficiency in a genetically diverse set of winter oilseed rape (Brassica napus L.). Crop Sci. 2012, 52, 2546–2554. [Google Scholar] [CrossRef]
Figure 1. Minirhizotron positioning in soil. Transparent tubes (105 cm length, 7 cm diameter) were inserted at 45° angle with a manual auger. Four observation zones were defined: A (0–12 cm), B (12–24 cm), C (24–36 cm) and D (36–48 cm).
Figure 1. Minirhizotron positioning in soil. Transparent tubes (105 cm length, 7 cm diameter) were inserted at 45° angle with a manual auger. Four observation zones were defined: A (0–12 cm), B (12–24 cm), C (24–36 cm) and D (36–48 cm).
Agronomy 10 00105 g001
Figure 2. Example of image analysis and quantification of root length surface density. The DK Exception (DKE) cultivar received N+ treatment. Root length surface density (RLSD) values (mm cm−2) are shown in different soil horizons and at different time points of the growing season. The soil horizons are indicated by different letters: A = 0–12 cm; B = 12–24 cm; C = 24–36 cm, and D = 36–48 cm. The BBCH-scale indicates phenological development stages [24,25]. GDD: growing degree days.
Figure 2. Example of image analysis and quantification of root length surface density. The DK Exception (DKE) cultivar received N+ treatment. Root length surface density (RLSD) values (mm cm−2) are shown in different soil horizons and at different time points of the growing season. The soil horizons are indicated by different letters: A = 0–12 cm; B = 12–24 cm; C = 24–36 cm, and D = 36–48 cm. The BBCH-scale indicates phenological development stages [24,25]. GDD: growing degree days.
Agronomy 10 00105 g002
Figure 3. Root length surface density of three oilseed rape cultivars in different soil depths and at different time points of the growing season. Measurements were done at BBCH 15 (a,f), BBCH 18 (b,g), BBCH 19 (c,h), BBCH 55 (d,i), and BBCH 65 (e,j). Mean results according to the fertilizer treatments (n = 9 ± s.e.) (ae) and to individual cultivars (n = 3 ± s.e) (fj). Stars indicate significant differences at level α = 0.05.
Figure 3. Root length surface density of three oilseed rape cultivars in different soil depths and at different time points of the growing season. Measurements were done at BBCH 15 (a,f), BBCH 18 (b,g), BBCH 19 (c,h), BBCH 55 (d,i), and BBCH 65 (e,j). Mean results according to the fertilizer treatments (n = 9 ± s.e.) (ae) and to individual cultivars (n = 3 ± s.e) (fj). Stars indicate significant differences at level α = 0.05.
Agronomy 10 00105 g003
Figure 4. Comparison of yield and harvest quality parameters. The spider plot shows the percentage variation of the normalized trait values relative to the average of the unfertilized plots (0% indicates no difference compared to the average of N− plots). Parameters are defined in Table 1.
Figure 4. Comparison of yield and harvest quality parameters. The spider plot shows the percentage variation of the normalized trait values relative to the average of the unfertilized plots (0% indicates no difference compared to the average of N− plots). Parameters are defined in Table 1.
Agronomy 10 00105 g004
Table 1. Definition of yield and harvest quality parameters.
Table 1. Definition of yield and harvest quality parameters.
Abbreviation.Description
SYSeed yield corrected to a standard water content of 9%
TSW1000 seed weight with moisture adjusted to 9%
ProtConcProtein concentration in dry seeds
NConcNitrogen concentration in dry seeds = ProtConc/6.25
OilConcOil concentration in dry seeds
SNUSeed nitrogen uptake = NConc × SY
ProteinYProtein yield = ProtConc × SY
OilYOil yield = OilConc × SY
GLSGlucosinolate concentration in dry seeds
Table 2. Above-ground biomass measurements and harvest parameters. Mean values according to the fertilizer treatments (n = 12 ± s.e.) and to individual cultivars (n = 4 ± s.e). Letters indicate statistical differences at significance level α = 0.05. Pooled standard errors (s.e.pooled) of N− and N+ populations are mentioned according to fertilizer treatments and cultivars. The effects of different factors as a percentage of the observed variation are given in Figure S2.
Table 2. Above-ground biomass measurements and harvest parameters. Mean values according to the fertilizer treatments (n = 12 ± s.e.) and to individual cultivars (n = 4 ± s.e). Letters indicate statistical differences at significance level α = 0.05. Pooled standard errors (s.e.pooled) of N− and N+ populations are mentioned according to fertilizer treatments and cultivars. The effects of different factors as a percentage of the observed variation are given in Figure S2.
BBCH 19BBCH 65
Dry BiomassN%shootNuptNUpEDry BiomassN%shootNuptNUpE
(t ha−1) (kg ha−1) (t ha−1) (kg ha−1)
Mean valuesN−2.31 b ± 0.172.63 b ± 0.0961.4 b ± 5.6 6.75 ± 0.422.50 b ± 0.10168.5 b ± 11.7
N+3.60 a ± 0.143.41 a ± 0.11122.2 a ± 5.01.01 ± 0.117.14 ± 0.393.20 a ± 0.08227.1 a ± 11.80.24 ± 0.06
s.e.pooled±0.22±0.14±7.4 ±0.57±0.12±16.6
DKEN−2.66 ± 0.382.50 ± 0.1767.9 ± 12.6 5.75 ± 0.762.53 ± 0.18144.3 ± 17.8
N+3.52 ± 0.143.42 ± 0.16120.4 ± 7.70.88 ± 0.235.92 ± 0.523.34 ± 0.13195.6 b ± 10.20.21 ± 0.08
s.e.pooled±0.41±0.23±14.7 ±0.92±0.22±20.5
DKAN−1.92 ± 0.172.67 ± 0.1551.8 ± 7.0 7.45 ± 0.842.63 ± 0.17194.8 ± 22.2
N+3.71 ± 0.283.36 ± 0.16124.0 ± 8.31.20 ± 0.207.62 ± 0.463.10 ± 0.11236.7 a ± 19.70.17 ± 0.14
s.e.pooled±0.33±0.22±10.8 ±0.96±0.21±27.9
DUAN−2.35 ± 0.202.72 ± 0.1864.6 ± 9.2 7.04 ± 0.332.35 ± 0.17166.4 ± 16.3
N+3.57 ± 0.343.45 ± 0.28122.1 ± 11.90.96 ± 0.147.88 ± 0.653.16 ± 0.15249.2 a ± 22.60.34 ± 0.13
s.e.pooled±0.39±0.33±15.0 ±0.73±0.23±27.9
Harvest Parameters
SYNConcOilConcGLSTSWSNUNUpE
(t ha−1)(%)(%)(µmol g−1)(g)(kg ha−1)
Mean valuesN−3.61 b ± 0.153.36 b ± 0.0346.6 a ± 0.316.2 ± 0.45.16 ± 0.11110 b ± 4
N+4.64 a ± 0.113.51 a ± 0.0245.3 b ± 0.217.0 ± 0.55.15 ± 0.11148 a ± 40.16 ± 0.02
s.e.pooled±0.18±0.04±0.3±0.7±0.15±6
DKEN−3.63 ± 0.343.40 ± 0.0345.6 bc ± 0.217.2 ± 0.25.46 ± 0.10112 ± 11
N+4.62 ± 0.213.50 ± 0.0144.7 c ± 0.217.9 ± 1.15.44 ± 0.08147 ± 60.14 ± 0.03
s.e.pooled±0.40±0.04±0.2±1.1±0.13±13
DKAN−3.40 ± 0.243.43 ± 0.0546.4 b ± 0.416.4 ± 0.95.30 ± 0.06106 ± 7
N+4.86 ± 0.133.55 ± 0.0345.7 bc ± 0.117.3 ± 0.85.32 ± 0.08157 ± 50.21 ± 0.02
s.e.pooled±0.27±0.06±0.4±1.2±0.10±9
DUAN−3.80 ± 0.203.24 ± 0.0447.7 a ± 0.115.0 ± 0.44.71 ± 0.09112 ± 5
N+4.43 ± 0.213.49 ± 0.0345.6 bc ± 0.215.7 ± 0.54.69 ± 0.03140 ± 50.12 ± 0.04
s.e.pooled±0.29±0.05±0.2±0.7±0.10±8
Table 3. Leaf optical parameters in response to nitrogen fertilization treatment. The Normalized Difference Vegetation Index (NDVI) and the Nitrogen Balance Index (NBI) were measured on canopy of control (N−) or fertilized (N+) plots respectively at BBCH 19 and BBCH 65. Mean values according to the fertilizer treatments (n = 12 ± std) and to individual cultivars (n = 4 ± std). Pooled standard deviations (stdpooled) of N− and N+ populations are mentioned according to the fertilizer treatment and cultivars. Letters indicate statistical differences at significance level α = 0.05.
Table 3. Leaf optical parameters in response to nitrogen fertilization treatment. The Normalized Difference Vegetation Index (NDVI) and the Nitrogen Balance Index (NBI) were measured on canopy of control (N−) or fertilized (N+) plots respectively at BBCH 19 and BBCH 65. Mean values according to the fertilizer treatments (n = 12 ± std) and to individual cultivars (n = 4 ± std). Pooled standard deviations (stdpooled) of N− and N+ populations are mentioned according to the fertilizer treatment and cultivars. Letters indicate statistical differences at significance level α = 0.05.
NDVINBI
(BBCH 19)(BBCH 65)
Mean valuesN−0.80 b ± 0.0329.8 b ± 2.5
N+0.84 a ± 0.0134.6 a ± 3.9
stdpooled±0.02±3.3
DKEN−0.79 ± 0.0329.4 ± 0.8
N+0.84 ± 0.0136.0 ± 5.5
stdpooled±0.02±3.9
DKAN−0.80 ± 0.0331.3 ± 3.9
N+0.84 ± 0.0136.0 ± 2.6
stdpooled±0.02±3.3
DUAN−0.82 ± 0.0328.8 ± 1.7
N+0.85 ± 0.0131.7 ± 1.8
stdpooled±0.02±1.7

Share and Cite

MDPI and ACS Style

Louvieaux, J.; Leclercq, A.; Haelterman, L.; Hermans, C. In-Field Observation of Root Growth and Nitrogen Uptake Efficiency of Winter Oilseed Rape. Agronomy 2020, 10, 105. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10010105

AMA Style

Louvieaux J, Leclercq A, Haelterman L, Hermans C. In-Field Observation of Root Growth and Nitrogen Uptake Efficiency of Winter Oilseed Rape. Agronomy. 2020; 10(1):105. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10010105

Chicago/Turabian Style

Louvieaux, Julien, Antoine Leclercq, Loïc Haelterman, and Christian Hermans. 2020. "In-Field Observation of Root Growth and Nitrogen Uptake Efficiency of Winter Oilseed Rape" Agronomy 10, no. 1: 105. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10010105

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop