The automated irrigation module consisted of a drainage tank, a nutrient solution tank, and a standard nutrient solution tank (Figure 1). In the drainage tank, the weight and EC of the collected drainage were measured by a load cell (JSB-20, CAS, South Korea) and an EC sensor (SCF-01A, DIK, South Korea). In the mixing tank, the nutrient solution’s mixing process was carried out according to the degree of nutrient recycling. The drainage tank’s storage capacity was 11.7 L, and the mixing tank capacity was 19.4 L. Ultrasonic-level sensors (UHA-300, Unics, South Korea) were installed at the top of the mixing tank to monitor the feed amount of drainage, raw water, and standard nutrient solution. Peristaltic pumps (M500, Verderflex, United Kingdom) were used to transfer standard nutrient solutions, drainage, and water from the pre-adjusted standard nutrient solution, drainage, and water tank to the mixing tank. For the irrigation of the mixed nutrient solution, a centrifugal pump (PUN-350M, Wilo Pump, South Korea) was used. A pyranometer (SP-110, Apogee, United States) was used for the integrated solar irradiance for automated irrigation control. A data logger (CR1000, Campbell Scientific, United States) and controller (SDM-CD16AC, Campbell Scientific, United States) were used for measurement, nutrient solution preparation, and irrigation control.

The cultivation experiment using the four automated irrigation modules was carried out in a Venlo-type greenhouse located at the experimental farm of Seoul National University (Suwon, South Korea, lat. 37.3°N, long. 127.0°E). The four automated irrigation modules consist of one OL, two SCL, and one CL soilless culture systems. Accordingly, the degree of nutrient recycling was manipulated by applying the OL, SCL, and CL soilless cultures, respectively. The OL soilless culture system uses no drainage in the irrigated solution. The SCL soilless culture system partially uses drainage according to the drainage dilution level. In the CL soilless culture system, most of the collected drainage is reused for the nutrient solution. Basically, we referred to the general process of automated nutrient solution preparation (Savvas and Manos, 1999; Savvas, 2002). The OL irrigation controls the mixing tank EC to the target value. Conventionally, the SCL irrigation reuses the drainage by diluting to a lower EC than the supply EC with water, and stock solutions are applied to the supply EC of the irrigation. However, during the experimental period, the target EC for the irrigated nutrient solution was maintained at 2.6 dS m^–1. Thus, in the greenhouse experiment, a pre-adjusted standard nutrient solution for the supply EC (2.6 dS m^–1) was stored in the standard nutrient solution tank, and this was used to replenish the standard nutrient solution. In the greenhouse experiment, the determination of the mixing ratio of drainage, water, and the standard nutrient solution was based on the calculation procedure described in section “Determination of the Mixing Ratio in the Soilless Culture System,” considering the amount of water already added in the pre-adjusted standard nutrient solution.

The four modules of the nutrient recycling degree treatment were composed of OL, CL, and SCL modules with two EC levels for a drainage dilution of 0.65 and 1.3 dS m^–1 (SCL 0.65 and SCL 1.3). An integrated solar irradiance-based irrigation control was applied. The modules automatically irrigated 150 ml of the nutrient solution prepared in the mixing tank per plant whenever the integrated solar irradiance reaches 100 J cm^–2. However, the transpiration capacity increases due to plant growth during the cultivation period. Thus, the irrigation frequency or irrigation amount must be adjusted to compensate for the transpiration change. Thus, the supply volume per irrigation event was adjusted by the meteorological condition to compensate for the transpiration capacity change and maintain a drainage ratio of 30%. Each module prepared the nutrient solutions according to the recycling degree before the irrigation event.

Compositions of the standard nutrient solution used in this experiment were 14.56 NO₃^–, 1.18 H₂PO₄^–, 3.20 SO₄^2–,5.96 K⁺, 9.56 Ca²⁺, 3.38 Mg²⁺, 0.21Na⁺, and 0.30Cl^– (in meq L^–1) as macro-elements and 18 Fe, 10 Zn, 2 Cu, 10 Mn, and 0.5 Mo (in μM) as micro-elements. The target EC of the nutrient solution mixing module was set to 2.6 dS m^–1, and the EC of the raw water was set to 0.15 dS m^–1. Daytime temperature, night time temperature, and relative humidity of the greenhouse remained between 25 and 30°C, 15 and 22°C, and 50 and 80%, respectively, by an environmental control system of the greenhouse. Sweet pepper (Capsicum annuum L. “Fiesta”) plants were used in the experiment and seeded on July 15, 2011. The sweet pepper plants were transplanted into 90 cm (L) × 15 cm (W) × 7 cm (H) rockwool slabs (Cultilene, Netherlands) on September 29, 2011. Three sweet pepper plants were planted per slab, and three slabs were used per experimental treatment. The planting density was 2.8 plants per square meter.

The sweet pepper plants were cultivated by OL nutrient supply until December 15, 2011, and after that, each treatment was applied, and the experiment was finished on March 9, 2012. On the day before the initiation of the treatments, a large amount of a standard nutrient solution was irrigated in an OL nutrient supply in order to make the initial condition of the nutrients balanced in the rockwool slab of all treatments close to the standard composition.

To observe the changes in the ratio between nutrients in the root zone of the soilless culture system, samples of nutrient solution in the rockwool slabs were extracted using a syringe. The collection points of the nutrient solution in the rockwool slab were randomly selected for the collection of representative samples of the overall concentration in the root zone. Then, 10 ml of a root zone nutrient solution was collected for each extraction, and this was performed five times to make a 50 ml sample. Six samples per treatment were collected each time every 2 weeks. K⁺, Ca²⁺, and Mg²⁺ were analyzed using an inductively coupled plasma optical emission spectrometer (ICP-730ES, Varian, Australia). The SAS system (version 9.2, SAS Institute, United States) was used for statistical analysis. In the experimental condition, three rockwool slabs were connected to the single automated module system. Thus, the replicate of the sample could have limited data variability. However, the soilless culture using slab-form substrates has an isolated system boundary of the root zone. Thus, even in the single system (i.e., a single greenhouse with a single irrigation system), each slab condition often shows significant variations in root zone nutrient and water according to the microclimate fluctuations of each slab’s location. Therefore, at least, it can be expected that the rockwool slab replicates could statistically rule out these effects in this experimental condition.

The models of the soilless culture system used in the simulation consist of water and nutrient transport, nutrient uptake, and nutrient solution mixing (Figure 2). This model’s basic structures were constructed based on previous models of the soilless culture system (Silberbush and Ben-Asher, 2001; Silberbush et al., 2005). The simulation model reflected the automated nutrient solution mixing process of the EC-based soilless culture system. For nutrients and water transport in the substrate, mass flow and diffusion models in porous media were applied (Shackelford and Daniel, 1991; Corwin et al., 1993; Snape et al., 1995). The nutrient uptake by plants follows the Michaelis–Menten equation. The rate of nutrient uptake is mainly driven by concentration. In this study, a modified Michaelis–Menten equation having an additional variable for the transpiration rate was used to reflect more stochastic changes in absorption parameter variations such as plant growth. This model considers reducing the nutrient depletion zone around the roots due to the accelerated mass flow generated by an increase in transpiration rate with plant growth (Sago et al., 2011; Nomiyama et al., 2013). Since transpiration is a variable that includes parameters for plant growth changes (Baille et al., 1994; Ta et al., 2012), the model was used to reflect the change in the nutrient absorption rate with plants’ growth. The definitions, values, units, and sources of the models’ parameters are summarized in Table 1.

Several assumptions were made to reduce the complexity of the simulation:

Water transport in a CL soilless culture system occurs between the substrate, the drainage tank, and the mixing tank, and these correspond to Eqs. 1–3, respectively.

where V(cm³) is the volume of water, subscript sub, n is the water content of each layer when the substrate is divided into n layers, drg is the drainage tank, and mix is the mixing tank. Q(cm³ min^–1) is the flow rate of the water; subscript ir is the irrigation flow rate to the top of the substrate (n = 1);subscripts dr,n-1 and dr,n are the flow rates of the drain from the n - 1 layer and the n + 1 layer, respectively; and tr,n represents the transpiration rate of the substrate layer n. Here, Q_dr,_n is the difference between the irrigation rate to the substrate and the transpiration rate (Q_ir – Q_tr,1or Q_dr,_n–1 – Q_tr,_n) (Shin et al., 2014). However, Q_dr,_n and Q_tr,_n are restricted by the field capacity (F, dimensionless) and difficult available water (W_DAW, dimensionless), respectively, for the volume of the layer of the substrate (S_sub, cm³). And Q_dr,_n flows only when V_sub,_n > S_subF, and Q_tr,_n flows only when V_sub,_n > S_subW_DAW. Q_drg, Q_stk, and Q_wtr correspond to the drain, stock solution, and raw water flow rates, respectively, introduced into the mixing tank when irrigation to the root zone is required. Q_disch is a variable that is applied to the SCL system, meaning the flow rate of the nutrient solution discharged to the outside when the volume of stored drainage (V_drg) exceeds the capacity of the drainage tank (S_drg, cm³).

In a soilless culture system, the movement of nutrients follows the same path as water, and the movement of nutrients between the substrate, drainage tank, and mixing tank can be expressed as Eqs. 4–6, respectively.

where C(meq cm^–3) represents the equivalent concentration of nutrients and superscript I is the type of ions (K⁺, Ca²⁺, and Mg²⁺). Subscripts sub, drg, mix, stk, and wtr represent the source of the nutrients and indicate substrate, drainage tank, mixing tank, stock solution, and raw water, respectively. J^I(meq plant^–1 min^–1) means the uptake rate of nutrients. The final term of Eq. 4 refers to the diffusion of solutes in a substrate, where D(10^–5 cm² min^–1) is the diffusion coefficient of an ion I, A(cm²) is the cross-sectional area of the substrate; θ_n is the water content of the n layer (V_sub,_n/S_sub, dimensionless); Z(cm) is the height of the n layer. In the last layer, the diffusion term from n + 1 is excluded because the diffusion occurs from the previous layer only.

The Michaelis–Menten equation was used to model the nutrient uptake rate of the plant (Eq. 7).

where JmaxI means the maximum absorption rate of the ion and KmI is the Michaelis–Menten constant. JnI is the absorption rate of individual nutrient I in layer n, and the equation applying the change in the rate of absorption of nutrients to the increase of transpiration rate (Q_tr,_n) is shown in Eq. 8.

In the simulation, the degree of nutrient recycling was manipulated by applying nutrient solution mixing processes of the OL, SCL, and CL soilless cultures. The nutrient solution preparation process for the OL, SCL, and CL soilless cultures of the above section of the greenhouse experiment of the OL, SCL, and CL soilless cultures was converted to this study’s soilless culture system model.

Q_stk, Q_drg, and Q_wtr of the SCL soilless culture system are determined by Eqs. 9–11, respectively.

where C_dil(meq cm^–3) is the equivalent target concentration for the dilution of the total nutrients in the drainage and C_tar is the target equivalent concentration for the total nutrients in the nutrient solution to be supplied to the plant. The OL system corresponds to the case when V_drg ≤ 0 in Eqs. 9–11.

Q_stk, Q_wtr, and Q_drg for the CL system are determined by Eqs. 12–14, respectively.

The CL soilless culture system discharges no drainage. Water replenishment into the system follows the transpiration amount in the system. Nutrients removed from the system correspond to the absorbed nutrients by plants. Thus, the Michaelis–Menten parameters of the nutrient uptake model were estimated by performing a numerical integration-based progress curve analysis using the water and nutrient inputs, transpiration, and nutrient concentration data measured in the CL soilless culture experiment. Simulation and progress curve analysis of these models were performed using Berkeley Madonna 8.3.23 (Berkeley Madonna, Inc., University of California, Berkeley).

The estimated nutrient uptake model parameters were applied to the OL, SCL, and CL soilless culture system models. We performed two simulation analyses as follows:

Under the simulation scenario, the transpiration rate, one of the main disturbance factors for nutrient uptake concentration, was manipulated by a random walk. The long-term feedback nutrient control was performed by changing the nutrient ratios in the stock solution for nutrient replenishment. The feedback period was set to 12 weeks, which is the total period of the simulation. After each simulation, the difference between the target nutrient ratio and average nutrient ratio (target average) was fed back to adjust the stock solution’s nutrient ratio (Figure 2). Thus, as an example, if the average nutrient ratio is higher than the target ratio, negative feedback is returned to the stock solution’s nutrient ratio.

Measured and simulated nutrient concentrations and ion ratios in the CL substrate showed close agreements with R² values of 0.81–0.88 and a root mean square error (RMSE) of 2.96–6.36 meq L^–1 (Figure 3A) and R² values of 0.72–0.94 and an RMSE of 0.0161–0.0188 (dimensionless) (Figure 3B), respectively. Measured and simulated values in all nutrient concentrations fluctuated with similar patterns. However, the nutrient ratio of the CL showed distinguishable directional trends. The cumulative nutrient supply to the simulated and experimental CL systems showed good agreement with the measured data with an R²-value of 0.95 (Figure 4).

Overall, the soilless culture system model closely simulated the nutrient concentration variations in the CL soilless culture. Furthermore, the nutrient supply into the CL soilless culture system’s boundary was simulated with a determination coefficient of 0.95. The replenished nutrients in a CL soilless culture system have a relation with the absorbed nutrient and water in the system (Savvas, 2002). Nutrient concentration variations observed in the CL soilless culture system represent the combined effect of nutrient and water absorption, nutrient replenishment, irrigation, and drainage recycling. Thus, the nutrient preparation process and nutrient variation aspects in the soilless culture system were verified here. However, it needs to be noted that the simulation model estimated a relatively high KmK. This might be derived from the current model’s structural limitations. The modified nutrient uptake model applies transpiration variations to the nutrient uptake parameter variations. Thus, the parameter estimation could be valid in the simulation condition.

The nutrient concentration in the substrate indicated similar trends throughout the experiment period. The concentration variations were not distinct from the treatment and showed little consistency in the treatments’ significant differences (Figure 5A). On the other hand, when the nutrients were converted to a ratio between the nutrients, the nutrient variation trends indicated more distinguishable trends in individual nutrients and each treatment (Figure 5B). The ratio of K⁺ in CL, which has the least nutrient supply, was significantly lower than that of the other treatments and deviated most from the initial value. OL, which has the highest nutrient input, was the closest to the initial ratios. In the CL treatment, Mg²⁺ and Ca²⁺ remained higher than those in other treatments after observing significant differences 4 weeks after treatment. Mg²⁺ and Ca²⁺were significantly lower in OL and tended to be the closest to the initial ratios (Figure 5B). Thus, unlike nutrient concentrations, the nutrient ratio conversion revealed increasing or decreasing individual nutrients under the nutrient concentration’s overall fluctuation.

Furthermore, those increasing or decreasing trends in nutrient ratios and K⁺ leveled off 6 weeks after treatment. These results suggest that the nutrient ratios have more deterministic aspects than the nutrient concentration. On the other hand, the nutrient concentrations were relatively difficult to distinguish between the nutrients. The concentration trends between the treatments also showed complicated fluctuations. Changes in nutrient concentrations in the root zone in soilless culture systems vary depending on the difference between irrigated nutrient concentration and nutrient uptake concentration (van Noordwijk, 1990). The nutrient uptake concentration represents the ratio of nutrients to water taken up over time (Le Bot et al., 1998). Thus, the change in the uptake concentration results from variations in the transpiration and nutrient uptake rate. The transpiration could affect mass flow near the root system (Marschner and Marschner, 2012). The transpiration could also partially drive the uptake rate of ions like Ca²⁺ (Adams and Ho, 1993; Mengel, 2009). Thus, an ion-specific interaction between transpiration and nutrient uptake could affect the substrate’s nutrient concentration variations.

However, in this study, the concentration of the substrate nutrient showed the synchronized fluctuations. As explained above, the nutrient concentration indicated similar trends and showed an overall gradual increase from the 8 weeks onwards. The drainage EC of CL, SCL, and OL systems followed similar tendencies (Figure 6). However, the CL system’s drainage EC showed more sensitive variations to the nutrient concentration change. In contrast, in SCL and OL, which have relatively low drainage reuse proportions, less sensitive EC changes were observed. The low drainage reuse leads to the accumulation of residual drainage in the tank, and subsequently, the EC in the drainage tank responds less sensitively to the substrate nutrient concentration changes. As a result, depending on each system’s storage conditions, the sensitivity of the EC fluctuation in the drainage tank was different, but the overall trends indicated similarity with variations in the substrate nutrient concentration. The overall increase in the nutrient concentration 8 weeks after treatment could mainly result from the increased transpiration rate as plant growth progresses. In general, a plant’s transpiration rate is accelerated as the plant expands its leaf area during growth (Bailey et al., 1993). In this study, the transpiration gradually increased to the end of the greenhouse experiment (Figure 8A). This may suggest that the change in transpiration has a more dominant global effect on the nutrient concentration changes in nutrient concentrations.

Changes in transpiration and consequent variations in nutrient absorption concentration mainly drove the macroscopic fluctuations in the overall nutrients and EC. On the other hand, nutrient absorption is mainly dominated by transporters’ arrangement located on the root surface (Arsova et al., 2020). Thus, converting the nutrient concentrations to nutrient ratios might have reduced the global effect of transpiration on the nutrient concentrations and revealed the distinguishable directional trends. Subsequently, it can be conjectured that the nutrient ratio changes in a soilless culture system mainly reflect the nutrient uptake phenomenon driven by the nutrient transporter. That is, the nutrient ratio changes could implicate the Michalis–Menten-driven nutrient uptake. Consequently, the nutrient ratio filtered the variations in nutrient uptake concentrations variations and showed the combined effect of nutrient recycling degree and nutrient uptake parameters.

Simulation (1), described in section “Parameter Estimation and Theoretical Analysis of the Nutrient Variations,” showed the nutrient ratio changes according to the degree of nutrient recycling, namely, OL, SCL, and CL irrigation models. The degree of nutrient recycling is also strongly related to the degree of nutrient supply into the system boundaries of OL, SCL, and CL systems. Thus, the nutrient ratio changes of the soilless culture systems were summarized on the basis of a nutrient supply (Figure 7). As a representative, when the K⁺ ratios in the soilless culture systems were summarized for nutrient inputs, specific trends were observed. As nutrient supply increases, i.e., as the degree of nutrient recycling decreases, the nutrient ratio’s deviation from the standard ratio approached zero. Furthermore, the K⁺ ratios in the experiment showed similar trends as the simulation. The deviation between the nutrient ratio measured in the experiment and the standard ratio decreased in the order of CL, SCL 1.3, SCL 0.65, and OL.

The nutrient variations in the CL soilless culture system are often considered a black box in some way for its dynamic behaviors in uptake concentration (van Noordwijk, 1990). Thus, in CL soilless culture, managing the nutrient uptake concentration effect is often approached as empirical methods. The literature on the EC-based CL soilless culture system investigated the appropriate interval or level of nutrient solution renewal (Ehret et al., 2005; Ko et al., 2014), drainage discharge (Massa et al., 2011), and EC control (Rouphael et al., 2016; Signore et al., 2016). The CL soilless culture system was featured due to uncertainties in nutrient management. However, Simulation (1) shows that although the OL, SCL, and CL soilless culture systems have been separated from each other, the systems’ nutrient ratios can be summarized continuously based on the systems’ nutrient inputs.

In Simulation (1) of section “Parameter Estimation and Theoretical Analysis of the Nutrient Variations,” a modified concept on the nutrient variations of OL, SCL, and CL soilless culture systems was presented. The nutrient ratio changes of each system were macroscopically integrated into nutrient inputs to the systems. In this respect, the soilless culture models could be simplified by the inflow and the outflow into the system boundary and the system’s uptake rate (Eq. 15). Moreover, with these critical parameters, the theoretical framework for the nutrient variations of OL, SCL, and CL soilless culture systems could be reconstructed.

where Q_stk and Q_wtr are the flow rates of the stock nutrient solution and raw water flowing into the boundary of the soilless culture system, respectively, and Q_out is the flow rate of the nutrient solution discharged outside the boundary. CsysI corresponds to the concentration of a nutrient in the system. CstkI and CwtrI are concentrations of a nutrient in the stock solution and raw water, respectively. V_sys is the volume of nutrient solution in the system. V_sys in Eq. 15 cannot exceed the system’s capacity, and the consumed water is replenished repeatedly and thus is assumed to be constant.

In the SCL system, the boundary of the system includes the mixing tank and the drainage tank because the drainage is partially circulated or discharged outside the system. Therefore, Q_stk and Q_wtr correspond to the supply rate of nutrients and water to the soilless culture system, respectively; and Q_out in Eq. 15 corresponds to the flow rate of the nutrient solution discharged from the drainage tank (Q_disch in Eq. 5). In the case of the CL system, the boundary is the same as that of the SCL system, but Q_out does not occur, and the steady-state solution, in this case, can be summarized as in Eq. 16.

This indicates that the point of convergence for CsysI is determined according to the concentrations of the supplied stock solution (CstkI) and raw water (CwtrI) and their inflow rates (Q_stk and Q_wtr, respectively) and means that when the inflow rates exceed the maximum nutrient uptake rate of the plant (Jm⁢a⁢xI), the system becomes unstable. However, in the EC-based CL soilless culture system, Q_stk is determined based on the total equivalent concentration of nutrients in the system; as a result, the feeding rate of the total nutrients follows the uptake rate of total nutrients at a certain level. Thus, too excessive increases or decreases can be controlled by general EC control practice. Under this condition, the nutrient ratio conversion could reveal the parameter-based nutrient variations. Unlike dynamic variations of nutrient concentration, the parameter-based nutrient variations could indicate deterministic changes. Thus, these deterministic changes could be controllable by long-term feedback.

Assuming that the CsysI is the total equivalent concentration of all nutrients in the system measured by EC and CsysI is controlled in steady-state through EC feedback, Eq. 16 theoretically predicts that CsysI converges to a constant value. This means that the individual nutrient variations could also level off to steady-state. Additionally, the terms including Q_wtr in Eq. 16 may contribute to the system fluctuation by replenishing the water consumed through transpiration. Thus, its influence is determined by the mineral concentration in the raw water. However, in most CL or SCL soilless culture studies, including the cultivation experiment of this study, and OL soilless culture, irregular fluctuations in the total nutrient concentrations relative to its initial values have been reported often (Hao and Papadopoulos, 2002; Massa et al., 2011; Shin and Son, 2016; Signore et al., 2016; Lee et al., 2017; Moon et al., 2018, 2019).

However, the recent literature theoretically deduced a modified nutrient replenishment method for steady-state management of EC and experimentally demonstrated its effect (Ahn and Son, 2019). Here, this study’s theoretical prediction on the effect of steady-state EC control was observed in the greenhouse experiment. When the total ion concentration was controlled to the initial target value, the individual nutrients of the EC-based CL soilless culture system rapidly converged into the average steady state.

Overall, Simulation (1) defines that the SCL and OL soilless culture systems are mainly driven by the nutrient input and discharge. On the other hand, the nutrient variation in a CL soilless culture system is mainly determined by the relationship between nutrient inputs and nutrient uptake parameters.

Simulation (2) shows the average changes of the nutrient ratios and their standard deviations over the simulation period when performing the nutrient balance control described in section “Parameter Estimation and Theoretical Analysis of the Nutrient Variations.” For each simulation, randomized transpiration rates were generated by the random walk method (Figure 8A). The mean values of K⁺, Ca²⁺, and Mg²⁺ ratios gradually converged to the control target values as the times of stock solution adjustment increased (Figure 8B). Also, the standard deviations gradually decreased. Even though the random walk generated arbitrary uptake concentration in every simulation, the nutrient ratios in the CL soilless culture system gradually approached the target ratios. Thus, every soilless culture system simulated in this scenario has experienced different nutrient uptake concentrations.

In this simulation, a different path of the transpiration by random walk was applied every time, within the range of maximum fluctuation during the cultivation period. Furthermore, the theoretical and experimental analyses were conducted under the substrate condition, which has a non-homogeneous root-zone condition. These can be regarded as extreme conditions. However, the average nutrient ratio approached the target value even in a long-term feedback period of 12 weeks under various nutrient uptake concentrations. This means that the discussion about Eq. 16 is also theoretically valid. Moreover, other soilless culture systems using no substrate such as nutrient film and deep-flow technique could have less spatial variations in nutrient and water distribution than a substrate-based system. Thus, this theoretical framework can be expected to be compatible with other soilless culture systems. Consequently, the results of the theoretical and experimental analyses suggest that the nutrient variations in EC-based soilless cultures could be managed even under the weekly based long-term nutrient analysis conditions.

In a previous study, nutrient control with a 2 weeks feedback period was conducted in an EC-based CL soilless culture system (Savvas, 2002). Here, the difference between the target concentration and the current concentration at every nutrient analysis event was applied as the feedback value. This corresponds to the proportional control from the viewpoint of control engineering. The proportional control can generate a steady-state error, and the literature reported similar trends in the nutrient concentrations. However, no discussion in this respect was conducted so far. Although the nutrient uptake concentration fluctuated by transpiration could be a short-term phenomenon, the cumulative difference between the supplied nutrient concentration and nutrient uptake concentration could result in a particular trend of nutrient concentration or nutrient ratio changes on a time series. Neutralizing the transpiration effects from the nutrient management in the recirculating system could reveal the kinetic domain of nutrient uptake. Thus, in this domain, nutrient variation complexity could be narrowed down to nutrient uptake parameter aspects.