Illinois River Nitrate-Nitrogen Concentrations and Loads: Long-term Variation and Association with Watershed Nitrogen Inputs

The Illinois River is a major contributor of nitrate-N to the Mississippi River and the Gulf of Mexico, where nitrate is a leading cause of summertime benthic hypoxia. Corn-soybean production on tile-drained land is a leading source of nitrate-N in this river system, in addition to municipal wastewater discharge. We calculated annual nitrate-N loads in the Illinois River at Valley City from 1976 to 2014 by linear interpolation. Although there was not a significant trend in annual loads during the entire study period, there was a significant downward trend in flow-weighted nitrate-N concentration after 1990 despite high concentrations in 2013 after the 2012 drought. Multivariate regression analysis revealed a statistically significant association between annual flow-weighted nitrate-N concentration and cumulative residual agricultural N inputs to the watershed during a 6-yr window. This suggests that declines in flow-weighted nitrate-N concentration may reflect increasing N use efficiency in agriculture and a depletion of legacy N stored in the watershed. The watershed appears to have transitioned from a state of stationarity in nitrate concentration to nonstationarity. The average annual nitrate-N load at Valley City from 2010 to 2014 was 10% less than the 1980-1996 average load, indicating recent progress toward Illinois' nutrient loss reduction milestone of 15% reduction by 2025 and ultimate target of 45% reduction.

Riverine Response of Sulfate to Declining Atmospheric Sulfur Deposition in Agricultural Watersheds

Sulfur received extensive study as an input to terrestrial ecosystems from acidic deposition during the 1980s. With declining S deposition inputs across the eastern United States, there have been many studies evaluating ecosystem response, with the exception of agricultural watersheds. We used long-term (22 and 18 yr) sulfate concentration data from two rivers and recent (6 yr) data from a third river to better understand cycling and transport of S in agricultural, tile-drained watersheds. Sulfate concentrations and yields steadily declined in the Embarras (from ∼10 to 6 mg S L) and Kaskaskia rivers (from 7 to 3.5 mg S L) during the sampling period, with an overall -23.1 and -12.8 kg S ha yr balance for the two watersheds. There was evidence of deep groundwater inputs of sulfate in the Salt Fork watershed, with a much smaller input to the Embarras and none to the Kaskaskia. Tiles in the watersheds had low sulfate concentrations (<10 mg S L), similar to the Kaskaskia River, unless the field had received some form of S fertilizer. A multiple regression model of runoff (cm) and S deposition explained much of the variation in Embarras River sulfate ( = 0.86 and 0.80 for concentrations and yields; = 46). Although atmospheric deposition was much less than outputs (grain harvest + stream export of sulfate), riverine transport of sulfate reflected the decline in inputs. Watershed S balances suggest a small annual depletion of soil organic S pools, and S fertilization will likely be needed at some future date to maintain crop yields.

Temperature and Substrate Control Woodchip Bioreactor Performance in Reducing Tile Nitrate Loads in East-Central Illinois

Tile drainage is the major source of nitrate in the upper Midwest, and end-of-tile removal techniques such as wood chip bioreactors have been installed that allow current farming practices to continue, with nitrate removed through denitrification. There have been few multiyear studies of bioreactors examining controls on nitrate removal rates. We evaluated the nitrate removal performance of two wood chip bioreactors during the first 3 yr of operation and examined the major factors that regulated nitrate removal. Bioreactor 2 was subject to river flooding, and performance was not assessed. Bioreactor 1 had average monthly nitrate removal rates of 23 to 44 g N m d in Year 1, which decreased to 1.2 to 11 g N m d in Years 2 and 3. The greater N removal rates in Year 1 and early in Year 2 were likely due to highly degradable C in the woodchips. Only late in Year 2 and in Year 3 was there a strong temperature response in the nitrate removal rate. Less than 1% of the nitrate removed was emitted as NO. Due to large tile inputs of nitrate (729-2127 kg N) at high concentrations (∼30 mg nitrate N L) in Years 2 and 3, overall removal efficiency was low (3 and 7% in Years 2 and 3, respectively). Based on a process-based bioreactor performance model, Bioreactor 1 would have needed to be 9 times as large as the current system to remove 50% of the nitrate load from this 20-ha field.

Denitrifying Bioreactors for Nitrate Removal: A Meta-Analysis

Meta-analysis approaches were used in this first quantitative synthesis of denitrifying woodchip bioreactors. Nitrate removal across environmental and design conditions was assessed from 26 published studies, representing 57 separate bioreactor units (i.e., walls, beds, and laboratory columns). Effect size calculations weighted the data based on variance and number of measurements for each bioreactor unit. Nitrate removal rates in bed and column studies were not significantly different, but both were significantly higher than wall studies. In denitrifying beds, wood source did not significantly affect nitrate removal rates. Nitrate removal (mass per volume) was significantly lower in beds with <6-h hydraulic retention times, which argues for ensuring that bed designs incorporate sufficient time for nitrate removal. Rates significantly declined after the first year of bed operation but then stabilized. Nitrogen limitation significantly affected bed nitrate removal. Categorical and linear assessments found significant nitrate removal effects with bed temperature; a of 2.15 was quite similar to other studies. Lessons from this meta-analysis can be incorporated into bed designs, especially extending hydraulic retention times to increase nitrate removal under low temperature and high flow conditions. Additional column studies are warranted for comparative assessments, as are field-based studies for assessing in situ conditions, especially in aging beds, with careful collection and reporting of design and environmental data. Future assessment of these systems might take a holistic view, reviewing nitrate removal in conjunction with other processes, including greenhouse gas and other unfavorable by-product production.

Chloride Sources and Losses in Two Tile-Drained Agricultural Watersheds

Chloride is a relatively unreactive plant nutrient that has long been used as a biogeochemical tracer but also can be a pollutant causing aquatic biology impacts when concentrations are high, typically from rock salt applications used for deicing roads. Chloride inputs to watersheds are most often from atmospheric deposition, road salt, or agricultural fertilizer, although studies on agricultural watersheds with large fertilizer inputs are few. We used long-term (21 and 17 yr) chloride water quality data in two rivers of east-central Illinois to better understand chloride biogeochemistry in two agricultural watersheds (Embarras and Kaskaskia), the former with a larger urban land use and both with extensive tile drainage. During our sampling period, the average chloride concentration was 23.7 and 20.9 mg L in the Embarras and Kaskaskia Rivers, respectively. Annual fluxes of chloride were 72.5 and 61.2 kg ha yr in the Embarras and Kaskaskia watersheds, respectively. In both watersheds, fertilizer chloride was the dominant input (∼49 kg ha yr), with road salt likely the other major source (23.2 and 7.2 kg ha yr for the Embarras and Kaskaskia watersheds, respectively). Combining our monitoring data with earlier published data on the Embarras River showed an increase in chloride concentrations as potash use increased in Illinois during the 1960s and 1970s with a lag of about 2 to 6 yr to changes in potash inputs based on a multiple-regression model. In these agricultural watersheds, riverine chloride responds relatively quickly to potash fertilization as a result of tile-drainage.

Characterizing the Performance of Denitrifying Bioreactors during Simulated Subsurface Drainage Events

The need to mitigate nitrate export from corn and soybean fields with subsurface (tile) drainage systems, a major environmental issue in the midwestern United States, has made the efficacy of field-edge, subsurface bioreactors an active subject of research. This study of three such bioreactors located on the University of Illinois South Farms during their first 6 mo of operation (July-Dec. 2012) focused on the interactions of seasonal temperature changes and hydraulic retention times (HRTs), which were subject to experimental manipulation. Changes in nitrate, phosphate, oxygen, and dissolved organic carbon were monitored in influent and effluent to assess the benefits and the potential harmful effects of bioreactors for nearby aquatic ecosystems. On average, bioreactors reduced nitrate loads by 63%, with minimum and maximum reductions of 20 and 98% at low and high HRTs, respectively. The removal rate per unit reactor volume averaged 11.6 g NO-N m d (range, 5-30 g NO-N m d). Multiple regression models with exponential dependencies on influent water temperature and on HRT explained 73% of the variance in NO-N load reduction and 43% of the variance in its removal rate. Although concentrations of dissolved reactive phosphorus and dissolved organic carbon in the bioreactor effluent increased relative to the influent by an order of magnitude during initial tests, within 1 mo of operation they stabilized at nearly equal values.

Nitrogen removal and greenhouse gas emissions from constructed wetlands receiving tile drainage water

Loss of nitrate from agricultural lands to surface waters is an important issue, especially in areas that are extensively tile drained. To reduce these losses, a wide range of in-field and edge-of-field practices have been proposed, including constructed wetlands. We re-evaluated constructed wetlands established in 1994 that were previously studied for their effectiveness in removing nitrate from tile drainage water. Along with this re-evaluation, we measured the production and flux of greenhouse gases (GHGs) (CO, NO, and CH). The tile inlets and outlets of two wetlands were monitored for flow and N during the 2012 and 2013 water years. In addition, seepage rates of water and nitrate under the berm and through the riparian buffer strip were measured. Greenhouse gas emissions from the wetlands were measured using floating chambers (inundated fluxes) or static chambers (terrestrial fluxes). During this 2-yr study, the wetlands removed 56% of the total inlet nitrate load, likely through denitrification in the wetland. Some additional removal of nitrate occurred in seepage water by the riparian buffer strip along each berm (6.1% of the total inlet load, for a total nitrate removal of 62%). The dominant GHG emitted from the wetlands was CO, which represented 75 and 96% of the total GHG emissions during the two water years. The flux of NO contributed between 3.7 and 13% of the total cumulative GHG flux. Emissions of NO were 3.2 and 1.3% of the total nitrate removed from wetlands A and B, respectively. These wetlands continue to remove nitrate at rates similar to those measured after construction, with relatively little GHG gas loss.

Navigating the socio-bio-geo-chemistry and engineering of nitrogen management in two illinois tile-drained watersheds

Reducing nitrate loads from corn and soybean, tile-drained, agricultural production systems in the Upper Mississippi River basin is a major challenge that has not been met. We evaluated a range of possible management practices from biophysical and social science perspectives that could reduce nitrate losses from tile-drained fields in the Upper Salt Fork and Embarras River watersheds of east-central Illinois. Long-term water quality monitoring on these watersheds showed that nitrate losses averaged 30.6 and 23.0 kg nitrate N ha yr (Embarras and Upper Salt Fork watersheds, respectively), with maximum nitrate concentrations between 14 and 18 mg N L. With a series of on-farm studies, we conducted tile monitoring to evaluate several possible nitrate reduction conservation practices. Fertilizer timing and cover crops reduced nitrate losses (30% reduction in a year with large nitrate losses), whereas drainage water management on one tile system demonstrated the problems with possible retrofit designs (water flowed laterally from the drainage water management tile to the free drainage system nearby). Tile woodchip bioreactors had good nitrate removal in 2012 (80% nitrate reduction), and wetlands had previously been shown to remove nitrate (45% reductions) in the Embarras watershed. Interviews and surveys indicated strong environmental concern and stewardship ethics among landowners and farmers, but the many financial and operational constraints that they operate under limited their willingness to adopt conservation practices that targeted nitrate reduction. Under the policy and production systems currently in place, large-scale reductions in nitrate losses from watersheds such as these in east-central Illinois will be difficult.

Variation in riverine nitrate flux and fall nitrogen fertilizer application in East-central illinois

In east-central Illinois, fertilizer sales during the past 20 yr suggest that approximately half of the fertilizer nitrogen (N) applied to corn ( L.) occurs in the fall; however, fall fertilizer N sales were greatly reduced in 2009 as wet soil conditions restricted fall fieldwork, including fertilizer N applications. In 2010, we observed unusually low flow-weighted nitrate concentrations (approximately 40% below the long-term average) in two east-central Illinois rivers (5.7 mg N L in the Embarras River and 5.6 mg N L in the Lake Fork of the Kaskaskia River). Using long-term river nitrate data sets (1993-2012 for the Embarras and 1997-2012 for the Kaskaskia), we examined nitrate concentrations and developed regression models to estimate the association between fall fertilizer N application on riverine nitrate yields in these tile-drained watersheds. During these periods of record, annual riverine nitrate yields ranged from 8 to 57 kg N ha yr (30 kg N ha yr average) for the Embarras River and 2.6 to 59 kg N ha yr (32 kg N ha yr average) for the Kaskaskia. Multivariate linear regression relationships with the current and previous year's annual water yields, previous year's corn yield, and nine-county fall fertilizer sales accounted for 96% of the annual variation in nitrate yield in both watersheds. Running the regression models with fall fertilizer sales set to the 2009 amount suggests that the average reduction in nitrate yield (for the period of record) would be 17 and 20% for the Embarras and Kaskaskia Rivers, respectively. These data suggest that shifting fertilizer N application to the spring can be detected in watersheds as large as 481 km.

Role of arthropod communities in bioenergy crop litter decomposition†

The extensive land use conversion expected to occur to meet demands for bioenergy feedstock production will likely have widespread impacts on agroecosystem biodiversity and ecosystem services, including carbon sequestration. Although arthropod detritivores are known to contribute to litter decomposition and thus energy flow and nutrient cycling in many plant communities, their importance in bioenergy feedstock communities has not yet been assessed. We undertook an experimental study quantifying rates of litter mass loss and nutrient cycling in the presence and absence of these organisms in three bioenergy feedstock crops-miscanthus (Miscanthus x giganteus), switchgrass (Panicum virgatum), and a planted prairie community. Overall arthropod abundance and litter decomposition rates were similar in all three communities. Despite effective reduction of arthropods in experimental plots via insecticide application, litter decomposition rates, inorganic nitrogen leaching, and carbon-nitrogen ratios did not differ significantly between control (with arthropods) and treatment (without arthropods) plots in any of the three community types. Our findings suggest that changes in arthropod faunal composition associated with widespread adoption of bioenergy feedstock crops may not be associated with profoundly altered arthropod-mediated litter decomposition and nutrient release.

Reduced nitrogen losses after conversion of row crop agriculture to perennial biofuel crops

Current biofuel feedstock crops such as corn lead to large environmental losses of N through nitrate leaching and NO emissions; second-generation cellulosic crops have the potential to reduce these N losses. We measured N losses and cycling in establishing miscanthus (), switchgrass ( L. fertilized with 56 kg N ha yr), and mixed prairie, along with a corn ( L.)-corn-soybean [ (L.) Merr.] rotation (corn fertilized at 168-202 kg N ha). Nitrous oxide emissions, soil N mineralization, mid-profile nitrate leaching, and tile flow and nitrate concentrations were measured. Perennial crops quickly reduced nitrate leaching at a 50-cm soil depth as well as concentrations and loads from the tile systems (year 1 tile nitrate concentrations of 10-15 mg N L declined significantly by year 4 in all perennial crops to <0.6 mg N L, with losses of <0.8 kg N ha yr). Nitrous oxide emissions were 2.2 to 7.7 kg N ha yr in the corn-corn-soybean rotation but were <1.0 kg N ha yr by year 4 in the perennial crops. Overall N balances (atmospheric deposition + fertilization + soybean N fixation - harvest, leaching losses, and NO emissions) were positive for corn and soybean (22 kg N ha yr) as well as switchgrass (9.7 kg N ha yr) but were -18 and -29 kg N ha yr for prairie and miscanthus, respectively. Our results demonstrate rapid tightening of the N cycle as perennial biofuel crops established on a rich Mollisol soil.

A spatial analysis of phosphorus in the Mississippi river basin

Phosphorus (P) in rivers in the Mississippi River basin (MRB) contributes to hypoxia in the Gulf of Mexico and impairs local water quality. We analyzed the spatial pattern of P in the MRB to determine the counties with the greatest January to June P riverine yields and the most critical factors related to this P loss. Using a database of P inputs and landscape characteristics from 1997 through 2006 for each county in the MRB, we created regression models relating riverine total P (TP), dissolved reactive P (DRP), and particulate P (PP) yields for watersheds within the MRB to these factors. Riverine yields of P were estimated from the average concentration of each form of P during January to June for the 10-yr period, multiplied by the average daily flow, and then summed for the 6-mo period. The fraction of land planted in crops, human consumption of P, and precipitation were found to best predict TP yields with a spatial error regression model ( = 0.48, = 101). Dissolved reactive P yields were predicted by fertilizer P inputs, human consumption of P, and precipitation in a multiple regression model ( = 0.42, = 73), whereas PP yields were explained by crop fraction, human consumption of P, and soil bulk density in a spatial error regression model ( = 0.49, = 61). Overall, the Upper Midwest's Cornbelt region and lower Mississippi basin had the counties with the greatest P yields. These results help to point out specific areas where agricultural conservation practices that reduce losses to streams and rivers and point source P removal might limit the intensity or spatial occurrence of Gulf of Mexico hypoxia and improve local water quality.

Sources of nitrate yields in the Mississippi River Basin

Riverine nitrate N in the Mississippi River leads to hypoxia in the Gulf of Mexico. Several recent modeling studies estimated major N inputs and suggested source areas that could be targeted for conservation programs. We conducted a similar analysis with more recent and extensive data that demonstrates the importance of hydrology in controlling the percentage of net N inputs (NNI) exported by rivers. The average fraction of annual riverine nitrate N export/NNI ranged from 0.05 for the lower Mississippi subbasin to 0.3 for the upper Mississippi River basin and as high as 1.4 (4.2 in a wet year) for the Embarras River watershed, a mostly tile-drained basin. Intensive corn (Zea mays L.) and soybean [Glycine max (L.) Merr.] watersheds on Mollisols had low NNI values and when combined with riverine N losses suggest a net depletion of soil organic N. We used county-level data to develop a nonlinear model ofN inputs and landscape factors that were related to winter-spring riverine nitrate yields for 153 watersheds within the basin. We found that river runoff times fertilizer N input was the major predictive term, explaining 76% of the variation in the model. Fertilizer inputs were highly correlated with fraction of land area in row crops. Tile drainage explained 17% of the spatial variation in winter-spring nitrate yield, whereas human consumption of N (i.e., sewage effluent) accounted for 7%. Net N inputs were not a good predictor of riverine nitrate N yields, nor were other N balances. We used this model to predict the expected nitrate N yield from each county in the Mississippi River basin; the greatest nitrate N yields corresponded to the highly productive, tile-drained cornbelt from southwest Minnesota across Iowa, Illinois, Indiana, and Ohio. This analysis can be used to guide decisions about where efforts to reduce nitrate N losses can be most effectively targeted to improve local water quality and reduce export to the Gulf of Mexico.

Miscanthus and switchgrass production in central Illinois: impacts on hydrology and inorganic nitrogen leaching

Biomass crops are being promoted as environmentally favorable alternatives to fossil fuels or ethanol production from maize (Zea mays L.), particularly across the Corn Belt of the United States. However, there are few if any empirical studies on inorganic N leaching losses from perennial grasses that are harvested on an annual basis, nor has there been empirical evaluation of the hydrologic consequences of perennial cropping systems. Here we report on the results of 4 yr of field measurements of soil moisture and inorganic N leaching from a conventional maize-soybean [Glycine max (L.) Merr.] system and two unfertilized perennial grasses harvested in winter for biomass: Miscanthus x giganteus and switchgrass (Panicum virgatum cv. Cave-in-Rock). All crops were grown on fertile Mollisols in east-central Illinois. Inorganic N leaching was measured with ion exchange resin lysimeters placed 50 cm below the soil surface. Maize--soybean nitrate leaching averaged 40.4 kg N ha(-1) yr(-1), whereas switchgrass and Miscanthus had values of 1.4 and 3.0 kg N ha(-1) yr(-1), respectively. Soil moisture monitoring (to a depth of 90 cm) indicated that both perennial grasses dried the soil out earlier in the growing season compared with maize-soybean. Later in the growing season, soil moisture under switchgrass tended to be greater than maize-soybean or Miscanthus, whereas the soil under Miscanthus was consistently drier than under maize--soybean. Water budget calculations indicated that evapotranspiration from Miscanthus was about 104 mm yr(-1) greater than under maize-soybean, which could reduce annual drainage water flows by 32% in central Illinois. Drainage water is a primary source of surface water flows in the region, and the impact ofextensive Miscanthus production on surface water supplies and aquatic ecosystems deserves further investigation.

Nitrogen mass balance of a tile-drained agricultural watershed in East-Central Illinois

Simple nitrogen (N) input/output balance calculations in agricultural systems are used to evaluate performance of nutrient management; however, they generally rely on extensive assumptions that do not consider leaching, denitrification, or annual depletion of soil N. We constructed a relatively complete N mass balance for the Big Ditch watershed, an extensively tile-drained agricultural watershed in east-central Illinois. We conducted direct measurements of a wide range of N pools and fluxes for a 2-yr period, including soil N mineralization, soybean N(2) fixation, tile and river N loads, and ground water and in-stream denitrification. Fertilizer N inputs were from a survey of the watershed and yield data from county estimates that were combined with estimated protein contents to obtain grain N. By using maize fertilizer recovery and soybean N(2) fixation to estimate total grain N derived from soil, we calculated the explicit change in soil N storage each year. Overall, fertilizer N and soybean N(2) fixation dominated inputs, and total grain export dominated outputs. Precipitation during 2001 was below average (78 cm), whereas precipitation in 2002 exceeded the 30-yr average of 97 cm; monthly rainfall was above average in April, May, and June of 2002, which flooded fields and produced large tile and riverine N loads. In 2001, watershed inputs were greater than outputs, suggesting that carryover of N to the subsequent year may occur. In 2002, total inputs were less than outputs due to large leaching losses and likely substantial field denitrification. The explicit change in soil storage (67 kg N ha(-1)) offsets this balance shortfall. Although 2002 was climatically unusual, with current production trends of greater maize grain yields with less fertilizer N, soil N depletion is likely to occur in maize/soybean rotations, especially in years with above-average precipitation or extremely wet spring periods.

Relationships between benthic sediments and water column phosphorus in Illinois streams

Sediments can be important in regulating stream water P concentrations, and this has implications for establishing nutrient standards that have not been fully investigated. We evaluated abiotic and biotic processes to better understand the role of sediments in determining stream water dissolved P concentrations. Sediment and stream water samples were collected during low discharge from 105 streams across Illinois and analyzed for equilibrium P concentration at zero release or retention (EPC(0)), P sorption characteristics, stream water P concentration, and sediment particle size. In addition, four east-central Illinois streams were repeatedly sampled to examine temporal patterns in sediment P retention and biotic processing of P. Median dissolved reactive P (DRP) and total P concentrations across the state were 0.081 and 0.168 mg L(-1), respectively. Sediment EPC(0) concentrations were related to stream water DRP concentrations (r(s) = 0.75). Sediment silt+clay (and co-correlated organic matter) was related to sorbed P (r(s) = -0.49) and the reactive sediment pool of P (r(s) = 0.76). However, for most sites this pool was small given the coarse textures present (median silt+clay was 5.7%). Repeated sampling at the four intensive sites showed little variation in EPC(0) values or alkaline phosphatase activity, suggesting overall stream conditions regulated the biotic processing. Biotic retention of P was 32% of short-term P removal. We conclude that sediments in Illinois streams are a reflection of and partially affected by stream water P concentrations through both abiotic and biotic processes. Sediments seem unlikely to alter annual stream P loads, but may affect concentrations at low discharge.

Long-term changes in mollisol organic carbon and nitrogen

Conversions of Mollisols from prairie to cropland and subsequent changes in crop production practices in the Midwestern USA have resulted in changes in soil organic matter. Few studies have used archived samples, long-term resampling of soils to a depth of 1 m, and space for time studies to document these changes. We resampled soils by depth (0-100 cm) in fields at 19 locations in central Illinois on poorly drained Mollisols that were in corn (Zea mays L.) and soybean (Glycine max L. Merr.) rotations, were tile drained, and had no known history of manure application in recent decades. Three fields were paired with virgin prairie remnants, two had grass borders that were sampled, and 16 had been previously sampled in 1901 to 1904 or 1957 under various land uses (virgin prairie, cultivation, grass cover). The soils had large amounts of C and N in the profile, with mean values of 175 [corrected] Mg C ha(-1) and 16.1 Mg N ha(-1) for the 18 cultivated fields sampled in 2001 and 2002. We confirmed a large reduction in organic C and total N pools from conversion of prairies to annual cultivation and artificial drainage and documented no change in these organic matter pools of cultivated soils during the period of synthetic fertilizer use (1957--2002). Cultivated fields had soil C and N concentrations typically 30 to 50% less than virgin prairie soils. Smaller but significant declines in C and N concentrations were found when comparing 1900s cultivated fields to concentrations in 2002, after another 100 yr of cultivation, and in comparing 1957 grass covered fields that had been converted to annual cultivation before 2002. The reduction in organic matter after cultivation of prairies occurred mostly in the top 50 cm of soil, with evidence of translocation of C and N from these upper layers to the 50- to 100-cm depth, possibly enhanced by tile drainage. For these Mollisols, declines in organic matter were likely completed by the 1950s, with organic matter pools in a steady state under the production practices in place from the late 1950s through 2002.

Assessment of chlorophyll-a as a criterion for establishing nutrient standards in the streams and rivers of Illinois

Nutrient enrichment is a frequently cited cause for biotic impairment of streams and rivers in the USA. Efforts are underway to develop nutrient standards in many states, but defensible nutrient standards require an empirical relationship between nitrogen (N) or phosphorus (P) concentrations and some criterion that relates nutrient levels to the attainment of designated uses. Algal biomass, measured as chlorophyll-a (chl-a), is a commonly proposed criterion, yet nutrient-chl-a relationships have not been well documented in Illinois at a state-wide scale. We used state-wide surveys of >100 stream and river sites to assess the applicability of chl-a as a criterion for establishing nutrient standards for Illinois. Among all sites, the median total P and total N concentrations were 0.185 and 5.6 mg L(-1), respectively, during high-discharge conditions. During low-discharge conditions, median total P concentration was 0.168 mg L(-1), with 25% of sites having a total P of > or =0.326 mg L(-1). Across the state, 90% of the sites had sestonic chl-a values of < or =35 microg L(-1), and watershed area was the best predictor of sestonic chl-a. During low discharge there was a significant correlation between sestonic chl-a and total P for those sites that had canopy cover < or =25% and total P of < or =0.2 mg L(-1). Results suggest sestonic chl-a may be an appropriate criterion for the larger rivers in Illinois but is inappropriate for small rivers and streams. Coarse substrate to support benthic chl-a occurred in <50% of the sites we examined; a study using artificial substrates did not reveal a relationship between chl-a accrual and N or P concentrations. For many streams and rivers in Illinois, nutrients may not be the limiting factor for algal biomass due to the generally high nutrient concentrations and the effects of other factors, such as substrate conditions and turbidity.

Relationships between water quality, habitat quality, and macroinvertebrate assemblages in Illinois streams

The influence of specific stressors, such as nutrient enrichment and physical habitat degradation, on biotic integrity requires further attention in Midwestern streams. We sampled 53 streams throughout Illinois and examined relationships between macroinvertebrate community structure and numerous physicochemical parameters. Streams were clustered into four major groups based on taxa dissimilarity. Habitat quality and dissolved nutrients were responsible for separating the major groups in a nonmetric multidimensional scaling ordination. Furthermore, the alignment of environmental factors in the ordination suggested there was a habitat quality-nutrient concentration gradient such that streams with high-quality habitats usually had low concentrations of nutrients. Discrimination by community measures further validated the major stream groups and indicated that forested streams had generally higher biological integrity than agricultural streams, which in turn had greater integrity than urban streams. Our results demonstrate that physical habitat degradation and nutrient pollution are important and often confounded determinants of biotic integrity in Illinois streams. In addition, we suggest that management of Midwestern streams could benefit from further implementation of multivariate data exploration and stream classification techniques.

Modeling riverine nitrate export from an East-Central Illinois watershed using SWAT

Reliable water quality models are needed to forecast the water quality consequences of different agricultural nutrient management scenarios. In this study, the Soil and Water Assessment Tool (SWAT), version 2000, was applied to simulate streamflow, riverine nitrate (NO(3)) export, crop yield, and watershed nitrogen (N) budgets in the upper Embarras River (UER) watershed in east-central Illinois, which has extensive maize-soybean cultivation, large N fertilizer input, and extensive tile drainage. During the calibration (1994-2002) and validation (1985-1993) periods, SWAT simulated monthly and annual stream flows with Nash-Sutcliffe coefficients (E) ranging from 0.67 to 0.94 and R(2) from 0.75 to 0.95. For monthly and annual NO(3) loads, E ranged from -0.16 to 0.45 and R(2) from 0.36 to 0.74. Annual maize and soybean yields were simulated with relative errors ranging from -10 to 6%. The model was then used to predict the changes in NO(3) output with N fertilizer application rates 10 to 50% lower than original application rates in UER. The calibrated SWAT predicted a 10 to 43% decrease in NO(3) export from UER and a 6 to 38% reduction in maize yield in response to the reduction in N fertilizer. The SWAT model markedly overestimated NO(3) export during major wet periods. Moreover, SWAT estimated soybean N fixation rates considerably greater than literature values, and some simulated changes in the N cycle in response to fertilizer reduction seemed to be unrealistic. Improving these aspects of SWAT could lead to more reliable predictions in the water quality outcomes of nutrient management practices in tile-drained watersheds.

Phosphorus transport pathways to streams in tile-drained agricultural watersheds

Agriculture is a major nonpoint source of phosphorus (P) in the Midwest, but how surface runoff and tile drainage interact to affect temporal concentrations and fluxes of both dissolved and particulate P remains unclear. Our objective was to determine the dominant form of P in streams (dissolved or particulate) and identify the mode of transport of this P from fields to streams in tile-drained agricultural watersheds. We measured dissolved reactive P (DRP) and total P (TP) concentrations and loads in stream and tile water in the upper reaches of three watersheds in east-central Illinois (Embarras River, Lake Fork of the Kaskaskia River, and Big Ditch of the Sangamon River). For all 16 water year by watershed combinations examined, annual flow-weighted mean TP concentrations were >0.1 mg L(-1), and seven water year by watershed combinations exceeded 0.2 mg L(-1). Concentrations of DRP and particulate P (PP) increased with stream discharge; however, particulate P was the dominant form during overland runoff events, which greatly affected annual TP loads. Concentrations of DRP and PP in tiles increased with discharge, indicating tiles were a source of P to streams. Across watersheds, the greatest DRP concentrations (as high as 1.25 mg L(-1)) were associated with a precipitation event that followed widespread application of P fertilizer on frozen soils. Although eliminating this practice would reduce the potential for overland runoff of P, soil erosion and tile drainage would continue to be important transport pathways of P to streams in east-central Illinois.

Denitrification and the nitrogen budget of a reservoir in an agricultural landscape

Denitrification is an important process in aquatic sediments, but its role has not been assessed in the N mass balance of upper-Midwestern (USA) reservoirs that receive large agricultural riverine N inputs. We used a 4400-ha reservoir to determine the role of denitrification in the N mass balance and effectiveness in reducing downstream transport of NO(3-)N. Sediment denitrification was (1) measured monthly (March 2002-March 2003) at eight sites in the Lake Shelbyville reservoir in central Illinois using the acetylene inhibition, chloramphenicol technique, (2) scaled to the overall reservoir and compared to N not accounted for in a mass balance, and (3) estimated indirectly using long-term (1981-2003) mass balances of N in the reservoir. Denitrification rates in the reservoir were high during spring and early summer of 2002, when maximum NO(3-)N concentrations were measured (10-14 mg NO(3-)N/L). We estimated that denitrification for the year was between 2580 and 5150 Mg N. Missing N from the mass balance was 3004 Mg N, suggesting that sediment denitrification was the sink. Areal rates of sediment denitrification in the reservoir ranged from 62 to 225 g N x m(-2) x yr(-1), with rates a function of both denitrification intensity (microg N x g dry mass x h(-1)) and the overall mass of sediment present. From 1981 to 2003 the average NO(3-)N inlet flux was 8900 Mg N/yr. About 58% of the total NO(3-)N input was removed, and annual NO(3-)N removed as a percentage of inputs was significantly related to reservoir retention time (average = 0.36 yr for the 23 years, range = 0.21-0.84 yr). By scaling denitrification in Lake Shelbyville to other reservoirs in Illinois, we estimated a sink of 48900 Mg N/yr. When combined with estimated in-stream denitrification, 60900 Mg N/yr was estimated to be removed by sediment denitrification. This reduces riverine export from Illinois to the Gulf of Mexico, where the flux during the 1990s was about 244000 Mg N/yr, and illustrates the importance of reservoir denitrification as an N sink in Midwestern agricultural landscapes.

Methods for measuring denitrification: diverse approaches to a difficult problem

Denitrification, the reduction of the nitrogen (N) oxides, nitrate (NO3-) and nitrite (NO2-), to the gases nitric oxide (NO), nitrous oxide (N2O), and dinitrogen (N2), is important to primary production, water quality, and the chemistry and physics of the atmosphere at ecosystem, landscape, regional, and global scales. Unfortunately, this process is very difficult to measure, and existing methods are problematic for different reasons in different places at different times. In this paper, we review the major approaches that have been taken to measure denitrification in terrestrial and aquatic environments and discuss the strengths, weaknesses, and future prospects for the different methods. Methodological approaches covered include (1) acetylene-based methods, (2) 15N tracers, (3) direct N2 quantification, (4) N2:Ar ratio quantification, (5) mass balance approaches, (6) stoichiometric approaches, (7) methods based on stable isotopes, (8) in situ gradients with atmospheric environmental tracers, and (9) molecular approaches. Our review makes it clear that the prospects for improved quantification of denitrification vary greatly in different environments and at different scales. While current methodology allows for the production of accurate estimates of denitrification at scales relevant to water and air quality and ecosystem fertility questions in some systems (e.g., aquatic sediments, well-defined aquifers), methodology for other systems, especially upland terrestrial areas, still needs development. Comparison of mass balance and stoichiometric approaches that constrain estimates of denitrification at large scales with point measurements (made using multiple methods), in multiple systems, is likely to propel more improvement in denitrification methods over the next few years.

Influence of geomorphological variability in channel characteristics on sediment denitrification in agricultural streams

Within fluvial systems, the spatial variability of geomorphological characteristics of stream channels and associated streambed properties can affect many biogeochemical processes. In agricultural streams of the midwestern USA, it is not known how geomorphological variability affects sediment denitrification rates, a potentially important loss mechanism for N. Sediment denitrification was measured at channelized and meandering headwater reaches in east-central Illinois, a region dominated by intensive agriculture and high NO(3)-N stream export, between June 2003 and February 2005 using the chloramphenicol-amended acetylene inhibition procedure. Sediment denitrification rates were greatest in separation zones, ranging from 0.6 to 76.4 mg N m(-2) h(-1), compared with riffles, point bars, pools, and a run ranging from 0 to 36.5 mg N m(-2) h(-1). Differences in benthic organic matter (r = 0.70) and the percentage of fine-grained sediments (r = 0.93) in the streambeds controlled much of the spatial variations in sediment denitrification among the geomorphological features. Although two meandering study reaches removed 390 and 99% more NO(3)-N by sediment denitrification than adjacent channelized reaches, NO(3)-N loss rates from all reaches were between 0.1 and 15.7% d(-1), except in late summer. Regardless of geomorphological characteristics, streams in east-central Illinois were not able to process the high NO(3)-N loads, making sediment denitrification in this region a limited sink for N.

Relationships among nutrients, chlorophyll-a, and dissolved oxygen in agricultural streams in Illinois

A better understanding of the controls on algae and dissolved O2 in agricultural streams of Illinois is needed to aid in development of nutrient standards. We investigated the relationships between dissolved nutrients, algal abundance, and dissolved O2 in five streams in east-central Illinois from March through November 2004. The streams drained watersheds from 25 to 777 km2 that were dominated by row crop agriculture. Three sites had open canopies and two were bordered by a narrow forest of deciduous trees. Algal abundance was measured as chlorophyll-a (chl-a) concentration in the water column (sestonic) and on the streambed (periphytic). Mean NO3-N concentrations ranged from 5.5 to 8.8 mg N L(-1) and did not relate to algal abundance. Sestonic chl-a values ranged from nearly zero to >15 mg m(-3) with no differences between open and shaded streams and only a weak correlation with dissolved reactive P (mean concentrations were 44-479 microg L(-1)). The results suggest that sestonic chl-a is a poor criterion for assessing nutrient-related problems in these streams. Greatest periphytic chl-a occurred during low flow from August through October, but periphyton occurred consistently in only two of the five streams. The abundance of filamentous algae explained 64% of the variation in diel O2 saturation, but was not correlated with nutrients. Currently it appears that hydrology and light, rather than nutrients, control algal abundance in these streams, and in the agricultural landscape of east-central Illinois, it may not be possible to reduce nutrient concentrations sufficiently to limit filamentous algal blooms.

Timing of riverine export of nitrate and phosphorus from agricultural watersheds in Illinois: implications for reducing nutrient loading to the Mississippi River

Agricultural watersheds in the upper Midwest are the major source of nutrients to the Mississippi River and Gulf of Mexico, but temporal patterns in nutrient export and the role of hydrology in controlling export remain unclear. Here we reporton NO3(-)-N, dissolved reactive phosphorus (DRP), and total P export from three watersheds in Illinois during the past 8-12 years. Our program of intensive, long-term monitoring allowed us to assess how nutrient export was distributed across the range of discharge that occurred at each site and to examine mechanistic differences between NO3(-)-N and DRP export from the watersheds. Last, we used simple simulations to evaluate how nutrient load reductions might affect NO3(-)-N and P export to the Mississippi River from the Illinois watersheds. Artificial drainage through under-field tiles was the primary mechanism for NO3(-)-N export from the watersheds. Tile drainage and overland flow contributed to DRP export, whereas export of particulate P was almost exclusively from overland flow. The analyses revealed that nearly all nutrient export occurred when discharge was > or = median discharge, and extreme discharges (> or = 90th percentile) were responsible for >50% of the NO3(-)-N export and >80% of the P export. Additionally, the export occurred annually during a period beginning in mid-January and continuing through June. These patterns characterized all sites, which spanned a 4-fold range in watershed area. The simulations showed that reducing in-stream nutrient loads by as much as 50% during periods of low discharge would not affect annual nutrient export from the watersheds.

Transport and Fate of Nitrate in Headwater Agricultural Streams in Illinois

Nitrogen inputs to the Gulf of Mexico have increased during recent decades and agricultural regions in the upper Midwest, such as those in Illinois, are a major source of N to the Mississippi River. How strongly denitrification affects the transport of nitrate (NO(3)-N) in Illinois streams has not been directly assessed. We used the nutrient spiraling model to assess the role of in-stream denitrification in affecting the concentration and downstream transport of NO(3)-N in five headwater streams in agricultural areas of east-central Illinois. Denitrification in stream sediments was measured approximately monthly from April 2001 through January 2002. Denitrification rates tended to be high (up to 15 mg N m(-2) h(-1)), but the concentration of NO(3)-N in the streams was also high (>7 mg N L(-1)). Uptake velocities for NO(3)-N (uptake rate/concentration) were lower than reported for undisturbed streams, indicating that denitrification was not an efficient N sink relative to the concentration of NO(3)-N in the water column. Denitrification uptake lengths (the average distance NO(3)-N travels before being denitrified) were long and indicated that denitrification in the streambed did not affect the transport of NO(3)-N. Loss rates for NO(3)-N in the streams were <5% d(-1) except during periods of low discharge and low NO(3)-N concentration, which occurred only in late summer and early autumn. Annually, most NO(3)-N in these headwater sites appeared to be exported to downstream water bodies rather than denitrified, suggesting previous estimates of N losses through in-stream denitrification may have been overestimated.

Stream transport of herbicides and metabolites in a tile-drained, agricultural watershed

The occurrence of metabolites of many commonly used herbicides in streams has not been studied extensively in tile-drained watersheds. We collected water samples throughout the Upper Embarras River watershed [92% corn, Zea mays L., and soybean, Glycine max (L.) Merr.] in east-central Illinois from March 1999 through September 2000 to study the occurrence of atrazine (2-chloro-4-ethylamino-6-isopropylamino-s-triazine), metolachlor 12-chloro-N-(2-ethyl-6-methylphenyl)-N-(methoxy-1-methylethyl) acetamide], alachlor [2-chloro-N-(2,6-diethylphenyl)-N-(methoxymethyl) acetamide], acetochlor [2-chloro-N-(ethoxymethyl)-N-(2-ethyl-6-methylphenyl) acetamide], and their metabolites. River water samples were collected from three subwatersheds of varying tile density (2.8-5.3 km tile km(-2)) and from the outlet (United States Geological Survey [USGS] gage site). Near-record-low totals for stream flow occurred during the study, and nearly all flow was from tiles. Concentrations of atrazine at the USGS gage site peaked at 15 and 17 microg L(-1) in 1999 and 2000, respectively, and metolachlor at 2.7 and 3.2 microg L(-1); this was during the first significant flow event following herbicide applications. Metabolites of the chloroacetanilide herbicides were detected more often than the parent compounds (evaluated during May to July each year, when tiles were flowing), with metolachlor ethanesulfonic acid [2-[(2-ethyl-6-methylphenyl)(2-methoxy-1-methylethyl)amino]-2-oxoethanesulfonic acid] detected most often (> 90% from all sites), and metolachlor oxanilic acid [2-[(2-ethyl-6-methylphenyl)(2-methoxy-1-methylethyl)amino]-2-oxoacetic acid] second (40-100% of samples at the four sites). When summed, the median concentration of the three chloroacetanilide parent compounds (acetochlor, alachlor, and metolachlor) at the USGS gage site was 3.4 microg L(-1), whereas it was 4.3 microg L(-1) for the six metabolites. These data confirm the importance of studying chloroacetanilide metabolites, along with parent compounds, in tile-drained watersheds.

Relating net nitrogen input in the Mississippi River basin to nitrate flux in the lower Mississippi River: a comparison of approaches

A quantitative understanding of the relationship between terrestrial N inputs and riverine N flux can help guide conservation, policy, and adaptive management efforts aimed at preserving or restoring water quality. The objective of this study was to compare recently published approaches for relating terrestrial N inputs to the Mississippi River basin (MRB) with measured nitrate flux in the lower Mississippi River. Nitrogen inputs to and outputs from the MRB (1951 to 1996) were estimated from state-level annual agricultural production statistics and NOy (inorganic oxides of N) deposition estimates for 20 states that comprise 90% of the MRB. A model with water yield and gross N inputs accounted for 85% of the variation in observed annual nitrate flux in the lower Mississippi River, from 1960 to 1998, but tended to underestimate high nitrate flux and overestimate low nitrate flux. A model that used water yield and net anthropogenic nitrogen inputs (NANI) accounted for 95% of the variation in riverine N flux. The NANI approach accounted for N harvested in crops and assumed that crop harvest in excess of the nutritional needs of the humans and livestock in the basin would be exported from the basin. The U.S. White House Committee on Natural Resources and Environment (CENR) developed a more comprehensive N budget that included estimates of ammonia volatilization, denitrification, and exchanges with soil organic matter. The residual N in the CENR budget was weakly and negatively correlated with observed riverine nitrate flux. The CENR estimates of soil N mineralization and immobilization suggested that there were large (2000 kg N ha-1) net losses of soil organic N between 1951 and 1996. When the CENR N budget was modified by assuming that soil organic N levels have been relatively constant after 1950, and ammonia volatilization losses are redeposited within the basin, the trend of residual N closely matched temporal variation in NANI and was positively correlated with riverine nitrate flux in the lower Mississippi River. Based on results from applying these three modeling approaches, we conclude that although the NANI approach does not address several processes that influence the N cycle, it appears to focus on the terms that can be estimated with reasonable certainty and that are correlated with riverine N flux.

Nitrate flux in the Mississippi River

Increased delivery of biologically available nitrogen to estuaries and coastal oceans in recent decades has been linked to eutrophication and seasonal hypoxia in the northern Gulf of Mexico and elsewhere. We have developed a model that accounts for 95% of annual variation in delivery of nitrate to the Gulf of Mexico by the Mississippi River in 1960-98. Retrospective analysis indicates that this nitrate flux could have been reduced by 33% if the use of nitrogen-containing fertilizer in the Mississippi River basin had been cut by 12%.

Estimated historical and current nitrogen balances for Illinois

The Midwest has large riverine exports of nitrogen (N), with the largest flux per unit area to the Mississippi River system coming from Iowa and Illinois. We used historic and current data to estimate N inputs, outputs, and transformations for Illinois where human activity (principally agriculture and associated landscape drainage) have had a dominant impact. Presently, approximately 800,000 Mg of N is added each year as fertilizer and another 420,000 Mg is biologically fixed, primarily by soybean (Glycine max L. Merr.). These annual inputs are greater than exports in grain, which results in surplus N throughout the landscape. Rivers within the state export approximately 50% of this surplus N, mostly as nitrate, and the remainder appears to be denitrified or temporarily incorporated into the soil organic matter pool. The magnitude of N losses for 1880, 1910, 1950, and 1990 are compared. Initial cultivation of the prairies released large quantities of N (approximately 500,000 Mg N year(-1)), and resulted in riverine N transport during the late 19th century that appears to have been on the same order of magnitude as contemporary N losses. Riverine flux was estimated to have been at a minimum in about 1950, due to diminished net mineralization and low fertilizer inputs. Residual fertilizer N from corn (Zea mays L.), biological N fixed by soybean, short-circuiting of soil water through artificial drainage, and decreased cropping-system diversity appear to be the primary sources for current N export.

Determination of the quench velocity and rewetting temperature of hot surfaces: Formulation of a nonisothermal microscale hydrodynamic model

A nonisothermal microscale model of the three-phase, solid-liquid-gas, contact zone is formulated in the context of rewetting phenomena. The model incorporates hydrodynamics, heat transfer, interfacial phenomena, and intermolecular long range forces, in a two-dimensional proximal region of the order of 1000 A in width and 100 A in thickness. The model comprises scaled mass, momentum, and energy balances, and their corresponding scaled boundary conditions. The small contact angles which are characteristic of rewetting situations facilitate the use of the lubrication approximation, and the dynamics of the liquid and gas phases is decoupled by applying the one-sided simplification. The microscale hydrodynamic model reflects the strong effect of the solid-liquid interactions on the film profile, and the attendant flow and thermal fields. Thinner films having smaller contact angles involve stronger solid-liquid attraction forces, and consequently they exhibit higher rewetting temperatures and lower evaporation and vapor recoil effects. Thermocapillary and evaporation and conduction effects are expressed by appropriate dimensionless numbers. A set of such numbers is defined in the context of the differential equations of the microscale model. This model covers the hydrodynamic aspect of rewetting phenomena, which are also controlled by thermodynamic and macroscale constraints. This calls for interfacing and appropriate combination between the microscale hydrodynamic model, thermodynamics, and other macroscale rewetting models, for the determination of rewetting temperatures and quench velocities of liquids on hot solid surfaces. This is addressed elsewhere, in subsequent papers that follow this work.

Soil Chemistry in a Loblolly/Longleaf Pine Forest with Interval Burning

We examined the 30-yr cumulative effects of prescribed fires at intervals of 1, 2, 3, and 4 yr in a loblolly and longleaf pine forest in the Coastal Plain of South Carolina. The fine fraction of the forest floor (Oe + Oa horizons) contained much more carbon and nitrogen per unit area in the control plots (1.7 and 0.05 kg/m² , respectively) than in the 1-yr burn interval plots (0.4 and 0.007 kg/m² , respectively). Mineral soils (0-0.2 m depth) were highly variable in chemistry, and showed only slight differences across the burning treatments for nitrogen and sulfur. No trends were apparent for phosphorus in the forest floor or mineral soil; differences in acidity and extractable cations were also slight. The nutrient content of foliage was generally low, with no differences across burning intervals. Our results are consistent with earlier studies that showed the biogeochemical effects of repeated surface fires in southern pine forests are generally limited to the forest floor, with the possible exception of overall reductions in nitrogen cycling.

Autumnal changes of sulfur fractions and the ratio of organic sulfur to total nitrogen in leaves and adjacent bark of eastern cottonwood, white basswood and actinorhizal black alder

Autumnal changes in organic-S, sulfate-S, total-S and the ratios of organic-S to total-N and sulfate-S to organic-S were followed in leaves and adjacent bark of actinorhizal (Frankia-nodulated) black alder (Alnus glutinosa (L.) Gaertn.) and eastern cottonwood (Populus deltoides Bartr. ex Marsh.) trees growing on a minespoil site high in extractable soil sulfate, and in black alder and white basswood (Tilia heterophylla Venten.) trees growing on a prairie-derived soil in Illinois. Organic-S concentrations decreased significantly (P < 0.05) during autumn only in foliage of trees growing on the prairie-derived soil where losses of leaf organic-S were 65% for black alder and 100% for white basswood. Leaf sulfate concentrations were relatively stable throughout autumn in white basswood growing on prairie-derived soil and in black alder at both sites. Sulfate-S concentrations in leaves were significantly (P < 0.05) higher in trees at the minespoil site than in trees growing in the prairie-derived soil (5.1 mg g(-1) for the minespoil site and 1.2 mg g(-1) for the prairie-derived soil), and in the non-actinorhizal species during late summer. During the autumn, the ratio of organic-S to total-N doubled in leaves of eastern cottonwood at the minespoil site, but in black alder and white basswood growing on the prarie-derived soil, it decreased by 60 and 74%, respectively. Organic-S concentrations in bark increased more during autumn in species unable to fix atmospheric N(2), than in black alder. The results suggest that patterns of autumnal translocation of leaf S can be site-dependent and that leaf S and leaf N are, at least in part, translocated independently in the fall. Black alder and eastern cottonwood seemed to incorporate sulfate-S readily into organic substances in leaves when grown in soils with a high sulfate content.