Effects of Phosphorus Ensembled Nanomaterials on Nutrient Uptake and Distribution in Glycine max L. under Simulated Precipitation

Nanoscale hydroxyapatite (nHA) was synthesized to investigate its potential as a phosphorus (P) ensembled nanofertilizer, using soybean (Glycine max L.) as a model plant. The conventional analogue phosphate (pi) was used for comparison with the synthesized nHA. Varied precipitation intensities (0%, 30%, 60%, and 100%) were simulated by adding selected volumes of the P fertilizers (nHA or pi) via foliar spray and soil amendment. The total amounts of added P were the same across all the treatments. The importance of a wash-off effect was investigated on foliar-treated seedlings by evaluating different watering heights (20, 120, and 240 cm above the seedlings). Fresh weight, pigment content, macro-, and micronutrient contents were measured in soybean tissues across all the treatments after 4 weeks of greenhouse cultivation. The synthesized nHA showed superior effects on plant nutrient content upon high precipitation intensities. For example, at 100% precipitation intensity, there was 32.6% more P and 33.2% more Ca in shoots, 40.6% more P and 45.4% more Ca in roots, and 37.9% more P and 82.3% more Ca in pods, as compared to those with pi treatment, respectively. No impact on soybean biomass was evident upon the application of nHA or pi. Further investigation into customizing nHA to enhance its affinity with crop leaves and to extend retention time on the leaf surface is warranted given that the present study did not show significant positive impacts of nHA on soybean growth under the effects of precipitation. Taken together, our findings increase understanding of the potential application of nHA as a nano-enabled fertilizer in sustainable agriculture.


Introduction
The ability to provide adequate and nutritional food worldwide will be challenged by the rapidly increasing global population, declining arable land, changing climate, and degradation of environmental resources [1][2][3]. The intensive application of chemical fertilizers has been a common and effective strategy in agriculture for decades [4,5]. However, fertilizer and pesticide efficiency of delivery and utilization is quite low, and as such, they have often been applied excessively to meet crop demands [6][7][8]. Importantly, over application of these agrochemicals has direct negative environmental consequences, such as toxicity to non-target species, contamination of soil and water, as well as disturbance to the plant diversity [9][10][11][12]. For example, the nutrient use efficiency of phosphorus is only 5-30%, and excessively applied amounts run-off into water bodies, serving as a primary driver of eutrophication and dangerous algal blooms [13][14][15]. Phosphorus has the added complication of being a finite resource, with declining sources of phosphorus rocks becoming an increasing concern [16][17][18][19]. Consequently, the environmental footprint of food production has become both massive and unsustainable. Given the challenges noted above, novel approaches and technologies to sustainably increase agricultural productivity are desperately needed [20][21][22][23].
There has been significant growing interest in the use of nanotechnology to enhance crop growth and agricultural production [24][25][26]. Potential applications of engineered nanomaterials in agriculture include, but are not limited to, nanofertilizers [8,27,28], nanopesticides [8,29], nano-enabled agrochemical carriers [30,31] and nanosensors [32][33][34]. Selected engineered nanomaterials have been applied as soil amendments, seed coatings, and foliar sprays [35,36]. To meet the goal of controlled release and delivery of active ingredients with reasonable environmental cost and energy input, nanomaterials need to be carefully designed and applied by sustainable methods [23,35,37]. For nanomaterials applied as fertilizers, there has been work on both nanoscale macro- [38][39][40][41] and micronutrients [42][43][44][45]. Importantly, material properties such as size, morphology, composition, and surface characterization all can be tuned to control release, accumulation, and translocation of nutrients in plants [46][47][48][49][50][51]. The benefits of using properly-controlled nanomaterials as fertilizers extends not only to precise and accurate delivery, but also to alleviating the negative environmental consequences derived from overapplication and run-off [27,52]. Given its extensive use in biomedical applications, nanoscale hydroxyapatite (nHA) has garnered significant interest as a potential phosphorus fertilizer [53,54]. The synthesis of rod-like hydroxyapatites can be readily achieved through one-step wet chemistry precipitation [55][56][57]. The synthesized nHA has shown positive effects on soybean, leading to 32.6% and 20.4% increases in growth rate and seed yield, respectively [38]. Exposure to 150 mg P/kg nHA with a negative charge increased sunflower biomass by twofold relative to plants treated with the same amounts of P in the conventional form of triple superphosphate [58]. Recently, nHA has been incorporated into a range of nanocomposites to achieve multifunctionality [40,59,60]. For example, hydroxyapatite-based nanostructured P fertilizers were achieved through an alkali-enhanced hydrothermal process; these materials exhibited significantly elevated phosphorus use efficiency (PUE) of 45.87% and 46.21% from two different precursors (CaHPO 4 and CaP 2 O 7 ) as compared to a chemical P fertilizer with a PUE of 23.44% [40]. Similarly, soil amendment of nHA functionalized with humic substances significantly improved early Zea mays growth, productivity, rhizosphere bacteria diversity, and tolerance to NaCl-induced abiotic stresses [60].
In the current study, nHA was synthesized using a wet chemistry precipitation method. The nanoscale fertilizer was applied to soybean foliage and roots simultaneously to simulate different precipitation intensities; comparisons to equivalent conventional bulk amendments were included. Additionally, the foliar retention ability of the synthesized nanofertilizer was evaluated through a wash-off test in which precipitation intensities were controlled at different watering heights (20-240 cm). The present study provides valuable information on the efficacy of nanoscale P materials as a component of sustainable nano-enabled agricultural practices. Moreover, discussion on the costs and benefits of both nanoscale and conventional fertilizers through the two application methods provides insight into the rational design of smart agrochemicals.

Synthesis and Characterization of Hydroxyapatite Nanoparticles
Hydroxyapatite nanoparticles (nHA) were synthesized using a wet chemistry precipitation method by mixing calcium (II) and phosphate (1.67:1, mole weight ratio) in an alkaline solution under controlled temperatures (80 • C). Briefly, 40 mL of 1.0 M Ca(NO 3 ) 2 solution was incubated in an water bath with continuous stirring. Then, 40 mL of 0.6 M NaH 2 PO 4 aqueous solution was added drop-wise into the cold Ca (NO 3 ) 2 solution. The solution pH was adjusted to approximately 12.0 by adding 2.0 M NaOH, and the solution was continuously stirred for 24 h. The resulting white precipitate was washed with deionized (DI) water to remove extraneous salts and alkaline materials. The solution was freeze dried to enable recovery of the hydroxyapatite nanoparticles. The resulting white precipitates were observed by transmission electron microscopy (TEM, JEOL JEM-2000FX, Tokyo, Japan), and functional group characterization was evaluated by Fourier Transform Infrared Spectroscopy (FT-IR spectrometer, PerkinElmer Spectrum, Waltham, MA, USA).

Greenhouse Cultivation
Soybean seeds were germinated in a mixture of vermiculite and PRO-MIX (1:3, v/v), and then seedlings were transplanted into soil (300 g/pot) after the two unifoliate leaves were fully developed. The soil was collected from the Crop and Animal Research and Education Farm at the University of Massachusetts, Amherst. Five replicate pots were used for each treatment. Plants were watered daily with 50 mL tap water through the trays under the pods. The greenhouse temperature was set at 22 • C for 16 h during the day and 18 • C for 8 h at night. Soybean seedlings were treated with nHA (P-nanofertilizer) or dicalcium phosphate (pi) through foliage (foliar spray) and root pathways (soil amendment) once per week over a 28-day growth period. The equivalent molar concentration of phosphorus (0.7 mmol/L) was used for both nHA and pi.

Simulated Precipitation Experiment
In order to simulate different precipitation intensities (0%, 30%, 60%, and 100%), 4 mL of prepared P-nanofertilizer or pi was sprayed onto the leaves and amended into the soil simultaneously. For the foliar spray, soil in the pots was covered with paper towels to avoid spraying onto the soil. For soil amendment, P-nanofertilizer or pi was pipetted through three pre-inserted tips. The pipette tips were inserted into soil (1 cm underneath) around the seedlings with equal distances. Detailed information for each treatment is listed in Table 1; A, B, C, and D are treatments with nHA, and A', B', C', and D' are treatments with pi, respectively. The simulated precipitation experiment was schematically illustrated in Figure 1A.

Wash-Off Test from Different Heights
In order to evaluate leaching upon precipitation, a separate experiment with varied precipitation heights (20 cm, 120 cm, and 240 cm) was conducted. A wash-off test was conducted one day prior to the next round of P fertilizer application. Five hundred millilitres of tap water was applied to each pot over the course of the experiment. Details of these treatments are listed in Table 2, in which E, F, and G represent the wash-off tests with nHA, and E', F', and G' represent the wash-off tests with pi. A and A' from the simulated precipitation experiment were used as no wash-off controls, respectively. The wash-off test is schematically illustrated in Figure 1B.  Table 2. Details of each treatment in the wash-off test.

Label Treatment Details
Control (CK) 4 mL water sprayed onto foliage A (no wash-off) 4 mL nHA (0.7 mmol/L) sprayed onto foliage A' (no wash-off) 4 mL pi (0.7 mmol/L) sprayed onto foliage E (20 cm) A with precipitation from 20 cm above seedlings F (120 cm) A with precipitation from 120 cm above seedlings G (240 cm) A with precipitation from 240 cm above seedlings E' (20 cm) A' with precipitation from 20 cm above seedlings F' (120 cm) A' with precipitation from 120 cm above seedlings G' (240 cm) A' with precipitation from 240 cm above seedlings

Plant Harvest
At harvest, soybean seedlings were washed with tap water, followed by deionized water rinsing three times. Intact leaves were washed with 1 mM HNO 3 and rinsed by ice cold 2 mM CaCl 2 and DI water to separate unabsorbed NPs from the leaves [61][62][63]. The fresh biomass of shoots and roots were separately weighed across all the treatments. Fresh shoots, roots, and pods were collected for pigment and elemental analysis.

Elemental Analysis in Plant Tissues
Approximately 100-400 mg of fresh tissues of pods, roots, or shoots of soybean were cut into pieces and transferred in vials, followed by freeze-drying in a lyophilizer (Free Zone 2.5 Liter-50C; Labcon co, Kansas, MI, USA) for 72 h. The mass difference between the fresh biomass and dry biomass was used to calculate water content. Freeze-dried tissues were ground into fine powder, and then approximately 50 mg of dry tissues were transferred into a digestion tube containing 3 mL of concentrated HNO 3 . The samples were pre-digested at the ambient temperature overnight. Then, the pre-digesta were heated at 105 • C on a hot block (DigiPREP System; SCP Science, Baie-d'Urfé, QC, Canada) for 40 min; after the digesta cooled to an ambient temperature, 500 µL of H 2 O 2 was then added prior to heating at 105 • C for an additional 20 min to complete the digestion. A volume of 300 µL of digesta was diluted to 30 mL with DI water prior to elemental analysis. The content of macro-(P, K, Ca, Mg, and S) and micronutrients (Cu, Zn, Mn, and Fe) was determined by inductively coupled plasma optical emission spectrometry (ICP-OES, iCAP 6500; Thermo Fisher Scientific, Waltham, MA, USA).

Statistical Analysis
Each treatment contained four or five replicates. All soybean pots were randomly arranged on greenhouse benches and rearranged on a weekly basis. SPSS statistics software (version 26, IBM Corp, Endicott, NY, USA) was used for the significance analysis. Two sets of one-way analyses of variance (one-yay ANOVA) followed by Waller-Duncan multiple comparison were used to perform statistical analyses for treatments in the simulated precipitation experiment and wash-off test separately. In figures and tables, values of each assay that are followed by different letters are significantly different at p < 0.05.

Characterization of the Synthesized Nano-Hydroxyapatites
Nanoscale HA particles exhibited a needle-like shape with dimensions of approximately 20 nm in diameter and 100-200 nm in length ( Figure 2A). In the FTIR spectrum ( Figure 2B), OH-bands appeared at around 3400 cm −1 . Bands that appeared at 1023 and 962 cm −1 represent asymmetrical and symmetrical stretching of P-O bond of the PO 4 3 , respectively; the peak at 860 cm −1 corresponds to the symmetrical stretching mode of the P−O(H) bond of the HPO 2 4− group [65,66]. The band at 1414 cm −1 indicates the presence of carbonate groups.

Physiological Responses of nHA-and pi-Treated Soybean under Simulated Precipitation and Wash-Off Conditions
As shown in Figure 3A,B and Figure S1A, externally applied P fertilizers had a slight but statistically insignificant enhancing effect on soybean biomass over the 28 days of treatment; differences between the nHA and pi treatments at the same precipitation intensities were not significant. This result is consist with the previous report by Liu and Lal, in which synthetic apatite nanoparticles only slightly improved soybean growth (biomass) relative to triple-phosphate [38]. In the wash-off test, the fresh biomass of nHAand pi-treated soybean were significantly decreased ( Figure S3) as compared to the control. Specifically, the root biomass of nHA-and pi-treated soybean was decreased by up to 27.6% and 21.0%, respectively, as compared to the control. A reasonable explanation for this effect could be the loss of surface-applied P fertilizer and successive leaching of nutrients from the soil. Foliar spray of the P fertilizers and subsequent wash-off could potentially cause phyllosphere dysbiosis and a subsequent decrease in crop biomass. It is also known that rain and/or high air humidity can induce phyllosphere dysbiosis and disease outbreaks in plants [67][68][69], which could potentially impede seedling growth and development. In addition, in comparison with foliar spray, soil irrigation could result in soil compaction and inhibition of root respiration [70], both of which could negatively affect plant health.  Pigment content is an important parameter for evaluating plant photosynthesis under a range of treatment conditions [71,72]. Leaf chlorophyll content was less than the untreated group under all pi treatments, and it decreased further along with increases in the precipitation intensity. For example, compared to the control, at 0% precipitation intensity, chlorophyll content decreased by 15.3% in pi-treated soybean shoots; at 100% precipitation intensity, chlorophyll contents decreased by 15.7% and 29.2% in nHA-and pi-treated shoots, respectively ( Figure 3C). Reduction effect on the chlorophyll content with pi treatments was stronger than those with nHA treatment under the same precipitation intensity. Specifically, as shown in Figure S1B,C, the slight decrease in chlorophyll contents in nHA-treated soybean shoots was primarily due to the reduction in the chlorophyll a content. A similar trend of decreasing chlorophyll content along with the increasing of precipitation intensity was also evident in the wash-off test ( Figure S3F). In this case, the decrease of chlorophyll content was mostly ascribed to the decrease of chlorophyll b levels.
For example, at a precipitation height of 240 cm, the chlorophyll b content in nHA-and pi-treated shoots decreased by 38.3 and 45.1%, respectively ( Figure S3E). Externally applied nHA and pi had no effects on carotenoid contents in both simulated precipitation and wash-off treated plants. It was reported that chlorophyll a content could be used as an indicator for abiotic stress [71,73,74]. At this point, there seemed to be no positive impact and some potentially negative impacts of external P fertilizers on the soybean photosynthesis.

Macronutrient Contents in Shoots, Roots, and Soil
The macronutrient content in soybean shoots, roots, and soil treated by nHA and pi in the simulated precipitation scenario are shown in Figures 4 and 5 and Figure S2. At 0% precipitation intensity, when all the external P fertilizers were applied through foliage, pi treatment resulted in a 36.9% increase in shoot P contents and a 26.8% increase in root P compared to the untreated group. Conversely, once all the external P fertilizers were amended into soil (100% precipitation intensity), nHA treatment had significantly enhanced P accumulation in the seedlings, while pi showed no effect with statistical difference. For example, the shoot and root P content in the nHA-treated soybeans at 100% precipitation intensity was 44.0% and 24.0% higher than the control, respectively ( Figure 4A,B). Calcium is an important structural element regulating cell wall stability and membrane permeability [75,76] and is another primary component of nHA and pi. Correspondingly, at 0% precipitation intensity, pi-treated soybean seedlings had 79.8% more Ca in shoots and 18.6% more Ca in roots compared to the controls, while nHA-treated soybean seedlings had 23.2% more Ca in shoots and no difference between roots and the control. At 100% precipitation intensity, pi-treated seedlings had 28.5% more Ca in the shoots and 20.3% less Ca in the roots compared with unamended ones, while nHA-treated seedlings had 71.1% more Ca in the shoots and 15.9% more Ca in the roots compared to the respective controls, which are statistically higher than those of the pi treatment ( Figure 4D,E).
nHA and pi treatments exhibited opposite trends on nutrient content as a function of precipitation intensity, both in shoots and roots. First, this could be a result of greater nHA availability in the soil ascribed to less binding to soil minerals that would form insoluble phosphates [77][78][79]; this subsequently results in greater P and Ca bioavailability for root-to-shoot transport. Second, plants require P in the phosphate form; thus, when nHA/pi were sprayed onto the foliar surfaces, longer retention time was required for the uptake of nHA than pi. Consequently, nHA may aggregate before uptake by the leaves, which could subsequently undermine P content [80,81].  Such trends were not limited only to the elements that were externally applied; the macro nutrients K, Mg, and S were increased in shoots and roots with increases in the soil P-nanofertilizer amendment ( Figure 5). Since there was no external K, Mg, and S source in the treated groups, the increase of these macronutrients was facilitated by the soil (Figure S2). Although there were no discernible patterns in the variations among the different treatments with different precipitation intensities, nutrients in the soil also contributed to the plant growth upon treatment with the external fertilizers. For example, at 100% precipitation intensity, 10.0% more P, 17.6% more Ca, 14.6% more Mg, and 54.8% more S were depleted from the soil in the nHA treatment compared to the untreated group. The K content in soil was not altered across all the nHA treatments. Since there were trays under each pot, there was no chance of the nutrient loss from the leachates, and since there was no visible leaching observed throughout the whole experiment period, the depleted soil nutrients were all being taken up by the seedlings. The input of external nanoscale P fertilizers to the soil could be regarded as a means to enhance the uptake of other soil nutrients. Therefore, when the P fertilizers were all added to the soil (at 100% precipitation intensity), 42.8% more K, 61.2% more Mg, and 22.7% more S in nHA-treated soybean shoots could be ascribed to root-to-shoot transport, as compared to the control.
For pi-treated soybeans, similar soil nutrient depletion was also evident. At 0% precipitation intensity, 9.7% P, 13.5% Ca, 14.1% Mg, and 28.3% S were depleted from soil and taken up by the seedlings ( Figure S2). The K content in soil was not altered across all the pi treatments. In contrast to the nHA treatment, shoot and root macronutrient content decreased with the increase in precipitation intensity ( Figures 3B,E and 4B,E,H). At 100% precipitation intensity, the root macronutrient content was notably less than the control (20.3% less Ca, 36.3% less Mg, and 61.1% less S); soil macronutrient content was equivalent with the control except for S, which was reduced by 42.6% compared to the control ( Figure S2E). There was no enhancement effect of externally applied conventional fertilizer on the uptake of other soil nutrients from soil due to inevitable loss of elements by soil binding. In addition, there might have been a potential root-to-shoot mobilization, which could have contributed to the accumulation of nutrients in the shoots. However, the growth of the seedlings could have been impeded due to the lack of root-shoot nutritional balance. Further investigation is needed for clarity. Sulfur is essential for pod formation, and its mobilization is critical to the plant's reproductive life cycle [82,83]. Moreover, it has been reported that N and S assimilation could be altered by changes in the P environment; N and S metabolism require rRNA for synthesis, and rRNA is the largest pool of organic P in leaf cells [84,85]. In the current study, the drastic alteration of S content in the plants could have been a response to the P environment change, as well as the development of pods.
In the wash-off test, soybean seedlings treated with nHA and pi responded differently ( Figure S4). Soybeans treated with pi accumulated the greatest amount of macronutrients in shoots and roots without wash-off. However, the macronutrient contents in both shoots and roots decreased significantly after wash-off, suggesting that much of the pi had been removed from the leaf surfaces. In addition, there were no statistically significant differences in shoot/root macronutrient content among different wash-off heights. Moreover, the mobilization of soil nutrients to the plants was not evident in the wash-off treated groups, with the exception being S, which was decreased by 25.5% at a wash-off height of 240 cm, as compared to the control ( Figure S6). With the treatment of wash-off, more K accumulated in the soil. For example, 24.2% more K accumulated in the soil at a wash-off height 120 cm. For nHA-treated seedlings, wash-off had a different impact. Other than no significant effects at a higher wash-off height (120 and 240 cm), the enhancement effect on the nutrients acquisition in soybean shoots was still evident at the wash-off height of 20 cm. Specifically, 11.9% more P, 37.0% more Ca, 35.5% more K, 41.7% more Mg, and 15.5% more S were accumulated in soybean shoots as compared to the untreated group ( Figure S4), which was even better than the no wash-off group.
Overall, nHA and pi could supply equivalent amounts of P and Ca to plants via root and foliar exposures, respectively. Given that complete breakdown of nHA, either on the leaf surface or in the soil, was not evaluated in this study, a longer-term full life cycle study is needed to fully understand the role of nHA in supplying P and Ca to crops. Furthermore, in a realistic environmental scenario, wash-off of the foliar-applied fertilizers is inevitable and is a key factor that lowers the nutrient availability; additional investigations on ways to increase of nutrient use efficiency through the careful design of nanofertilizers are needed.

Micronutrient Contents in Shoots, Roots and Soil
The content of Cu, Zn, Fe, Mn, and B in soybean tissues across all the treatments varied to different extents as compared to the respective controls ( Figure 6). Unlike the macronutrients, there were few discernible trends in the soybean tissue micronutrients. For the nHA-treated seedlings, the enhancing effect of nHA application on micronutrient uptake into the shoots/roots was evident at 100% precipitation intensity on select elements. For example, there was 69.7% more Zn and 60.1% more Mn accumulated in nHA-treated soybean shoots, and 76.7% more Cu accumulated in nHA-treated soybean roots, as compared to the control ( Figure 6B). For the pi-treated seedlings, the enhancing effect of pi application on micronutrient uptake into the shoots/roots was only evident at 0% precipitation intensity on selected elements. For example, there was 94.7% more Zn, 25.9% more Fe, and 51.1% more Mn that accumulated in pi-treated soybean shoots and 83.7% more Cu in pi-treated roots, as compared to the control. With the increase of the soil pi supply and the corresponding decrease of total available pi to the plants, root micronutrient contents were notably decreased compared to the control. For example, Cu, Zn, Fe, Mn, and B content at 100% precipitation intensity were reduced by 45.3%, 84.4%, 14.7%, 23.9%, and 25.1%, respectively, as compared to the ones at 0% precipitation intensity ( Figure 6B). Moreover, Zn, Fe, Mn, and B contents in pi-treated soybean roots were unexpectedly less than the control. The mechanism by which the soil amendment of pi impeded the root accumulation of these micronutrients is unknown. In this case, further investigations optimizing the fertilizer application route, pattern, and techniques are needed.
Similar to the above-mentioned macronutrients that were not externally applied, increased micronutrient contents in the soybean tissues were all due to soil uptake ( Figure S2). For example, at 100% precipitation intensity, when nHA/pi was supplied through soil amendment, Cu, Zn, and Fe content in soil were depleted by 18.6%, 18.2%, and 12.5%, respectively, with nHA treatment as compared to the control; while there was no significant alteration with pi treatment. Conversely at 0% precipitation intensity, when nHA/pi was supplied through foliar spray, Cu, Zn, and Fe content in soil were depleted by 15.7%, 16.9%, and 11.1%, respectively, with pi treatment as compared to the control; while there was no significant alteration with nHA treatment. nHA is more suitable for soil amendment to enhance soybean shoot nutrient acquisition, while pi is better for foliar spray. At 0% precipitation intensity, pi-treated seedlings had 15.5% more Cu, 84.0% more Zn, 19.2% more Fe, 52.3% more Mn, and 10.1% more B in shoots than nHA-treated plants; while at 60% precipitation intensity, nHA-treated soybean seedlings had 91.3% more Cu, 137.5% more Zn, 25.5% more Fe, 63.4% more Mn, and 21.2% more B in shoots than pi-treated ones ( Figure 6).
In the wash-off test, the enhancement effect of the nHA/pi application on micronutrient uptake into the shoots/roots was gradually lost with increasing of the wash-off height ( Figure S5), with the exception being Cu, which showed increased accumulation in soybean roots upon the wash-off.

Nutrient Contents in Edible Tissues
The content of both micro-and macronutrients in the immature soybean pods were increased upon exposure to both P fertilizers across all treatments (Figures 4C,F, 5C,F,I  and 6C,F,I,L,O), although the effects were not dose-dependent. The highest P and Ca accumulation occurred in the nHA treatment with 0% precipitation intensity and the pi treatment with 30% precipitation intensity. nHA-treated soybean pods had 250.6% more P and 140.1% more Ca than the control ( Figure 4C, at 0% precipitation intensity), and pitreated soybean pods had 518.7% more P and 167.7% more Ca than the control ( Figure 4F, at 30% precipitation intensity). Upon the application of nHA/pi, the accumulation of other macro/micronutrients in the soybean pods was also stimulated. For example, upon the nHA treatment at 100% precipitation intensity, as compared to the controls, soybean pods had 342.3% more K, 100.7% more Mg, 85.1% more S, 85.3% more Cu, 1387.0% more Zn, 343.2% more Mn, and 87.5% more B; upon the pi treatment at 30% precipitation intensity, soybean pods had 130.8% more K, 51.5% more S, 47.0% more Cu, 667.1% more Zn, and 364.6% more Mn. The enhancing effect of nHA on the nutrients acquisition in the soybean pods was significantly greater than pi across nearly all treatments. The nutrient composition and content in the early stage of pods is highly dynamic, which likely explains the high variability. In addition, it is noteworthy that the reproductive growth period of a crop can be altered upon environmental condition changes, such as P deficiency and warmer temperatures, among others [86,87]. The slow-release property of the externally applied P-nanofertilizer may potentially impact the time required for pod development [88]; this effect needs additional investigation. Regarding the effects of wash-off on the pods, both the macro-and micronutrients content in pods were notably decreased across most of the wash off-treatments, likely due to the loss of nutrients from the leaves (Figures S4 and S5). It is noteworthy that after wash-off, P content in pi-treated soybean pods at wash-off heights of 20 and 120 cm were still equivalent to the plants without wash-off; this could be ascribed to the readily uptake of pi by the above-ground portion of the soybean.

Conclusions
This study evaluated the effects of nanoscale and conventional P fertilizers on the development and growth of soybean seedlings. A precipitation intensity variable was incorporated into the experiments, and the costs and benefits of both P fertilizers under foliar and soil application scenarios were evaluated. For foliar spray of the P fertilizers, maintaining high retention on the leaf surface to facilitate direct uptake is important, while for soil amendment, maintaining high bioavailability is the most important. Our findings suggest that a P-nanofertilizer (nHA) could be a better option as a soil amendment compared to a conventional phosphate fertilizer, while pi may be the better option for foliar fertilization. However, encapsulating the nHA into a carrier material, such as a biopolymer, could significantly promote efficacy of the foliar application approach. Importantly, the P-nanofertilizer was superior in stimulating P nutrient accumulation in soybean pods during early development. The findings of this study add to the growing body of literature demonstrating that nanoscale fertilizer strategies may offer increased efficiency of delivery and utilization, while simultaneously alleviating some of the negative impacts associated with agriculture.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/agronomy11061086/s1, Physiological effects of simulated precipitation on the growth of treated ( Figure S1); Soil nutrient contents under the simulated precipitation scenario ( Figure S2); Physiological effects of wash-off on the growth of treated soybeans ( Figure S3); Macronutrient contents under the wash-off test ( Figure S4); Micronutrient contents under the wash-off test ( Figure S5); Soil nutrient contents under wash-off test ( Figure S6).

Data Availability Statement:
The data presented in this study are available upon reasonable request from the corresponding author.