Introduction

A vigorous, adaptive root system is critical for improved water and nutrient uptake, which in turn increases crop grain yield (GY), especially under stress conditions (Fang et al. 2017). Therefore, the measurement of root-related traits is just as important in agriculture as those of the shoot (Fageria 2013). As it is practically impossible to extract an intact root system from the field soil, non-destructive techniques need to be applied for root characterization, despite their inherent limitations (Ehosioke et al. 2020; Cabal et al. 2021).

The electrical capacitance method is considered a promising in situ phenotyping tool in plant breeding because it allows numerous plants to be measured repeatedly during the growing season, with the subsequent harvest of seeds for further reproduction (Chloupek et al. 2010; Středa et al. 2020). Unlike some other methods (root core, rhizotron), the capacitance technique provides information on individual plants due to the electrical separation of adjacent plants (Ellis et al. 2013). The electrical capacitance of the root–soil system (CR), measured between a ground and a plant electrode with alternating current (AC), was found to correlate with the root system size (RSS), including both root mass and length (Chloupek et al. 2010; Carlson and Smart 2016). The first model by Dalton (1995) assumed that roots are leaky cylindrical capacitors with membrane dielectrics storing electric charges. The CR measured is related to the polarized membrane surface area. Various experiments verified that the root system was electrically connected along most of its length, but the contribution of distal fine roots to CR was uncertain (Ozier-Lafontaine and Bajazet 2005; Ellis et al. 2013). Mary et al. (2018) found that AC reached the majority of the grapevine root system, penetrating to a depth of 30–40 cm into the soil. Conversely, others questioned Dalton’s model and the efficacy of the CR method, observing localized current leakage in the proximal root part (Urban et al. 2011; Peruzzo et al. 2020). Dietrich et al. (2012) proposed a revised model, showing that CR was dominated by the stem between the plant electrode and the substrate surface, and was related to the stem cross-sectional area (CR–RSS correlations were attributed merely to root–shoot allometry). Nevertheless, Cseresnyés et al. (2020b) experimentally demonstrated a strong dependence of CR on the mass of root branches in the soil, though the stem-base capacitance was considerable as well. Recent work revealed the potential of CR to directly quantify root traits in dry soil (which is more electrically resistive than the roots) rather than under wet conditions, i.e. at field capacity (Gu et al. 2021), a result which is relevant mainly for field studies. It is commonly agreed that, as root current pathways are affected by root water content, tissue density and suberization, the CR values represent not only the geometrical size but also the functional activity of the root system (Dalton 1995; Ellis et al. 2013; Peruzzo et al. 2020).

Besides numerous pot studies (summarized by Ehosioke et al. 2020), the CR method was efficiently used for evaluating RSS in field-grown wheat (Středa et al. 2012; Heřmanská et al. 2015; Cseresnyés et al. 2021), barley (Chloupek et al. 2006, 2010; Svačina et al. 2014), canola (Wu and Ma 2016), maize and soybean (Cseresnyés et al. 2018). Owing to the sensitivity of CR to external factors, capacitance data can only be compared when the same species is grown in the same substrate and measured with the same electrode type and placement (Chloupek et al. 2006; Středa et al. 2012). Variation in soil water content (SWC) has a marked effect on CR due to changes in soil–ground electrode contact (Ellis et al. 2013) and in the root to soil conductance ratio (Gu et al. 2021). A close positive relationship between CR and the SWC of the root zone was previously shown for several species (Cseresnyés et al. 2018, 2020a). Using a species-specific exponential function, a saturation capacitance, CR* (which can be detected in water-saturated soil) was calculated from the measured CR and the corresponding SWC. This transformation allowed a comparison of the field data recorded under different SWC values at different dates (for a detailed description, see Cseresnyés et al. 2018). In this manner, the ontogenic pattern of CR* was revealed, as a proxy for root uptake activity for field-grown maize and soybean with a peak at flowering, in accordance with previous findings obtained using conventional root methods (Fageria 2013).

Nevertheless, to date there has been no field research aimed at using the capacitance method to monitor root dynamics in situ in crops grown under different cultivation and environmental conditions. Due to global climate change, many experiments have aimed to study the plant response to increased atmospheric carbon dioxide concentration ([CO2]) over the past few decades, focusing on the major food crops (Kirkham et al. 2011; Wang et al. 2013; Broberg et al. 2019). Root dynamics and production are reported to have key roles in plant adaptation to rising [CO2] (Benlloch-Gonzalez et al. 2014). Several studies have shown that elevated [CO2] resulted in an increased root growth rate in cereals, including wheat (Wechsung et al. 1999; Uddin et al. 2018). Nowadays, free-air CO2 enrichment (FACE) experiments have been prioritized worldwide for evaluating crop responses to the future climate (Kirkham 2011), as they provide more realistic field conditions compared to enclosed systems such as open-top chambers (Broberg et al. 2019). As the small plot sizes and the huge investment and maintenance costs arising in FACE systems necessitate minimizing soil disturbance and destructive plant sampling during the growing season, the use of high-throughput, in situ root methods is definitely preferable.

Therefore, a FACE experiment (“MartonFACE”) was chosen to test the efficiency of monitoring CR* to evaluate the response of wheat to an altered environment, such as different N fertilization rates and elevated [CO2], non-intrusively. The specific aims were: (i) to determine the temporal pattern of CR* in relation to plant phenology stages in two winter wheat cultivars from different maturity groups over two growing seasons; (ii) to investigate the effect of N level, elevated [CO2] and their interaction on CR* during the crop growth period; and (iii) to study the relationship between CR* and leaf area index (LAI) detected non-destructively on several occasions up till wheat flowering. Considering that CR* is likely influenced by shoot traits, and LAI is a determinant of root activity (water use) via whole-plant transpiration, CR*–LAI regressions served to indirectly validate the capacitance method. Furthermore, correlations were determined between: (iv) the maximum CR* measured at anthesis and the total aboveground biomass (TAB) and grain mass (GM) per plant measured at maturity for both cultivars at the plant scale; and (v) the maximum CR* and shoot dry mass (SDM; at the end of flowering) and GY (at maturity) across treatments and years at stand scale. Overall, the aim was to provide a useful tool to assess plant performance and yield in the field for, among other things, the selection of crop genotypes with improved adaptability to changing conditions.

Materials and methods

Site description

The field trial was carried out in the nursery of the Centre for Agricultural Research, Martonvásár, Central Hungary (N 47°19’, E 18°47’, 110 m asl.) over the 2018–19 and 2020–21 growing seasons (2019 and 2021 hereafter). The soil is classified by FAO-WRB (IUSS Working Group 2015) as a Haplic Chernozem (34% sand, 42% silt and 24% clay in the 0–25 cm layer), with a pHH2O of 7.59, 1.84% CaCO3, 3.39% humus, 1799/374/429 mg kg–1 total N/P/K, and 0.322 and 0.476 cm3 cm–3 water content at field capacity and saturation point, respectively.

The climate is continental with a mean (1988–2017; recorded by an on-site weather station) annual temperature of 10.9 °C (January: –1.2 °C, June: 21.2 °C) and total rainfall of 552 mm with 193 mm falling during the main winter wheat growing season (March–June; Fig. 1).

Fig. 1.
figure 1

Monthly rainfall (mm; columns) and mean air temperature (°C; lines) in the nursery of the Centre for Agricultural Research (Martonvásár, Hungary) during the winter wheat growing seasons. Long-term (1988–2017) average is displayed for reference.

Crop cultivation and FACE system

A factorial experiment was set up in three replicates with (1) two winter wheat (Triticum aestivum L.) cultivars: the early maturing Mv Nemere and the medium-early Mv Dandár, (2) two levels of N supply: low (80 kg ha–1) and high (160 kg ha–1), and (3) two levels of [CO2]: ambient and elevated to ~600 ppm. A split-plot arrangement was used, where the [CO2] treatment was in the main plot and the wheat cultivar by N combination in the subplots. A moldboard plow was used for tillage to a depth of 25 cm. The soil was manually fertilized with ammonium nitrate a week before planting (30% of the total doses) and in early spring (70%). Wheat was sown in mid-October 2018 and 2020 at a density of 500 seedlings m–2 with a row spacing of 12.5 cm, using a standard 3 m wide seeder to ensure uniform row and plant spaces for the whole area. The crop stand was sprayed with a pesticide combination at the 3-leaf stage in April.

The FACE system was engineered by the Institute for Biometeorology, Italian National Research Council, Florence (for details, see Miglietta et al. 2001). The facility consisted of three 18-m diameter octagonal rings (~250 m2) of horizontal tubes releasing pure CO2 on the upwind side to a targeted [CO2] of 600 ppm. The fumigation of the small amount of CO2 at high velocity through a large number of small gas jets allows open-air elevation of [CO2] without altering the microclimate (Ainsworth and Long 2020). A GMP343 type sensor (Vaisala Co. Ltd., Helsinki, Finland) was installed centrally in each ring to monitor [CO2] and control the venting. The FACE rings were installed right after wheat planting. The plots were fumigated throughout the whole vegetation season (from crop emergence to the fully ripe stage) every time when the plants were photosynthetically active, i.e. during the daylight hours with ambient temperature above 0 °C. The tubes were kept at a height of 0.1–0.2 m above the wheat canopy. The average [CO2] in the FACE rings during the treatment period was 597 and 587 ppm in 2019 and 2021, respectively, and was within 600 ppm ± 10% for 76.2% and 69.5% and within 600 ppm ± 20% for 94.2% and 90.8% of the operational time in the given years. Three further rings without CO2 enrichment were established as control (ambient) plots.

Monitoring of root electrical capacitance and leaf area index

CR was measured eight times between the 2-node stage (late March) and the over-ripe stage of wheat (early July) to monitor the seasonal root dynamics. The phenology stages were documented using the BBCH scale (Meier 2001). On each occasion, twelve plants were randomly selected from various rows in the centre of each plot (36 plants per treatment, 288 plants in total), and new plants were chosen on the next sampling date. SWC was measured at a depth of 0–12 cm in the root zone of each sample plant with a calibrated CS620 portable TDR meter (Campbell Sci. Ltd., Loughborough, UK), inserting the sensor vertically in the interrow, 6 cm away from the plant.

Thereafter, CR was measured for each plant in a parallel circuit at 1 kHz, 1 V AC with a U1733C handheld LCR instrument (Agilent Co. Ltd., Penang, Malaysia). One kHz is widely considered to be the optimum operating frequency, which induces an efficient electrostatic energy storage, and eliminates electrode polarization effects and stray capacitances (Ozier-Lafontaine and Bajazet 2005; Středa et al. 2020). A stainless steel rod, 15 cm in length and 6 mm in diameter, was used as a ground electrode, pushed vertically into the soil to a depth of 12 cm, 6 cm distance from the plant (in the place of the TDR probe). The plant electrode was clamped to all basal parts of the plant 15 mm above the soil (Svačina et al. 2014), after smearing them with conductivity gel to ensure good electric contact. The measured volumetric SWC values were divided by the saturation water content to obtain the relative water saturation (θrel). The saturation capacitance was calculated for each of the 288 θrel–CR data pairs as CR* = CR × 5.807e–1.115θrel to eliminate the SWC effect. The procedure for calculating the exponential function was previously reported in detail by Cseresnyés et al. (2018).

LAI was measured non-destructively in each plot five times in each growing season, on the same dates as the CR measurements (or 1–2 days later due for weather reasons). The monitoring terminated at the wheat flowering stage (mid-May), when LAI peaked but the reading was not yet affected by leaf senescence (Pokovai and Fodor 2019). LAI was detected during clear midday hours (from 10 a.m. to 2 p.m.) using an LP-80 handheld ceptometer (Meter Group Inc., Pullmann, WA, USA). Each LAI value, calculated as a mean of 22 readings, represented a ~0.6 m2 area under the canopy. The 80 cm long probe was placed parallel and perpendicular to the crop rows, and was read twice in each position.

Plant sampling and harvest

In 2019 destructive plant sampling had to be postponed from anthesis to the end of flowering owing to bad weather conditions, so sampling was also carried out in this stage in 2021 to ensure comparability. A 0.5 m long row was manually cut just above the soil surface in each plot. After drying the samples at 70 °C, SDM was determined (±0.01 g) and expressed as t ha–1 values.

Plants from the central 2 m × 6 m area of each plot were harvested and threshed at grain maturity (mid-July 2019 and 2021) using a plot combine. GY was determined in t ha–1 on a dry weight basis.

Considering that only one GY value was obtained per plot (though they represented a relatively large area), a plant scale investigation was also carried out to confirm the relationship of CR* to TAB and GM. On the day in 2019 when CR* measurements were taken at the flowering stage (day of year, DOY 140), the 36 sample plants in the “NL(–)” and “DL(–)” treatments were individually tagged at the stem base. These plants were cut separately at maturity, after which TAB was weighed (±0.001 g) after drying at 70 °C. The spikes were hand threshed to record GM.

Statistical analysis

The data were analysed in R programming language (R Core Team 2021). At first, two-way ANOVA was used to examine the effect of N (two levels: low and high) and CO2 (two levels: ambient and elevated), as categorical variables, and their interactions N × CO2 on quantitative variables of CR*, LAI, SDM and GY within each year for both cultivars. Thereafter, multivariate ANOVA was performed to evaluate the effect of year, N, CO2, as categorical variables, and their interactions on the maximum CR*, maximum LAI, SDM and GY, as quantitative variables for both cultivars. The normality, and the equality of variances in the data groups were examined with the Shapiro–Wilk test and Levene test, respectively. Statistical significance was assessed at p < 0.05 in each case.

Linear regression models were compiled: (i) to evaluate the relationships between CR* and LAI (across the N and CO2 treatments), and between the maximum CR*, SDM and GY (across the N and CO2 treatments and years) at a stand scale; and (ii) to relate CR* to TAB and GM at a plant scale for the cultivars. The line slopes and y-intercepts were compared using ANOVA. The reason for applying the maximum CR* (measured at flowering) was that it previously proved to be the best predictor of yield (Chloupek et al. 2006).

Results

Weather conditions

2019 was 1–1.5 °C warmer and 55% drier from January to April than the long-term mean. In contrast, May was very wet with 139 mm (239%) precipitation, followed by a relatively normal June. In 2021, April and May were very cold (3 °C below the average), and March and June were extremely dry with only 3 mm (10%) and 4 mm (6%) total rainfall, respectively (Fig.1).

Root electrical capacitance

In 2019 the mean CR* increased sharply during vegetative growth until the flag-leaf stage (BBCH 37–39; DOY 113) and then moderately until flowering (BBCH 61–69; DOY 140), when it reached a maximum of 10.6–12.7 nF, finally decreasing during the maturity stages (Fig. 2). The higher N level caused a significant increase in CR* from the flag-leaf stage to the milk stages (BBCH 73–75; DOY 154). At flowering this increase was 2–8% for Mv Nemere and 4% for Mv Dandár (Table 1). The effect of elevated [CO2] proved to be significant from the 3-node stage (BBCH 33; DOY 88) to the dough stages (BBCH 83–87; DOY 169), enhancing CR* at flowering by 12–18% and 12–16% for Mv Nemere and Mv Dandár, respectively. A significant N × CO2 interaction was only found for Mv Nemere at DOY 113.

Fig. 2.
figure 2

Changes in the saturation root electrical capacitance (CR* in nanofarads, nF) over time (DOY – day of year) for wheat cultivars Mv Nemere (˝N˝) and Mv Dandár (˝D˝) grown under low (˝L˝) or high (˝H˝) nitrogen supply under ambient (˝–˝) or elevated (˝+˝) [CO2] in two years. Vertical bars represent SDs (n = 36).

Table 1. Two-way ANOVA results for the effect of nitrogen (N), CO2 and N × CO2 on plant parameters for wheat cultivars Mv Nemere and Mv Dandár in two years. DOY – Day of year, CR* – Saturation root electrical capacitance, LAI – Leaf area index, SDM – Shoot dry mass, GY – Grain yield. (*)p < 0.1, *p < 0.05, **p < 0.01, ***p < 0.001, NS not significant.

In 2021 a severe drought in March and a prolonged cold period led to a considerable delay in wheat development and postponed anthesis by 8–10 days. This was also indicated by the capacitance data: initially, CR* showed a slower increase during vegetative growth, but thereafter it reached a maximum of 10.5–13.5 nF for Mv Nemere and 10.8–13.0 nF for Mv Dandár at flowering (DOY 145). Thus, there was no significant year effect on the maximum CR* for the cultivars (Table 2). The higher N rate positively influenced the root activity from the 3-node stage (DOY 99) to the late milk stage (BBCH 77; DOY 174). In the case of Mv Nemere, a 9–10% increase in CR* was observed at flowering. For Mv Dandár, however, the maximum CR* remained unchanged (0%) or increased by 20% under ambient and elevated [CO2], respectively. Elevated [CO2] significantly influenced CR* during most of the growth cycle for Mv Nemere with a 17–18% increase at flowering. The effect was weaker for Mv Dandár with no change (0%) and a 19% increase in CR* in the low and high N treatments, respectively. This cultivar showed a significant N × CO2 effect at several measurement dates.

Table 2. Multivariate ANOVA results for the effect of year (Y), nitrogen (N) and CO2 and their interactions on the maximum saturation root electrical capacitance (CR*), maximum leaf area index (LAI), shoot dry mass (SDM) and grain yield (GY) for wheat cultivars Mv Nemere and Mv Dandár. (*)p < 0.1, *p < 0.05, **p < 0.01, ***p < 0.001, NS not significant.

Leaf area index

In 2019 LAI steadily increased over time until wheat flowering, up to 5.6–7.7 m2 m–2 for Mv Nemere and 4.7–7.3 m2 m–2 for Mv Dandár (Fig. 3). The higher N supply increased the maximum LAI by 11–18% and 16–22% for these cultivars. However, differences in LAI between the N treatments were only marginally significant (p < 0.1) in some cases, likely partly due to the small number of replicates (n = 3; Table 1). The elevated [CO2] significantly increased LAI during the vegetative growth (DOY 88 and 99) for both cultivars, and also at the heading (DOY 128) and the flowering stages for Mv Dandár. The maximum LAI proved to be 18–25% higher for Mv Nemere and 27–33% higher for Mv Dandár under [CO2] enrichment.

Fig. 3.
figure 3

Changes in the leaf area index (LAI) over time (DOY – day of year) for wheat cultivars Mv Nemere (˝N˝) and Mv Dandár (˝D˝) grown under low (˝L˝) or high (˝H˝) nitrogen supply under ambient (˝–˝) or elevated (˝+˝) [CO2] in two years. Vertical bars represent SDs (n = 3).

Very small LAI was measured early in the season in 2021 due to the slow shoot growth rate. Thereafter, this was partly compensated for by more vigorous development up to an LAI of 3.0–4.3 m2 m–2 for Mv Nemere and 3.1–6.2 m2 m–2 for Mv Dandár. The maximum LAI was significantly lower in 2021 than in 2019, by 44–49% and 16–46% for Mv Nemere and Mv Dandár, respectively (Table 2). Enhanced N supply resulted in a significant increase in LAI in most phenophases, with the relative effect of 2% (Mv Nemere) and 24% (Mv Dandár) under ambient, and 23% and 79% under elevated [CO2] at flowering. The positive impact of CO2 treatment on LAI proved to be significant in most cases. The effect on the maximum LAI was 18% (Mv Nemere) and 9% (Mv Dandár) in the low N, and 42% and 60% in the high N treatments. A significant N × CO2 effect was observed chiefly for Mv Dandár.

A strongly significant positive linear relationship was found between CR* and the corresponding LAI for each cultivar and year (R2: 0.897–0.962; p < 0.001; n = 20) when the data were pooled across the N and CO2 treatments (Fig. 4). The regression parameters were statistically equal for the cultivars in 2019, whereas Mv Nemere had a steeper line (p < 0.001) in 2021. The slope of the regressions differed significantly between years for Mv Nemere (p < 0.001) and the y-intercept for both cultivars (p < 0.05).

Fig. 4.
figure 4

Relationships between the saturation root electrical capacitance (CR*) and leaf area index (LAI) for wheat cultivars Mv Nemere (˝N˝) and Mv Dandár (˝D˝) in two years. For each cultivar, data were pooled across low and high nitrogen and ambient and elevated [CO2] treatments (n = 20). Regressions were significant at the p < 0.001 level.

Shoot dry mass and grain yield

Wheat SDM ranged from 14.9 to 23.5 t ha–1 for Mv Nemere and from 14.4 to 20.6 t ha–1 for Mv Dandár in 2019, and from 11.7 to 17.3 t ha–1 and 8.2 to 14.0 t ha–1 in 2021 (Fig. 5). SDM was significantly lower in 2021 than in 2019 (Table 2), by 5–26% and 18–56% for Mv Nemere and Mv Dandár, respectively. There was no significant N effect on SDM, except for a marginal (p < 0.1) increase for Mv Dandár in 2021 (Table 1). Elevated [CO2] led to a significant increase in SDM (but only marginally for Mv Dandár in 2019): the relative change ranged from –1% to 58%, with higher responses under low N conditions. A significant N × CO2 effect was shown for the cultivars in 2019.

Fig. 5.
figure 5

Shoot dry mass (SDM) and grain yield (GY) of wheat cultivars Mv Nemere (˝N˝) and Mv Dandár (˝D˝) grown under low (˝L˝) or high (˝H˝) nitrogen supply under ambient (˝–˝) or elevated (˝+˝) [CO2] in two years. Vertical bars represent SDs (n = 3).

A GY of 4.75–7.47 t ha–1 and 4.81–6.80 t ha–1 was produced by Mv Nemere and Mv Dandár, respectively, in 2019, and 4.22–6.98 t ha–1 and 5.24–7.66 t ha–1 in 2021 (Fig. 5). GY for Mv Nemere was marginally significantly (p = 0.05) lower in 2021 than in 2019; the difference was insignificant for Mv Dandár. The N effect was significant only for Mv Dandár in 2021, when 2% and 46% yield increases were observed under ambient and elevated [CO2], respectively. The [CO2] enrichment significantly increased GY for the two cultivars in 2019, and for Mv Nemere in 2021, with more pronounced differences in the low N (41–57%) than in the high N treatment (12–37%). As for Mv Dandár in 2021, GY was reduced by 9% but enhanced by 30% in the CO2 treatments under low and high N supply, respectively. Significant N × CO2 effects were found for Mv Nemere in 2019 and for Mv Dandár in 2021.

The relationship between the maximum CR* and SDM proved to be insignificant (Mv Nemere: R2: 0.423; Mv Dandár: R2: 0.187), when the data were pooled across N and CO2 treatments and years (Fig. 6). In contrast, the linear correlations between the maximum CR* and GY were significant both for Mv Nemere (R2: 0.805; p < 0.01) and Mv Dandár (R2: 0.867; p < 0.01), with equal regression parameters for the cultivars.

Fig. 6.
figure 6

Relationships between the maximum saturation root electrical capacitance (CR*) and shoot dry mass (SDM) and grain yield (GY) for wheat cultivars Mv Nemere (˝N˝) and Mv Dandár (˝D˝). For each cultivar, data were pooled across low and high nitrogen and ambient and elevated [CO2] treatments and years (n = 8). Horizontal and vertical bars represent SDs. The CR*–SDM regressions were insignificant, the CR*–GY regressions were significant at the p < 0.01 level.

Plant-scale correlations

Considering the individual wheat plants, CR* measured at anthesis was significantly correlated with TAB (R2: 0.715 and 0.727; p < 0.001) and GY (R2: 0.648 and 0.661; p < 0.001) for each cultivar (Fig. 7). No significant differences in line parameters were found between the cultivars for the CR*–TAB and CR*–GY regressions.

Fig. 7.
figure 7

Relationships between the saturation root electrical capacitance (CR*) measured at the flowering stage and the total aboveground biomass (TAB) and grain mass (GM) per plant for wheat cultivars Mv Nemere (˝N˝) and Mv Dandár (˝D˝) grown under low nitrogen level and ambient [CO2]. Regressions were significant at the p < 0.001 level (n = 36).

Discussion

Seasonal patterns

The monitoring of CR* demonstrated marked seasonal dynamics in wheat root development with a peak at the flowering stage, as well as delayed early vegetative growth in 2021, linked to substantially smaller LAI. Similar temporal variations in wheat root biomass, root length and water use were previously shown in the field using the minirhizotron (Asseng et al. 1998; Uddin et al. 2018) or soil core methods (Xue et al. 2003; Jha et al. 2017). The root pattern is related to ontogenic changes in LAI and whole-plant transpiration (Yang et al. 2018). In cereals this trend is caused by strong resource allocation to the roots until flowering and subsequent nutrient translocation for reproductive growth (Fageria 2013). The sudden decrease in CR* after the flowering stages can be explained by reduced relative permittivity in aging roots owing to tissue maturation, lignification and senescence (Dalton 1995; Ellis et al. 2013).

The shorter growth period (6–8 days earlier flowering) of Mv Nemere compared to Mv Dandár was indicated by a moderate increase in CR* between DOY 128 and 140 in 2019 (Fig. 2). Elevated [CO2] hastened the flowering of the cultivars by 3–4 days in 2019 but not in 2021. This is in line with reports that CO2 treatment generally had a small and inconsistent effect on wheat phenology (Kirkham 2011; Wang et al. 2013).

Responses to higher N level and elevated [CO2]

In the present study, the positive effect of elevated N and [CO2] on crop growth and their interaction (for Mv Dandár in 2019) were clearly indicated by the higher CR* values, confirmed by increased LAI, SDM and GY. Rising [CO2] stimulates photosynthesis and inhibits photorespiration in C3 plants, and also improves water use efficiency (WUE) by reducing stomatal conductance (Wang et al. 2013; Dubey et al. 2015; Manderscheid et al. 2018). These physiological responses often lead to higher above- and belowground productivity, plant height, LAI and GY, usually to a greater extent under adequate N supply and drought conditions (Amthor 2001, Cai et al. 2016). Increased leaf area and reduced transpirational cooling, however, may counterbalance the improved WUE, increasing the total transpiration and in turn root uptake activity (Manderscheid et al. 2018). Wheat root growth may be stimulated even more than shoot growth by rising [CO2], resulting in a higher root/shoot ratio (Benlloch-Gonzalez et al. 2014).

The capacitance technique offered an opportunity to evaluate the treatment effects dynamically. When applying the root core method in a FACE system, Wechsung et al. (1999) revealed that the response of wheat root growth to elevated [CO2] was initiated at the 3-leaf stage, and was the largest at stem elongation. These findings agree with the present non-intrusive CR* monitoring and the concurrent crop LAI data. It is worth noting that although CR* was closely related to LAI for each cultivar and year (90–96% explained variance; Fig. 4), this could be at least partially explained by the wide range of phenology stages involved in data evaluation, and indicated the strong allometry between root and shoot traits in wheat (Bektas et al. 2016). Furthermore, ceptometer-based LAI may be somewhat overestimated at very large LAI values (Casa et al. 2019), such as those observed in 2019, perhaps contributing to an enhanced difference in regressions between years.

Relationship between electrical capacitance and yield

The significant relationships observed between CR* (at anthesis) and TAB or GM (65–73% explained variance) demonstrated the predictive power of capacitance for the aboveground biomass production and yield at a plant scale (Fig. 7). Chloupek et al. (2006) found a significant correlation between CR directly measured and shoot biomass (R2: 0.42) for barley. The correlations previously shown between the capacitance-based RSS and GM values of field-grown cereals were significant in dry environments (R2: 0.21–0.63), but not under optimal water conditions (R2: 0.11–0.14) in many cases (Chloupek et al. 2010; Středa et al. 2012; Svačina et al. 2014; Heřmanská et al. 2015). In the present study, a significant relationship was also obtained between CR* and GY at a stand scale using two years’ data (R2: 0.805 and 0.867). A weaker correlation (R2: 0.61) was established previously between directly measured CR and GY for barley across four varieties and over four years (Chloupek et al. 2010). In the present case, the higher R2 values implied that the incorporation of SWC into the measured capacitance improved the predictive accuracy of the method.

Several former studies reported that GY was positively related to wheat RSS traits and the water uptake rate measured at flowering, but these findings were based on the use of destructive techniques (Muñoz-Romero et al. 2010; Man et al. 2016) and minirhizotrons (Postic et al. 2019), or were performed in rhizoboxes (Fang et al. 2017). Although Mv Nemere and Mv Dandár presented statistically similar regressions, this is always so when cultivars possessing different phenotypic characteristics and root architecture are compared (Středa et al. 2020; Cseresnyés et al. 2021).

Overall evaluation

In a given year, the responses of wheat to the N and CO2 treatments were clearly indicated by CR*, in accordance with changes in biomass and yield. Comparing the two years, however, the maximum CR* values were equivalent, and GY decreased only slightly for Mv Dandár despite the substantially reduced maximum LAI and SDM for both cultivars in 2021. These findings were supported by the different CR*–LAI relationships between years, as well as by the insignificant and significant CR*–SDM and CR*–GY correlations, respectively, established for the two-year data (Fig. 6). As destructive root quantification was unfeasible in the present FACE study to further validate the CR* data, this result is hypothesized to be due to an enhanced root/shoot ratio in response to the severe spring drought in 2021. The increased biomass allocation to the roots was later able to compensate for the yield loss under the more favourable growing conditions that prevailed around flowering. Furthermore, CR* was unaffected by the reduction in SDM, implying that the capacitance is more or less determined by the root traits and not merely by the root-neck and stem-base properties. This corroborates previous findings by Gu et al. (2021), who reported poorer correlations between stem diameter and CR under relatively dry soil conditions when the current used for measurement penetrated deep into the root system. In a pot experiment, Wu et al. (2016) also observed a stronger relationship between CR and RSS traits before than after surface irrigation. It should be noted that though CR* values were calculated for water-saturated soil, the capacitance measurements in the present study were taken under much drier conditions (θrel from 0.2 to 0.5) in most cases.

Conclusions

The present study convincingly demonstrated the potential of measuring CR* to reveal the seasonal pattern of root development and uptake activity, and to assess the response of root and shoot growth and seed production to various treatments in the field. The non-destructive character of the method enables it to evaluate environmental effects over time using the same plant population (or even the same individual plants) for repeated measurements. Moreover, the approach proved to be useful for the early estimation of changes in GY across years. Nevertheless, more research is required to strengthen the link between CR* measured at flowering and GY at a stand scale over many years, and to further validate the capacitance results via the concurrent investigation of changes in root biomass and root/shoot ratio (even destructively outside the FACE area).

While it is true that neither the underlying biophysical principles nor the current pathways between the electrodes are as yet fully understood, the capacitance method should be considered to be a rapid, inexpensive approach for the dynamic evaluation of root responses to environmental changes without damaging the plant or disturbing the soil. As CR* represents a functional root size, it may provide an adequate prediction of GY. One current challenge facing breeding programs is the selection of crop genotypes having higher physiological plasticity and responsiveness to changing climatic events, such as rising [CO2], extreme temperatures and drought. Capacitance measurement, as a simple, in situ, high-throughput root phenotyping tool could partially replace invasive routine field techniques, and may thus be of interest for future application.