Soil was collected from agriculturally managed areas of southwestern Missouri and western and eastern Arkansas for use in a laboratory rainfall-runoff simulation study. Bulk soil was collected from the 0- to 15-cm depth at all locations. A Roxana fine sandy loam (coarse-silty, mixed, superactive, nonacid, thermic Typic Udifluvents) [31] was collected from a row-crop-cultivated field at the University of Arkansas, Division of Agriculture’s Vegetable Research Station near Kibler, AR (35˚22'50"N; 94˚14'01"W) [32] . A Calloway silt loam (fine-silty, mixed, active, thermic Aquic Fraglossudalfs) [33] was collected from within an approximate 3-m² area from the edge of a row-crop-cultivated field at the University of Arkansas, Division of Agriculture’s Pine Tree Research Station near Colt, AR (35˚07’23”N; 90˚55'46''W). The Roxana soil was collected in December 2017 and the Calloway soil was collected in early spring 2021. Two soils were also collected from managed pasturelands at the University of Missouri’s Southwest Research Center near Mount Vernon, MO in June 2021. One soil was a Creldon silt loam (fine, mixed, active, mesic Oxyaquic Fragiudalfs; 37˚04'45''N; 93˚52'17''W) [34] and the second soil was a Dapue silt loam (fine-silty, mixed, active, mesic Fluventic Hapludolls; 37˚05'07''N; 93˚52'17''W) [35] . All four soils represented various typical agricultural areas in the mid-Southern US.

After collection, all moist soil was gently, manually pushed through a 6-mm sieve to remove large plant or/and rock debris and to unify the aggregate-size distribution to facilitate use in the rainfall-runoff simulation tests. Sieved soil was then air-dried in a greenhouse for approximately one week at approximately 37.8˚C and stored in 20-L plastic buckets. Three soil sub-samples were gathered from each of the four air-dried bulk soils. Soil sub-samples were oven-dried at 70˚C for 48 hours to determine gravimetric water content, manually crushed using a mortar and pestle, and sieved through a 2-mm sieve before performing chemical analyses.

Soil pH and electrical conductivity (EC) were measured potentiometrically in a 2:1 water volume: soil mass slurry. Soil organic matter (SOM) concentration was determined gravimetrically through weight-loss-on ignition after 2 hours of combustion at 360˚C in a muffle furnace [36] . Total N (TN) and C (TC) were determined by high-temperature combustion with a VarioMAX CN analyzer (Elementar Americas Inc., Mt. Laurel, NJ). Particle-size analyses were conducted using a modified 12-hr hydrometer method [37] to determine sand, silt, and clay fractions and confirm soil textural classes.

Extractable soil nutrient concentrations were also determined for each of the four soils. A water extraction was performed with a 1:10 soil mass: water volume ratio, where the soil suspensions were agitated for 1 hour, filtered through a 0.45-μm filter, and analyzed by inductively coupled, argon-plasma spectrometry (ICAPS; Spectro Arcos ICP, Spectro Analytical Instruments, Inc., Kleve, Germany) [38] to determine WS element (i.e., P, K, Ca, Mg, S, Na, Fe, Mn, Zn, and Cu) concentrations. A M3 extraction [39] was conducted with a 1:10 (mass: volume) soil: extractant solution ratio to determine weak-acid-extractable nutrient (i.e., P, K, Ca, Mg, S, Na, Fe, Mn, Zn, and Cu) concentrations. A strong-acid digestion was conducted using the Environmental Protection Agency (EPA) method 3050B [40] and analyzed by ICAPS to determine total-recoverable element (i.e., P, K, Ca, Mg, S, Na, Fe, Mn, Zn, and Cu) concentrations. All measured initial soil physical and chemical properties are summarized in Table 1.

The fertilizer-P treatments used in the rainfall-runoff simulation were CPST (i.e., Crystal Green) [17] , MAP, ECSTreal, ECSTsyn, and an unamended control (UC). The CPST material is pelletized and has nutrient concentrations of 11.7% P, 6% N, and 8.3% Mg [32] . Monoammonium phosphate is a pelletized, commonly used, commercially available P and N fertilizer that has a reported fertilizer grade of 11-52-0 and contains 20.9% P, 11% N, and 1.5% Mg [27] . Both ECSTsyn and ECSTreal were electrochemically precipitated struvite materials, but ECSTsyn was precipitated from a synthetic wastewater source, produced to have a similar average P and N composition as typical municipal wastewater [19] . The nutrient composition of ECSTsyn was 18.5% P, 3.3% N, and 13.3% Mg [13] . Using similar methods as for the creation of the ECSTsyn material, the ECSTreal material was created from an actual municipal wastewater source collected from the West Side Wastewater Treatment Facility in Fayetteville, AR and contained 15.4% P, 3.3% N, and 13.6% Mg.

Chemical analyses had been previously performed on the ECSTsyn, CPST, and MAP fertilizer-P sources by Anderson et al. [13] , and similar procedures were used to chemically characterize the ECSTreal fertilizer-P source. Since both the ECSTsyn and ECSTreal materials were electrochemically precipitated, the material produced was a crystalline-flake powder. Since CPST and MAP were originally in pellet form, both were finely ground to match the powder consistency of the ECSTsyn and ECSTreal materials for more valid fertilizer comparisons.

Three water sources were used for the rainfall-runoff simulations. Rainfall was collected in July 2021 from a single rainfall event. The rainwater used in the rainfall-runoff simulation was to mimic natural rainfall events on the soil. Groundwater was obtained from an existing well in July 2021 from the Savoy Research Center west of Fayetteville, AR (36˚07'52''N; 94˚18'54''W) to represent an irrigation-water source. The third water source was struvite-removed wastewater (SRW) produced in July 2021 as a result of the precipitation of the ECSTreal material described above. The purpose of using the SRW was to evaluate its effectiveness as a potential irrigation-water source after struvite removal. All water sources were refrigerated at 4˚C until used. Chemical characterization of the three water sources occurred from the analyses of blanks as part of the rainfall-runoff simulations described below.

^† Means in a row with different letters are different at P < 0.05.

Due to the experimental nature of ECSTsyn and ECSTreal, only small quantities of both ECST materials were able to be produced for testing. Consequently, despite abundant runoff-water-quality studies using simulated rainfall onto field plots > 1 m², to accomplish the objective of the current study required a much smaller scale in terms of soil and fertilizer mass used, treatment area exposed to stimulated rainfall, and the rainfall-simulation apparatus itself.

To simulate rainfall onto fertilizer-soil mixtures and generate runoff, a wooden rainfall-runoff simulator was constructed measuring 91.1-cm tall by 87.0-cm wide by 73.7-cm deep. A frame containing four sets of seven, 31.5-mL∙min⁻¹ drip emitters (MLD-AAD, Mister Landscaper Inc., Dundee, FL), separated by approximately 15.9 cm, rested parallel to one another on top of the simulator. For each of the four rows of emitters, three drip emitters were installed in a row 10 cm apart starting at 10.8 cm from the water-source end of the series of emitters and four additional drip emitters were installed 5 cm apart 10.5 cm from the third emitter. All emitters were connected with 4.8-mm-inside-diameter plastic tubing. Tapered, plastic trays (i.e., short sections of rain gutters for houses), 62-cm long by 12-cm wide at the top by 7-cm wide at the base, were manually placed free-standing at a slope of 22.1% below each row of drip emitters to contain the soil-fertilizer mixtures and collect runoff.

On a day of a set of rainfall-runoff simulations, water stored in a 20-L plastic carboy was removed from the refrigerator approximately 1 hr prior to use. The carboy was set on top of a ladder 176.5 cm above and immediately adjacent to the rainfall simulator. The carboy’s spigot was connected to 6.4-mm-inside-diameter plastic tubing to gravity-flow to and through the drip emitters. The water delivery rate was approximately 3.5 cm∙hr⁻¹ to represent a substantial, high-intensity storm event.

Tea bags (TamBee, B07TCDT76Q), 15-cm wide by 20-cm long, made from a synthetic fiber, were filled with 175 g of air-dried soil that had previously been sieved through the 6-mm mesh screen. An agriculturally relevant fertilizer-P rate (56 kg P₂O₅ ha⁻¹) was determined based on the fertilizer-P recommendation for soybeans (Glycine max L. [Merr.]) grown on a loamy Arkansas soil [41] . The appropriate mass of each fertilizer material (i.e., CPST, ECSTsyn, ECSTreal, and MAP) was added to the soil-filled tea bags to deliver the fertilizer-P rate of 56 kg P₂O₅ ha⁻¹ based on the air-dry mass of soil added to the tea bag, but the nitrogen rate was not balanced among fertilizer treatments since the objective of the study did not include evaluating any plant response. The tea bags were then gently, manually massaged to mix the fertilizer with the soil. Since air-dried soil was initially used, once fertilizers were mixed with the soil, each tea bag was pre-wet using a squirt bottle with 20 mL of the appropriate water source (i.e., rainwater, groundwater, or SRW) one week before conducting the rainfall simulation to increase the soil moisture content and to minimize potential soil and fertilizer loss through the tea bags upon handling. Three replications of each soil-fertilizer combination were prepared for a total of 60 tea bags for each of the three water sources, totaling 180 tea bags.

Rainfall simulations were performed on 21, 22, and 23 July 2021 (Trial 1). Each day, a single water source was used for rainfall simulations. The sequence of tea bags containing each soil-fertilizer combination was randomized within each water source. Prior to any rainfall simulations, the plastic tubing was primed with the water source to be used for raining on treatment-filled tea bags. Prior to raining on each set of four treatment-filled tea bags, rainfall was simulated onto the four empty trays and then rainfall was simulated into the four trays with an empty, dry tea bag in the tray to serve as two sets of blanks. Treatment-filled tea bags were placed four at a time directly under the four drip emitters that were closest together and near the down-slope end of the tray to simulate raindrops impacting the soil and water moving underneath and through the soil from up-slope. Each set of four rainfall simulations were conducted for 6 minutes. Pre-trial tests were conducted to determine that 6 minutes of rainfall would yield an estimated 220 mL of runoff water for conducting all planned chemical analyses.

Runoff water from each tray was collected in 2400-mL, rectangular containers at the end of each tray and poured into 250-mL plastic cups for immediate processing for water quality attributes. The process was repeated until each sample was rained on and all runoff water samples were collected.

After raining on each set of four tea bags, tea bags were laid out on plastic bags to air-dry at room temperature (~22˚C) for one week. Each day the tea bags were gently, manually turned over to facilitate uniform drying on each side of the tea bag. Once air-dried, the three replications of each treatment-filled tea bag were placed into its own separate plastic bag for incubation for approximately six months when a second rainfall-runoff simulation was performed. The incubated samples for the rainwater and groundwater water sources were wetted using a beaker with 5 to 10 mL of each appropriate water source once a month.

Following the procedures described above, the second rainfall-runoff simulation experiment was conducted on 15 and 16 December 2021 (Trial 2). Only the rainwater and groundwater sources were used in the second rainfall-runoff simulation experiment due to insufficient quantity of the SRW source.

Runoff-sample pH and EC measurement, filtering, and acidification were conducted immediately after sample collection. For pH and EC measurements, approximately 25 mL of each sample were poured into a 50-mL glass beaker. A combination pH/EC probe (Orion Star^TM A215 pH/Conductivity Benchtop Multiparameter Meter, Thermo Scientific, Beverly, MA) was placed in the runoff sample for 30 seconds, after which pH and EC readings were recorded. The pH/EC probe was rinsed with distilled water between each individual runoff sample. The pH and EC for all runoff samples were measured within seven minutes of sample collection.

Following pH/EC measurements, approximately 100 mL of each runoff water sample were poured into a small plastic vessel (UC475, Ultimate Clean Environmental Express, Charleston, SC) and filtered through a 0.45-µm filter (SC0409, FilterMate, Environmental Express, Charleston, SC) attached to the bottom of a plunger. Approximately 20 mL of filtered runoff sample were poured into four, 20-mL, plastic scintillation vials. Two of the four vials were acidified with two drops of concentrated hydrochloric acid for preservation and analysis for total soluble elemental concentrations [i.e., total P (TP), Ca, Mg, and Fe], while the other two vials were left unacidified for nitrate + nitrite (NO₃ + NO₂), NH 4 + , and soluble-reactive P (SRP; primarily orthophosphates) concentration analyses. Nitrate + nitrite was measured using the salicylate method (EPA 351.2) [42] , NH 4 + was measured using the cadmium-reduction method (EPA 353.2) [43] , and SRP was measured using the acid-molybdate method (EPA 365.1) [44] . Blanks were also analyzed for total soluble S, Na, Mn, Zn, and Cu in acidified runoff samples. When a measured concentration was below the instrument’s detection limit (DL), the value used was one-half of the instrument’s DL.

To prepare the resulting runoff-water data for statistical analyses, separately for each water source, the average of all eight measured blank values for a water quality parameter was subtracted from the raw water quality parameter measurement to create a blank-corrected value. Since each of the four soils used had somewhat differing initial soil properties, to eliminate potential data basis from dissimilar initial properties, the average was calculated for the UC replications for each water-source-soil-fertilizer-treatment combination and was also subtracted from the blank-corrected value. Correcting the raw water quality parameter measurements for blanks and subtracting the average UC value resulted in a final water quality parameter data set. The final data set varied from negative to positive values for all measured water quality parameters, which is hereafter referred to as the change in water quality data, to better address actual soil-fertilizer interactions without being skewed by dissimilar initial soil properties.

Based on a completely random design, a one-factor analysis of variance (ANOVA) was conducted in SAS (version 9.4, SAS Institute, Inc., Cary, NC) using PROC GLIMMIX to evaluate the differences in initial soil properties among the four soils used and the differences in initial water properties among the three water sources used. A gamma distribution was used for analyses of initial soil and water data.

It was impractical to randomize the water sources for the rainfall-runoff simulations, but the soil-fertilizer treatment combinations were randomized within a water source. Therefore, a randomized block design was assumed for the purposes of statistical analyses of runoff-water-quality data. Consequently, a three-factor ANOVA was performed using PROC GLIMMIX in SAS to evaluate the effects of water source, soil, fertilizer-P source, and their interactions on the change in runoff water pH, EC, total soluble P, K, Ca, Mg, Fe, NO₃ + NO₂, NH₄, and SRP separately for the first rainfall-runoff trial. Since the change in water quality parameters varied between negative and positive values, a normal distribution was used for statistical analyses of water-quality-change data.

To evaluate treatment effects for Trial 2 and changes over time (i.e., results from the first compared to the second rainfall-runoff trial) on the soil-fertilizer treatment combinations for the rainwater and groundwater sources only, a four-factor ANOVA was performed using PROC GLIMMIX in SAS evaluate the effects of water source, soil, fertilizer-P source, time (i.e., trial), and their interactions on the change in runoff water pH, EC, total soluble P, K, Ca, Mg, Fe, NO₃ + NO₂, NH₄, and SRP. Since the change in water-quality parameters over time varied between negative and positive values, a normal distribution was also used for analyses of water-quality data. For all data analyses, significance was judged at P < 0.05. When appropriate, means were separated by the least significant difference at the 0.05 level.

All initial soil properties evaluated differed (P < 0.01) among soils (Table 1). Percent clay was largest in the Calloway soil, smallest in the Dapue soil, and was intermediate for the Roxana and Creldon soils, which did not differ (Table 1). Similar to percent clay, soil pH was most alkaline in the Calloway soil, most acidic in the Dapue soil, and was intermediate for the Roxana and Creldon soils (Table 1). Total C, TN, SOM, and NO₃-N were largest in the Dapue soil and smallest in the Roxana soil (Table 1).

Water-soluble P, Ca, Mg, and Fe concentrations differed (P < 0.01) among all soils (Table 1). Water-soluble P, Mg, and Fe were numerically largest in the Roxana soil and numerically smallest in the Calloway soil (Table 1). However, unlike P and Fe, Mg did not differ between the Calloway and Creldon soils. Water-soluble Ca was largest in the Dapue soil and smallest in the Roxana soil (Table 1).

Mehlich-3-extracable P, Mg, Ca, and Fe concentrations differed (P < 0.01) among all soils (Table 1). Mehlich-3 Ca and Fe were largest in the Calloway soil, but M3 Ca was smallest in the Roxana soil, while M3 Fe was smallest in the Creldon soil (Table 1). Mehlich-3 P was largest in the Roxana soil and smallest in the Calloway soil, while M3 Mg was largest in the Creldon soil and smallest in the Roxana soil (Table 1).

Total-recoverable P, Ca, Mg, and Fe concentrations differed (P < 0.01) among all soils (Table 1). Total-recoverable Fe and Ca were largest in the Calloway soil, but TR Fe was smallest in only Roxana, whereas TR Ca was smallest in both Roxana and Creldon, which did not differ (Table 1). Total-recoverable P was largest in the Dapue soil and smallest in the Creldon soil, while TR Mg was largest in the Roxana soil and smallest in the Dapue soil (Table 1). Differences in initial soil pH and nutrient concentrations necessitated inclusion of the UC treatment and subtraction of the mean water-quality response for the UC treatment from the water-quality response from each replicate soil-fertilizer treatment combination.

All initial water properties evaluated (i.e., pH, EC, and TP, K, Ca, Mg, S, Na, Fe, Mn, Zn, Cu, B, NO₃ + NO₂, NH₄, and SRP concentrations), except Fe, differed (P < 0.01) among water types (Table 2). Initial water pH was most alkaline for the SRW and least alkaline for rainwater (Table 2). Of the initial nutrient concentrations evaluated, most were largest in the SRW, except for Ca and NO₃ + NO₂, which were largest in groundwater (Table 2). Initial Ca, Mg, and NO₃ + NO₂ were numerically smallest in rainwater, while TP, NH₄, and SRP were numerically smallest in groundwater (Table 2). The alkalinity differences among initial water types likely affected elemental and compound solubilities, particularly after 6 months of incubation for Trial 2.

With the exception of runoff EC, every other measured water quality parameter (i.e., pH, TP, Ca, Mg, Fe, NO₃ + NO₂, NH₄, and SRP) change was affected (P < 0.05) by a combination of two or more treatments (i.e., water type, soil, and/or fertilizer-P source; Table 3). Runoff water quality changes are being presented where responses from blanks without soil and the UC were both subtracted so that the change results isolate the actual effect of the fertilizer addition rather than being complicated by the inherent differences among initial soil (Table 1) and water (Table 2) properties.

^† Means in a row with different letters are different at P < 0.05.

^† Electrical conductivity, EC; total phosphorus, TP; calcium, Ca; magnesium, Mg; iron, Fe; nitrate + nitrite, NO₃ + NO₂; ammonium, NH₄; soluble-reactive phosphorus, SRP.

Averaged across fertilizer-P sources, runoff-water pH, NH₄, Ca, Mg, and Fe concentration changes differed (P < 0.05) among water type-soil combinations (Table 3). Water pH influences the fraction of soluble nutrients compared to what might remain sorbed to soil solids. Of the three water types, runoff-water pH changes in response to rainwater differed from zero in all soils, except for Creldon where runoff pH decreased in the Dapue and Calloway soils, but increased in the Roxana soil (Figure 1). In addition, runoff-pH changes decreased from zero in the Dapue soil-wastewater treatment combination (Figure 1). Runoff-pH changes from all other water type-soil combinations did not differ from zero or from each other (Figure 1). The runoff-pH decrease among the Dapue- and Calloway-rainwater and Dapue-wastewater combinations, which did not differ, was lower than the runoff pH increase from rainwater in the Roxana soil (Figure 1). It is likely that the runoff pH response in the Dapue soil was dominated by the lowest initial soil pH and largest initial soil TC and SOM concentrations that imparted an acidifying effect to the runoff water. The runoff pH increase from rainwater in the Roxana soil was likely related to low buffer capacity from the coarsest soil texture so that the alkaline rainwater dominated the response. For the Calloway soil, the decreased runoff pH from rainwater was likely due to the lower initial rainwater pH compared to the initial soil pH, whereas the initial pH for the groundwater and wastewater were both greater than the initial soil pH.

Excluding the groundwater-Calloway and rainwater-Roxana combinations, runoff-water NH₄ concentration changes for all other water type-soil combinations differed (P < 0.05) from zero (Figure 1). The greatest runoff-water NH₄ concentration decrease occurred from the wastewater treatment within each soil (Figure 1), in which all runoff-water NH₄ concentration changes for the wastewater-soil combinations also differed from each other (Figure 1). The wastewater-soil combinations had the greatest decrease in NH₄. The large initial NH₄ concentration in the wastewater was likely retained in the soil during rainfall and was responsible for the runoff NH₄ decrease.

In addition, averaged across water types and soils, runoff-water NH₄ concentration changes differed (P = 0.02) among fertilizer treatments (Table 3). Runoff-water NH₄ concentration changes for all four fertilizer treatments (i.e., CPST, MAP, ECSTreal, and ECSTsyn) differed from zero (i.e., −1.35, −1.14, −1.28, and −1.25 mg∙L⁻¹, respectively). The MAP treatment had the numerically smallest runoff-water NH₄ concentration change, which differed from ECSTreal and CPST, while CPST had the numerically largest runoff-water NH₄ concentration change (Table 3). Runoff-water NH₄ concentration differences from ECSTsyn did not differ from any other fertilizer treatment (Table 3). Added NH₄ from the fertilizers was retained by the soils, likely due to cation exchange sites from clays, TC, and/or SOM.

In contrast to NH₄, runoff-water Ca concentration changes differed from zero in three of four groundwater-soil combinations (i.e., Creldon, Dapue, and Roxana), but did not differ from zero for any of the rainwater- or wastewater-soil combinations (Figure 1). However, runoff-water Ca concentration changes did not differ from each other in the groundwater-Calloway combination or among any other rainwater-soil combination (Figure 1). In addition, in contrast to NH₄, runoff-water Ca concentrations decreased the most in the groundwater-soil combinations, which was a 5 to 30 times greater decrease in runoff-water Ca concentration compared to all other water type-soil combinations (Figure 1). Groundwater had a large initial pH (7.9) in addition to having the largest initial Ca concentration (Table 2). The combination of large initial pH and large initial Ca concentration likely facilitated the precipitation of a portion of the Ca as Ca-P compounds to retain Ca and result in lower runoff Ca for the Creldon, Dapue, and Roxana soils. However, the Calloway soil was unaffected by all three water types, possibly due to its large initial pH leading to no additional precipitation of Ca.

Similar to NH₄, runoff-water Mg concentration changes differed from zero in all wastewater-soil combinations and in the groundwater-Roxana combination (Figure 2). Runoff-water Mg concentration changes in all rainwater- and groundwater-soil combinations did not differ from each other. Runoff-water Mg concentration changes in the wastewater treatment across the four soils decreased 5 to 25 times more than any other water type-soil combination (Figure 2). Wastewater had the largest initial pH of the three water types (9.8; Table 2),

which likely caused a fraction of the Mg to precipitate as Mg-P and/or Ca-Mg-P compounds in all four soils, resulting in the runoff-Mg decrease from zero. The runoff-Mg concentration change was also positive in the groundwater-Roxana combination, but the positive change was small.

Runoff-water Fe concentration changes differed from zero in three of four rainwater-soil combinations, excluding the Dapue soil, and was greater than zero for the groundwater-Calloway combination (Figure 2). Runoff-water Fe concentration from all water type-soil combinations, except for those that differed from zero and the wastewater-Dapue combination, did not differ from each other (Figure 2). Of the water type-soil combination whose runoff-water Fe concentration change differed from zero, the rainwater-Calloway and -Roxana combinations did not differ from each other and both changes decreased from zero (Figure 2). Similarly, runoff-water Fe concentration change for the rainwater-Creldon and groundwater-Dapue combinations did not differ from each other, but both changes were greater than zero (Figure 2). The decrease in runoff Fe in the Calloway- and Roxana-rainwater combinations was likely due to the relatively large initial soil and water pHs, as Fe solubility in soil decreases as pH increases. For similar reasons, in combination with the low initial Fe in the Calloway and Creldon soils (Table 1), there was an increase in runoff Fe in the Creldon-rainwater and Calloway-groundwater combinations.

Averaged across soils, runoff-water Fe concentration changes also differed (P = 0.04) among water type-fertilizer-P source combinations (Figure 3). Runoff-water Fe concentration changes differed from zero in three of four rainwater-fertilizer-P source combinations, excluding CPST (Figure 3). Runoff-water Fe concentration changes were less than zero in all wastewater-fertilizer-P source combinations, but were greater than zero in all groundwater-fertilizer-P source combinations (Figure 3). In general, the rainwater treatment caused the greatest runoff-water Fe concentration change from zero among all water sources (Figure 3). In all rainwater-fertilizer-P source combinations, except for ECSTreal, the runoff-water Fe concentration change was less than zero (Figure 3). Similarly, the runoff-water Fe concentration change from rainfall for CPST and ECSTsyn was similar to MAP, except for the rainwater-ECSTreal combination (Figure 3). All treatment combinations did not differ from each other, except for the rainwater-CPST, -MAP, and -ECSTsyn and wastewater-ECSTsyn combinations, which did not differ from each other (Figure 3). The increase in runoff Fe from only ECSTreal in the rainwater treatment was likely due to the original conditions under which the ECSTreal was formed. Unlike ECSTsyn, ECSTreal was precipitated from a real municipal wastewater source, while CPST and MAP were also created from a source with known composition. There could

have been unknown or unexpected elemental concentrations within the wastewater from which ECSTreal was originally created. Because Fe is more available at more acidic pHs, the non-significant runoff Fe change that occurred in the groundwater and wastewater treatments was likely due to the groundwater’s and wastewater’s initial pHs > 7.8 (Table 2).

Runoff-water TP, NO₃ + NO₂, and SRP concentration changes differed (P < 0.05) among water type-soil-fertilizer-P source treatment combinations (Table 4). There were few consistent trends among water type-soil-fertilizer-P source treatment combinations, and most were complex. Most runoff-water TP concentration changes differed from zero among treatment combinations, particularly within the wastewater treatment (Table 4). All runoff-water TP concentration changes from wastewater were negative in three of four soil-fertilizer-P source combinations, excluding the Roxana soil (Table 4). The numerically largest runoff-water TP concentration changes were in the rainwater-Roxana combination, excluding the CPST fertilizer-P source. The numerically smallest runoff-water TP concentration changes occurred in the wastewater-Calloway combinations (Table 4). Few runoff-water TP concentration changes differed from each other among fertilizer-P sources within water type-soil combinations (Table 4). In most treatment combinations, runoff-water TP concentration changes from the struvite-P sources (i.e., CPST, ECSTsyn, and ECSTreal) were similar to MAP (Table 4). However, within the rainwater-Creldon and groundwater-Dapue combinations, runoff-water TP concentration changes from ECSTsyn differed from all other fertilizer treatments. Within the rainwater-Roxana combination, runoff-water TP concentration changes did not differ among MAP, ECSTreal, or ECSTsyn (Table 4).

The negative runoff-TP concentration changes in the wastewater-soil-fertilizer treatment combinations could have been caused by the large initial wastewater Ca and Mg concentrations binding with P due to the wastewater’s alkaline pH (Table 2). Phosphorus solubility in the soil is largely dictated by pH, where phosphate forms poorly soluble Fe³⁺ and Al³⁺ compounds at low pH, readily soluble Ca²⁺ and Mg²⁺ compounds at pH near neutrality, and poorly soluble Ca²⁺ compounds at greater pH [45] . However, in the wastewater-Roxana combination, the TP concentration change was positive. The positive TP concentration change, in contrast to the other three soils’ negative TP concentration change, was likely due to the large initial TP, but small initial Ca (Table 1), causing little to no precipitation of Ca-P compounds. In every water-soil-fertilizer combination but one (i.e., wastewater-Roxana), either ECSTsyn or ECSTreal consistently had the largest runoff-TP concentration (Table 4). The large runoff-TP concentration of ECSTreal and ECSTsyn was unexpected since the initial P concentration of the ECST fertilizers was the intermediate among the four fertilizer treatments, but the initial Mg concentration for ECSTsyn and ECSTreal was the largest (i.e., 13.3% and 13.6%, respectively). The large runoff-TP concentration in the ECSTsyn and ECSTreal treatments, in contrast to CPST and MAP, may

^† Means in a column with different letters are different at P < 0.05. * An asterisk (*) indicates the mean change differs from zero at P < 0.05. ^‡ Chemically precipitated struvite, CPST; monoammonium phosphate, MAP; real-wastewater-derived electro chemically precipitated struvite, ECSTreal; synthetic electrochemically precipitated struvite, ECSTsyn; total phosphorus, TP; nitrate + nitrite, NO₃ + NO₂; soluble-reactive phosphorus, SRP.

have been due to the timing of when each fertilizer material solubilized and subsequently ran off during the rainfall event. Being readily soluble, MAP may have dissolved within the first few days of pre-wetting and had sufficient time for Ca-Mg-P compounds to form, lowering the runoff-TP concentration. In contrast, ECST’s slow-release nature, as reported by Anderson et al. (2021c), could have led to the simulated rainfall promoting ECST solubilization, increasing runoff-TP concentrations.

Runoff-water SRP concentration changes were similar to runoff-water TP concentration changes in terms of which treatments differed from zero and which treatments differed from each other (Table 4). Similar to runoff-water TP, runoff-water SRP concentration changes differed from zero in almost all water type-soil-fertilizer-P source combinations, except for the groundwater-Calloway combination (Table 4). Also similar to runoff-water TP, the numerically largest runoff-water SRP concentration changes occurred for the rainwater-Roxana treatment combination, excluding CPST (Table 4). The numerically smallest runoff-water SRP concentration changes occurred for the wastewater-Calloway combinations (Table 4). Few runoff-water SRP concentration changes differed from each other among fertilizer-P sources within water type-soil combinations (Table 4). In general, except for the rainwater-Calloway and -Creldon and wastewater-Dapue combinations, runoff-water SRP concentration changes were numerically smaller than runoff-water TP concentration changes. In most treatment combinations, runoff-water SRP concentration changes from the struvite-P sources were similar to MAP (Table 4).

Runoff-SRP concentration responses were likely the result of reasons similar to that for runoff-TP responses. The negative runoff-SRP concentration changes in the wastewater-soil-fertilizer combinations could have been caused by the large initial wastewater Ca and Mg concentrations binding with P due to the wastewater’s alkaline pH. However, in the wastewater-Roxana combination, the SRP concentration change was generally positive. The positive TP concentration change, in contrast to the other three soils’ negative SRP concentration change, was likely due to the large initial TP, but small initial Ca causing no precipitation of Ca-P compounds. Runoff-SRP concentration changes were smaller than for runoff TP in almost all instances, which was consistent with TP measuring all P forms and SRP only measuring plant-available P forms.

In contrast to runoff-water TP and SRP, most runoff-water NO₃ + NO₂ concentration changes did not differ from zero (Table 4). Runoff-water NO₃ + NO₂ concentration changes differed from zero in three of four rainwater-Callo-way-fertilizer-P source combinations, excluding ECSTreal (Table 4). Runoff-water NO₃ + NO₂ concentration changes for most treatment combinations also did not differ from each other, except the numerically largest and smallest runoff-water NO₃ + NO₂ concentration changes (Table 4). The largest runoff-water NO₃ + NO₂ concentration changes occurred for the wastewater-Calloway-CPST and the rainwater-Calloway-CPST and -MAP combinations (Table 4). The smallest runoff-water NO₃ + NO₂ concentration change occurred for the groundwater-Dapue-ECSTreal combination (Table 4). In most treatment combinations, runoff-water NO₃ + NO₂ concentration changes for the struvite-P sources were similar to MAP (Table 4). Complex interactions between soil, water type, and fertilizer properties made it challenging to explain runoff NO₃ + NO₂ concentration changes, as there were no consistent patterns among treatment combinations. In contrast to all other measured water-quality parameters, runoff-EC changes were unaffected by water type, soil, or fertilizer-P source (Table 3) and averaged 4.32 dS∙m⁻¹ across all treatments.

After a 6-month period of incubation between rainfall simulations, with monthly wetting and in contrast to Trial 1, every measured water quality parameter, with the exception of runoff-Mg concentration changes, differed (P < 0.05) over time (Table 5). For Trial 2, SRW was not used, thus the water types only consisted of rainwater and groundwater. The remainder of the results in this section will focus on changes over time and differences among treatment combinations from Trial 2 results only.

^†Electrical conductivity, EC; total phosphorus, TP; calcium, Ca; magnesium, Mg; iron, Fe; nitrate + nitrite, NO₃ + NO₂; ammonium, NH₄; soluble-reactive phosphorus, SRP.

Changes in runoff-pH, which partially dictates the fractionation between WS and sorbed nutrient concentrations, and runoff-TP, -Ca, and -SRP concentration changes differed (P < 0.02; Table 5) among water type-soil combinations over time (Table 5). Averaged across fertilizer-P sources, runoff-pH changes decreased over time in all water type-soil combinations, where the decrease was larger for rainwater-compared to groundwater-soil combinations (Table 6). In response to rainwater and groundwater, after 6 months of incubation (Trial 2), all runoff-pH changes differed from zero (Table 6). Runoff-pH changes in groundwater-soil combinations did not differ from each other in Trial 2 (Table 6). In contrast, runoff-pH changes in all rainwater-soil combinations differed from each other in Trial 2, except for the Roxana soil which did not differ from the Dapue or Creldon soils (Table 6). The decrease in runoff pH over time was likely due to the nitrifying processes that occurred during incubation, producing H⁺ ions, thus lowering runoff pH [46] . The general temporal decrease in runoff pH was similar to the results of Anderson et al. [27] , who reported a decrease in soil pH over time in a silty clay loam soil with large initial SOM and clay concentration.

^† Means in a column with different letters are different at P < 0.05. * An asterisk (*) indicates the mean value differs from zero at P < 0.05. ^‡ Total phosphorus, TP; calcium, Ca; soluble-reactive phosphorus, SRP.

Similar to runoff-pH, runoff-TP concentration changes also decreased over time in all water type-soil combinations, with the larger decreases over time generally occurring in the rainwater-soil combinations (Table 6). In six of the eight water type-soil combinations in Trial 2, with the exception of the groundwater- and rainwater-Roxana soil combinations, runoff-TP concentration changes were small and did not differ from zero and did not differ among themselves (Table 6). However, runoff-TP concentration changes the groundwater- and rainwater-Roxana combinations were larger than for the other six water-type combinations and differed from zero, but did not differ from one another. Soluble P may have reacted with Ca and/or Mg to form Ca- and/or Mg-P compounds, partially explaining the decrease in runoff TP over the 6 months from Trial 1 to Trial 2. However, due to the decrease in runoff pH, it was more likely that soluble P forms were adsorbed onto soil particles over the 6-month incubation period. Shigaki et al. [29] and Smith et al. [30] reported a similar decrease in WS P over time in a loam and silt loam soil, respectively.

In contrast to runoff-pH and -TP, runoff-water Ca concentrations at least numerically increased over time in all groundwater-soil and two of four rainwater-soil combinations, excluding the rainwater-Creldon, which did not change over time, and the rainwater-Calloway combination, which significantly decreased over time (Table 6). For Trial 2 only, the largest runoff-water Ca concentration increase occurred in the rainwater-Dapue, while the largest decrease occurred in the rainwater-Calloway combination (Table 6). The increase in runoff Ca concentration over time was likely due to the decrease in runoff pH and Ca being released from cation exchange sites. The runoff pH of the rainwater-Dapue combination was already low, likely resulting in no change of the cation exchange site components and explaining the decrease in runoff Ca in the rainwater-Dapue combination.

Similar to runoff-water TP, runoff-water SRP concentration changes also decreased significantly over time in all water type-soil combinations, with the larger decreases generally occurring in response to rainwater compared to groundwater combinations (Table 6). Similar to runoff-water TP, in Trial 2, runoff-water SRP concentration changes were generally small, where most did not differ from a change of zero and many did not differ from one another (Table 6). For Trial 2 only, the largest runoff-water SRP concentration increases occurred in the groundwater-Roxana and rainwater-Creldon combinations, which differed from zero, but did not differ from one another (Table 6). The explanation for changes in runoff SRP is likely similar to that of runoff TP, where it was likely that SRP adsorbed onto soil particles over time. Runoff-SRP was lower than runoff-TP concentration changes in almost all instances, which was consistent with TP measuring all P forms and SRP only measuring plant-available P forms.

Averaged across water type, runoff-water TP and SRP concentration changes also differed (P < 0.01; Table 5) among soil-fertilizer-P source combinations over time (Table 7). Treatment combination effects on runoff-TP and -SRP concentration changes over time were the same (Table 7). Runoff-TP and -SRP concentration changes were generally small for Trial 2, where, with the exception of the same two treatment combinations, most changes did not differ from zero, and decreased over time in all soil-fertilizer-P source treatment combinations (Table 7). The largest runoff-TP and -SRP concentration changes from Trial 2 occurred in the Roxana-CPST and- ECSTsyn combinations, which were both larger than zero, but did not differ from one another, and both were larger than the changes measured in all other soil-fertilizer-P source combinations (Table 7). Runoff-TP and -SRP concentration changes for Trial 2 did not differ between the three struvite-P sources in the Creldon, Dapue, or Calloway soils and behaved similarly to MAP (Table 7). However, runoff-TP and -SRP concentration changes for Trial 2 from the Roxana-CPST and -ECSTsyn were both greater than changes from Roxana-ECSTreal and -MAP combinations, which did not differ (Table 7).

The explanation for changes in runoff-SRP and -TP concentration changes was likely similar, as both SRP and TP measure soil-P pools that were likely adsorbed to soil particles over time. In most instances, runoff-TP and -SRP concentration changes were largest in the ECSTsyn or ECSTreal combinations in Trial 1. The large runoff-TP and -SRP concentration changes were unexpected

^† Means in a column with different letters are different at P < 0.05. * An asterisk (*) indicates the mean value differs from zero at P < 0.05. ^‡ Chemically precipitated struvite, CPST; monoammonium phosphate, MAP; real-wastewater-derived electrochemically precipitated struvite, ECSTreal; synthetic electrochemically precipitated struvite, ECSTsyn; total phosphorus, TP; soluble-reactive phosphorus, SRP.

since the initial P concentration of the ECST fertilizers was intermediate among the four fertilizer treatments, but the initial Mg concentration for ECSTsyn and ECSTreal was the largest (i.e., 13.3% and 13.6%, respectively). As previously, the large runoff-TP concentration change in the ECSTsyn and ECSTreal treatments, in contrast to CPST and MAP, may have been due to the timing of when each fertilizer material solubilized and subsequently ran off during the rainfall event. Being readily soluble, MAP may have dissolved within the first few days of pre-wetting and had sufficient time to form Ca-Mg-P compounds, lowering the runoff-TP concentration. The increase in runoff-TP concentration for ECST may have been due to the slow-release nature of ECST (Anderson et al., 2021c) and the promotion of ECST dissolution by rainfall. After the 6-month incubation period, there was likely sufficient time for Ca-Mg-P compounds to form in all fertilizer treatment combinations. Anderson et al. (2021c) reported that CPST, not ECST, had the greatest increase in WS-P concentration over nine months of plant-less soil incubation. However, Anderson et al. [27] also reported that, in general, there was a decrease in WS-P concentration after the initial wetting event, which was similar to the results of the current study.

In contrast to the other soluble nutrients, averaged across water type and fertilizer-P source, runoff-Fe and -NH₄ concentration changes differed (P < 0.01; Table 5) among soils over time (Figure 4). For three of the four soils (i.e., Creldon,

Dapue, and Roxana), runoff-Fe concentration changes did not differ over time, changes for Trial 2 did not differ from zero, and changes were all similar to one another (Figure 4). However, for the Calloway soil, runoff-Fe concentration changes decreased > 0.5 mg∙L⁻¹ over the six months between trials, in which the decrease differed from zero and differed from all other soil-time combinations (Figure 4). The decrease in runoff Fe over time was likely due to the large initial pH of the Calloway soil (Table 1), where Fe availability decreased at the large pH.

In contrast to Fe, averaged across water type and fertilizer-P sources, runoff-NH₄ concentration changes decreased over time for the Creldon, Dapue, and Roxana soils, where the change for Trial 2 only was greater than zero for the Dapue soil, but did not differ from the change in the Creldon or Roxana soils (Figure 4). However, for the Calloway soil, the runoff-NH₄ concentration change remained similar over time, while the change for Trial 2 only was less than zero, which also did not differ from the change in the Roxana soil, but was smaller than the change in the Creldon and Dapue soils (Figure 4). The decrease in runoff-NH₄ concentrations over time was likely due to nitrification occurring during the 6-month incubation period. There was no change over time in runoff-NH₄ concentration in the Calloway soil, but the Calloway soil had the smallest initial NH₄ concentration and largest initial pH (Table 1), which may have contributed to an increased nitrification rate and the runoff NH₄ decrease from zero. In contrast to all other measured water-quality parameters, runoff-EC and -NO₃ + NO₂ concentration changes were complex and differed (P < 0.03; Table 5) among water type-soil-fertilizer-P source combinations over time (Table 8). Runoff-EC changes for Trial 2 did not differ from zero, except for the groundwater-Calloway-MAP and -ECSTsyn, rainwater-Dapue-MAP, -ECSTreal, and -ECSTsyn combinations, where all changes were greater than zero (i.e., an increase in runoff EC; Table 8). Similarly, most runoff-EC changes did not change over time, except for the groundwater-Calloway-MAP, rainwater-Dapue-MAP, and -ECSTsyn combinations that increased over time, where the Trial 2 change did not differ among themselves, and the rainwater-Calloway-CPST and -MAP combinations that decreased over time, where the Trial 2 changes did not differ between them (Table 8). Runoff-EC changes for most other water type-soil-fertilizer-P source combinations did not differ among themselves (Table 8). Most of the changes and interactions in runoff EC among treatment combinations were likely due to the initial EC differences among the soils (Table 1) and water types (Table 3).

Similar to EC, runoff-NO₃ + NO₂ concentration changes for Trial 2 did not differ from zero, except for the groundwater-Calloway-MAP, -ECSTreal, and -ECSTsyn, and rainwater-Dapue-MAP, ECSTreal, and -ECSTsyn combinations, where all Trial 2 changes were greater than zero (i.e., an increase in runoff-NO₃ + NO₂ concentration; Table 8). Similarly, most runoff-NO₃ + NO₂ concentration changes did not change over time, except for the groundwater-Calloway-MAP and rainwater-Dapue-MAP and -ECSTsyn combinations that increased over

^† Means in a column with different letters are different at P < 0.05. * An asterisk (*) indicates the mean value differs from zero at P < 0.05. ^‡ Chemically precipitated struvite, CPST; monoammonium phosphate, MAP; real-wastewater-derived electrochemically precipitated struvite, ECSTreal; synthetic electrochemically precipitated struvite, ECSTsyn; electrical conductivity, EC; nitrate + nitrite, NO₃ + NO₂.

time, where the Trial 2 change did not differ among themselves, and the rainwater-Calloway-MAP and -ECSTsyn combinations that decreased over time, where the Trial 2 changes did not differ between them (Table 8). Runoff-NO₃ + NO₂ concentration changes for most other water type-soil-fertilizer-P source combinations did not differ among themselves (Table 8). The change in runoff-NO₃ + NO₂ concentrations over time was likely due to the mineralization of organic matter, nitrification of NH 4 + , and/or the denitrification of fertilizer-N additions, especially in the Calloway soil, which had the largest pH and largest runoff-NO₃ + NO₂ changes among the four soils.

On a global scale, capturing and reusing P from wastewater instead of converting mined PR into fertilizer could potentially lead to an increase in global food security due to lowered reliance on a decreasing, finite resource. Struvite also has the potential to maintain or increase crop yields, potentially further increasing food production and security. On a regional scale, struvite could lower nutrient loads to already sensitive aquatic environments and resources. Recovery of struvite from wastewater also leads to a reduction of nutrients in wastewater effluent and sewage sludge. Furthermore, as results of the current study show, nutrient runoff concentrations are not substantially different from runoff concentrations from a widely used, commercially available fertilizer-P source (i.e., MAP). However, as there is not yet a pelletized form of ECST, this study showed that the timing of ECST application may influence runoff water quality. Furthermore, as there were few differences between runoff nutrient concentrations from MAP and ECST, the application of and potential runoff from ECST would not appear to exacerbate current environmental, nutrient-input issues. Because of reduced nutrient loads to receiving waters from nutrient recovery and having similar nutrient concentrations as MAP, struvite appears to be an alternative fertilizer-P source that may help minimize eutrophication of surface receiving waters, without negatively affecting runoff water quality in areas where struvite (i.e., ECSTreal or ECSTsyn) may be land-applied and incorporated.

Currently, it is too expensive to produce ECST on a commercial scale. Struvite is more costly than traditional fertilizers to produce, partially due to the experimental state of ECST materials and the lack of commercial-scale cost information [47] . However, struvite generally contains more N than many traditional fertilizer-P-only sources, such as TSP, which can reduce the need to purchase and apply additional fertilizer-N inputs in some cropping systems [48] . Furthermore, struvite also has the potential to outperform traditional fertilizer-P sources in terms of crop yield [48] [49] . Continued research is needed to compare the cost of struvite, particularly ECST, production to the production of traditional, commercially available fertilizer-P sources, but, as more research is conducted, there will be the opportunity to increase the efficiency and decrease the cost of struvite production.