Efficient and Practical Approaches to Ground‐Water Right Transfers Under the Prior Appropriation Doctrine and the Snake River Example1

Abstract: Water right transfers are one of the basic means of implementing changes in water use in the highly appropriated water resource systems of the western United States. Many of these systems are governed by the Prior Appropriation Doctrine, which was not originally intended for application to ground‐water pumping and the conjunctive management of ground water and surface water, and thus creates an administrative challenge. That challenge results from the fact that ground‐water pumping can affect all interconnected surface‐water bodies and the effects may be immeasurably small relative to surface water discharge and greatly attenuated in time. Although we may have the ability to calculate the effects of ground‐water pumping and transfers of pumping location on surface‐water bodies, mitigating for all the impacts of each individual transfer is sufficiently inefficient that it impedes the transfer process, frustrates water users, and consequently inhibits economic development. A more holistic approach to ground‐water right transfers, such as a ground‐water accounting or banking scheme, may adequately control transfer third‐party effects while reducing mitigation requirements on individual transfers. Acceptance of an accounting scheme can accelerate the transfer process, and possibly reduce the administrative burden.     

Sensitivity and Uncertainty of Ground‐Water Discharge Estimates for Semiarid Shrublands1

Abstract: We proposed a step‐by‐step approach to quantify the sensitivity of ground‐water discharge by evapotranspiration (ET) to three categories of independent input variables. To illustrate the approach, we adopt a basic ground‐water discharge estimation model, in which the volume of ground water lost to ET was computed as the product of the ground‐water discharge rate and the associated area. The ground‐water discharge rate was assumed to equal the ET rate minus local precipitation. The objective of this study is to outline a step‐by‐step procedure to quantify the contributions from individual independent variable uncertainties to the uncertainty of total ground‐water discharge estimates; the independent variables include ET rates of individual ET units, areas associated with the ET units, and precipitation in each subbasin. The specific goal is to guide future characterization efforts by better targeting data collection for those variables most responsible for uncertainty in ground‐water discharge estimates. The influential independent variables to be included in the sensitivity analysis are first selected based on the physical characteristics and model structure. Both regression coefficients and standardized regression coefficients for the selected independent variables are calculated using the results from sampling‐based Monte Carlo simulations. Results illustrate that, while as many as 630 independent variables potentially contribute to the calculation of the total annual ground‐water discharge for the case study area, a selection of seven independent variables could be used to develop an accurate regression model, accounting for more than 96% of the total variance in ground‐water discharge. Results indicate that the variability of ET rate for moderately dense desert shrubland contributes to about 75% of the variance in the total ground‐water discharge estimates. These results point to a need to better quantify ET rates for moderately dense shrubland to reduce overall uncertainty in estimates of ground‐water discharge. While the approach proposed here uses a basic ground‐water discharge model taken from an earlier study, the procedure of quantifying uncertainty and sensitivity can be generalized to handle other types of environmental models involving large numbers of independent variables.     

Evaluation of Ground‐Water and Surface‐Water Exchanges Using Streamflow Difference Analyses1

Abstract: In the karstic lower Flint River Basin, limestone fracturing, jointing, and subsequent dissolution have resulted in the development of extensive secondary permeability and created a system of major conduits that facilitate the exchange of water between the Upper Floridan aquifer and Flint River. Historical streamflow data from U.S. Geological Survey gaging stations located in Albany and Newton, Georgia, were used to quantify ground‐water and surface‐water exchanges within a 55.3 km section of the Flint River. Using data from 2001, we compared estimates of ground‐water flux using a time adjustment method to a water balance equation and found that these independent approaches yielded similar results. The associated error was relatively large during high streamflow when unsteady conditions prevail, but much lower during droughts. Flow reversals were identified by negative streamflow differences and verified with in situ data from temperature sensors placed inside large spring conduits. Long‐term (13 years) analysis showed negative streamflow differentials (i.e., a losing stream condition) coincided with high river stages and indicated that streamflow intrusion into the aquifer could potentially exceed 150 m³/s. Although frequent negative flow differentials were evident, the Flint River was typically a gaining stream and showed a large net increase in flow between the two gages when examined over the period 1989‐2003. Ground‐water contributions to this stream section averaged 2‐42 m³/s with a mean of 13 m³/s. The highest rate of ground‐water discharge to the Flint River occurred during the spring when regional ground‐water levels peaked following heavy winter and spring rains and corresponding rates of evapotranspiration were low. During periods of extreme drought, ground‐water contributions to the Flint River declined.     

Use of Streamflow Data to Estimate Base Flow/Ground‐Water Recharge For Wisconsin1

Abstract: The average annual base flow/recharge was determined for streamflow‐gaging stations throughout Wisconsin by base‐flow separation. A map of the State was prepared that shows the average annual base flow for the period 1970‐99 for watersheds at 118 gaging stations. Trend analysis was performed on 22 of the 118 streamflow‐gaging stations that had long‐term records, unregulated flow, and provided aerial coverage of the State. The analysis found that a statistically significant increasing trend was occurring for watersheds where the primary land use was agriculture. Most gaging stations where the land cover was forest had no significant trend. A method to estimate the average annual base flow at ungaged sites was developed by multiple‐regression analysis using basin characteristics. The equation with the lowest standard error of estimate, 9.5%, has drainage area, soil infiltration and base flow factor as independent variables. To determine the average annual base flow for smaller watersheds, estimates were made at low‐flow partial‐record stations in 3 of the 12 major river basins in Wisconsin. Regression equations were developed for each of the three major river basins using basin characteristics. Drainage area, soil infiltration, basin storage and base‐flow factor were the independent variables in the regression equations with the lowest standard error of estimate. The standard error of estimate ranged from 17% to 52% for the three river basins.     

Riparian Ground‐Water Flow Patterns Using Flownet Analysis: Evapotranspiration‐Induced Upwelling and Implications for N Removal1

Abstract: Ground‐water flow paths constrain the extent of nitrogen (N) sinks in deep, stratified soils of riparian wetlands. We examined ground‐water flow paths at four forested riparian wetlands in deep, low gradient, stratified deposits subjected to Southern New England’s temperate, humid climate. Mid‐day piezometric heads were recorded during the high water table period in April/May and again in late November at one site. Coupling field data with a two‐dimensional steady‐state ground‐water flow model, flow paths and fluxes were derived to 3 m depths. April/May evapotranspiration (ET) dominated total outflux (44‐100%) while flux to the stream was <10% of total outflux. ET exerted upward ground‐water flux through shallow carbon‐rich soils, increasing opportunities for N transformations and diverting flow from the stream. Dormant season results showed a marked increase in flux to the stream (27% of the total flux). Riparian sites with deep water tables (naturally or because of increased urbanization or other hydrologic modifications) or shallow root zones may not generate ground‐water upwelling to meet evaporative demand, thereby increasing the risk of N movement to streams. As water managers balance issues of water quality with water quantity, they will be faced with decisions regarding riparian management. Further work towards refining our understanding of ET mediation of N and water flux at the catchment scale will serve to inform these decisions.     

Evaluation of Rio Grande Management Alternatives Using a Surface‐Water/Ground‐Water Model1

Abstract: Previous investigations observed significant seepage losses from the Rio Grande to the shallow aquifer between Socorro and San Antonio, New Mexico. High‐resolution telescopic modeling was used along a 10‐km reach of the Rio Grande and associated drains and canals to evaluate several management alternatives aimed at improving river conveyance efficiency. Observed data consisted of ground‐water and surface‐water elevations, seepage rates along the Rio Grande and associated canals and drains, and borehole geology. Model calibration was achieved by adjusting hydraulic conductivity and specific storage until the output matched observed data. Sensitivity analyses indicated that the system was responsive to changes in hydrogeologic properties, especially when such alterations increased vertical connectivity between layers. The calibrated model predicted that removal of the low flow conveyance channel, a major channel draining the valley, would not only decrease river seepage by 67%, but also decrease total flow through the reach by 75%. The decreased flow through the reach would result in increased water logging and an average increase in ground‐water elevations of 1.21 meter. Simulations of the system with reduced riparian evapotranspiration rates or a relocated river channel also predicted decreased river seepage, but to a much lesser degree.     

Drought Resilience of the California Central Valley Surface‐Ground‐Water‐Conveyance System1

Abstract: A series of drought simulations were performed for the California Central Valley using computer applications developed by the California Department of Water Resources and historical datasets representing a range of droughts from mild to severe for time periods lasting up to 60 years. Land use, agricultural cropping patterns, and water demand were held fixed at the 2003 level and water supply was decreased by amounts ranging between 25 and 50%, representing light to severe drought types. Impacts were examined for four hydrologic subbasins, the Sacramento Basin, the San Joaquin Basin, the Tulare Basin, and the Eastside Drainage. Results suggest the greatest impacts are in the San Joaquin and Tulare Basins, regions that are heavily irrigated and are presently overdrafted in most years. Regional surface water diversions decrease by as much as 70%. Stream‐to‐aquifer flows and aquifer storage declines were proportional to drought severity. Most significant was the decline in ground water head for the severe drought cases, where results suggest that under these scenarios the water table is unlikely to recover within the 30‐year model‐simulated future. However, the overall response to such droughts is not as severe as anticipated and the Sacramento Basin may act as ground‐water insurance to sustain California during extended dry periods.     

Applying Zonation Methods and Tabu Search to Improve the Ground‐Water Modeling1

Abstract: An optimization procedure combining zonation methods with the Tabu Search method is proposed to identify the spatial distribution of hydraulic conductivity field. Three zonation methods, Voronoi diagram (VD), multiplicatively weighted Voronoi diagram (MWVD), and pattern zonation are adopted for the parameterization purposes. With the consideration of the homogeneity and the heterogeneity, there are four spatial distributions of hydraulic conductivity designed to test whether the parameter structure can be successfully identified. The fitting residual error is first considered to determine an adequate number of zones without over parameterization. Then, the parameter uncertainty is evaluated the decision of the number of zones. The results indicate that the MWVD performs better than other two methods because the MWVD has better flexibility in describing the zonal boundaries with small number of decision variables.     

Surface Water Seepage Effects on Shallow Ground‐Water Quality Along the Rio Grande in Northern New Mexico1

Abstract: Interactions between surface irrigation water, shallow ground water, and river water may have effects on water quality that are important for both drinking water supplies and the ecological function of rivers and floodplains. We investigated water quality in surface water and ground water, and how water quality is influenced by surface water inputs from an unlined irrigation system in the Alcalde Valley of the Rio Grande in northern New Mexico. From August 2005 to July 2006, we sampled ground water and surface water monthly and analyzed for concentrations of major cations and anions, specific conductance, pH, dissolved oxygen, and water levels. Results indicate that irrigation ditch seepage caused an increase in ground water levels and that the Rio Grande is a gaining stream in this region. Temporal and spatial differences were found in ion concentrations in shallow ground water as it flowed from under the ditch toward the river. Ground‐water ion concentrations were higher when the ditch was not flowing compared with periods during peak irrigation season when the ditch was flowing. Ditch inputs diluted ion concentrations in shallow ground water at well positions near the ditch. Specifically, lower ion concentrations were detected in ground water at well positions located near the ditch and river compared with well positions located in the middle of an agricultural field. Results from this project showed that ditch inputs influenced ion concentrations and were associated with ground‐water recharge. In arid region river valleys, careful consideration should be given to management scenarios that change seepage from irrigation systems, because in some situations reduced seepage could negatively affect ground‐water recharge and water quality.     

Effects of Rainfall and Ground‐Water Pumping on Streamflow in Mākaha,O’ahu,Hawai’i1

Abstract: Land‐use/land‐cover changes in Mākaha valley have included the development of agriculture, residential dwellings, golf courses, potable water supply facilities, and the introduction of alien species. The impact of these changes on surface water and ground water resources in the valley is of concern. In this study, streamflow, rainfall, and ground‐water pumping data for the upper part of the Mākaha valley watershed were evaluated to identify corresponding trends and relationships. The results of this study indicate that streamflow declined during the ground‐water pumping period. Mean and median annual streamflow have declined by 42% (135 mm) and 56% (175 mm), respectively, and the mean number of dry stream days per year has increased from 8 to 125. Rainfall across the study area appears to have also declined though it is not clear whether the reduction in rainfall is responsible for all or part of the observed streamflow decline. Mean annual rainfall at one location in the study area declined by 14% (179 mm) and increased by 2% (48 mm) at a second location. Further study is needed to assess the effect of ground‐water pumping and to characterize the hydrologic cycle with respect to rainfall, infiltration, ground‐water recharge and flow in the study area, and stream base flow and storm flow.     

Managing the Commons Texas Style: Wildlife Management and Ground‐Water Associations on Private Lands1

Abstract: As nearly all of Texas’ rural lands are privately owned, landowner associations for the management of white‐tailed deer and ground‐water have become increasingly popular. Deer are a common‐pool resource with transboundary characteristics, requiring landowner cooperation for effective management. Ground‐water reserves are economically important to landowners, but are governed by the “rule of capture” whereby property rights are not defined. One ground‐water association and four wildlife management associations (WMAs) were surveyed to characterize their member demographics, land use priorities, attitudes, and social capital. Members of the ground‐water cooperative were part of a much larger, more heterogeneous, and more recently formed group than members of WMAs. They also placed greater importance on utilitarian aspects of their properties, as opposed to land stewardship for conservation as practiced by members of WMAs. If ground‐water association members could be more locally organized with more frequent meetings, social capital and information sharing may be enhanced and lead to land stewardship practices for improved hydrologic functions and sustained ground‐water supply. This, coupled with pumping rules assigned by the local ground‐water district, could yield an effective strategy that is ecologically and hydrologicaly sound, and that allows rural provision of water supply to urban consumers.     

Field‐Scale Application of Three Types of Neural Networks to Predict Ground‐Water Levels1

Abstract: In this paper, a field‐scale applicability of three forms of artificial neural network algorithms in forecasting short‐term ground‐water levels at specific control points is presented. These algorithms are the feed‐forward back propagation (FFBP), radial basis networks (RBN), and generalized regression networks (GRN). Ground‐water level predictions from these algorithms are in turn to be used in an Optimized Regional Operations Plan that prescribes scheduled wellfield production for the coming four weeks. These models are up against each other for their accuracy of ground‐water level predictions on lead times ranging from a week to four weeks, ease of implementation, and execution times (mainly training time). In total, 208 networks of each of the three algorithms were developed for the study. It is shown that although learning algorithms have emerged as a viable solution at field scale much larger than previously studied, no single algorithm performs consistently better than others on all the criteria. On average, FFBP networks are 20 and 26%, respectively, more accurate than RBN and GRN in forecasting one week ahead water levels and this advantage drops to 5 and 9% accuracy in forecasting four weeks ahead water levels, whereas GRN posted a training time that is only 5% of the training time taken by that of FFBP networks. This may suggest that in field‐scale applications one may have to trade between the type of algorithm to be used and the degree to which a given objective is honored.     

Estimating Evapotranspiration and Seepage for a Sinkhole Wetland From Diurnal Surface‐Water Cycles1

Abstract: This study used measured diurnal surface‐water cycles to estimate daily evapotranspiration (ET) and seepage for a seasonally flooded sinkhole wetland. Diurnal surface‐water cycles were classified into five categories based on the relationship between the surface‐water body and the surrounding ground‐water system (i.e., recharge/discharge). Only one class of diurnal cycles was found to be suitable for application of this method. This subset of diurnal cycles was used to estimate ET and seepage and the relative importance of each transfer process to the overall water budget. The method has limited utility for wetlands with erratic hydrologic regimes (e.g., wetlands in urban environments). This is due to violation of the critical assumption that the inflow/outflow rate remains constant throughout the day. For application to surface‐water systems, the method is typically applied with an assumed specific yield of 1.0. This assumption was found to be invalid for application to surface‐water systems with a noncylindrical pond geometry. An overestimation of ET by as much as 60% was found to occur under conditions of low pond stage and high water loss. The results demonstrate the high ET rates that can occur in isolated wetlands due to contrasting roughness and moisture conditions (oasis and clothesline effects). Estimated ET rates ranged from 4.1 to 18.7 mm/day during the growing season. Despite these large ET rates, seepage (recharge) was found to be the dominant water loss mechanism for the wetland.     

Differences in Dissolved Organic Carbon and Nitrogen Responses to Storm‐Event and Ground‐Water Conditions in a Forested,Glaciated Watershed in Western New York1

Abstract: Differences in the storm‐event responses of dissolved organic carbon (DOC) and nitrogen (DON) in streamflow and ground water were evaluated for a glaciated forested watershed in western New York. Eight to ten storm events with varying rainfall amounts, intensities, and antecedent moisture conditions were studied for three catchments (1.6, 3.4, and 696 ha) over a three‐year period (2003‐2005). Concentrations of DOC in streamflow exiting the catchments were significantly higher for storm events following a dry period, whereas no similar response was observed for DON. Highest DON concentrations in streamflow were typically associated with storm events following wet antecedent moisture conditions. In addition to antecedent moisture conditions, DOC concentrations were also positively correlated with precipitation amounts, while DON did not reveal a consistent pattern. Streamwater and ground‐water concentrations of DOC during storm events were also strongly correlated with riparian ground‐water depths but a similar relationship was not observed for DON. Ground‐water DON concentrations were also more variable than DOC. We hypothesized that the differences in DOC and DON responses stemmed from the differences in catchment sources of these solutes. This study suggests that while DOC and DON are intrinsically linked as components of dissolved organic matter, their dynamics and exports from watersheds may be regulated by a different set of mechanisms and factors. Identifying these differences is critical for developing more reliable and robust models for transport of dissolved organic matter.     

四川省环境饮水的^3H水平,分布及评价

本文报道了四川省20市、地、州的环境饮水(自来水、井水及部分江河、湖库、泉水)中之~3H辐射水平、变化规律与分布特点,对其影响因素进行了简要探讨和评价。 3H(氚)是氢的唯一放射性同位素,发射低能炖β粒子,在许多生产和研究领域得到了日益广泛的应用。它是核燃料循环中可能对环境产生较大影响的核素之一,更是未来的能源——可控热核反应堆可能对环境产生影响的主要核素。因此,氚是人们进行核环境评价时感兴趣的主要放射性核素之一。事实上氚,是水体中含量最高的放射性核素。因此,氚在环境水中的现状,以及对氚进行有效的监测监督,已成为人们十分关注的问题。     

Hydrologic Model Framework for Water Resources Planning in the Santa Cruz River,Southern Arizona1

Abstract: The authors develop a model framework that includes a set of hydrologic modules as a water resources management and planning tool for the upper Santa Cruz River near the Mexican border, Southern Arizona. The modules consist of: (1) stochastic generation of hourly precipitation scenarios that maintain the characteristics and variability of a 45‐year hourly precipitation record from a nearby rain gauge; (2) conceptual transformation of generated precipitation into daily streamflow using varied infiltration rates and estimates of the basin antecedent moisture conditions; and (3) surface‐water to ground‐water interaction for four downstream microbasins that accounts for alluvial ground‐water recharge, and ET and pumping losses. To maintain the large inter‐annual variability of streamflow as prevails in Southern Arizona, the model framework is constructed to produce three types of seasonal winter and summer categories of streamflow (i.e., wet, medium, or dry). Long‐term (i.e., 100 years) realizations (ensembles) are generated by the above described model framework that reflects two different regimes of inter annual variability. The first regime is that of the historic streamflow gauge record. The second regime is that of the tree ring reconstructed precipitation, which was derived for the study location. Generated flow ensembles for these two regimes are used to evaluate the risk that the regional four ground‐water microbasins decline below a preset storage threshold under different operational water utilization scenarios.     

Clogging of an Effluent Dominated Semiarid River: A Conceptual Model of Stream‐Aquifer Interactions1

Abstract: Water managers in arid and semiarid regions increasingly view treated wastewater (effluent) as an important water resource. Artificial recharge basins allow effluent to seep into the ground relieving stressed aquifers, however these basins frequently clog due to physical, chemical, and biological processes. Likewise effluent is increasingly used to maintain perennial base flow for dry streambeds, however, little is known about the impact of effluent on streambed hydraulic conductivity and stream‐aquifer interactions. We address this issue by investigating: if a clogging layer forms, how the formation of a clogging layer alters stream‐aquifer connections, and what hydrologic factors control the formation and removal of clogging layers. We focused on the Upper Santa Cruz River, Arizona where effluent from the Nogales International Waste Water Treatment Plant sustains perennial flow. Monthly sampling, along a 30 km river reach, was done with two foci: physical streambed transformations and water source identification using chemical composition. Historical dataset were included to provide a larger context for the work. Results show that localized clogging occurs in the Upper Santa Cruz River. The clogging layers perch the stream and shallow streambed causing desaturation below the streambed. With these results, a conceptual model of clogging is established in the context of a semiarid hydrologic cycle: formation during the hot premonsoon months when flow is nearly constant and removal by large flood flows (>10 m³/s) during the monsoon season. However, if the intensity of flooding during the semiarid hydrologic cycle is lessened, the dependent riparian area can experience a die off. This conceptual model leads us to the conclusion that effluent dominated riparian systems are inherently unstable due to the clogging process. Further understanding of this process could lead to improved ecosystem restoration and management.     

Assessing Satellite and Ground‐Based Potential Evapotranspiration for Hydrologic Applications in the Colorado River Basin

This study evaluates a remotely sensed and two ground‐based potential evapotranspiration (PET) products for hydrologic application in the Upper Colorado River Basin (UCRB). The remotely sensed Moderate Resolution Imaging Spectroradiometer product (MODIS‐PET) is a continuous, daily time series with 250 m resolution derived using the Priestley‐Taylor (P‐T) equation. The MODIS‐PET is evaluated against regional flux tower data as well as a synthetic pan product (Epan; 0.125°, daily) derived from the North American Land Data Assimilation System (NLDAS) and a Hargreaves PET derived from DAYMET variables (DAYMET‐PET; 1 km, daily). Compared to point‐scale PET computed using regional flux tower data, the MODIS‐PET had lower errors, with RMSE values ranging from 2.24 to 2.85 mm/day. Epan RMSE values ranged from 3.70 to 3.76 mm/day and DAYMET‐PET RMSE values ranged from 3.55 to 4.58 mm/day. Further investigation showed biases in temperature and radiation data contribute to uncertainty in the MODIS‐PET values, while bias in NLDAS temperature, downward shortwave (SW↓), and downward longwave (LW↓) propagate in the Epan estimates. Larger discrepancies between methods were observed in the warmer, drier regions of the UCRB, however, the MODIS‐PET was more responsive to landcover transitions and better captured basin heterogeneity. Results indicate the satellite‐based MODIS product can serve as a viable option for obtaining spatial PET values across the UCRB.     

Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data1

Nishat, Bushra and S.M. Mahbubur Rahman, 2009. Water Resources Modeling of the Ganges‐Brahmaputra‐Meghna River Basins Using Satellite Remote Sensing Data. Journal of the American Water Resources Association (JAWRA) 45(6):1313‐1327. Abstract: Large‐scale water resources modeling can provide useful insights on future water availability scenarios for downstream nations in anticipation of proposed upstream water resources projects in large international river basins (IRBs). However, model set up can be challenging due to the large amounts of data requirement on both static states (soils, vegetation, topography, drainage network, etc.) and dynamic variables (rainfall, streamflow, soil moisture, evapotranspiration, etc.) over the basin from multiple nations and data collection agencies. Under such circumstances, satellite remote sensing provides a more pragmatic and convenient alternative because of the vantage of space and easy availability from a single data platform. In this paper, we demonstrate a modeling effort to set up a water resources management model, MIKE BASIN, over the Ganges, Brahmaputra, and Meghna (GBM) river basins. The model is set up with the objective of providing Bangladesh, the lowermost riparian nation in the GBM basins, a framework for assessing proposed water diversion scenarios in the upstream transboundary regions of India and deriving quantitative impacts on water availability. Using an array of satellite remote sensing data on topography, vegetation, and rainfall from the transboundary regions, we demonstrate that it is possible to calibrate MIKE BASIN to a satisfactory level and predict streamflow in the Ganges and Brahmaputra rivers at the entry points of Bangladesh at relevant scales of water resources management. Simulated runoff for the Ganges and Brahmaputra rivers follow the trends in the rated discharge for the calibration period. However, monthly flow volume differs from the actual rated flow by (?) 8% to (+) 20% in the Ganges basin, by (?) 15 to (+) 12% in the Brahmaputra basin, and by (?) 15 to (+) 19% in the Meghna basin. Our large‐scale modeling initiative is generic enough for other downstream nations in IRBs to adopt for their own modeling needs.     

Estimating Nitrogen Loading to Ground Water and Assessing Vulnerability to Nitrate Contamination in a Large Karstic Springs Basin,Florida1

Abstract: A nitrogen (N) mass‐balance budget was developed to assess the sources of N affecting increasing ground‐water nitrate concentrations in the 960‐km² karstic Ichetucknee Springs basin. This budget included direct measurements of N species in rainfall, ground water, and spring waters, along with estimates of N loading from fertilizers, septic tanks, animal wastes, and the land application of treated municipal wastewater and residual solids. Based on a range of N leaching estimates, N loads to ground water ranged from 262,000 to 1.3 million kg/year; and were similar to N export from the basin in spring waters (266,000 kg/year) when 80‐90% N losses were assumed. Fertilizers applied to cropland, lawns, and pine stands contributed about 51% of the estimated total annual N load to ground water in the basin. Other sources contributed the following percentages of total N load to ground water: animal wastes, 27%; septic tanks, 12%; atmospheric deposition, 8%; and the land application of treated wastewater and biosolids, 2%. Due to below normal rainfall (97.3 cm) during the 12‐month rainfall collection period, N inputs from rainfall likely were about 30% lower than estimates for normal annual rainfall (136 cm). Low N‐isotope values for six spring waters (δ¹⁵N‐NO₃ = 3.3 to 6.3‰) and elevated potassium concentrations in ground water and spring waters were consistent with the large N contribution from fertilizers. Given ground‐water residence times on the order of decades for spring waters, possible sinks for excess N inputs to the basin include N storage in the unsaturated zone and parts of the aquifer with relatively sluggish ground‐water movement and denitrification. A geographical‐based model of spatial loading from fertilizers indicated that areas most vulnerable to nitrate contamination were located in closed depressions containing sinkholes and other dissolution features in the southern half of the basin.