Characterization of preferential flow in undisturbed, structured soil columns using a vertical TDR probe

Rapid movement of agricultural chemicals through soil to groundwater via preferential flow pathways is one cause of water contamination. Previous studies have shown that time domain reflectometry (TDR) could be used to characterize solute transport in soil. However, previous studies have only scarcely addressed preferential flow. This study presents an extended application of TDR for determining preferential flow properties. A TDR method was tested in carefully controlled laboratory experiments using 20-cm long and 12-cm diameter undisturbed, structured soil columns. The method used a vertically installed TDR probe and a short pulse of tracer application to obtain residual mass (RM) breakthrough curves (BTC). The RM BTC obtained from TDR were used to estimate mobile/immobile model (MIM) parameters that were compared to the parameter estimates from effluent data. A conventional inverse curve fitting method (CXTFIT) was used to estimate parameters. The TDR-determined parameters were then used to generate calculated effluent BTC for comparison with observed effluent BTC for the same soil columns. Time moments of the calculated and observed BTC were calculated to quantitatively evaluate the calculated BTC. Overall, the RM BTC obtained from TDR were similar to the RM BTC obtained from effluent data. The TDR-determined parameters corresponded well to the parameters obtained from the effluent data, although they were not within the 95% confidence intervals. Correlation coefficients between the parameters obtained from TDR and from effluent data for the immobile water fraction (theta im/theta), mass exchange coefficient (alpha), and dispersion coefficient (Dm) were 0.95, 0.95, and 0.99, respectively. For three of the four soil cores, theta im/theta ranged from 0.42 to 0.82, indicating considerable preferential flow. The TDR-calculated effluent BTC also were similar to the observed effluent BTC having an average coefficient of determination of 0.94. Time moments obtained from calculated BTC were representative of those obtained from observed BTC. The vertical TDR probe method was simple and minimally destructive and provided representative preferential flow properties that enabled the characterization of solute transport in soil.     

Soil-test N recommendations augmented with PEST-optimized RZWQM simulations

Improved understanding of year-to-year late-spring soil nitrate test (LSNT) variability could help make it more attractive to producers. We test the ability of the Root Zone Water Quality Model (RZWQM) to simulate watershed-scale variability due to the LSNT, and we use the optimized model to simulate long-term field N dynamics under related conditions. Autoregressive techniques and the automatic parameter calibration program PEST were used to show that RZWQM simulates significantly lower nitrate concentration in discharge from LSNT treatments compared with areas receiving fall N fertilizer applications within the tile-drained Walnut Creek, Iowa, watershed (>5 mg NL(-1) difference for the third year of the treatment, 1999). This result is similar to field-measured data from a paired watershed experiment. A statistical model we developed using RZWQM simulations from 1970 to 2005 shows that early-season precipitation and early-season temperature account for 90% of the interannual variation in LSNT-based fertilizer N rates. Long-term simulations with similar average N application rates for corn (Zea mays L.) (151 kg N ha(-1)) show annual average N loss in tile flow of 20.4, 22.2, and 27.3 kg N ha(-1) for LSNT, single spring, and single fall N applications. These results suggest that (i) RZWQM is a promising tool to accurately estimate the water quality effects of LSNT; (ii) the majority of N loss difference between LSNT and fall applications is because more N remains in the root zone for crop uptake; and (iii) year-to-year LSNT-based N rate differences are mainly due to variation in early-season precipitation and temperature.     

Evaluation of nitrate nitrogen fluxes from a tile-drained watershed in central Iowa

Nitrate N fluxes from tile-drained watersheds have been implicated in water quality studies of the Mississippi River basin, but actual NO3-N loads from small watersheds during long periods are poorly documented. We evaluated discharge and NO3-N fluxes passing the outlet of an Iowa watershed (5134 ha) and two of its tile-drained subbasins (493 and 863 ha) from mid-1992 through 2000. The cumulative NO3-N load from the catchment was 168 kg ha(-1), and 176 and 229 kg ha(-1) from the subbasins. The outlet had greater total discharge (1831 mm) and smaller flow-weighted mean NO3-N concentration (9.2 mg L(-1)) than the subbasins, while the larger subbasin had greater discharge (1712 vs. 1559 mm) and mean NO3-N concentration (13.4 vs. 11.3 mg L(-1)) than the smaller subbasin. Concentrations exceeding 10 mg L(-1) were common, but least frequent at the outlet. Nitrate N was generally not diluted by large flows, except during 1993 flooding. The outlet showed smaller NO3-N concentrations at low flows. Relationships between discharge and NO3-N flux showed log-log slopes near 1.0 for the subbasins, and 1.2 for the outlet, considering autocorrelation and measurement-error effects. We estimated denitrification of subbasin NO3-N fluxes in a hypothetical wetland using published data. Assuming that temperature and NO3-N supply could limit denitrification, then about 20% of the NO3-N would have been denitrified by a wetland constructed to meet USDA-approved criteria. The low efficiency results from the seasonal timing and NO3-N content of large flows. Therefore, agricultural and wetland best management practices (BMPs) are needed to achieve water quality goals in tile-drained watersheds.     

Nitrate loss in subsurface drainage as affected by nitrogen fertilizer rate

The relationships between N fertilizer rate, yield, and NO3 leaching need to be quantified to develop soil and crop management practices that are economically and environmentally sustainable. From 1996 through 1999, we measured yield and NO3 loss from a subsurface drained field in central Iowa at three N fertilizer rates: a low (L) rate of 67 kg ha(-1) in 1996 and 57 kg ha(-1) in 1998, a medium (M) rate of 135 kg ha(-1) in 1996 and 114 kg ha(-1) in 1998, and a high (H) rate of 202 kg ha(-1) in 1996 and 172 kg ha(-1) in 1998. Corn (Zea mays L.) and soybean [Glycine max (L.) Merr.] were grown in rotation with N fertilizer applied in the spring to corn only. For the L treatment, NO3 concentrations in the drainage water exceeded the 10 mg N L(-1) maximum contaminant level (MCL) established by the USEPA for drinking water only during the years that corn was grown. For the M and H treatments, NO3 concentrations exceeded the MCL in all years, regardless of crop grown. For all years, the NO3 mass loss in tile drainage water from the H treatment (48 kg N ha(-1)) was significantly greater than the mass losses from the M (35 kg N ha(-1)) and L (29 kg N ha(-1)) treatments, which were not significantly different. The economically optimum N fertilizer rate for corn was between 67 and 135 kg ha(-1) in 1996 and 114 and 172 kg ha(-1) in 1998, but the net N mass balance indicated that N was being mined from the soil at these N fertilizer levels and that the system would not be sustainable.