Phosphorus Transport in a Lowland Stream Derived from a Tracer Test with 32 P

: Small streams in urbanized rural areas receive loads of P from various, often episodic, sources. This paper addresses, through a tracer test with 32 P, retention and transport of a pulse input of phosphorus in a 2.6 km long stretch of a channelized second-order lowland stream. Tritiated water was introduced alongside the 32 P-labelled ortophosphate in order to isolate the inﬂuence of the hydrodynamic factors on P behavior. Tracer concentrations in unﬁltered water samples were measured by liquid scintillation counting. Four in-stream and ﬁve hyporheic breakthrough curves were collected at four points along the stream, two of which encompass a beaver dam impoundment. The overall retention efﬁciency of 32 P along the studied reach was 46%. The transient storage transport model OTIS-P provided reasonable ﬁts for in-stream breakthrough curves (BTCs) but failed at reproducing the hyporheic BTCs. The overall small effect of transient storage on solute transport was higher in the stretch with a more pronounced surface storage. Transient storage and phosphorus retention were not enhanced in the beaver dam impoundment.


Introduction
Anthropogenic phosphorus enrichment is one of the principal factors in the deterioration of surface water quality, with rivers being an important component in human-affected phosphorus cycling. Rivers receive excessive phosphorus loads from urban and agricultural sources and carry this nutrient to lakes and coastal seas. Therefore, the understanding of phosphorus cycling and retention in streams is important for the improvement of eutrophication management tools in river catchments and in their recipient seas, which may not rely on load reduction only [1][2][3]. Riverine phosphorus cycling is complex and involves transformations between the inorganic-organic and dissolved-particulate forms of different bio-availability [4,5]. Furthermore, the riverine ecosystems and their internal phosphorus processing are constrained and controlled by hydrological and hydromorphological factors [6]. The coupled action of biotic and abiotic phosphorus cycling and downstream transport of dissolved and particulate phosphorus forms is conceptualized as spiralling, with the spiralling length reflecting the efficiency of nutrient utilization in stream ecosystems [7][8][9].
Due to the complexity of riverine phosphorus cycling the quantification of uptake rates and lengths is not straightforward and is best addressed through phosphorus addition experiments. Tracing of intentionally introduced substances (artificial tracers) at field and laboratory scales in order to understand the functioning and status of hydrological systems is an essential tool in water resources management [10][11][12][13][14]. In this regard, whole-stream additions of radioactive phosphorus present advantages in quantification of phosphorus cycling over stable phosphorus additions and batch experiments [15,16]. Phosphorus has two radioactive isotopes ( 32 P, T 1/2 = 14.268 d; 33 P, T 1/2 = 25.34 d) manufactured artificially by nuclear transmutation and used widely as radiopharmaceuticals and tracers in biochemical and biomolecular research. Both isotopes are also naturally produced in the atmosphere through interaction of cosmic radiation with argon nuclei [17,18], but their applications as the environmental tracers are limited to studies of atmospheric and oceanic mixing. Artificially produced 32 P and 33 P are used at field-and laboratory-scale to trace and quantify phosphorus uptake, speciation and transport in ecosystems and soils [19][20][21]. Applications related to the aquatic environment were initially concerned with the pathways of phosphorus uptake by biota in lakes [22,23] and rivers [24][25][26]. Results of tracer tests with radio-phosphorus were instrumental in the identification of the key features of phosphorus cycling in streams, such as uptake rates [7,8,27,28], food-web pathways [29] and the influence of transient storage [30]. An important advantage of radioactive phosphorus injections, comparing to stable phosphorus additions, is that they introduce a negligible amount of phosphorus, which does not affect stream uptake capacity. Despite the high relevance of radio-phosphorus tracing to deconvolute phosphorus cycling in rivers the application of this method at field scale is hindered, and in many countries practically impossible, due to the restrictions imposed by radiation protection regulations. Consequently, methods based on stable phosphorus addition to quantify the relative contribution of abiotic and biotic processes [16] and to quantify uptake lengths [31] have been refined [16,31] and methods based on tracing stable isotopic composition of oxygen in phosphate [32][33][34] have been developed.
Phosphorus retention in small streams in urbanized rural areas is an important, and probably globally relevant, but not well quantified aspect of phosphorus cycling in catchments. These streams receive significant loads of phosphorus from various, often episodic, sources that collectively constitute a semi-diffuse source not accounted for in pollution assessments [1,35,36]. Among such sources are farmyard, urban and road runoff, sewage and wastewater disposal, fish ponds. In rural Poland, where only 42% of the population is connected to sewerage systems, the improper wastewater management may be a significant source of nutrient pollution [37]. Wastewater effluent phosphorus may be effectively sequestered in stream biota and sediment and remobilized under storm flows [1,16]. However, large inputs of nutrients, whether continuous or episodic, may exceed the retention capacity of streams [38][39][40].
This work presents and discusses results of a tracer test with 32 P performed to characterize phosphorus retention and transport in a 2.6 km long stretch of a second-order rural stream subjected to multiple sources of phosphorus contamination. Tritium as tritiated water was injected simultaneously with 32 P-labelled ortophosphate in order to assess in-stream transport and transient storage exchange characteristics of a conservative tracer. Previous tracer tests with 32 P which aimed at quantification of phosphorus uptake rates/lengths in streams [7,8,27,28,30] were performed in short (<200 m) reaches of streams characterized by low mean discharges (<20 L/s), with tracers released at constant rates during 30-60 min periods. In comparison, the instantaneous injection performed in this study was traced during 6 days along the 2.6 km long stream stretch at the discharges varying between 100-200 L/s. While most phosphorus addition experiments aimed at quantification of phosphorus cycling in streams were performed with constant injection rates, the instantaneous injections provide equally reliable estimates and may be more appropriate at variable backgrounds [41], which might be a common feature of small streams in urbanized rural areas. Furthermore, pulse injections are more adequate for the analysis of long-term behavior of conservative solutes as affected by the exchange with the hyporheic and other transient storage zones [42,43]. Finally, the pulse injection of radioactive tracers simplifies the injection procedure and reduces exposure times of the personnel to ionizing radiation.
The main goal of this work is to characterize, in the timescale of the carried out tracer test (6 days), three subsequent stretches of the stream with respect to: (1) overall phosphorus retention capacity and (2) influence of transient storage on phosphorus retention and transport. The first goal is achieved by a comparison of tracer recoveries between the sampling points. The influence of transient storage is assessed through the metrics derived from parameters of the transient storage transport model OTIS-P [44]. The results of this study contribute to the understanding of the fate of the episodic inputs of inorganic phosphorus in channelized streams during the non-growing season. The relative importance of such semi-diffuse sources of phosphorus to small streams is expected to increase due to urban sprawl [45]. On the other hand, in the Northern Hemisphere the degraded stream ecosystems are now impacted by the expansion of beavers [46,47]. Accordingly, the presented study addresses the influence of a beaver dam impoundment on phosphorus retention and transport.

Study Site
The Kocinka (length 40.2 km; catchment area 257.8 km 2 ; average discharge at the outlet 1 m 3 /s) is a lowland river in Southern Poland, in the Baltic Sea basin. The river is polluted by nitrate of agricultural and wastewater origin and as such was a subject of an interdisciplinary study concerned with the development of novel methods for the governance of agricultural water pollution [48,49]. The river is groundwater-fed [50] with the streamflow and water quality being controlled by nitrate-polluted groundwaters discharging from a large fissured-karstic aquifer in Jurassic limestones and, in the headwater part, from shallow porous aquifers in Quaternary sands and gravels of fluvioglacial origin. Total dissolved phosphorus (TDP, 0.45 µm filtered) concentrations were observed occasionally in the period 2015-2017 at two sampling points encompassing the stream stretch studied in this tracer test ( Figure 1). Mean TDP concentrations at the upstream and downstream points were, respectively 1.810 mg/L and 0.948 mg/L. The high and variable (within one order of magnitude) TDP concentrations reflected domestic wastewater inputs from the partly urbanized catchment, including incidental sewer overflows and dumping of untreated sewage. A consistent drop in TDP between the two sampling points might be, at least partly, explained by dilution by lateral inflows received through ditches. The stream may also receive groundwater discharges directly through the bed, as evidenced in the more downstream parts of the Kocinka [50]. Indeed, stream discharge showed a significant downstream increase during the experiment (Table 1).  The artificially straightened stream channel is incised in a flat and wide valley bottom covered by meadows and wetlands ( Figure 2). Sampling points were selected to encompass river stretches with different hydromorphological conditions ( Figure 1). The distance between the injection point I and the first sampling point R1 (305 m) was large enough to ensure good transversal mixing of tracers, according to even the most conservative estimates [51,52]. Stream channel in stretches R1-R2 and R3-R4 has regular, parallel edges with the width and cross-section area fluctuating around 3.5 m and 1.3 m 2 , respectively. The R2-R3 stretch overlaps with a backwater zone upstream a beaver dam. Channel width and cross section area increase gradually along this stretch as water flow is slowed down.

Tracer Test Design
The usefulness of a tracer test in characterization of solute transport in rivers depends on the quality of tracer breakthrough curves (BTCs). The reliability of the estimates of different transport parameters depends on how well particular regions of a BTC are represented by sampling [53]. Additionally, given the short half-life of 32 P (14.268 days) the application of this tracer requires a careful design of sampling and analysis schemes. In order to reduce uncertainties in the analysis of low activity samples their measurements have to be performed as quickly as possible after sampling. Therefore the time frames for the passage of the low activity parts of the BTCs at particular sampling locations need to be predicted before tracer injection. Application of auxiliary dye tracers for tracing of the injected plume was inappropriate because of possible interferences in liquid scintillation counting. Therefore, the conservative tracer BTCs were predicted for the sampling points by the OTIS-P model (without transient storage) using discharge and water velocity values measured on the day preceding tracer injection and the dispersion coefficient derived from an independent tracer experiment. The BTCs obtained during the nitrate addition experiment performed in October 2017 in the stretch I-R1 (Q = 123 L/s) were used to estimate the dispersion coefficient. Based on these predictions, sampling schedule was designed for the period of 72 h after the injection. Sampling was extended to 147 h because samples from the tail regions of BTCs collected at points R3 and R4 appeared to have measurable 32 P activities.
Tracer injection was performed on 5 April 2018, 1 pm at point I ( Figure 1). The injected activities were estimated by the supplier as 30 GBq of 3 H (tritiated water) and 10 GBq of 32 P (aqueous solution of phosphoric acid) with 10% relative uncertainties. Potential doses of ionising radiation to the personnel, general public and aquatic organisms were estimated as not harmful and approved by the Polish Atomic Agency through their permit to introduce the above activities of both tracers into the Kocinka River. Tracers were delivered to the injection site by a licensed transport company. The glass ampules with tracers were crushed in a plastic container filled with 5 L of water. Contents of the container were emptied in a riffle area of the stream, which facilitated cross-sectional mixing of tracers. The weather was stable during the experiment, with no precipitation and air temperatures fluctuating between −2 to + 10 • C.

Tracer Sampling and Analysis
Water samples were collected from the central part of the stream channel by a plastic bottle attached to a sampling pole. Samples were immediately transferred to 22 mL HDPE liquid scintillation vials containing 0.1 mL of 0.1 N of H 3 PO 4 to prevent tracer adsorption on vial surfaces. Dilution by the added acid volume was taken into account in calculations of 32 P specific activities. Sample activities were later measured in the same vials. Contrary to the previous 32 P tracer studies conducted in streams our water samples were not filtered, primarily because this study was not concerned with the quantification of phosphorus pathways and pools in stream ecosystem. The main objective of this work was to assess an overall effect of the abiotic and biotic processes on the retention and transport of phosphorus in the stream. It must be also noted that filtration through the commonly used 0.45 µm membranes does not separate the finest particulates and the filtering process itself may affect dissolved P concentration in the filtrate [54]. Indeed, tests performed on 35 stream water samples showed a large variability (5%-76%) in the fraction of 32 P lost during filtration. Furthermore, a large fraction of phosphorus in surface waters is associated with iron-rich colloids [55,56] most of which pass through the 0.45 µm membranes. The truly dissolved, easily bioavailable ortophosphate ion is probably not the dominant movable form of phosphorus in rivers.
At points R3 and R4 water samples were collected also from the hyporheic zone by use of a dedicated sampler consisting of an array of stainless steel tubings pushed into bottom sediment. Each of the tubings had at its lower end two orifices of 0.1 mm diameter through which pore water was sucked out. Samples were collected from the depths of 5, 10 and 15 cm in sediment, however, at site R4 water could not be retrieved from the depth of 15 cm. Pore water samples were slowly sucked into syringes connected to the sampler through a 1.5 m long small diameter teflon tubing. Before sampling the dead volume of the teflon tubing was removed. Water collected in syringes was immediatelly transfered to vials. Each syringe was used only once to avoid cross-contamination of samples. The sampling arrays remained in sediment between samplings, however they were removed after the first 72 h and placed in the same locations for the extended part of the experiment. Activity concentrations of 3 H and 32 P in stream water samples were measured by liquid scintillation counting method with the Quantulus 1220 liquid scintillator beta/alpha spectrometer. Two different scintillators were used for 32 P (Insta-Gel Plus) and 3 H (Ultima Gold LLT) measurements. The ratios of scintillator to water sample volumes optimized to obtain the highest efficiency of the detection were 10 mL/10 mL for 32 P and 12 mL/8 mL for 3 H. After scintillators were added the vials were kept in the apparatus for 24 h before the measurement in order to reduce chemiluminescence. Measurements were performed for batches containing 18 samples, a reference activity solution and a blank sample. The activities of the reference solutions were used to determine counting efficiency for each batch of samples. Dead water was used as the 3 H blank and stream water collected before tracer injection was used as the 32 P blank. Therefore, the measured 3 H activities include the environmental component. Activity concentrations of the environmental 3 H measured in the Kocinka at point I in the years 2015-2016 varied from 7.2 × 10 −4 kBq/L to 8.6 × 10 −4 kBq/L. Due to the partial overlap of the 3 H and 32 P spectra the activity of each sample was measured twice. The first series of measurements was performed within 10 weeks after collection of samples, beginning with samples with the lowest expected activities. Measurements performed for 3 H and 32 P reference solutions allowed to identify the spectral overlap range and to evaluate a correction factor. The proportion of 32 P counts in the overlap range was 5% of total counts. The second series of measurements began 4 months later, after the 32 P activities of samples decreased to around 0.3% of the initial values and practically only 3 H was present in the samples, which improved the accuracy of 3 H activity determination. The measured activity concentrations of the tracers are presented in Appendix A.

Data Treatment and Otis-P Modeling
Because sampling extended for 43% of 32 P half-life, which corresponds to the radioactive decay of 26% of the initial activity, measured activities were corrected for the decay. While the transport model predicts the actual concentrations (specific activities) of the tracer at the moment of sampling, the estimates of tracer recovery at sampling points have to be based on the activities at the time of injection. Accordingly, the measured 32 P activities were recalculated for each sample for the moments of sampling and injection.
For the estimation of tracer recovery and of the retention and transport metrics the discharge, water velocity and channel dimensions have to be estimated. Particularly, the OTIS-P model requires quantification of discharge, lateral inflows and channel crosssection areas along the stream. Channel widths and cross-section areas were measured on the day before tracer injection at 20 cross-sections roughly evenly distributed between points I and R4. Discharge was measured a few hours before tracer injection at points R1 and R4 and 2 days after the injection at points R1, R3 and R4 (Table 1). Discharge was calculated by the area-averaged integration of water velocities measured by a 801 Valeport electromagnetic velocity meter. While discharge increased downstream due to lateral inflows, it also significantly decreased at each point during the course of the experiment. The input data for OTIS-P model were interpolated for different cross sections and times based on the measured discharges, assuming that discharges did not change after the second round of measurements performed on 7 April. The BTCs obtained at each point were treated as tracer input functions for subsequent stretches. Parameters of the transient storage model OTIS-P were then fitted to in-stream and hyporheic BTCs obtained at points R2-R4. The BTCs were integrated using the trapezoidal rule, taking into account variations of the discharge interpolated in time and space, to obtain total activities of tracers recovered at sampling points ( Table 2). Figure A1 presents processes and the corresponding model parameters used in OTIS-P to describe the advective-dispersive-reactive transport of solutes in the main channel of a stream, exchange of solutes between the main channel and the storage zones, and the reactive processes in the latter [44]. Model parameters were fitted stretch-wise treating the BTCs obtained at each sampling point as tracer input functions for subsequent stretches. The model was run at the spatial resolution of 1 m and the temporal resolution of 0.01 h. The downstream increase of discharge was implemented as lateral inflow, with a uniform distribution along each stretch. Tracer concentrations in the lateral inflow were assumed to be 8.6 × 10 −4 kBq/L for 3 H and 0 Bq/L for 32 P. Parameters of the conservative tracer model (D, A, A s , α) were fitted automatically to main channel BTCs, assuming constant discharge, by the use of the Nonlinear Least Squares technique implemented in OTIS-P with weights emphasizing the tail concentrations [44,57]. The A values were in the next step adjusted manually to account for the observed spatially and temporary variable discharge. The final conservative tracer model values of D, A, A s and α were sought by trial-and-error approach to obtain satisfactory fits to the hyporheic 3 H BTCs. The above values were used as the initial values for the 32 P model in which, again, automatic fitting was used for the main channel BTCs and trial-and-error fitting for the hyporheic BTCs.  Table 2 presents tracer activities recovered at subsequent sampling points. The BTCs were integrated using the trapezoidal rule, taking into account variations of the discharge interpolated in time and space, to obtain total activities of tracers recovered at sampling points ( Table 2). Assuming that tritium behaved like a conservative tracer, the downstream decrease of the ratio of recovered activities (the fifth column) reflected the gradual retention of 32 P along the stream. Because of the large uncertainty of the injected activities these ratios were recalculated relative to activities recovered at point R1 (the sixth column). Finally, the last column presents fractions of 32 P retained in the subsequent stretches of the stream relative to point R1. Table 3 presents metrics of 32 P tracer retention for the three subsequent stretches and for the whole studied stretch. Assuming that tracer removal rate along the stream is proportional to its amount remaining in water, the activities that pass through subsequent stream cross-sections must decrease exponentially with the distance [8]. The uptake length S w was thus calculated as:

Whole-Stretch Retention of Phosphorus
where x is the distance of a cross section from the injection point (j being further downstream than i) and A is total tracer activity that passed the cross section. The first-order uptake coefficient λ r : where v is water velocity, reflects the overall uptake capacity of the stream, including both the in-stream and transient storage uptake of the tracer through abiotic and biotic processes. It is not equivalent to either of uptake coefficients of the OTIS-P model, which describe separately uptake in the main channel and the storage zone [44]. The λ r values calculated from Equation (2) were corrected for radioactive decay by subtracting the decay constant of 32 P (λ 32P = 5.6 × 10 −7 [1/s]). The uptake velocity v f was used to compare retention capacity at varying stream depth [58]: where Q is discharge and W is stream width. The shortest uptake length S w , found for the beaver dam impoundment, did not translate into the highest uptake velocity v f , which had nearly the same value for stretches R1-R2 and R2-R3. Apparently, the sedimentation of coarser particulates that might be expected to occur in the impoundment [59] did not affect the movable fraction of 32 P. The stretch R3-R4, where water velocity was the highest among the three stretches, had a significantly lower uptake velocity.
The values of uptake velocity v f obtained in this study (Table 3) fall in the lower range of the estimates derived from stable phosphorus additions [16,41,[60][61][62]. Spiralling lengths derived from results of 32 P tracer tests for two contrasting streams with small and large volumes of transient storage zone were estimated at 643 m and 111 m, respectively [30]. Based on the data provided in [30] the respective uptake velocities could be calculated at 0.18 mm/min and 1.07 mm/min. A direct comparison of v f estimates obtained for rivers of different size, hydraulics, bed substrate and P loads and for different modes of P additions (pulse-continuous, stable-radioactive) was of limited value. Nevertheless, the consistency of our results with other studies indicated that tracing of all movable P fractions, without water sample filtration, did not bias the estimates of P retention metrics. Since colloidal organic and inorganic particles, majority of which is 0.45 µm filterable, may be the main carrier for phosphorus transport in streams [55,63] the unfiltered stream water samples properly represented the movement of a reactive tracer plume. The OTIS-P model was capable to fit reasonably well the in-stream 3 H BTCs for all sampling points, but the hyporheic BTCs could not be satisfactorily reproduced. At point R3 ( Figure 3) the model overestimated the hyporheic 3 H for the first 2 days (48 h) of the experiment and underestimated it for the later period. The spatial distribution of 3 H in sediment observed at R3 during the first 2 days (C 5cm < C 15cm < C 10cm ) is the result of competition between two effects that the depth in sediment has on hyporheic exchange: the hydrodynamically induced flux decreases with depth [64], while the hyporheic travel times are longer for the deeper penetrating flow paths. The apparent lag times between the in-stream and hyporheic BTCs at different depths arose due to not only delayed tracer penetration into sediment but also due to tracer dilution with the pre-experiment interstitial water. Therefore, the 15 cm BTC seemed to lag behind the stream BTC less than the 10 cm BTC. The same model explained the convergence of concentrations at 10 cm and 15 cm after the first 2 days and the excess of 3 H at 15 cm seen after 60 h. The linear plot inserted in Figure 3 confirms the above interpretation by showing the fastest response of the 5 cm BTC to in-stream tracer variations and the decrease of concentrations at 15 cm relative to 10 cm. These observations clearly confirmed that the well-mixed transient storage zone of the OTIS-P model did not properly represent actual distribution of tracers in the hyporheic zone [64,65]. Similar relationships between the 5 cm and 10 cm BTCs can be seen at point R4 ( Figure 4). In this case, however, the hyporheic concentrations increased between hours 60 to 120. Only a tentative explanation can be provided for the stabilization and increase of hyporheic concentrations seen during the second half of the experiment at R3 and R4. The decrease in stream discharge observed during the experiment, through its influence on stream water velocity and depth, could have resulted in the reduction of hyporheic exchange. The reduction of water flux into sediment might have retarded flushing of tracer from the hyporheic zone at R3 and facilitated release of the delayed tracer plume from the less-connected regions of the hyporheic zone at R4. Table 4 presents parameters of the OTIS-P model obtained through the inverse modeling (α-storage zone exchange coefficient, A-cross sectional area of the main channel, A s -cross-sectional area of the storage zone, D-dispersion coefficient [44]) as well various metrics used to characterize transient storage, derived from the above model parameters and the advective velocity u. Of these metrics only F 200 med is based on all three quantities (α, A s , v) that influence the effect of transient storage on tracer transport [66]. The F 200 medthe fraction of median travel time due to transient storage normalized to the 200 m long steam stretch-was evaluated as:

Characteristics of 3 H Transport
where L = 200 m.     Table 4 are relatively low, indicative of a small influence of transient storage on conservative solute transport. Our values of α are low, but still within the range reported by a meta-analysis based on results of 67 papers [67]. The modeled BTCs obtained for higher α did not reproduced well the observed concentrations underestimating them in the tail region. The attempts to fit higher α values did not significantly affect the A s values. Figure 5 shows an example of an alternative fit for point R2. On the other hand, lower α values led to the narrowing of BTC peaks and to the elevated concentrations in the tail. Our results showed that the α values can be overestimated when too much weight is placed on fitting of the peak of the BTC or when the tail region is not well represented by the measurements.
Transient storage in the studied stream may be associated with the hyporheic zone, narrow strips of dead vegetation extending along the otherwise straight and regular banks, patches of submerged vegetation and debris accumulations on the streambed. The metrics of transient storage reflected the overall small effect of these transient storage zones. While the metrics related to travel times in the main channel and storage zone as well as their sum (t m mean , t s mean , t mean ) were largest for the beaver dam impoundment (R2-R3), the F 200 med values indicated that transient storage accounted for the largest fraction of total travel time in stretch R1-R2 and for the smallest fraction in stretch R3-R4. Apparently, effects of transient storage on fluxes of conservative and non-conservative solutes were not enhanced in the beaver dam impoundment. Stretch R3-R4 differed from the more upstream part of the stream by higher water velocity, coarser bottom sediments and less frequent occurrence of submerged vegetation. The higher F 200 med values between R1 to R3 reflected a higher role of surface storage comparing to R3-R4, where the hyporheic storage became relatively more significant.

Characteristics of 32 P Transport
Figures 7-9 present semi-log plots of the observed and modeled BTCs for points R2-R4. As for tritium, the OTIS-P model was capable to fit reasonably well the in-stream 32 P BTCs for all sampling points. Differently than for 3 H the highest hyporheic concentrations of 32 P occurred at 5 cm depth. Only the 5 cm BTCs showed a roughly monotonic decay after peak concentrations which could be reasonably reproduced by the model for R4, except for the three highest concentrations. Concentrations at 10 and 15 cm fluctuated showing at R3 distinct minima during the second 24 h of the experiment. Transient storage of 32 P in the hyporheic zone occurred mostly in the top layer of sediment, which, when compared with the deeper penetration of 3 H, reflected filtration of particulate P and/or adsorption of dissolved P. The irregular behavior of 32 P at greater depths might have been partly related to the changing patterns of hyporheic flow, as discussed above for 3 H.
A comparison of OTIS-P parameters and various metrics of transient storage presented in Table 4 for 3 H and in Table 5 for 32 P shows similar values of parameters and metrics connected with tracer transport in the main channel (A, u, D, P D , t m mean ). This provides more evidence that non-filtered water samples adequately represented in-stream transport of 32 P. Additionally, the storage zone exchange coefficient α was similar for both tracers, except for R1-R2 where it was almost twice higher for 32 P. Consequently, L s -the average distance a molecule travels downstream within the main channel prior to entering the storage zone-was roughly twice lower for 32 P in R1-R2. The higher than in other stretches capacity of transient storage zone in R1-R2 to capture tracers from the main channel was additionally enhanced for 32 P. As could be expected the parameters and metrics characterizing transient storage were higher for reactive 32 P than for non-reactive 3 H. Interestingly, the F 200 med values obtained from 32 P and 3 H BTCs varied between stretches almost proportionally, the former being 6.5-6.9 times higher. The effect of transient storage on 32 P flux was thus proportional to the median of the travel time it spent in storage zones. The F 200 med metric correlated well (R 2 = 0.997) with the effective uptake coefficient λ r (Equation (2), Table 3). This correlation confirmed that the enhanced 32 P retention in R1-R2 was not due to a higher retention capacity of sediments or periphyton, nor to filtering of the coarser 32 P-bearing particles but was related to the fraction of median travel time due to storage being the highest in that stretch. A direct comparison between the uptake coefficient λ r derived from overall tracer recoveries (Table 3) and the first order uptake coefficient λ obtained by inverse modeling of 32 P BTCs (Table 5) was meaningless because λ r represented the apparent uptake in the timescale of the experiment, not distinguishing between the unidirectional uptake (represented in OTIS-P by λ and λ s ) and the reversible sorption process. Nevertheless, the lambdas derived by both methods differed for the same stretches by not more than a factor of 2.

Conclusions
The 32 P BTCs obtained for unfiltered stream water samples provided a consistent picture of the retention and transport of a pulse input of phosphorus in the 2.6 km long stream reach comprising stretches of diverse hydromorphological conditions. The overall retention of phosphorus tracer in the studied stream reach in the timescale of 6 days was 46%, with the least efficient retention in the stretch differing from others by higher water velocity, coarser bottom sediments and less developed surface storage. The OTIS-P model of transient storage provided reasonable fits to in-stream BTCs of the simultaneously injected 3

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: