Evaluating Status Change of Soil Potassium from Path Model

The purpose of this study is to determine critical environmental parameters of soil K availability and to quantify those contributors by using a proposed path model. In this study, plot experiments were designed into different treatments, and soil samples were collected and further analyzed in laboratory to investigate soil properties influence on soil potassium forms (water soluble K, exchangeable K, non-exchangeable K). Furthermore, path analysis based on proposed path model was carried out to evaluate the relationship between potassium forms and soil properties. Research findings were achieved as followings. Firstly, key direct factors were soil S, ratio of sodium-potassium (Na/K), the chemical index of alteration (CIA), Soil Organic Matter in soil solution (SOM), Na and total nitrogen in soil solution (TN), and key indirect factors were Carbonate (CO3), Mg, pH, Na, S, and SOM. Secondly, path model can effectively determine direction and quantities of potassium status changes between Exchangeable potassium (eK), Non-exchangeable potassium (neK) and water-soluble potassium (wsK) under influences of specific environmental parameters. In reversible equilibrium state of , K balance state was inclined to be moved into β and χ directions in treatments of potassium shortage. However in reversible equilibrium of , K balance state was inclined to be moved into θ and λ directions in treatments of water shortage. Results showed that the proposed path model was able to quantitatively disclose moving direction of K status and quantify its equilibrium threshold. It provided a theoretical and practical basis for scientific and effective fertilization in agricultural plants growth.


Introduction
The status and transformation of potassium in the soil is of great significance to crop growth [1,2]. Investigations on potassium status gained a numerous achievements [3,4]. The effectiveness of potassium in soil is controlled by four forms, e.g., mineral potassium, non-exchangeable potassium (neK), exchangeable potassium (eK) and water-soluble potassium (wsK) which can be transformed into each other [5][6][7][8]. Absorption capacity of crops and application of farming fertilizer have influence on the potassium forms, thus impact on the release and fixation of potassium in the soil [9]. Within these different forms of soil potassium, there exists a complex and dynamic chemical balance, and it greatly depends on the situation of each form and its environment conditions [10][11][12][13]. However, there are few experimental studies on ternary systems involving K + in spite of the plant nutritional importance of K + in soils [3,5,14]. Research indicated that the amount of K ions retained on the solid phase of soil was much larger than that dissolved in the soil solution. It suggested that a root segment become easily permeable to Ca 2+ followed by Na + and K + with aging [15,16]. Thus, it is indispensable to take into account the cation exchange processes in modeling the response of the soil solution to fertilizer application and transport of cationic solutes in soils.
The transform of soil non-exchangeable potassium to soil exchangeable potassium and soil soluble potassium is a slowly process. It is hard to be determined by using the routine method of ion exchange, although it shows that sodium tetraphenylborate method is more accurate than conventional methods to reflect changes in soil potassium through long-term field fixed experiments of potash fertilization [17,18]. Path analysis is a statistical technique that distinguishes coefficient and causation by partitioning coefficients into direct and indirect effects. Researchers used path analysis to analysis phosphorus retention capacity in allophanic and non-allophanic andisols and soil organic matter effects on phosphorus sorption is successful [19,20]. The method of path coefficients, which was proposed by , was effective in disclosing relationships between variables by diagrams [21][22][23][24][25]. However, path coefficients approach has not yet applied to the study on soil potassium. The potassium absorb efficiency of plant root is affected by multiple factors [26], so it is not suitable to make estimation by use a single parameter [27][28][29]. The above mentioned studies focused on qualitatively descriptive analysis of single factor or multifactor, and failed to quantitatively disclose the dynamic process of potassium status [30].
In this study, we conducted soil K ions exchange experiments for soil samples and further calculated their selectivity coefficients, and explored potassium status on basis of path model to determine dynamic changes of potassium status and specific environmental parameters. The model assumptions were related to the adsorption mechanism at molecular level. The first step, we attempted to examine different indicators of potassium in the soils, and then to analyze direct and indirect correlations between water-soluble potassium (wsK), non-exchangeable potassium (neK) and exchangeable potassium (eK), and further to calculate path coefficients of different potassium status, finally to construct path model for describing the changes of potassium status. The model was calibrated and validated by using 48 cotton rhizosphere/ nonrhizosphere soils samples to test its accuracy. We quantitatively investigated on interaction process of plant-soil-microorganisms and relationship between environmental parameters and potassium release through path coefficient model, and the model was used to predict the multi-component ion exchange equilibrium in soil.

Materials and Methods
Hereby, I, along with coauthor, confirm that no specific permissions were required for our experiment locations/activities since this experiment field belongs to our institute and for scientific research only. And the field studies did not involve endangered or protected species.
In this study, meadow soil was selected as experimental agents in order to investigate potassium nutrient supply and physiological mechanism (Table 1). Rhizosphere soil samples were collected from 90 days cotton plant. When collecting soil samples, we firstly loosed root zone soil to collect rhizosphere soil, and uprooted whole plant cotton, and then gently shook off root zone soil, and finally get root surface adhesion soil. The non-rhizosphere soil samples were collected 10-15 cm depth from surface.
For purpose of this study, we proposed a detail research scheme to evaluate status changes of soil potassium from path model (Fig. 1). We conducted soil K ions exchange experiments for soil samples, and further calculated their selectivity coefficients, and explored potassium status on basis of path model to determine dynamic changes of potassium status and specific environmental parameters. Soil samples were Non-Rhizosphere soil in OPTimum of K of High efficiency genotype cotton (NROPTH); Non-Rhizosphere soil in OPTimum of K of Low efficiency genotype cotton (NROPTL); Non-Rhizosphere soil in Shortage of K of High efficiency genotype cotton (NRSKH); Non-Rhizosphere soil in Shortage of K of Low efficiency genotype cotton (NRSKL); Non-Rhizosphere soil in Shortage of Water of High efficiency genotype cotton (NRSWH); Non-Rhizosphere soil in Shortage of Water of Low efficiency genotype cotton (NRSWL); Non-Rhizosphere soil in Shortage of Water and K of High efficiency genotype cotton (NRSWKH); Non-Rhizosphere soil in Shortage of Water and K of Low efficiency genotype cotton (NRSWKL); Rhizosphere soil in OPTimum of K of High efficiency genotype cotton (ROPTH); Rhizosphere soil in OPTimum of K of Low efficiency genotype cotton (ROPTL); Rhizosphere soil in Shortage of K of High efficiency genotype cotton (RSKH); Rhizosphere soil in Shortage of K of Low efficiency genotype cotton (RSKL); Rhizosphere soil in Shortage of Water of High efficiency genotype cotton (RSWH); Rhizosphere soil in Shortage of Water of Low efficiency genotype cotton (RSWL); Rhizosphere soil in Shortage of Water and K of High efficiency genotype cotton (RSWKH); Rhizosphere soil in Shortage of Water and K of Low efficiency genotype cotton (RSWKL). Firstly, we divided experiment soil into 8 treatments, and then implemented 256 soil samples test (Fig. 2). We measured 33 parameters of the different attributes of 256 soil samples. Those parameters were selected based on standards of Hashimoto and Kang research methods of phosphorus [19,31]. Humic acid (NHA) in soil which can be extracted by NaOH solution (0.1 mol?L 21 , pH = 3), humic acid (PHA) in soil which can be extracted by Na 2 P 2 O 4 solution (0.05 mol?L 21 , pH = 9.2), humus linked to iron (HMi) and humic linked to clay (HMc). Key experimental sampling test datasets were listed in Table 2-Table 6. Secondly, statistical analysis was further implemented to investigate relationship between environmental parameters and potassium status. Datasets were normalized, followed by correlation analysis (Table 7), 8 out of 16 parameters were selected thus to determine the most important parameters. Path coefficients of absolute value of represented the size effect on potassium morphology change. The size of the ''+'' meant the same as the arrow direction with arrow, ''2'' represents in contrast to the arrow direction with arrow. Finally, path model of potassium status was constructed on basis of direct and indirect correlated parameters. In this model, ''e1'' represents unknown variable and its impact factor, the straight line with arrows stands for direct impact factor, and double arc arrow is direction of interaction between parameters.

Results
Path analysis of water-soluble potassium (wsK) For purpose of path analysis of water-soluble potassium, direct and indirect impact factors of the path analysis were derived from multiple linear regressions coefficients of wsK, and correlation coefficients between soil properties. Direct coefficient variables of wsK include neK, eK, PHA, SOM, S, Na and Na/K; while indirect coefficient variables of wsK include CIA, Na, Mg, CO 3 and pH. Impact factor of path size of water-soluble potassium was determined by pH, CIA, Oxidation Reduction Potential (ORP), S, Na ( Fig. 3 and Table 8). In the path model of wsK, direct effects of soil properties on the normalized value of wsK (ZwsK) were represented by single-headed arrows, while coefficients between soil properties were represented by double-headed arrows. Direct and indirect effects were indicated by value and marked with ''6''.   The direct effects of soil properties on the ZwsK were termed path coefficients and were standardized partial regression coefficients for each of the soil properties in the multiple linear regression against the ZwsK. Indirect effects of soil properties on the ZwsK were determined from the product of the simple coefficient between soil properties and the path coefficient (i.e., one two-headed arrow and one single-headed arrow). The coefficient between the ZwsK and soil property was the sum of the entire path connecting two variables, as described by r neK,wsK~PwsK,neK zr neK,PHA : P wsK,PHA z r neK,S : P wsK,S zr neK,Na=K : P wsK,Na=K z r neK,CIA : P wsK,CIA zr neK,Na : P wsK,Na z r neK,SOM : P wsK,SOM zr neK,eK : P wsK,eK  where, r ij is the simple correlation between the ZwsK and a parameter of soil property, P ij is the path coefficient between the ZwsK and a parameter of soil property, and r ij P ij is the indirect effect of a parameter of on the ZwsK. An uncorrelated residue (e wsK ) that represented the unexplained part of an observed variable in the path model was where, R wsk 2 is the coefficient of determination of the multiple regression equation between the ZwsK and the eight soil properties. Those parameters of soil properties were neK, PHA, S, Na/K, CIA, Na, SOM, eK. Backward regression analysis was performed with the ZwsK as the dependent variable and the eight soil properties were regarded as independent variables. Backward regression was a multiple regression procedure in which all the independent variables were entered into the regression equation at the beginning, and variables that did not contribute significantly to the fit of regression model were successively eliminated until only statistically significant variables remain. Fig. 4 and Fig. 5 were similar expressions.
In the path of water-soluble potassium, path contribution coefficients were quite different. Coefficient values were 2.10, 21.52 and 0.20 for wsK to neK, neK to wsK, and other variables to wsK respectively. It disclosed that most part of wsK was changed into neK (with values of 2.10). In the path of neK to wsK, coefficient was negative (with value of 21.52); it indicated that there was an inverse change from wsK to neK, a lot of wsK was changed into neK. Those dynamic changes indicated that wsK and eK reached equilibrium in soil solution, and this equilibrium represented its quick reaction process.
These changes impacted greatly on nutrients dynamic equilibrium of rhizosphere soil. This changed nutrients recycling process of the soil environment since continuous-release of humic acid had direct effects on absorption and uptakes of those variables including neK, eK and PHA, SOM, S, Na and Na/K. It was clear that rhizosphere effects were likely to result in changes in root exudates, pH and ORP (eV). Consequently, it led to changes in pH values and population, quality and activity of microorganisms while pH change was caused by root exuded organic acids. Under conditions of different pH values, humic acid significantly correlated to adsorption and desorption of potassium. Humic acid of potassium adsorption and desorption presented an upward trend when it increased in initial concentration (pH was 4.0-8.0), however, desorption rate declined. The migration distance also showed very significant linear relationship within ranges of migrations between water-soluble potassium and exchangeable potassium. Organic matter nutrients and chemical efficiency of humic acid altered different variables such as physical, chemical, and biotic in root-soil interface.

Path analysis of non-exchangeable potassium (neK)
For path analysis of non-exchangeable potassium, we further investigated its direct and indirect impact factors. Results showed that direct coefficient variables of neK include eK, wsK, SOM, S, Na, Na/K, TN, CIA; and indirect coefficient variables of neK include Mg, CO 3 . Impact factor of path size of neK was determined by eK, CO3, pH ( Fig. 4 Table 9). Table 9 showed there was significant correlation between eK and neK (with coefficient of 0.81). In the Path of non-exchangeable potassium, contribution coefficient was 0.01 for wsK to neK, although wsK was determined by general effects of N, Na, S and CIA. The SOM was also an important contributor to eK and wsK of rhizosphere soil. Coefficients were 0.69 and 0.55 between SOM and eK, SOM and wsK, respectively.
The correlation between the ZneK and a soil property was the sum of the entire path connecting two variables, as described by where, r ij is the simple correlation coefficient between the ZneK and a soil property, P ij is the path coefficient between the ZneK and a soil property, and r ij P ij is the indirect effect of a soil property on the ZneK. An uncorrelated residue (e) that represents the unexplained part of an observed variable in the path model was where R nek 2 is the coefficient of determination for the multiple regression equation between the ZneK and those soil properties. Those soil properties were wsK, SOM, S, Na/K, TN, Na, SOM, eK.
It disclosed that plant nutrition was indirectly influenced by nutrients release of microbial decomposition and organic matter. It was particularly obvious in nutrition release of low molecular weight organics in rhizosphere soil. High contents of associated ions (Na + ) were limit factors that control crop yields and soil nutrients use. Ammonium nitrogen was likely to reduce fixation of potassium ions of soil, and increase risk of potassium leaching. NH 4 + and K + were competitive ions for absorption of ORP (eV) of eK and neK in rhizophere soil, and different status of nitrogen effects on efficient use of potassium. Therefore, suitable proportion of N and K is important for plant. Soil nitrogen availability determines potassium absorption of root surface and cell membrane. Root membrane could directly affected by NH 4 + , NO 3 2 , K + when it absorbs potassium ions, and indirect affected by potential balance from assimilation of NH 4 + . Strong affinity of K + , high efficiency and fast speed of absorption of potassium are likely to improve efficient absorption of potassium.

Path analysis of exchangeable potassium (eK)
As to exchangeable potassium, it was quite different situation in comparison with those of water soluble and non-exchangeable potassium. It showed that direct coefficient variables of eK include neK, wsK, SOM, S, Na, Na/K, TN, CIA, HMi, PHA; and indirect coefficient variables of eK include Mg, CO 3 . Impact factor of path size of eK was determined by neK, CO 3 , pH, ORP, Na (Table 10).
From Fig. 5, it showed that neK was closely related to wsK (with correlation coefficients of 0.72), and there was an inverse proportion between ZneK and ZwsK (the normalized value of neK and wsK) in condition of changes of micro-bioorganic and geochemical environmental factors in rhizosphere soil. Path contribution coefficients were 1.10 and 0.04 for neK to eK, and wsK to eK. It suggests that neK is more likely to be changed into eK while the possibility of wsK to eK is relatively low.
The correlation between the normalized value of eK (ZeK) and a soil property is the sum of the entire path connecting two variables, as described by r neK,eK~PeK,neK zr neK,SOM : P eK,SOM z r neK,Na=K : P eK,Na=K zr neK,Na : P eK,Na z r neK,HMi : P eK,HMi zr neK,TN : P eK,TN z r neK,wsK : P eK,wsK zr neK,S : P eK,S where r ij is the simple correlation coefficient between the ZeK and a soil property, P ij is the path coefficient between the ZeK and a soil property, and r ij P ij is the indirect effect of a soil property on the ZeK. An uncorrelated residue (e) that represents the unexplained part of an observed variable in the path model was Where, R eK 2 is the coefficient of determination of the multiple regression equation between the ZwsK and the eight soil properties. Those soil properties are neK, SOM, Na/K, Na, HMi, TN, wsK, S. Backward regression analysis was performed with the ZeK as the dependent variable and the eight soil properties as independent variables. Backward regression is a multiple regression procedure in which all the independent variables are entered into the regression equation at the beginning, and variables that do not contribute significantly to the fit of the regression model are successively eliminated until only statistically significant variables remain.
Cation exchange capacity (CEC) or the number of exchangeable cation depends on the type and number of clay and organic matter in the soil. To a certain extent, TN, SOM, HMi and the PHA and the CIA are able to well reflect type and number of clay and organic matte in rhizosphere soil. Proper pH value of solution, hot water pressure, concentration increases of magnesium carbonate and ion are useful to exchangeable potassium (neK) when it is absorbed by crops. A plenty of CO 2 could be produced by root respiration and soil microbe respiration, and is fixed in form of CO 3 22 . For stable pH values, it is beneficial for organic acid to reach in equilibrium of absorption in soil solution. The exchangeable potassium and fixed potassium which comes from layers of secondary clay mineral are released when there occurs hydrolysis and cation substitution resulted by chemical element H + . Those variables, such soil organic matter (ZSOM), ironbinding Humic (ZHMi), and water-soluble nitrogen (ZTN), not only provide nutrition for cotton root, but also changed soil chemical compositions and absorption characteristics of potassium in rhizosphere soil. Additionally, concomitant phenomenon (ZS) appears between S (the lost chemical element in soil) and potassium. Antagonism occurred between Na and K, and it implied that sodium salt was able to have significant effect on plant (ZNa/K, ZNa).
From above analysis, we conclude that key direct factors of wsK, neK, and eK were S, Na/K, the CIA, SOM, Na, and TN, and key indirect factors were CO 3 , Mg, pH, Na, S, and SOM. Whereas, direct impact factors of water soluble potassium (wsK) were PHA, S, Na, Na/K, K, the CIA, SOM, eK and neK.

Path model of soil potassium status changes of cotton
It is of great importance to further investigate on the equilibrium movement in soil. The following part will discuss the shift of dynamic balance among wsK, eK and neK, and their flow path. To test the reliability of the models, we collect another 48 cotton soil samples in the same experiment, using stepwise regression method, observation of potassium mobility and morphological change.
We carried out a plot experiment of potash application with high K-efficiency genotype (HEG) and low K-efficiency genotype (LEG) cotton. The soil in each treatment weights 8.5 kg. Detail fertilizer schemes are listed in Table 11. LEG and HEG cottons' planting effects is shown in Fig. 6.
(1) Equilibrium shift in the path for non-exchangeable potassium Supposed that V wsK R nek is path for wsK to neK, V ek R nek is path for eK to neK, therefore, equilibrium shift in path for nonexchangeable potassium can be expressed by contribution coefficient functions as below,   samples RSWKH, RSKH, RSKL, ROPTH, ROPTL, NRSWKH, NRSWH, NROPTH, NROPTL, K balance cycle were moved in b and c direction movement; soil samples NRSKH, NRSKL, K balance cycle were moved in b, x direction movement. In Fig. 6, neK reaches its equilibrium in following soil types which includes LEG in water stress treatment (SW), HEG after potassium stress treatment (SK), LEG after optimum fertilization treatment (OPT), and LEG after water stress treatment (SW). However, all the value in those soils with HEG in SW, LEG in SK, HEG in SW and SK, the values were all greater than the equilibrium, wsK reached its saturation. For other treatment soils, when wsK was lower than the equilibrium and PHA and eK was greater than equilibrium, neK changed into wsK. It proves that eK equales to total K in matured crops. When eK is below a certain level in soil, plants can no longer get eK. When eK was at low levels, it was absorbed by other ions with stronger absorption. This reduces chance of potassium to enter into the soil solution. In this condition, neK begins to release through low concentrations of wsK and eK. However, the release capacity of neK depends on situations of crop types (2) Equilibrium shift in the path for exchangeable potassium Supposed that V wsK R ek is path for wsK to eK, V nek R ek is path for neK to eK, therefore, equilibrium shift in path for exchangeable potassium can be expressed by contribution coefficient functions as below, Equilibrium occurs when x wsK , x PHA , x neK satisfy with Functions 11, x wsk~1 :64 : x neK z0:61 : Here, y ek is the equilibrium threshold. And the value of x wsK , x PHA , x neK is determined by soil characteristics. Calculated results from Equations 9 and 10 are listed in Table 12. In this table,  column 5  As is shown in Fig. 7, neK reaches its equilibrium in following soil plots, including LEG (SW), HEG (SK), LEG (OPT), and LEG (SW). And the exchangeable K reaches its equilibrium in soil types of LEG (SW). For other treatment soils, such as LEG (SW and , neK changes to eK more or less. In soil of HEG (SW and SK), eK changes to wsK when wsK was lower than equilibrium. In soil of LEG (SW and SK), neK changes to eK when wsK of soil was greater than that of crops. Soils of both genotypes cotton (HEG and LEG) can efficiently absorb exchangeable potassium.
During a lot of vermiculite in the soil, those negative charges caused by vermiculite isomorphous substitution are near p-site, the electrostatic attraction of potassium ion became bigger, the adsorption capacity much more than other type (2:1) minerals, therefore, wsK and neK are fixed and the soil changed into deficiency. In general, drying accelerated the fixation of wsK adsorption, but neK still moved in the direction of eK in the reversible equilibrium in HEG rhizosphere soil because rhizosphere soil exchangeable potassium content was low to meet the needs of the cotton growth in the SW and SWK plant. When moisture was adequate, redox reaction of soil became strengthened and oxidation reduction potential was reduced. For example, Fe 2+ and Mn 2+ ions were increased rapidly, Fe 2+ and Mn 2+ replaced p-site potassium of soil colloids and part of i-site potassium, and this improved the content of rapidly-available potassium and biological effective parts. In addition, the dissolved iron and manganese mineral alteration and also increased the release of mineral K.

Discussion and Conclusions
In this study, changes of potassium status in soil were investigated through our proposed path model on basis of laboratory experimental datasets. Through this investigation, following findings were achieved.
Firstly, changes of potassium status in the rhizosphere soil were controlled by different environmental variables. wsK, neK and eK were controlled by some common factors and have their speculiar characteristics respectively. Those parameters, such as soil organic matter (SOM), Na and Na/K and S, were significant direct coefficient factors which control dynamic process of wsK, neK and eK. CO 3 and Mg are significant indirect coefficient factors which control dynamic of wsK, neK and eK. The pH is the primary factor which controls path of the wsK, neK and eK. Organic acids of HEG root system secretion increased soil acidity; the amount of fixed potassium was reduced. PHA surface is porous structure, PHA increased the adsorption of wsK in the experiment, and HMi could promote the moving to wsK, eK in the K form exchange. Both of them were negative correlation significantly negative correlation. First, under acid condition, hydroxyl aluminum and aluminum ions can occupy the potassium fixation point stronger. At the same time, the larger diameter of hydroxy aluminium ions entered into the interlayer of mineral, formed as the ''island'' and played a supporting role, provided the potassium ion diffusion of mobile channel in the interlayer. These reduced the soil potassium fixation. Second, water protonated in acidic condition, the radius of H 3 O + 12-13.3 nm were similar to K + , H 3 O + can produced competitive adsorption in soil where fixed on K + parts, so potassium reduced fixation. Those dynamic processes, such as eK to wsK, low concentration wsK and eK to neK, are likely to occur, in condition that wsK is lower than equilibrium value and PHA and eK are greater than equilibrium values. eK is changed neK when wsK is greater than equilibrium value. In soil of HEG (SW and SK), eK is changed into wsK when plant wsK is lower than soil wsK. In soil of LEG (SW and SK), neK is changed into eK when plant wsK is higher than soil wsK. Both genotype cottons can effectively absorb eK in treatment of stress of potassium (SK).
Secondly, it discloses that SOM, S, TN, HMi and PHA in rhizosphere soil determined dynamic balance of potassium status which includes adsorption and desorption, precipitation and dissolution, complexion and chelating. Special characteristics of humus, e.g., acidic and hydrophilic, cationic exchange, complexing capacity and high absorption capacity, could improves the sustained-release effect of potassium fertilizer. Many variables, pH, ORP, etc., are closely related to transformation and utilization of potassium fertilizer. The pH change in soil is mainly caused by coupling effects of nutrient uptake by plant roots and the secretion of organic acid and absorption imbalance of positive ion and negative ion leads to rhizosphere pH change, however, different factors affect imbalance of absorption of positive ion and negative ion. Dynamic equilibrium of potassium is affected by equilibrium constants, temperature, and products of this equilibrium system (wsK). Therefore, it is useful for cotton planting to improve to appropriate temperature, to alter SOM, TN and pH in soil, and to reduce the equilibrium constant to positive reaction. Meanwhile, it is good to choose HEG cotton for better use soil wsK. It is important to understand ways to slow down the fast response and accelerate speed of slow response in potassium balance system. Finally, path model can effectively determine direction and quantities of potassium status changes between eK, neK and wsK under influences of specific environmental parameters. The proposed path model was able to quantitatively disclose moving  direction of K status and quantify its equilibrium threshold. It provided a theoretical and practical basis for scientific and effective fertilization in agricultural plants growth. The significance of this study is that we are able to implement investigation on the equilibrium movement of potassium status through path model analysis. It discloses the gradual dynamic process of potassium status, and to decouple the sophisticated interactions between different variables. It is useful to guide the use of potassium fertilizer for cotton crops in practical. The model provides a biogeochemical theory basis of controlling moisture in cotton plant, applying potash fertilizer and organic fertilizer, and appropriately reducing nitrogen fertilizer. However, environmental parameters such as soil moisture and temperature are some of the most important factors that influenced crop growth. Observation of dynamic change of soil moisture and temperature in every other time-step is quite sophisticated for us to design the control experiments while it involves in more than 4 times treatments, which means that we have more than 500 pots. By limitation, we had to set the soil moisture into two states, water sufficient (.35%) and deficient (,25%), and temperature were set to indoor temperature.  Anyway, we were very aware of the importance of soil moisture and temperature impacts on control experiment potassium status. The next step of our work is to further investigate efficient utilization of potassium and the influence of soil moisture and temperature on exchange of potassium forms. We hope this study is able to gain knowledge of promoting high efficient use of potassium for crops.

Supporting Information
Appendix S1 Explanation of symbols. (DOC)