Different Responses in Root Water Uptake of Summer Maize to Planting Density and Nitrogen Fertilization

Modifying farming practices combined with breeding has the potential to improve water and nutrient use efficiency by regulating root growth, but achieving this goal requires phenotyping the roots, including their architecture and ability to take up water and nutrients from different soil layers. This is challenging due to the difficulty of in situ root measurement and opaqueness of the soil. Using stable isotopes and soil coring, we calculated the change in root water uptake of summer maize in response to planting density and nitrogen fertilization in a 2-year field experiment. We periodically measured root-length density, soil moisture content, and stable isotopes δ¹⁸O and δD in the plant stem, soil water, and precipitation concurrently and calculated the root water uptake based on the mass balance of the isotopes and the Bayesian inference method coupled with the Markov Chain Monte Carlo simulation. The results show that the root water uptake increased asymptotically with root-length density and that nitrogen application affected the locations in soil from which the roots acquired water more significantly than planting density. In particular, we find that reducing nitrogen application promoted root penetration to access subsoil nutrients and consequently enhanced their water uptake from the subsoil, while increasing planting density benefited water uptake of the roots in the topsoil. These findings reveal that it is possible to manipulate plant density and fertilization to improve water and nutrient use efficiency of the summer maize and the results thus have imperative implications for agricultural production.

Introduction

Maize (Zea mays L.) is a staple crop that is economically significant and most produced ahead of wheat (Hodgkinson et al., 2017). Since the available croplands and water resources for agriculture have both been dwindling, maize production is facing an unprecedented challenge due to its high demand for water and nutrients. Improving soil resource use efficiency is hence critical to sustaining maize production to feed an increasing global population projected to reach 9 billion in 2050 (Giri et al., 2018).

Crop management has substantial effects on spatiotemporal availability of soil resources, as well as root growth and its acquisition of water and nutrients (Chapman et al., 2012; Peake et al., 2013; Guan et al., 2014). Modifying farming practices to coordinate roots to grow in a way that facilitates water and nutrient acquisition is thus attractive and has been studied intensively over the past two decades (Widdicombe and Thelen, 2002; Liu et al., 2010, 2021; Ramezani et al., 2011; Li et al., 2015; Prechsl et al., 2015; Wu et al., 2015; Piao et al., 2016; Wang et al., 2017b). Although the experimental results are not conclusive, a consensus is that cultivars with “deep, steep, and cheap” root traits are effective for accessing resources in the subsoil thereby achieving high yield under terminal or intermittent drought conditions, whereas under conditions without water and nitrogen stresses, enhancing proliferation of thin lateral roots in the topsoil benefits resource acquisition and yield (Lynch, 2013). This is corroborated by experiments over the past decade. For example, Wang et al. (2019) found that reducing nitrogen application proliferated fine roots of maize in the subsoil, thereby enhancing its water use efficiency when the crop was under moderate water stress, and Ma and Song (2016) reported that adjusting fertilizer application reshaped root-length distribution of maize and impacted its water uptake from the soil profile as a result.

Along with fertilization, planting density and tillage also alter root growth and change the way crops take up water and nutrients from soil (Li et al., 2017b; Fiorini et al., 2018; Shao et al., 2018; Zhang et al., 2018). The impact of farming practice on soil resource acquisition has been well-documented (Majdi and Andersson, 2005; Kou et al., 2017; Li et al., 2017a), but the role played by root adaptation and the associated uptake of water and nutrient is not well-established due to the difficulty of in situ root phenotyping in opaque soils (Whalley et al., 2017). Traditional root measurement involving soil excavation is not only labor-intensive and time-consuming, but it is also insufficient to quantify water uptake as root-length density is not necessarily proportional to soil resource acquisition (Meinzer et al., 2001). Measuring soil-dying is another method to quantify root uptake (Whalley et al., 2017), but it is reliable only when soil is relatively dry where water dynamics is predominantly induced by root water uptake. It is inadequate when soil is wet following irrigation or rainfall in which soil-water redistribution becomes significant.

The impact of planting density and fertilization on crop yield has been fairly studied, but there is a lack of understanding of their impact on root growth and the consequence for water and nutrient uptake (Testa et al., 2016; Li et al., 2021; Liu et al., 2021). For example, it remains elusive that to what extent a change in fertilization and planting pattern alters root growth and water and nutrient uptake from different soil layers (Goebel et al., 2015; Du et al., 2018; Penna et al., 2020). Irrigation is known to influence root water uptake patterns significantly (Yang et al., 2015a; Ma and Song, 2016; Wu et al., 2016, 2018b; Zheng et al., 2019), while in semi-arid regions, such as northern China, maize is often grown rain-fed. Improving planting density and fertilization and understanding their combined impact on root growth and water acquisition is paramount to safeguard maize production in these regions but has been overlooked. The purpose of this article is to bridge this knowledge gap. A 2-year field experiment has been conducted with different combinations of nitrogen (N) fertilization and plant patterns. During the experiment, we periodically measured soil moisture content, root-length density, as well as the concentration of δ¹⁸O and δD in soil water, plant stem, and precipitation, concurrently. Root water uptake from different soil layers was calculated based on the intersection method and mass balance of the isotopes coupled with the Bayesian inference method (Eggemeyer et al., 2009; Brooks et al., 2010; Zhang et al., 2011; West et al., 2012; Schwendenmann et al., 2015; Yang et al., 2015b; Ma and Song, 2016). We analyzed how N fertilization and planting patterns combined to modulate root water uptake from different soil layers, as well as the underlying mechanisms.

Materials and Methods

Experimental Site and Planting

The 2-year field experiment was conducted in 2016 and 2017 at the Experimental Station of Institute of Farmland Irrigation, Chinese Academy of Agricultural Science (CAAS), located at Qiliying (35⁰08′ N, 113⁰45′ E and 81 m altitude) in Xinxiang, Henan province, China. The mean annual rainfall and temperature (1951–2014) at the station are 578 mm and 14.3°C, respectively, and the average precipitation over the growing season of the maize (June to September) is 418 mm. The physical and chemical properties of the soil are given in Table 1. During the experiment, the groundwater table was >5 m below the ground surface and had negligible impacts on the crop.

TABLE 1

Table 1. Physical and hydraulic parameters of the soil at the experimental station.

The maize variety Denghai-605 (Zea mays L.) was used as the model plant and the experiment consisted of three planting patterns and two nitrogen applications, each having three replicates, on a number of 6 m× 18 m plots designed using the randomized complete block method. The seeds were sown on June 13, 2016, and June 10, 2017, and the associated harvest was on October 8, 2016, and October 4, 2017, respectively. The intra-row planting distance was 30 cm in all treatments, and the treatments differed in the inter-row spacing and pattern as shown in Figure 1A: constant 60 cm spacing (P1); alternate 40 cm and 70 cm spacing (P2); zig-zag pattern with alternate 40 cm and 70 cm spacing (P3). The planting densities associated with P1, P2, and P3 were 5.6 × 10⁴, 6.1 × 10⁴, and 7.3 × 10⁴ plants ha⁻¹, respectively. Each planting pattern had two N applications: 240 kg N ha⁻¹ (N₂₄₀) and 120 kg N ha⁻¹ (N₁₂₀) to cover the varying N applications used by the local farmers (Zhang et al., 2021); in each treatment, 50% of N fertilizer was broadcasted with 90 kg ha⁻¹ of P₂O₅ and 75 kg ha⁻¹ of K₂O as basal fertilizer over each plot, followed by immediate plowing to avoid nitrogen volatilization; the remaining 50% of N fertilizer was applied as top-dressing at early jointing stage in the middle of July. All plots were irrigated with 60 mm of water on June 18 in 2016 and June 15 in 2017, using sprinkler irrigation to facilitate seed germination and emergence.

FIGURE 1

Figure 1. (A) Sketch of the three planting patterns: P1 - Constant 60 cm inter-row spacing; P2 – alternate 40 cm and 70 cm inter-row spacing, and P3 - zig-zag pattern with alternate 40 cm and 70 cm inter-row spacing; the intra-row spacing is 30 cm in all treatments. (B) Changes in air temperature, reference evapotranspiration (ET₀) and precipitation in the 2-year experimental period.

Sampling

Plant and soil in each plot were concurrently sampled on June 24, July 14, August 3 and 22, September 12 and 30 in 2016, and on June 23, July 13 and 30, August 22, and September 10 and 27 in 2017. During each sampling day, three plants in each plot were randomly sampled first and the stem epidermis of each sampled plant was then immediately removed to avoid potential contamination by isotopically enriched water (Querejeta et al., 2007). The stem was cut into 4–5 cm segments before being frozen in a sealed vial for isotope analysis. Following the plant sampling, soil samples at the depths of 5, 10, 20, 40, 60, 80, and 100 cm were taken from three sites close (about 5 cm) to each sampled plant using a hand-driven auger. Soils taken from the same depths from the three sites were pooled, a half of the pooled soils was used for isotope analysis, and the remaining was used for soil moisture measurement using the gravimetric method.

The soil and plant samples for stable isotope analysis were kept in a freezer at −20°C. Prior to the analysis, water in the stem and the soil water was extracted using a vacuum extraction system (LI-2000, LICA, China) by applying a suction pressure for 0.5–1.5 h. The results revealed this had extracted more than 99.0% of the water in both soil and stem and was therefore deemed adequate to obtain the unfractionated water from the soil and stem for isotope analysis (West et al., 2006; Meißner et al., 2014).

The following day after the soil and plant sampling, further soil cores were taken down to 100 cm deep at 10-cm intervals close to each sampled plant using an auger with an internal diameter of 6.91 cm and an external diameter of 7 cm to measure root-length density. Each cored sample was transferred to a sieve (0.5 mm mesh), and the sieve was then suspended in a water-filled trough to wash the soil away at low pressure to leave the root segments only. The root-length density (RLD) in each sample was analyzed using the WinRHIZO Reg. 2007d (Regent Instrument Inc.) based on the images of the roots scanned using the EPSON PREFECTION_TM V700 Photo Flatbed Scanner at resolution of 6,400 dpi × 9,600 dpi. The RLD was calculated as the ratio of the lengths of all root segments in each core to the volume of the core (375.0 cm³). The root-length densities for samples taken from each soil profile were normalized by the total root lengths in all samples taken from the soil profile.

Precipitation was sampled from a funnel filled with a ping-pong ball inside to prevent evaporation. A polyethylene bottle was connected to the bottom of the funnel to collect raindrops. It was difficult to collect sufficient rainwater, as well as strong fractionation, when precipitation was <5 mm, and we hence only considered precipitation events which were >5 mm. The irrigation water was sampled for helping data interpretation. To prevent isotope fractionation caused by evaporation, both rain and irrigation water samples were stored in an air-tight container at 4°C prior to isotope analysis.

Isotope Analysis

All water samples were analyzed for deuterium (δ²H) and oxygen (δ¹⁸O). We compared the isotopes in the water samples by injecting the water into a high-temperature conversion/elemental analyzer (TC/EA) coupled with a Con-Flo III interface on a Delta V Advantage isotope ratio mass spectrometer (Thermo-Electron Corporation, Bremen, Germany). All water samples were calibrated and normalized to an internal laboratory water standard previously calibrated against the Vienna Standard Mean Ocean water (VSMOW, 0‰). The results were expressed as δ-value related to the VSMOW:

$\begin{array}{ll}\delta (\%0)=\left(\frac{R_{sample}}{R_{VSMOW}}-1\right)1000& (1)\end{array}$

where R is the ¹⁸O/¹⁶O or ²H/¹H ratio with the subscript “sample” and “VSMOW” referring to the value of the samples and of the Vienna Standard Mean Ocean water, respectively. The error of the measurements was ±1.4‰ for δ²H and ±0.1‰ for δ¹⁸O.

Evapotranspiration

Daily precipitation, average temperature, relative humidity, sunshine, and wind speed were recorded from an automatic weather station installed on the site (RX3000, HOBO, USA). The reference evapotranspiration, ET_o (mm d⁻¹), was calculated from the Penman–Monteith equation (Allen et al., 1998):

$\begin{array}{ll}E{T}_o=\frac{0.408\Delta (R_n-G)\gamma \frac{900}{T+273}{u}_2(e_s-e_a)}{\Delta +\gamma (1+0.34u_2)}& (2)\end{array}$

where R_n is the net radiation (MJ m⁻² d⁻¹), G is the soil heat flux density (MJ m⁻² d⁻¹), T (^oC) and u₂ (m s⁻¹) are the mean daily air temperature and wind speed at 2 m above the ground surface, respectively, e_s and e_a are the saturated and actual vapor pressure (kPa), respectively, Δ is the slope of saturation vapor pressure curve (kPa ^oC⁻¹), and γ is the psychrometric constant (kPa ^oC⁻¹).

Root Water Uptake

There is no isotopic fraction when water flows from soil into roots (Dawson and Ehleringer, 1991; Ehleringer and Dawson, 1992), and we hence use mass balance to track the origin of the isotopes in the plant stem. We divided the soil profile into N layers, assuming all water taken up by the roots came from these layers. During a time period of Δt, if the transpiration is S and the water taken up by the roots from i_th soil layer isS_i, from the mass balance we have

$\begin{array}{l}\begin{array}{c}\begin{array}{l}\delta D_{stem}={{{\displaystyle \sum }}^{\text{}}}_{i=1}^N{f}_i·\delta D_i,\\ {\delta }^{18}{O}_{stem}={{{\displaystyle \sum }}^{\text{}}}_{i=1}^N{f}_i·{\delta }^{18}{O}_i,\end{array}\end{array}\\ {{{\displaystyle \sum }}^{\text{}}}_{i=1}^N{f}_i=1,\\ f_i=S_i/S,& (3)\end{array}$

where δD_stemand ${\delta }^{18}{O}_{stem}$are the concentration of H² and O¹⁸ measured from the stem water, respectively, δD_iand ${\delta }^{18}{O}_i$ are their counterparts in the ith soil layer, and f_i is the fraction of the water taken up by roots in the ith soil layer. If the two isotopes are not correlated and the soil profile is divided into three layers, Equation (3) can be solved exactly. However, given the irregular distribution of isotopes over the soil profile, to be consistent with the root density measurement, we divided the soil profile into four uneven layers: 0–5 cm, 5–20 cm, 20–60 cm, and 80–100 cm, calculating the water uptake of roots in each layer. There are more variables than the number of equations, and we hence solved Equation (3) using the Bayesian inference method, which gives the most likelihood rather than the exact amount of water taken up by roots in each layer.

Given a set of measurements Y, the Bayesian formalism postulates the problem as follows:

$\begin{array}{ll}p(\text{f}|\text{Y})=\frac{p(\text{f})p(\text{Y}|\text{f})}{p(\text{Y})}\propto p(\text{f})L(\text{f}|\text{Y}),& (4)\end{array}$

wherep(f|Y)andp(f)represent the posterior and prior distributions of the normalized root uptake from each soil layer, p(Y|f)andp(Y)are the posterior and priori distributions of the measured isotopes, andL(f|Y) is the likelihood function, respectively; L(f|Y) was estimated using the measured data as follows assuming that the measured isotopes were normally distributed

$\begin{array}{ll}L(\text{f}|\text{Y})={\displaystyle \prod _{t=1}^N}\frac{1}{\sqrt{2\pi {\sigma }_t^2}}exp\left[\frac{1}{2{\sigma }_t^2}{(\frac{y_t-y_t(\text{f})}{{\sigma }_t^{}})}^2\right],& (5)\end{array}$

wherey_tandy_t(f)are the measured and calculated isotopes in the plant stem using Equation (4), respectively, and${\sigma }_t^2$ is the variance. The posterior distribution of the root uptake from all soil layers defined in Equation (5) was calculated numerically based on the Monte Carlo Markov chain simulation, using the open-source SIAR software [CRAN—Package siar (r-project.org)] (Parnell et al., 2013).

Statistical Analysis

The Kolmogorov–Smirnov (K-S) test showed that the isotopes measured for the plant stem and soil water were normally distributed and these distributions were used to calculate the above likelihood function. We used hierarchical cluster analysis to classify the soil water and the analysis of variance (ANOVA) to investigate the difference in root water uptake, δ¹⁸O and δD in water, as well as the difference between treatments with p < 0.05 deemed significance. Multiple comparisons were made using the least significant difference (LSD) to quantify significant variation in δ¹⁸O and δD in the soil profile. All statistical analyses were performed using the SPSS 21.0 program.

Results

Environmental Factors and Isotopes in Precipitation

The changes in precipitation, air temperature, and ET_o during the experimental period are given in Figure 1B. The daily average ET_o and temperature were 3.6 mm d⁻¹ and 26.4°C in 2016 and 3.5 mm d⁻¹ and 25.4°C in 2017, respectively. The seasonal rainfall in 2016 and 2017 was 438.9 mm and 319.6 mm, respectively, during the experimental period; there were 15 and 17 precipitation events with rainfall >5 mm in 2016 and 2017, accounting for 95% and 92% of the precipitation in each year, respectively. During the experiment, the δD measured from the precipitation ranged from −64.34 to −16.83% with a mean of −45.37%, and the δ¹⁸O varied from −9.39% to −3.25% with a mean of −6.92% (Figure 2). The local meteoric water line (LMWL) measured from the precipitations was approximately linear: δD = 6.7935δ¹⁸O+2.0593 (R² = 0.88) for 2016 and δD = 6.7818δ¹⁸O+1.2127 (R² = 0.91) for 2017.

FIGURE 2

Figure 2. The relationship between δ¹⁸O and δD for data measured from all water samples in comparison with the global meteoric water line (GMWL) in 2016 (A), and 2017 (B), as well as the linear fittings to stem water and water in the 10–40 cm of soil in all treatments respectively.

Isotopes in Soil Water and Stem Water

The isotopes measured from the stem water varied with time significantly in all treatments (Figure 2 and Supplementary Figures 1, 2), implying an impact on the environment. The relationship between δ¹⁸O and δD was fitted to δD = 6.4326δ¹⁸O-5.7838 (R² = 0.95) and δD = 5.0586δ¹⁸O-17.289 (R² = 0.89) for 2016 and 2017, respectively. Precipitation and irrigation altered isotopes in soil water, with the measured δ¹⁸O declining along the soil depth (Supplementary Figures 1, 2). The δ¹⁸O and δD in soil water fall into two groups, each following its own linear relationship (Figure 2). Samples taken from the 10–40-cm soil layer with δ¹⁸O>-10% were consistent with the δ¹⁸O -δD relationship for the stem water, suggesting that most water taken up by the roots originated from this soil layer. In contrast, data with δ¹⁸O <-10% for samples taken from 60- to 100-cm soil layer deviated from the δ¹⁸O -δD relationship for the stem water, indicating that they are unlikely the main source of the water respired by the crop. As a comparison, the linear fitting curve for the stem water and the water in the 10–40-cm soil layer is also plotted in Figure 2.

Root-Length Density and Root Water Uptake

The region in soil where the roots took up water can be estimated by plotting δ¹⁸O (or δD) measured from the stem water along the soil profile. Its intersections with the distribution of δ¹⁸O (or δD) in the soil water are the points around which the roots took up the water (Figures 3A,B, Supplementary Figures 1, 2). While uncertainties might be raised when there are more than one such intersection, these did not appear in our experiment. The roots took most of the water from soil around the depth of 10 cm at the seedling stage, and they then progressively moved downward to the depth of 40 cm during the flowering and harvesting stages (Figures 3A,B, and Supplementary Figures 1, 2).

FIGURE 3

Figure 3. Spatiotemporal changes in root uptake depth estimated from the intersection method in 2016 (A), and 2017 (B). The impact of different combinations of N fertilization and planting density on root water uptake depth estimated using the intersection method for all treatments over the 2-year experiment (C).

The intersection method illustrates the approximate regions from which the roots take up water. To quantify the contribution of each soil layer to the transpiration of the crop at different growing stages, we estimated the posterior distribution of the root water uptake from each layer. The posterior distribution for each soil layer is a probability distribution of the root water uptake, and we present the mean only as this is the most likely water that the roots in this layer had taken up (Figure 4).

FIGURE 4

Figure 4. Spatiotemporal changes in the normalized root water uptake rate under different treatments (see Figure 1A) calculated based on the mass balance of isotopes and the Bayesian inference method; 2016 (A), and 2017 (B).

In the early stage before the N top-dressing (13 July in 2016 and 15 July in 2017), about 60% of the transpired water emanated from the top 10 cm of soil, while following the top-dressing, the water uptake from the subsoil increased steadily. The timing in 2016 at which the water uptake from the subsoil (20–60-cm layer) exceeded that from the topsoil was in the late August for N240, while for N120 this was delayed to the early August (Figure 4A). For P1 and P2 in 2017, the difference in root water uptake between the topsoil and the subsoil was similar to that in 2016. For P3, root water uptake from the topsoil was higher than that from the subsoil during the whole growing period (Figure 4B); even at the late stage, the roots in the top 10 cm of soil contributed more than 30% of the total transpiration. The relative significance of water uptake from different soil layers varied between the 2 years, but the trend was the same: N top-dressing increased the water uptake from the subsoil, and when the planting density (P1 and P2) and N fertilization were both low, the water uptake rate of roots in the 20–60-cm soil layer was higher than that in the top 10-cm soil layer.

Water uptake of a root segment is proportional to the difference between water potential on the soil-root surface and in its xylem network. For a region where roots are sparse, the roots are unlikely to interfere with each other and water uptake from this region is hence proportional to root-length density. In contrast, for regions where the roots are dense, continuous water uptake drives distant water moving into the rhizosphere; when the influencing zone of the roots meets, the roots start competition for water. Therefore, in general, root water uptake in a soil profile increases with the root-length density asymptotically rather than linearly. At different growth stages, the root-length density in all treatments decreases approximately exponentially with the soil depth (Figure 5). Pooling the root water uptake calculated using the statistical method and the measured root-length density at different growing stages in all treatments, Figure 6 shows that the root water uptake rate (y) calculated using the statistical methods increased asymptotically with the normalized root-length density (S). The increase is fitted toy = kS/(A+S), with k = 234.24 and A = 207.72. The slight deviation from linear increase when root-length density is high (in the topsoil) indicates the existence of competition between roots in the topsoil for water (Figure 5).

FIGURE 5

Figure 5. Spatiotemporal variation in the normalized root density under different treatments (see Figure 1A) in 2016 (A), and 2017 (B).

FIGURE 6

Figure 6. The increase in the normalized root water uptake rate with the normalized root-length density for all treatments in the 2-year experiment.

N Fertilization and Planting Pattern Effects

To elucidate the combined impact of fertilization and planting density on root water uptake, we plotted in Figure 3C the depths estimated from the cross-points between the isotopes in soil profile and the stem water as shown in Supplementary Figures 1, 2 (called root water uptake depth hereafter) under different treatments. It is manifest that, except for a few points measured following the N top-dressing, the majority of the depths are below the 1:1 line, indicating that reducing N fertilization from 240 to 120 kg N ha⁻¹ increased water uptake from the subsoil. Linearly fitting the data for the same planting pattern reveals that the fitting for P3 (most densely planted) deviates from the 1:1 line more significantly, followed by P2 though the difference between P1 and P2 is not significant. These phenomena suggest that the impact of N fertilization and planting density on root water uptake is confounded. This is corroborated by the results calculated from the statistical method (Figure 4).

The impact of basal fertilization and planting density on root water uptake from different soil layers appeared to be minor in the early stage, and significant differences started to emerge after the N top-dressing (Figure 4). Reducing N application promoted root penetration to take more water from the subsoil (Figure 4). For example, for P1 in the middle of September 2016, the roots in the 80–100-cm soil layer contributed 30% of the respired water under N₁₂₀, while under N₂₄₀ the same soil layer contributed 25% of the transpired water (Figure 4A). For P3 (the highest planting density) combined with N₁₂₀, the critical depth below which the roots contributed <10% of the transpired water increased from 40 cm in the early August to around 100 cm at the end of September. In contrast, during the same period but under the combination of P3 and N₂₄₀, the critical depth below which the roots contributed <10% of the transpired water increased from 40 to 60 cm only (Figure 4A). The results in 2017 did not replicate those in 2016, but the trend of the response of root water uptake to planting pattern and N fertilization was the same (Figure 4B).

Discussion

Changes in Isotopes

The slope of the LMWL was smaller than that of the global meteoric water line δD = 8δ¹⁸O+10 (Craig, 1961), indicating that humidity change and secondary evaporation might have enriched the O¹⁸ in the precipitation in our experimental site (Araguás-Araguás et al., 1998). Most δ¹⁸O and δD in soil water were plotted beneath the LMWL (Figure 2), implying that soil water with its origin from the rainfall was likely to have undergone evaporation which enriched δ¹⁸O (Wang et al., 2017a). The significant difference was found in isotopes in soil water between the treatments (Supplementary Figures 1, 2) due to the effects of rainfall and root water uptake. The heavy rain (123 mm) on July 9, 2016, (Figure 1B) enriched the isotopes, with δ¹⁸O measured from soil water varied approximately 1-fold along the soil profile (Supplementary Figure 1). Soil evaporation enriched δ¹⁸O in the vicinity of the soil surface, and δ¹⁸O in the topsoil fluctuated seasonally due to the periodic precipitation (Tang and Feng, 2001; McCole and Stern, 2007; Dai et al., 2015; Sprenger et al., 2016; Wu et al., 2018a). Our results agreed with those found in Tang and Feng (2001) and Wu et al. (2018a), all showing that hydrological processes modulated the seasonal change of isotopes in the topsoil water.

The δD and δ¹⁸O in the precipitation and stem water were linearly correlated, but this correlation does not apply to soil water (Figure 2). The δD and δ¹⁸O in soil water were categorized into two groups: one for δ¹⁸O < −10‰ and the one for otherwise, with the data in each group following its own linear relationship. The data in the above-part of the graph are for soil samples taken mainly from the 0–60 cm of soil, which can be fitted to a linear relationship, while the data for the 0–5 cm of soil in the graph fall below the δD - δ¹⁸O line for the stem water due to δ¹⁸O enrichment by evaporation. We linearly regressed the data taken from the 10–40-cm soil layer and plotted the results in Figure 2; it is close to the regression curve for the stem water.

Soil evaporation enriched δ¹⁸O, thereby reducing the δD/δ¹⁸O ratio for soil samples taken from the regions proximal to the soil surface (Figure 2). In our results, soil water with low δ¹⁸O concentration was mainly for samples taken from the subsoil (60–100 cm). The likely mechanism is the attenuation of δ¹⁸O by molecular and hydrodynamic diffusion when it moves downward with the rainfall and irrigation water, which spreads the δ¹⁸O along the soil profile. Winter wheat–maize rotation is the main cropping system in the region, and before wheat harvest, the soil normally endures a prolonged drought which was likely to have enriched δ¹⁸O in soils proximal to the surface. In 2016, the 60 mm of irrigation after the seed drilling and the 123 mm of rainfall on July 9 (Figure 1B) combined to have leached the enriched δ¹⁸O below the depth of 100 cm as δ¹⁸O measured on July 14 was >-9‰ in most treatments, corroborated by the spatiotemporal changes in soil-water content (Figure 7). Although the enriched δ¹⁸O was likely to have moved upward since, driven by evaporation and root uptake, it remained in the soil deeper than 40 cm in most treatments (Supplementary Figures 1, 2), consistent with the spatiotemporal change in soil-water content (Figure 7) and the root water uptake (Figure 4). This is corroborated by the results in 2017 in which there was less rainfall and hence less water infiltration (Figure 1B). As a result, δ¹⁸O concentration in the subsoil is low and the δD/δ¹⁸O ratio is small compared to that in 2016. The δD - δ¹⁸O curve hence further deviates from the GMWL line (Figure 2B).

FIGURE 7

Figure 7. Spatiotemporal variations in soil moisture content under different treatments (see Figure 1A) in 2016 (A), and 2017 (B).

Planting Pattern Effects

The planting pattern modulated the local environment for crops to grow both above-ground and below-ground (Dodd et al., 1998; Gani et al., 2002; Dass et al., 2015). Root-length density and root water uptake interactively affect each other (Coleman, 2007; Zhao et al., 2018), both varying with planting pattern and fertilization (Figure 4). Increasing planting density appeared to have enhanced water uptake by roots in the topsoil, especially under high N fertilization after the jointing stage (Figure 4). For example, in 2016, roots under combination of N₁₂₀ and P1 took 25% of the transpired water from the soil around the depth of 60 cm, while under combination of N₂₄₀ and P3, this soil layer supplied only 18% of the transpired water (Figure 4). Such variation in root water uptake with planting density was also found in other treatments in 2017, although the significance in the variation varies (Figure 4).

Roots regulate their water uptake from different soil layers as a response to change in soil water and other environmental factors (Figure 7). In P3, for example, the roots took 63% of the transpired water from the top 0–20 cm of soil on June 26, 2016. After the rainfall on July 19, the roots in the same soil layer contributed 67% to the transpired water despite the significant decrease in fraction of the root lengths in this soil layer (Figure 5). This increase in water uptake accompanied by a decrease in relative root length indicates that the roots in the top 0–20-cm soil layer were likely to have been water-stressed before the rainfall.

Compared to P1 and P2, P3 increased the root-length density in the top 0–10 cm of soil by 1.10-fold and 1.18-fold, respectively (Figure 5), similar to those found by others (Guan et al., 2007; Loades et al., 2010; Li et al., 2011). The consequence of the increased root density in the topsoil due to the increased planting density is that it increased the root water uptake and competition between roots for water in the topsoil (Tan et al., 2010; Li et al., 2011). Although the increased root water uptake depleted the topsoil water quickly, we did not find noticeable compensation by deep roots to increase their water uptake from the subsoil (Figure 4), indicating that the topsoil water remained sufficient for roots to take up and that the increased planting density did not enhance root penetration (Figure 5). The likely reason is that high planting density increased leaf area index thereby reducing water loss from evaporation (Hodge, 2004), as manifest from the soil-water change (Figure 7); this is consistent with the findings of others (Rossatto et al., 2013; Ma and Song, 2016).

Maize cultivars with shallow-root traits proliferate roots in the topsoil (Yi et al., 2009; Ma and Song, 2016), thereby increasing their water uptake from the topsoil (Schenk and Jackson, 2005). Zhao et al. (2018) suggested to express root water uptake as a function of dry-root weight rather than root-length distribution, claiming that the latter was inadequate to explain the spatiotemporal variation in root water uptake as often observed in the field (Ehleringer and Dawson, 1992). Physically, water acquisition by a root is driven by the difference between water potential at the root–soil surface and in the xylem network, and it depends on both root architecture and soil properties. Root traits that do not contain information on soil are hence insufficient to quantity root water uptake. Water needs to pass through the root surface prior to moving into the xylems, and the root-length density is hence one of the best proxies of root architecture to quantify root water uptake although root diameter and root age might also play an important part (Hodge, 2004). This is also corroborated by our results that the root water uptake increased asymptotically, rather than linearly, with the root-length density (Figure 6). The mechanisms underlying, such asymptotic increases, are manifold, including the variation in root diameter (hence the root–soil interfaces), which is not accounted for by the root-length density. The most likely reason, however, is the competition between roots in the topsoil for water, which reduces the uptake of individual roots (Wu et al., 2014). This is also consistent with previous studies on water acquisition of winter wheat where the root uptake rate is proportional to root-length density only when root-length density is <1 cm/cm³ (Gregory et al., 1978; Zhang et al., 2020).

Effects of N Fertilization

Top-dressing N relieved root competition for N and hence enhanced root proliferation in the topsoil (Figure 5). On average, the roots in the top 20 cm of soil were longer under N₂₄₀ than under N₁₂₀ (Figure 5). The increased length of the shallow roots increased their water uptake from the topsoil (Figure 4). Depending on water availability in the topsoil, the roots regulate the ways they take up soil water. At the early growth stage, since the demand for water and nutrients was low and water and nutrients in the topsoil were sufficient, there was no noticeable difference in root water uptake between treatments (Figure 4). As crops grew and their demands for water and nutrients increased, the differences in root water uptake between treatments emerged (Figure 4).

As the planting density increases, the demand for water and nutrients increases. If one or two of them becomes limiting, the roots penetrate to access water and nutrients in the subsoil (Rogers and Benfey, 2015). Our results showed that when the planting density was the same, the roots took more water from the subsoil under low N fertilization (Figure 4). Denser planting density associates with high root-length density and needs more water, especially in the topsoil. The fact that the root water uptake from the subsoil is affected by N more than by the planting density implies that N was more likely the limiting factor in our experiment (Figure 4). For example, roots in the combination of P3 and N₂₄₀ took approximately 35% of the transpired water from the top 0–15 cm of soil (Figure 4), while the roots in the same soil layer took only 30% of the transpired after halving the N application (Figure 4). This is corroborated by the recent finding of Ma and Song (2016) where the roots tended to exploit water and nutrients from the deep soil when fertilizer was in deficiency. The results of Ma and Song (2016) showed that the root water uptake increased linearly with root-length density, indicating the absence of water stress and root competition. This differs from our experiment where there was a slight water stress in the topsoil, and the root water uptake hence increases asymptotically with the root-length density (Figure 6).

Developing sustainable agriculture in arid and semi-arid regions, such as northern China (Guan et al., 2015), requires improving the use efficiency of soil water and rainfalls (Fang et al., 2010). While various efforts have been made, the efficacy of manipulating planting density and fertilization to help achieve this goal has been overlooked. On average, maize planted in high density in the zig-zag pattern took more water from the topsoil than from the subsoil due to the proliferation of shallow roots, especially when combined with high N application (Figure 5). It thus improves rainfall use efficiency as it preferentially uses the topsoil water, consistent with other studies (Hatfield et al., 2001; Zhang et al., 2005).

Conclusion

Spatiotemporal change in root water uptake of summer maize was studied in a 2-year field experiment comprising three planting patterns and two N fertilizations. Water uptake by roots at different soil layers was calculated using stable isotopes δD and δ¹⁸O measured concurrently from soil water and crop stem. The results showed that the crop took most of the required water from the 0–60 cm of soil, but the uptake pattern varied with treatments. Regardless of the planting patterns, reducing N fertilization boosted root penetration to access nutrients in the subsoil and consequently increased the root water uptake from the subsoil. Increasing planting density and uniformity enhanced proliferation of shallow roots and their water uptake from the topsoil due to the increased leaf index area which reduced water loss from evaporation, especially in the early growth stage. In soil profile, because of the competition of shallow roots for water, root water uptake increased asymptotically rather than linearly with root-length density. In summary, high planting density combined with high N fertilization improves the preferential use of the topsoil water which is prone to evaporation and is hence more water-use-efficient for rain-fed maize production in semi-arid regions.

Data Availability Statement

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

Author Contributions

YG: funding acquisition and writing—original draft preparation. JC: investigation. GW: visualization. ZL: data curation and software. YZ and WS: investigation and data curation. XZ: methodology and writing—reviewing and editing. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the China Agriculture Research System (CARS-02) and the National Natural Science Foundation of China (51879267). The work at Rothamsted Research was supported by the United Kingdom Biotechnology and Biological Science Research Council (BBSRC)-funded Soil to Nutrition strategic program (BBS/E/C/000I0310).

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.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.918043/full#supplementary-material

References

Allen, R. G., Pereira, L. S., Raes, D., and Smith, M. (1998). Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements. United Nations FAO irrigation and drainage paper no.56, Rome, Italy.

PubMed Abstract | Google Scholar

Araguás-Araguás, L., Froehlich, K., and Rozanski, K. (1998). Stable isotope composition of precipitation over southeast Asia. J. Geophys. Res. Atmos. 103, 28721–28742. doi: 10.1029/98JD02582

CrossRef Full Text | Google Scholar

Brooks, J. R., Barnard, H. R., Coulombe, R., and McDonnell, J. J. (2010). Ecohydrologic separation of water between trees and streams in a Mediterranean climate. Nat. Geosci. 3, 100. doi: 10.1038/ngeo722

CrossRef Full Text | Google Scholar

Chapman, N., Miller, A. J., Lindsey, K., and Whalley, W. R. (2012). Roots, water, and nutrient acquisition: let's get physical. Trends Plant Sci. 17, 701–710. doi: 10.1016/j.tplants.2012.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, M. (2007). Spatial and temporal patterns of root distribution in developing stands of four woody crop species grown with drip irrigation and fertilization. Plant Soil. 299, 195–213. doi: 10.1007/s11104-007-9375-5

CrossRef Full Text | Google Scholar

Craig, H. (1961). Isotopic Variations in Meteoric Waters. Science. 133, 1702–1703. doi: 10.1126/science.133.3465.1702

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Zheng, X.-J., Tang, L.-S., and Li, Y. (2015). Stable oxygen isotopes reveal distinct water use patterns of two Haloxylon species in the Gurbantonggut desert. Plant Soil. 389, 73–87. doi: 10.1007/s11104-014-2342-z

CrossRef Full Text | Google Scholar

Dass, A., Chandra, S., Choudhary, A. K., Singh, G., and Sudhishri, S. (2015). Influence of field re-ponding pattern and plant spacing on rice root-shoot charactgeristics, yield and water productivity of two modern cultivars under SRI management in Indian Mollisols. Paddy Water Environ. 14, 45–59. doi: 10.1007/s10333-015-0477-z

CrossRef Full Text | Google Scholar

Dawson, T. E., and Ehleringer, J. R. (1991). Streamside trees that do not use stream water. Nature. 350, 335. doi: 10.1038/350335a0

CrossRef Full Text | Google Scholar

Dodd, M. B., Lauenroth, W. K., and Welker, J. M. (1998). Differential water resource use by herbaceous and woody plant life-forms in a shortgrass steppe community. Oecologia. 117, 504–512. doi: 10.1007/s004420050686

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, J., Ma, Y., Hu, X., Tong, J., Zhang, B., Sun, N., et al. (2018). Applying dual stable isotopes and a MixSIAR model to determine root water uptake of winter wheat. Acta Ecol. Sin. 38, 6611–6622. doi: 10.5846/stxb201802250384

CrossRef Full Text | Google Scholar

Eggemeyer, K. D., Awada, T., Harvey, F. E., Wedin, D. A., Zhou, X., and Zanner, C. W. (2009). Seasonal changes in depth of water uptake for encroaching trees Juniperus virginiana and Pinus ponderosa and two dominant C4 grasses in a semiarid grassland. Tree Physiol. 29, 157–169. doi: 10.1093/treephys/tpn019

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehleringer, J. R., and Dawson, T. E. (1992). Water uptake by plants: perspectives from stable isotope composition. Plant Cell Environ. 15, 1073–1082. doi: 10.1111/j.1365-3040.1992.tb01657.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Q. X., Ma, L., Green, T. R., Yu, Q., Wang, T. D., and Ahuja, L. R. (2010). Water resources and water use efficiency in the North China Plain: current status and agronomic management options. Agric. Water Manage. 97, 1102–1116. doi: 10.1016/j.agwat.2010.01.008

CrossRef Full Text | Google Scholar

Fiorini, A., Boselli, R., Amaducci, S., and Tabaglio, V. (2018). Effects of no-till on root architecture and root-soil interactions in a three-year crop rotation. Eur. J. Agron. 99, 156–166. doi: 10.1016/j.eja.2018.07.009

CrossRef Full Text | Google Scholar

Gani, A., Rahman, A., Rustam, D., and Hengsdijk, H. (2002). “Synopsis of water management experiments in Indonesia,” in Water-Wise Rice Production, eds B. A. M. Bouman, H. Hengsdijk, B. Hardy, P. S. Bindraban, T. P. Tuong, and J. K. Ladha (Philippines: IRRI), 29–37.

Google Scholar

Giri, J., Bhosale, R., Huang, G. Q., Pandey, B. K., Parker, H., Zappala, S., et al. (2018). Rice auxin influx carrier OsAUX1 facilitates root hair elongation in response to low external phosphate. Nat. Commun. 9, 7. doi: 10.1038/s41467-018-04280-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Goebel, T. S., Lascano, R. J., Paxton, P. R., and Mahan, J. R. (2015). Rainwater use by irrigated cotton measured with stable isotopes of water. Agric. Water Manage. 158, 17–25. doi: 10.1016/j.agwat.2015.04.005

CrossRef Full Text | Google Scholar

Gregory, P. J., McGowan, M., and Biscoe, P. V. (1978). Water relations of winter wheat: 2 soil water relations. J. Agric. Sci. 91, 103–116 doi: 10.1017/S0021859600056665

CrossRef Full Text | Google Scholar

Guan, D., Al-Kaisi, M. M., Zhang, Y., Duan, L., Tan, W., Zhang, M., et al. (2014). Tillage practices affect biomass and grain yield through regulating root growth, root-bleeding sap and nutrients uptake in summer maize. Field Crop. Res. 157, 89–97. doi: 10.1016/j.fcr.2013.12.015

CrossRef Full Text | Google Scholar

Guan, D., Zhang, Y., Al-Kaisi, M. M., Wang, Q., Zhang, M., and Li, Z. (2015). Tillage practices effects on root distribution and water use efficiency of winter wheat under rain-fed condition in the North China Plain. Soil Tillage Res. 146, 286–295. doi: 10.1016/j.still.2014.09.016

CrossRef Full Text | Google Scholar

Guan, J., XInxu, G., Yang, L., Keli, L., Jihua, W., and Xiaodong, G. (2007). Study on dynamic variation of root dry weight space distribution on different densities of maize. J. Maize Sci. 15, 105–108. doi: 10.3969/j.issn.1005-0906.2007.04.028

CrossRef Full Text | Google Scholar

Hatfield, J. L., Sauer, T. J., and Prueger, J. H. (2001). Managing soils to achieve greater water use efficiency. Agron. J. 93, 271–280. doi: 10.2134/agronj2001.932271x

CrossRef Full Text | Google Scholar

Hodge, A. (2004). The plastic plant: root responses to heterogeneous supplies of nutrients. New Phytol. 162, 9–24. doi: 10.1111/j.1469-8137.2004.01015.x

CrossRef Full Text | Google Scholar

Hodgkinson, L., Dodd, I. C., Binley, A., Ashton, R. W., White, R. P., Watts, C. W., et al. (2017). Root growth in field-grown winter wheat: Some effects of soil conditions, season and genotype. Eur. J. Agron. 91, 74–83. doi: 10.1016/j.eja.2017.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kou, L., McCormack, M. L., Chen, W., Guo, D., Wang, H., Gao, W., et al. (2017). Nitrogen ion form and spatio-temporal variation in root distribution mediate nitrogen effects on lifespan of ectomycorrhizal roots. Plant Soil. 411, 261–273. doi: 10.1007/s11104-016-3018-7

CrossRef Full Text | Google Scholar

Li, H., Mei, X., Nangia, V., Guo, R., Liu, Y., Hao, W., et al. (2021). Effects of different nitrogen fertilizers on the yield, water-and nitrogen-use efficiencies of drip-fertigated wheat and maize in the North China Plain. Agric. Water Manage. 243, 106474. doi: 10.1016/j.agwat.2020.106474

CrossRef Full Text | Google Scholar

Li, H., Mollier, A., Ziadi, N., Shi, Y., Parent, L.-É., and Morel, C. (2017a). The long-term effects of tillage practice and phosphorus fertilization on the distribution and morphology of corn root. Plant Soil. 412, 97–114. doi: 10.1007/s11104-016-2925-y

CrossRef Full Text | Google Scholar

Li, H. X., Mollier, A., Ziadi, N., Shi, Y. C., Parent, L. E., and Morel, C. (2017b). Soybean root traits after 24 years of different soil tillage and mineral phosphorus fertilization management. Soil Tillage Res. 165, 258–267. doi: 10.1016/j.still.2016.09.002

CrossRef Full Text | Google Scholar

Li, J., Xie, R. Z., Wang, K. R., Ming, B., Guo, Y. Q., Zhang, G. Q., et al. (2015). Variations in maize dry matter, harvest index, and grain yield with plant density. Agron. J. 107, 829–834. doi: 10.2134/agronj14.0522

CrossRef Full Text | Google Scholar

Li, Y., Zhang, P., Wu, G., Li, J., Bi, J., and Liu, J. (2011). Research on dynamic variation of dry matter accumation in maize cultivated by “double-row interlaced planting” method. S. Agr. Sci. 6, 35–38. doi: 10.14083/j.issn.1001-4942.2011.06.021

CrossRef Full Text | Google Scholar

Liu, J., Si, Z., Wu, L., Chen, J., Gao, Y., and Duan, A. (2021). Using stable isotopes to quantify root water uptake under a new planting pattern of high-low seed beds cultivation in winter wheat. Soil Tillage Res. 205, 104816. doi: 10.1016/j.still.2020.104816

CrossRef Full Text | Google Scholar

Liu, W., Liu, W., Li, P., Duan, W., and Li, H. (2010). Dry season water uptake by two dominant canopy tree species in a tropical seasonal rainforest of Xishuangbanna, SW China. Agric. Forest Meteorol. 150, 380–388. doi: 10.1016/j.agrformet.2009.12.006

CrossRef Full Text | Google Scholar

Loades, K. W., Bengough, A. G., Bransby, M. F., and Hallett, P. D. (2010). Planting density influence on fibrous root reinforcement of soils. Ecol. Eng. 36, 276–284. doi: 10.1016/j.ecoleng.2009.02.005

CrossRef Full Text | Google Scholar

Lynch, J. P. (2013). Steep, cheap and deep: an ideotype to optimize water and N acquisition by maize root systems. Ann. Bot. 112, 347–357. doi: 10.1093/aob/mcs293

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Y., and Song, X. (2016). Using stable isotopes to determine seasonal variations in water uptake of summer maize under different fertilization treatments. Sci. Total Environ. 550, 471–483. doi: 10.1016/j.scitotenv.2016.01.148

PubMed Abstract | CrossRef Full Text | Google Scholar

Majdi, H., and Andersson, P. (2005). Fine root production and turnover in a norway spruce stand in northern sweden: effects of nitrogen and water manipulation. Ecosystems. 8, 191–199. doi: 10.1007/s10021-004-0246-0

CrossRef Full Text | Google Scholar

McCole, A. A., and Stern, L. A. (2007). Seasonal water use patterns of Juniperus ashei on the Edwards Plateau, Texas, based on stable isotopes in water. J. Hydrol. 342, 238–248. doi: 10.1016/j.jhydrol.2007.05.024

CrossRef Full Text | Google Scholar

Meinzer, F. C., Clearwater, M. J., and Goldstein, G. (2001). Water transport in trees: current perspectives, new insights and some controversies. Environ. Exp. Bot. 45, 239–262. doi: 10.1016/S0098-8472(01)00074-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Meißner, M., Köhler, M., Schwendenmann, L., Hölscher, D., and Dyckmans, J. (2014). Soil water uptake by trees using water stable isotopes (δ2H and δ18O)–a method test regarding soil moisture, texture and carbonate. Plant Soil. 376, 327–335. doi: 10.1007/s11104-013-1970-z

CrossRef Full Text | Google Scholar

Parnell, A. C., Phillips, D. L., Bearhop, S., Semmens, B. X., Ward, E. J., Moore, J. W., et al. (2013). Bayesian stable isotope mixing models. Environmetrics 24, 387–399. doi: 10.1002/env.2221

CrossRef Full Text | Google Scholar

Peake, A. S., Huth, N. I., Kelly, A. M., and Bell, K. L. (2013). Variation in water extraction with maize plant density and its impact on model application. Field Crop. Res. 146, 31–37. doi: 10.1016/j.fcr.2013.02.012

CrossRef Full Text | Google Scholar

Penna, D., Geris, J., Hopp, L., and Scandellar, F. (2020). Water sources for root water uptake: Using stable isotopes of hydrogen and oxygen as a research tool in agricultural and agroforestry systems. Agric. Ecosys. Environ. 291, 106790. doi: 10.1016/j.agee.2019.106790

CrossRef Full Text | Google Scholar

Piao, L., Qi, H., Li, C., and Zhao, M. (2016). Optimized tillage practices and row spacing to improve grain yield and matter transport efficiency in intensive spring maize. Field Crop. Res. 198, 258–268. doi: 10.1016/j.fcr.2016.08.012

CrossRef Full Text | Google Scholar

Prechsl, U. E., Burri, S., Gilgen, A. K., Kahmen, A., and Buchmann, N. (2015). No shift to a deeper water uptake depth in response to summer drought of two lowland and sub-alpine C3-grasslands in Switzerland. Oecologia. 177, 97–111. doi: 10.1007/s00442-014-3092-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Querejeta, J. I., Estrada-Medina, H., Allen, M. F., and Jiménez-Osornio, J. J. (2007). Water source partitioning among trees growing on shallow karst soils in a seasonally dry tropical climate. Oecologia. 152, 26–36. doi: 10.1007/s00442-006-0629-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramezani, M., Ab, R. R. S., Mobasser, H., and Amiri, E. (2011). Effects of row spacing and plant density on silage yield of corn (Zea mays L. cv. sc704) in two plant pattern in North of Iran. Afr. J. Agric. Res. 6, 1128–1133. doi: 10.5897/AJAR10.898

CrossRef Full Text | Google Scholar

Rogers, E. D., and Benfey, P. N. (2015). Regulation of plant root system architecture: implications for crop advancement. Curr. Opin. Biotech. 32, 93–98. doi: 10.1016/j.copbio.2014.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossatto, D. R., da Silveira Lobo Sternberg, L., and Franco, A. C. (2013). The partitioning of water uptake between growth forms in a Neotropical savanna: do herbs exploit a third water source niche? Plant Biol. 15, 84–92. doi: 10.1111/j.1438-8677.2012.00618.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Schenk, H. J., and Jackson, R. B. (2005). Mapping the global distribution of deep roots in relation to climate and soil characteristics. Geoderma 126, 129–140. doi: 10.1016/j.geoderma.2004.11.018

CrossRef Full Text | Google Scholar

Schwendenmann, L., Pendall, E., Sanchez-Bragado, R., Kunert, N., and Hölscher, D. (2015). Tree water uptake in a tropical plantation varying in tree diversity: interspecific differences, seasonal shifts and complementarity. Ecohydrology. 8, 1–12. doi: 10.1002/eco.1479

CrossRef Full Text | Google Scholar

Shao, H., Xia, T. T., Wu, D. L., Chen, F. J., and Mi, G. H. (2018). Root growth and root system architecture of field-grown maize in response to high planting density. Plant Soil. 430, 395–411. doi: 10.1007/s11104-018-3720-8

CrossRef Full Text | Google Scholar

Sprenger, M., Leistert, H., Gimbel, K., and Weiler, M. (2016). Illuminating hydrological processes at the soil-vegetation-atmosphere interface with water stable isotopes. Rev. Geophys. 54, 674–704. doi: 10.1002/2015RG000515

CrossRef Full Text | Google Scholar

Tan, X., Bi, J., Liu, J., and Ye, B. (2010). Development trends of maize planting patterns. S. Agr. Sci. 5, 57–59. doi: 10.14083/j.issn.1001-4942.2010.05.014

CrossRef Full Text | Google Scholar

Tang, K., and Feng, X. (2001). The effect of soil hydrology on the oxygen and hydrogen isotopic compositions of plants' source water. Earth Planet. Sci. Lett. 185, 355–367. doi: 10.1016/S0012-821X(00)00385-X

CrossRef Full Text | Google Scholar

Testa, G., Reyneri, A., and Blandino, M. (2016). Maize grain yield enhancement through high plant density cultivation with different inter-row and intra-row spacings. Eur. J. Agron. 72, 28–37. doi: 10.1016/j.eja.2015.09.006

CrossRef Full Text | Google Scholar

Wang, J., Fu, B., Lu, N., and Zhang, L. (2017a). Seasonal variation in water uptake patterns of three plant species based on stable isotopes in the semi-arid Loess Plateau. Sci. Total Environ. 609, 27–37. doi: 10.1016/j.scitotenv.2017.07.133

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., He, H., Zhang, Z., Dai, Z., Tan, Y., and Chang, P. (2017b). Effects of wide-narrow row interlaced thick-planting pattern on photosynthetic characterisric and yield in summer maize. J. Maize Sci. 25, 65–72. doi: 10.13597/j.cnki.maize.science.20170312

CrossRef Full Text | Google Scholar

Wang, Y., Zhang, X. Y., Chen, J., Chen, A. J., Wang, L. Y., Guo, X. Y., et al. (2019). Reducing basal nitrogen rate to improve maize seedling growth, water and nitrogen use efficiencies under drought stress by optimizing root morphology and distribution. Agric. Water Manage. 212, 328–337. doi: 10.1016/j.agwat.2018.09.010

CrossRef Full Text | Google Scholar

West, A. G., Dawson, T. E., February, E. C., Midgley, G. F., Bond, W. J., and Aston, T. L. (2012). Diverse functional responses to drought in a Mediterranean-type shrubland in South Africa. New Phytol. 195, 396–407. doi: 10.1111/j.1469-8137.2012.04170.x

PubMed Abstract | CrossRef Full Text | Google Scholar

West, A. G., Patrickson, S. J., and Ehleringer, J. R. (2006). Water extraction times for plant and soil materials used in stable isotope analysis. Rapid Commun. Mass Sp. 20, 1317–1321. doi: 10.1002/rcm.2456

PubMed Abstract | CrossRef Full Text | Google Scholar

Whalley, W. R., Binley, A., Watts, C. W., Shanahan, P., Dodd, I. C., Ober, E. S., et al. (2017). Methods to estimate changes in soil water for phenotyping root activity in the field. Plant Soil. 415, 407–422. doi: 10.1007/s11104-016-3161-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Widdicombe, W. D., and Thelen, K. D. (2002). Row width and plant density effects on corn grain production in the northern corn belt. Agron. J. 94, 1020–1023. doi: 10.2134/agronj2002.1020

CrossRef Full Text | Google Scholar

Wu, H., Li, J., Zhang, C., He, B., Zhang, H., Wu, X., et al. (2018a). Determining root water uptake of two alpine crops in a rainfed cropland in the Qinghai Lake watershed: First assessment using stable isotopes analysis. Field Crop. Res. 215, 113–121. doi: 10.1016/j.fcr.2017.10.011

CrossRef Full Text | Google Scholar

Wu, X., Chen, Y., Sui, P., Gao, W., Yan, P., and Tao, Z. (2015). Effect of planting geometries on canopy structure of spring maize under high-density condition in North China Plain. Chin. J. Ecol. 34, 1–8. doi: 10.13292/j.1000-4890.2015.0004

CrossRef Full Text | Google Scholar

Wu, Y., Du, T., Li, F., Li, S., Ding, R., and Tong, L. (2016). Quantification of maize water uptake from different layers and root zones under alternate furrow irrigation using stable oxygen isotope. Agric. Water Manage. 168, 35–44. doi: 10.1016/j.agwat.2016.01.013

CrossRef Full Text | Google Scholar

Wu, Y., Du, T., Yuan, Y., and Shukla, M. K. (2018b). Stable isotope measurements show increases in corn water use efficiency under deficit irrigation. Sci. Rep. 8, 14113. doi: 10.1038/s41598-018-32368-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Zhou, H., Zheng, X.-J., Li, Y., and Tang, L.-S. (2014). Seasonal changes in the water use strategies of three co-occurring desert shrubs. Hydrol. Process. 28, 6265–6275. doi: 10.1002/hyp.10114

CrossRef Full Text | Google Scholar

Yang, B., Wen, X., and Sun, X. (2015a). Irrigation depth far exceeds water uptake depth in an oasis cropland in the middle reaches of Heihe River Basin. Sci. Rep. 5, 15206. doi: 10.1038/srep15206

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, B., Wen, X., and Sun, X. (2015b). Seasonal variations in depth of water uptake for a subtropical coniferous plantation subjected to drought in an East Asian monsoon region. Agric. Forest Meteorol. 201, 218–228. doi: 10.1016/j.agrformet.2014.11.020

CrossRef Full Text | Google Scholar

Yi, Z., Wang, P., and Tu, N. (2009). Responses of the distribution and nitrogen content in roots of summer maize to nitrogen fertilization types and amounts. Plant Nutr. Fert. Sci. 15, 91–98. doi: 10.3321/j.issn:1008-505X.2009.01.013

CrossRef Full Text | Google Scholar

Zhang, C., Zhang, J., Zhao, B., Zhu, A., Zhang, H., Huang, P., et al. (2011). Coupling a two-tip linear mixing model with a δD–δ18O plot to determine water sources consumed by maize during different growth stages. Field Crop. Res. 123, 196–205. doi: 10.1016/j.fcr.2011.04.018

CrossRef Full Text | Google Scholar

Zhang, X., Chen, S., Liu, M., Pei, D., and Sun, H. (2005). Improved water use efficiency associated with cultivars and agronomic management in the North China Plain. Agron. J. 97, 783–790. doi: 10.2134/agronj2004.0194

CrossRef Full Text | Google Scholar

Zhang, X. F., Zhu, A. N., Xin, X. L., Yang, W. L., Zhang, J. B., and Ding, S. J. (2018). Tillage and residue management for long-term wheat-maize cropping in the North China plain: I. crop yield and integrated soil fertility index. Field Crop. Res. 221, 157–165. doi: 10.1016/j.fcr.2018.02.025

CrossRef Full Text | Google Scholar

Zhang, X. X., Whalley, P. A., Ashton, R. W., Evans, J., Hawkesford, M. J., Griffiths, S., et al. (2020). A comparison between water uptake and root length density in winter wheat: effects of root density and rhizosphere properties. Plant Soil. 451, 345–356. doi: 10.1007/s11104-020-04530-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yu, Z., Zhang, Y., and Shi, Y. (2021). Split nitrogen fertilizer application improved grain yield in winter wheat (Triticum aestivum L.) via modulating antioxidant capacity and ¹³C photosynthate mobilization under water-saving irrigation conditions. Ecol. Process. 10, 1–13. doi: 10.1186/s13717-021-00290-9

CrossRef Full Text | Google Scholar

Zhao, X., Li, F., Ai, Z., Li, J., and Gu, C. (2018). Stable isotope evidences for identifying crop water uptake in a typical winter wheat–summer maize rotation field in the North China Plain. Sci. Total Environ. 618, 121–131. doi: 10.1016/j.scitotenv.2017.10.315

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, L., Ma, J., Sun, X., Guo, X., Li, Y., Rong, R., et al. (2019). Effective root growth zone of apple tree under water storage pit irrigaiton using stable isotope methodology. Arch. Agron. Soil Sci. 65, 1521–1535. doi: 10.1080/03650340.2019.1566718

CrossRef Full Text | Google Scholar