Vegetation and Residence Time Interact to Inﬂuence Metabolism and Net Nutrient Uptake in Experimental Agricultural Drainage Systems

: Agricultural drainage networks within the Lower Mississippi River Basin (LMRB) have potential to attenuate nutrient loading to downstream aquatic ecosystems through best management practices. Nutrient uptake (nitrogen, phosphorus), gross primary production (GPP), ecosystem respiration (ER), and denitriﬁcation rates were estimated using a combination of sensor measurements and hourly discrete samples for dissolved nutrients and gases at three hydraulic residence times (2, 4, and 6 h) in three vegetated and three unvegetated ditches. We also measured vegetation and soil nutrient content. GPP and ER were signiﬁcantly higher in vegetated drainages and increasing hydraulic residence time increased respiration rates. Shorter hydraulic residence times were associated with increased uptake rates for both N and P, and vegetation increased N uptake rates in all hydraulic residence time (HRT) treatments. Vegetation and sediment assimilated N and P over the course of the experiment. Overall, our experimental results demonstrate the strong role of emergent vegetation in nutrient retention and removal processes in agricultural drainage ditch networks. Author Contributions: Conceptualization, R.L.N. and J.M.T.; methodology, R.L.N. and J.M.T.; soft-ware, R.L.N. and J.M.T.; validation, R.L.N. and J.M.T.; formal analysis, R.L.N.; investigation, R.L.N. J.M.T.; resources, J.M.T.; data curation, R.L.N.; writing—original draft preparation, R.L.N.; writing—review and editing, J.M.T.; visualization, R.L.N.; supervision, J.M.T.; project administration,


Introduction
It is estimated that 70% of the nitrogen (N) and phosphorus (P) loads to the Gulf of Mexico via the Mississippi River originate from agricultural sources within the basin [1]. N and P runoff may result in a higher incidence of harmful algal blooms, hypoxia-induced degraded habitat for freshwater and marine biodiversity and commercial fisheries, and impacts to other economies that rely on the availability of clean water, all phenomena associated with the outflow of the Mississippi River [2,3]. Reducing ecological impacts of nutrient runoff will require improving nutrient use efficiencies on farms but also require land managers to mitigate nutrient loading through best management practices (BMPs).
Agricultural areas within the Lower Mississippi River Basin (LMRB) are drained by networks of low-gradient, open-canopy ditches and streams that serve as important intermediate aquatic ecosystems between agricultural fields and downstream water bodies that support a variety of uses. These drainage networks are dominated by ditches that traditionally have been viewed as simple conduits to route excess moisture away from farm fields [4]. In fact, the productive farmland of the Mississippi Alluvial Plain (MAP) was not utilized until the early 20th century when much of the region was drained via ditches [5]. While this drainage improved production on existing agricultural lands; actions also encouraged further draining, resulting in the encroachment of agriculture into previously undisturbed lands [5]. Increased and sustained profitability of large-scale row-crop agriculture within the MAP has required increasingly intensive fertilizer inputs of N and P, and yields within the region have doubled in response to the four-fold increase in N fertilizer use since 1950 [6]. Approximately 40% of applied N is utilized by plants [7] with the remainder running off agricultural land into adjacent ditches. Today, drainage shallow, low-gradient ditches or streams which drain farmed lands of the lower MAP. Each flume had a 0.34 m 3 head pool followed by a sequence of run pool sections, where each section had a volume of 0.91 m 3 , resulting in a total volume of 3.96 m 3 per flume. Flumes were filled with 1.54-1.74 m 3 of a mixture of Dundee silty loam and Sharkey clay soil typical of the nearby MAP region to provide a relatively consistent sediment depth of~5 cm across all sections. Water was supplied by a pond located~50 m from our experimental setup. This water had stable nutrient concentrations (0.091-0.224 NO 3 -N mg L −1 , 0.004-0.012 PO 4 −3 -P mg L −1 ) and entered each flume through a PVC pipe at the main headpool and flowed through the run, pool, run, pool sequence before exiting through an outlet standpipe in the final pool. Outlet water flowed into a series of small wooded and vegetated wetlands. Half of the ditches (# 2, 3, and 6) were vegetated with Leersia oryzoides while the other half (# 1, 4, and 5) had bare soil ( Figure 1B). Three months prior to the start of the experiment, L. oryzoides was planted at a density of approximately 132 stems m −2 (µ ± SE, 132 ± 27) from a nearby floodplain wetland complex area located on the field station property. We chose L. oryzoides as our experimental plant because our previous small-scale lab and mesocosm studies have demonstrated its potential to enhance nutrient removal [36]. The ditch complex was maintained with flowing water gravity fed by the nearby pond and monocultures in vegetated ditches were maintained by hand weeding.

Study Site
A set of 6 artificial flumes at the University of Mississippi Field Station in Abbeville, MS was employed for this experiment ( Figure 1A). Flumes were designed to simulate shallow, low-gradient ditches or streams which drain farmed lands of the lower MAP. Each flume had a 0.34 m 3 head pool followed by a sequence of run pool sections, where each section had a volume of 0.91 m 3 , resulting in a total volume of 3.96 m 3 per flume. Flumes were filled with 1.54-1.74 m 3 of a mixture of Dundee silty loam and Sharkey clay soil typical of the nearby MAP region to provide a relatively consistent sediment depth of ~5 cm across all sections. Water was supplied by a pond located ~50 m from our experimental setup. This water had stable nutrient concentrations (0.091-0.224 NO3-N mg L −1 , 0.004-0.012 PO4 −3 -P mg L −1 ) and entered each flume through a PVC pipe at the main headpool and flowed through the run, pool, run, pool sequence before exiting through an outlet standpipe in the final pool. Outlet water flowed into a series of small wooded and vegetated wetlands. Half of the ditches (# 2, 3, and 6) were vegetated with Leersia oryzoides while the other half (# 1, 4, and 5) had bare soil ( Figure 1B). Three months prior to the start of the experiment, L. oryzoides was planted at a density of approximately 132 stems m −2 (μ ± SE, 132 ± 27) from a nearby floodplain wetland complex area located on the field station property. We chose L. oryzoides as our experimental plant because our previous smallscale lab and mesocosm studies have demonstrated its potential to enhance nutrient removal [36]. The ditch complex was maintained with flowing water gravity fed by the nearby pond and monocultures in vegetated ditches were maintained by hand weeding. We conducted three experimental loading events across all ditches on 3 different dates in a 15 d period. Each ditch was randomly assigned to one of 3 HRT so that over the We conducted three experimental loading events across all ditches on 3 different dates in a 15 d period. Each ditch was randomly assigned to one of 3 HRT so that over the 3 sampling events, each ditch experienced a 2, 4, and 6 h residence time [43]. HRTs were set manually with PVC spigot structures and verified with conductivity tracer tests [43]. During each nutrient loading event we dripped a solution of 17.92 and 8.3 g L −1 of NaNO 3 and KPO 4 into each ditch at a constant rate with a peristaltic pump. Pumps were calibrated to deliver nutrient stock solution at 5 mL min −1 into the headpool of each drainage, where the solution could mix thoroughly with the inflow water. We purposely loaded all ditches (15.5 ± 0.6 and 9.9 ± 0.4 g of NO 3 -N and PO 4 -P, respectively; Table 1) at the same rate regardless of HRT treatment to better separate the effects of HRT on reduction in a specific load. We acknowledge that this creates different nutrient concentrations experienced by experimental ditches under each HRT treatment. The estimated expected nutrient concentrations for each HRT treatment were 1, 2 and 4 mg L −1 NO 3 -N and 0.5, 1 and 2 mg L −1 PO 4 −3 -P for HRT 2, 4 and 6, respectively.

Dissolved Nutrient Sampling
We conducted diel dissolved nutrient sampling at upstream and downstream locations within each flume to estimate nutrient uptake rates. At each hour for a 24 h period, dissolved nutrient samples were collected and analyzed for nitrate (NO 3 − -N), and phosphate (PO 4 −3 -P). Samples were collected in 60 mL plastic syringes followed by pushing water through a 0.45 mm syringe filter into collection vials. All dissolved nutrient samples were frozen after collection until analysis. Dissolved nutrient samples were analyzed using the Cd reduction, and molybdate methods [44] on a Lachat QuickChem 8500 Series auto analyzer (Lachat Instruments) to measure NO 3 − -N and PO 4 −3 -P concentrations, respectively.
We estimated reach-scale nutrient uptake rates on an hourly basis [45]. Travel time for a given parcel of water between the upstream and downstream sampling stations was accounted for using the following equation: where U a,t is uptake at hour (t) (mg N or P m −2 h −1 ), Q is the flow rate (L h −1 ), A is benthic area in m 2 , N upstream,t is the measured nutrient concentration at hour t, and N downstream, t + i is the measured nutrient concentration at the hour plus travel time (mg L −1 ). Given the timing of sampling, each drainage ditch had between 18 and 22 hourly uptake rate estimates.

Reach-Scale Metabolism Estimates
Reach-scale metabolism rates were estimated based on changes in dissolved oxygen (O 2 ) through time with a two-station method [46,47]. The two-station method was appropriate due to expected low reaeration coefficients (K) and a short distance to source water which can confound dissolved gas flux estimates if experimental systems are not treated discretely [48,49]. We used 2 HOBO ® (U26-001, Onset ® Computer Corporation, Bourne, ME, USA) loggers placed in the middle of each ditch at upstream and downstream ends to measure dissolved O 2 and temperature at 15 min intervals during each 24 h nutrient loading event. A Bayesian inverse modeling procedure was used to estimate GPP, ER and K from dissolved O 2 data using the following equation [50]: where O 2up(t) is the upstream O 2 concentration (g O 2 /m 3 ) and O 2down(t + τ ) is the downstream O 2 concentration of that same parcel of water following travel time, τ. GPP and ER are expressed in g O 2 m -2 d -1 and represent positive and negative rates of O 2 production and consumption, respectively. z is mean depth (m) calculated as stream volume (m 3 ) divided by surface area (m 2 ), O 2,satup and O 2,satdown are the saturation concentrations upstream and downstream (g O 2 m −3 ), K (d -1 ) is the gas exchange rate, and PPFD is the sum of the photosynthetic photon flux density (µmol m -2 s -1 ) accumulated in the time interval from t to (t + τ) divided by the daily total of PPFD (PPFD total ) for any given parcel of water.
We used the function metrop in the mcmc package [51] in R to sample the posterior probability distribution P (µ|D). Each chain was run for 20,000 iterations following a 1000 iteration burn-in, with no chain thinning, and probability distributions adjusted to achieve an acceptance rate of~20%. We used minimally informative prior probability distributions for GPP, ER, and K (GPP~N [µ = 5, SD = 10]; ER~N µ = -5, SD = 10]; K[ µ = 5, SD = 5]). The model can have difficulties estimating all three parameters due to equifinality [52], but our previous work on the influence of measured verses modelled K conducted in study flumes demonstrates that the models are robust without measured K estimates in systems with high GPP and ER rates and low K [43].

Denitrification Estimates
Hourly triplicate dissolved gas samples were collected for dinitrogen (N 2 ) and argon (Ar) [53,54]. Using a polyvinyl chloride sampler, exetainer vials were overfilled 3x, preserved with zinc chloride (ZnCl 2 ), capped, and stored until analysis. We measured dissolved N 2 :Ar in dissolved gas samples with a Membrane Inlet Mass Spectrometer (MIMS) (Bay Instruments, Easton, MD, USA) [49,55]. Sample N 2 :Ar for each sample was converted to N 2 -N concentrations based on the following equation [56]: where DG/Ar sample is the measured dissolved gas sample signal and DG/Ar standard is the measured dissolved gas signal for well-mixed, deionized water open to the atmosphere at the same temperature as the sample. The terms [Ar] expected and [DG]/[Ar] expected are the theoretical saturated concentration and ratio, respectively, calculated for each field sample using gas solubility tables [57]. Diel patterns of dissolved gas data indicated the strong patterns in O 2 concentrations throughout the 24 h period which corresponded with diel N 2 patterns [43]. We developed and applied a correction factor to eliminate instrument specific O 2 interference with N 2 :Ar measurements which may have contributed to below equilibrium N 2 :Ar ratios observed during daytime hours [43,58]. N 2 :Ar values were then used to model daily fluxes in a Bayesian framework [50,53]. Specifically, we used a single-station approach which allows for the possibility of N 2 consumption to estimate posterior distributions for denitrification and K using the following model [49]: where N 2(t) is the dissolved N 2 at time t (g m −3 ); DN is the denitrification rate (g m −2 d −1 ), z is depth of water (m), N 2,equil(t) is the equilibrium concentration of N 2 at time t (g m −3 ); and K t is the gas exchange rate (d −1 ) for N 2 at time t calculated from K 600 . K 600 is utilized to normalize the model daily estimates of K from the metabolism model (Equation (2)) as K 600 is comparable across temperatures and gases [59]. N consume is the N 2 consumption rate (g m −2 d −1 ) from ebullition or fixation. For any parcel of water, the fraction of light it accumulates is photosynthetic photon flux density (PPFD t , µmol m −2 s −1 ) at time t divided by the daily total of PPFD (PPFD total , µmol m −2 s −1 d −1 ). DN was transformed from a daily rate to an average hourly rate of daily DN flux in mg m −2 h −1 . We acknowledge the application of the single station model may violate minimum reach assumptions [60] and thus the rates should be used with caution; however, the relative magnitude of denitrification and uptake as well as qualitative differences in rates across the vegetation and hydraulic residence times may still be informative to land managers because any effects from the upstream water source would be equal across the treatments.

Vegetation and Sediment Data Analysis
Vegetation samples were collected from each drainage ditch to analyze differences between pre-experimental conditions and post-experimental conditions (15 d). Triplicate shallow cores were taken from each vegetated drainage. Samples were transported to the lab, where sediment was washed away, and vegetation gently cleaned using deionized water. Each sample was divided into above-and belowground tissues, placed in paper bags, and dried at 60 • C until mass stabilized. We measured the carbon (C) and N content of plant tissues using a Vario Max CNS elemental analyzer (Elementar Analysensysteme GmbH, Hanau, Germany). The phosphorus (P) content of plant material was measured by combusting a 4 ± 1 mg sample at 500 • C for 2 h followed by digestion of material in 1 N hydrochloric acid for 30 min at 85 • C [61]. Samples were diluted with deionized water and analyzed with the molybdate method [44] using a Turner Designs Trilogy laboratory fluorometer (Model # 7200-000) equipped with the phosphate absorbance module (Model # 7200-070; Turner Designs, Sunnyvale, CA, USA).
Sediment samples were collected from each drainage ditch to analyze differences among C and N content between pre-experiment and post-experiment sampling periods. Samples were oven dried at 50 • C to a constant mass, weighed to estimate dry mass (DM), and ground with a Thomas-Wiley Intermediate Mill (Thomas Scientific). We measured C and N content of oven-dried material in a Vario Max CNS elemental analyzer (Elementar) using ground peach leaf standard (National Institute of Standard and Technology [NIST] Standard Reference Material [SRM] 1547) to correct all elemental analyses based on recovery efficiencies.

Statistical Analyses
Linear mixed-effects models that included random effects (~1|drainage) were used to account for multiple measurements taken across the study from the same drainage. All models were fit with the restricted maximum likelihood criterion. We assessed assumptions of all models visually with normality plots (qqnorm) and standardized residual plots across treatments [62]. When variances in error differed across treatments, variance was modeled separately among treatments with the VarIdent command [62]. The nlme package was used to run linear mixed-effects models [63]. We used emmeans within the emmeans package to assess interactions between vegetation and hydraulic residence time when significant [64]. All analyses were run in R [65].

Agroecosystem Gas Fluxes
Mean rates of GPP were approximately 2-fold higher and ER rates were more than doubled in vegetated versus unvegetated ditches (Figure 2). Gross primary production was significantly higher in vegetated drainages than unvegetated drainages across all hydraulic residence times (F 1,4 = 14.81, p = 0.02; Figure 2A). Vegetated ditch GPP was highest in the 6 h HRT (8.84 ± 1.56 g O 2 m −2 d −1 , µ ± se), followed by the 4 h (7.91 ± 1.56) and 2 h treatment (7.08 ± 1.16). Gross primary production was not significantly different among HRT treatments (p = 0.23) nor were there significant interactions between vegetation and HRT (p = 0.91).

Agroecosystem Gas Fluxes
Mean rates of GPP were approximately 2-fold higher and ER rates were more than doubled in vegetated versus unvegetated ditches (Figure 2). Gross primary production was significantly higher in vegetated drainages than unvegetated drainages across all hydraulic residence times (F1,4 = 14.81, p = 0.02; Figure 2A). Vegetated ditch GPP was highest in the 6 h HRT (8.84 ± 1.56 g O2 m −2 d −1 , μ ± se), followed by the 4 h (7.91 ± 1.56) and 2 h treatment (7.08 ± 1.16). Gross primary production was not significantly different among HRT treatments (p = 0.23) nor were there significant interactions between vegetation and HRT (p = 0.91).  Vegetation also significantly increased rates of ER (F 1,4 = 25.95, p = 0.01) and HRT had a marginally significant influence on ER (p = 0.06). Again, the highest rates of ER were observed in the 6 h HRT treatment within vegetated streams (−7.20 ± 1.47 g O 2 m −2 d −1 , µ ± se) followed by the 4 h (−5.69 ± 1.30) and 2 h (−4.72 ± 1.34) treatments. A similar pattern was observed for unvegetated drainages, with 6 h HRT treatment exhibiting the highest respiration values (−2.56 ± 1.09 g O 2 m −2 d −1 ) followed by 4 h and 2 h HRT treatments. Though fluxes in unvegetated ditches were smaller than vegetated ditch fluxes, no significant interaction between vegetation and HRT (p = 0.78) was observed ( Figure 2B).
Rates of denitrification were higher in vegetated ditches across all HRTs though differences were not statistically significant (

Nutrient Uptake
Nitrate uptake rates were more than 2-fold greater in vegetated ditches but decreased with increasing HRT (Figure 3, Table 1). HRT had overall significant effects on nitrate net uptake(F 1,351 = 45.35, p < 0.001), but was also influenced by a significant interaction between HRT and vegetation with stronger effects of HRT observed in vegetated compared to unvegetated ditches (F 1,351 = 4.29, p= 0.04; Figure 3A). Initially, observed concentrations were similar across HRTs (2 h = 2.43 ± 0.46 mg L −1 ; 4 h = 2.2 ± 0.47; 6 h = 2.41 ± 0.53). Shorter residence times (2 h) displayed significantly higher uptake rates and removal (32% in unvegetated ditches) than longer residence times (4 h and 6 h, 21% and 17% for unvegetated ditches) with the magnitude of fluxes higher in vegetated drainage ditches (2 h: 63%, 4 h: 44%, 6 h: 32%). For P net uptake, HRT treatment was found to significantly influence rates overall (F 1,351 = 4.23, p = 0.04) but differences between 2 h, 4 h, and 6 h were not significant at the p < 0.05 level ( Figure 3B). There was no significant overall effect of vegetation as similar patterns between unvegetated and vegetated drainages were observed in the 2 h treatment (14% vs. 16% P removal) and rates exceeded 7 mg P m −2 h −1 (F 1,4 = 2.41, p = 0.19; Table 1). The interaction of vegetation and HRT was not significant (F 1,351 = 0.17, p = 0.68); rates of P uptake in vegetated drainages remained high while rates in unvegetated drainages declined in the 4 h and 6 h treatments ( Figure 3B).

Vegetation and Sediment
Stem density measured in the post-experiment sampling period (µ = 48) was higher than the pre-experimental sampling period (µ = 29) though differences were marginally significant (t = −2.06, df = 11.13, p = 0.06). However, we observed differences in nutrient

Vegetation and Sediment
Stem density measured in the post-experiment sampling period (µ = 48) was higher than the pre-experimental sampling period (µ = 29) though differences were marginally significant (t = −2.06, df = 11.13, p = 0.06). However, we observed differences in nutrient content and ratios for aboveground vegetation metrics between pre-experiment and postexperiment sampling periods (Figure 4). Aboveground vegetation percent N, percent P were significantly greater, molar C:N and C:P ratios were lower, and N:P ratios were greater in post-experiment samples ( Figure 4B-F), indicating that uptake of N and P in plant tissue occurred during the experiment. Aboveground percent carbon was not significantly different between sampling periods ( Figure 4A). No significant differences in belowground tissue metrics was observed including percent C, percent N, percent P, C:N, C:P, nor N:P ( Table 2).   There were no significant differences between pre-and post-sampling periods for soil percent C and percent N ( Figure 4G,H). However, molar C:N ratios in soils were found to be significantly higher during the pre-experiment sampling period than the postexperiment sampling period ( Figure 4I; 12.74 ± 0.29 pre, 11.84 ± 0.27 (µ ± se) post, t = 2.23,  There were no significant differences between pre-and post-sampling periods for soil percent C and percent N ( Figure 4G,H). However, molar C:N ratios in soils were found to be significantly higher during the pre-experiment sampling period than the post-experiment sampling period ( Figure 4I; 12.74 ± 0.29 pre, 11.84 ± 0.27 (µ ± se) post, t = 2.23, df = 105.42, p = 0.03). This was likely driven by accumulation of N in sediments as evidenced by the higher though non-significant %N values during the post-experiment sampling period (0.06 ± 0.0018 post vs. 0.056 ± 0.0019 pre, µ ± se).

Discussion
Drainage ditches are a ubiquitous feature within the MAP and other alluvial plain agroecosystems throughout the world. Our results provide much needed experimental but reach-scale support for the hypothesis that the presence of wetland vegetation in drainage networks will increase rates of biogeochemical processing including GPP, ER, and nutrient removal. Surprisingly, we also found that while longer hydraulic residence times did increase metabolism rates, net nutrient uptake rates actually decreased, contrary to our second hypothesis. We also expected that the additive effects of vegetation and increasing HRT would result in higher metabolism and nutrient retention rates in vegetated ditches with longer HRTs. This third hypothesis was also supported for metabolism, but we found that vegetated ditches combined with 2 h HRTs had the highest N retention rates. These confounding results may be due to unique discharge-nutrient concentration relationships created by our experimental design but may also be attributable to more physical mixing at higher flows and are discussed further below. Overall, our results provide experimental evidence that reach-scale biogeochemical processes that govern nutrient uptake are enhanced by managing wetland vegetation within agricultural ditches.

Vegetation
The presence of vegetation significantly increased gross primary production, ecosystem respiration, as well as rates of N and P uptake across all hydraulic residence times (Table 1). Increasing vegetative nutrient content combined with decreasing C:N and C:P ratios indicates effective short-term nutrient removal occurred within ditches via luxury uptake [66][67][68]. As no significant belowground tissue changes were observed, it appears L. oryzoides roots are able to effectively transfer nutrients to aboveground tissues in this environment [69] Plant nutrient uptake and storage may be underestimated by our data because we rinsed plant material, potentially dislodging epiphytic algae, bacteria, and associated nutrient storage. Unaccounted for uptake and storage may help to explain strong patterns of enhanced GPP and ER coupled with weaker nutrient storage patterns within vegetated systems. However, macrophyte cover has been shown to enhance GPP and ER rates and contribute proportionally more to ecosystem metabolism than other plant forms; moreover, plant biomass and chlorophyll a have been found to explain the majority of variation in daily GPP (80%) and ER (63%) in agricultural streams [70,71]. Additionally, C assimilation and processing via metabolic processes require nutrient inputs; plant and associated microbial uptake can be a significant nutrient removal mechanism in vegetated aquatic systems [42]. In smaller-scale studies, mesocosms planted with L. oryzoides removed >90% of the retained NO 3 − -N and >70% of PO 4 −3 -P [36,39,72].
Higher GPP observed in systems vegetated with L. oryzoides likely fuels higher ER rates, which may enhance conditions favoring denitrification by reducing O 2 and increasing C [43,73]. Previous investigations have found that models that predict rates based on diel N 2 patterns were best fit by including terms for both N 2 consumption (i.e., N fixation, ebullition) and production (i.e., denitrification), particularly in highly productive, lower flow reaches; estimates of denitrification among reaches of various flows were further improved by using a two-station approach [49]. The use of a single-station approach, which introduces the potential for upstream influence on rates, necessitates denitrification rate values reported herein be used with caution. However, the observation that reach-scale denitrification increased in ditches planted with L. oryzoides is reasonable, as all estimates were equally influenced by the upstream water source, yet rates were higher in vegetated ditches. Additionally, a previous comparison of one and two-station models in larger experimental ditches found minimal differences in estimated denitrification rates, suggesting that biological activity within short reaches can have a larger effect on N metabolism than potential upstream influences, providing some support for our estimates [49]. Further, our previous smaller scale mesocosm studies have demonstrated that a significant proportion of N retention can be attributed to high denitrification efficiencies in L. oryzoides beds [36,37,39].
While the idea that aquatic vegetation increases metabolism and nutrient uptake in streams and ditches is not necessarily novel, this study fills a critical intermediate gap between small-scale experiments and field observations. Our results scale up previous mesocosm studies which are often constrained by the small temporal and spatial scale of the experiments and demonstrate approaches based on diel data collection can be used to make ecosystem estimates of C, N and P cycling processes in agricultural ditch environments. Our experiment demonstrates the benefits of managing L. oryzoides in agricultural drainage networks by quantifying whole system metabolism and nutrient uptake processes in a replicated experimental design conducive to testing hypotheses. Our approach represents a balance between larger scale, unreplicated field studies [13,42] and small-scale, mesocosm studies [36,72] that have all demonstrated enhanced nutrient retention through managing vegetation in agricultural ditches. Furthermore, our results demonstrate that vegetated linear systems can achieve rates of nutrient removal (up to 63%) comparable to or approaching those of constructed wetlands (56-79%) [74,75].

Residence Time
Our experimental manipulations provided support for our hypothesis that higher rates of gross primary production and ecosystem respiration would be observed in longer HRT treatments. Field studies in streams have found that GPP and ER rates increase at lower discharge, likely because of greater accumulation and development of autotrophic and heterotrophic communities at low flows compared to high flow that may disturb biological controls on metabolism [76,77]. All flow rates were extremely low so it is more likely that higher rates of metabolism would be observed with greater assimilatory capacity at lower flows.
However, the highest rates of net nutrient uptake were observed in the 2 h HRT treatment which was contrary to expectations. In the MAP, storm events often create a combination of high nutrient loads and high flow rates can lower concentrations and allow a significant portion of the load to escape the uptake capacity of different processes in streams [78,79]. In contrast, within this study there were no rainfall events, and we supplied equal nutrient loads to all streams and varied discharge to manipulate HRT but not load. This should have caused a discharge × nutrient concentration relationship where inflow nutrient concentrations were greater in ditches with higher HRTs. However, examination of inflow NO 3 -N concentrations across HRT treatments revealed this was not the case (2 h = 2.43 ± 0.46 mg L −1 ; 4 h = 2.2 ± 0.47; 6 h = 2.41 ± 0.53) and highlights the difficulty of managing experimental dosing regimens across extremely low flows which are prone to error. However, all inflow NO 3 --N concentrations were above values previously reported as associated with maximum NO 3 --N retention (1 mg L −1 ) in sediment cores vegetated with L. oryzoides [37], so observed differences in HRT are likely driven by some other factor than potential discharge effects on uptake efficiency. A study of vegetated ditches found lower flows to physically limit N removal [80]; here, the highest flows in the 2 h HRT could experience the highest mixing activity due to increased turbulence associated with slightly higher velocities at the sediment water interface thus enhancing N uptake. Epiphytic algae, while visibly negligible, was not captured in these rates and may account for lower observed uptake rates with higher respiration rates in lower HRT treatments. Longer residence times may also create preferential flow paths and actually decrease contact between reactive sites and nutrients as longitudinal dispersion becomes dominated by inefficient exchange in dense canopies with shallow submergence [81]. Finally, these results are somewhat contrary to the literature suggesting that manipulating for longer HRTs through water control structures in agricultural ditches may enhance nutrient processing and retention [19,20]; our results are more like the result of an unexpected experimental artifact or a potential HRT range that exceeds thresholds for positive impacts on nutrient processing. Unfortunately, facilitating shorter residence times, in order to examine a wider range of HRT, while still maintaining sensitivity in our uptake measurements will likely require longer experimental reaches. More research on nutrient uptake patterns across a wider gradient of residence time in agricultural ditches is needed.

Conclusions
Drainage ditches represent a diffuse interface between agricultural fields and downstream freshwater networks that may be manipulated to enhance nutrient processing and retention through management of vegetation and HRT. This study confirms findings from numerous field-and small-scale mesocosm studies by providing evidence at the reach scale that vegetation enhances metabolism and net nutrient uptake under controlled, whole-ditch conditions. Vegetated ditches and increasing hydraulic residence time have been investigated separately, but have not been experimentally tested together. This study attempted to examine the efficacy of combining conservation practices in improving within ditch N and P uptake rates. Unfortunately, our results were complicated by considerable uncertainty in manipulating dosing and HRT at such low flows. Despite our inconclusive results, previous studies have demonstrated that hydrologic management is needed to maximize nutrient retention and removal in vegetated ditches [30,82,83], and we still recommend stacking conservation practices by managing vegetation and HRTs within ditches. Low-grade weirs can increase the depth and duration of wetted areas, time to peak flow during storm events, and nutrient uptake in poorly drained landscapes of the Lower Mississippi River Basin [82,84,85]. Despite this effort, more research needs to be conducted to better understand potential additive benefits of managing both vegetation and HRT in agricultural drainage networks.
Nutrient retention efficiency is likely overwhelmed by high concentrations during stormwater-driven runoff events [82,86], and managing vegetation and hydrology within drainage networks likely only provides nutrient uptake benefits under certain conditions. Therefore, it is also important for watershed managers to couple field-based best management practices such as cover crops, vegetated buffer strips, and increased nutrient use efficiency in agricultural production areas to enhance agroecosystem sustainability through reducing mobilization and runoff of excess nutrients into ditch systems [86,87]. The role of vegetated ditches as important modulators of N export from agricultural watersheds will likely be enhanced through hydrologic and vegetation management practice within ditches, in conjunction with field management practices (e.g., cover crops). More research is needed to understand cumulative benefits of integrated management practices within drainage networks, agricultural areas, and intervening habitats for N and P retention in agroecosystems.