Spatiotemporal Variability of Modeled Watershed Scale Surface‐Depression Storage and Runoff for the Conterminous United States

This study explores the viability of using simulated monthly runoff as a proxy for landscape‐scale surface‐depression storage processes simulated by the United States Geological Survey’s National Hydrologic Model (NHM) infrastructure across the conterminous United States (CONUS). Two different temporal resolution model codes (daily and monthly) were run in the NHM with the same spatial discretization. Simulated values of daily surface‐depression storage (treated as a decimal fraction of maximum volume) as computed by the daily Precipitation‐Runoff Modeling System (NHM‐PRMS) and normalized runoff (0 to 1) as computed by the Monthly Water Balance Model (NHM‐MWBM) were aggregated to monthly and annual values for each hydrologic response unit (HRU) in the CONUS geospatial fabric (HRU; n = 109,951) and analyzed using Spearman’s rank correlation test. Correlations between simulated runoff and surface‐depression storage aggregated to monthly and annual values were compared to identify where which time scale had relatively higher correlation values across the CONUS. Results show Spearman’s rank values >0.75 (highly correlated) for the monthly time scale in 28,279 HRUs (53.35%) compared to the annual time scale in 41,655 HRUs (78.58%). The geographic distribution of HRUs with highly correlated monthly values show areas where surface‐depression storage features are known to be common (e.g., Prairie Pothole Region, Florida).     

Surface‐Water Nutrient Conditions and Sources in the United States Pacific Northwest1

Wise, Daniel R. and Henry M. Johnson, 2011. Surface‐Water Nutrient Conditions and Sources in the United States Pacific Northwest. Journal of the American Water Resources Association (JAWRA) 47(5):1110‐1135. DOI: 10.1111/j.1752‐1688.2011.00580.x Abstract: The SPAtially Referenced Regressions On Watershed attributes (SPARROW) model was used to perform an assessment of surface‐water nutrient conditions and to identify important nutrient sources in watersheds of the Pacific Northwest region of the United States (U.S.) for the year 2002. Our models included variables representing nutrient sources as well as landscape characteristics that affect nutrient delivery to streams. Annual nutrient yields were higher in watersheds on the wetter, west side of the Cascade Range compared to watersheds on the drier, east side. High nutrient enrichment (relative to the U.S. Environmental Protection Agency’s recommended nutrient criteria) was estimated in watersheds throughout the region. Forest land was generally the largest source of total nitrogen stream load and geologic material was generally the largest source of total phosphorus stream load generated within the 12,039 modeled watersheds. These results reflected the prevalence of these two natural sources and the low input from other nutrient sources across the region. However, the combined input from agriculture, point sources, and developed land, rather than natural nutrient sources, was responsible for most of the nutrient load discharged from many of the largest watersheds. Our results provided an understanding of the regional patterns in surface‐water nutrient conditions and should be useful to environmental managers in future water‐quality planning efforts.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.     

Identification of Putative Geographically Isolated Wetlands of the Conterminous United States

Geographically isolated wetlands (GIWs) are wetlands completely surrounded by uplands. While common throughout the United States (U.S.), there have heretofore been no nationally available, spatially explicit estimates of GIW extent, complicating efforts to understand the myriad biogeochemical, hydrological, and habitat functions of GIWs and hampering conservation and management efforts at local, state, and national scales. We used a 10‐m geospatial buffer as a proxy for hydrological or ecological connectivity of National Wetlands Inventory palustrine and lacustrine wetland systems to nationally mapped and available stream, river, and lake data. We identified over 8.3 million putative GIWs across the conterminous U.S., encompassing nearly 6.5 million hectares of wetland resources (average size 0.79 ± 4.81 ha, median size 0.19 ha). Putative GIWs thus represent approximately 16% of the freshwater wetlands of the conterminous U.S. The water regime for the majority of the putative GIWs was temporarily or seasonally flooded, suggesting a vulnerability to ditching or hydrologic abstraction, sedimentation, or alterations in precipitation patterns. Additional analytical applications of this readily available geospatially explicit mapping product (e.g., hydrological modeling, amphibian metapopulation, or landscape ecological analyses) will improve our understanding of the abundance and extent, effect, connectivity, and relative importance of GIWs to other aquatic systems of the conterminous U.S.     

Testing and Improving Temperature Thresholds for Snow and Rain Prediction in the Western United States

The phase of precipitation at the land surface is critical to determine the timing and amount of water available for hydrological and ecological systems. However, there are few techniques to directly observe the precipitation phase and many prediction tools apply a single temperature threshold (e.g., 0°C) to determine phase. In this paper, we asked two questions: (1) what is the accuracy of default and station optimized daily temperature thresholds for predicting precipitation phase and (2) what are the regions and conditions in which typical temperature‐based precipitation phase predictions are most suited. We developed a ground truth dataset of rain vs. snow using an expert decision‐making system based on precipitation, snow depth, and snow water equivalent observations. This dataset was used to evaluate the accuracy of three temperature‐threshold‐based techniques of phase classification. Optimizing the temperature threshold improved the prediction of precipitation phase by 34% compared to using 0°C threshold. Developing a temperature threshold based on station elevation improved the error by 12% compared with using the 0°C temperature threshold. We also found the probability of snow as a function of temperature differed among ecoregions, which suggests a varied response to future climate change. These results highlight a current weakness in our ability to predict the effects of regional warming that could have uneven impacts on water and ecological resources.     

Estimating the Variability and Confidence of Land Use and Imperviousness Relationships at a Regional Scale1

Abstract: Impervious cover is a commonly used metric to help explain or predict anthropogenic impacts on aquatic resources; often it is used as a surrogate for intensity of human impacts when evaluating effects on aquatic resources. The most common way to estimate imperviousness is based on relationships with land use. Few studies have evaluated how the relationship between impervious surface and land use varies among geographies with different levels of development and between types of imagery used to assign land use type. In this study, we assess variability in estimates of imperviousness based on two locally available land use datasets: one based on aerial imagery (2‐m resolution) and another based on satellite imagery (30‐m resolution). The ranges and variability in imperviousness within land use categories were assessed at several spatial scales, including within counties, between counties, and between watersheds. Results indicate that there was considerable variability for all developed land use types. Estimated impervious cover often varied over a range of 20‐40% points within a land use category. Furthermore, there were clear spatial patterns both between and within counties, with impervious cover for a given land use type being higher near the urban centers and lower at the margins of development. Estimates of imperviousness for 12 study watersheds indicated that variability increased with increasing watershed development, making it difficult to confidently set management or regulatory targets based on impervious cover. This study suggests that locally derived, high resolution satellite or aerial imagery should be used to estimate imperviousness when a high level of accuracy and precision is required for regulatory or management decisions. Furthermore, the error associated with impervious land use relationships should be accounted for when using impervious cover in runoff or water quality models, or when making management decisions regarding stream health.     

A New Open‐Access HUC‐8 Based Downscaled CMIP‐5 Climate Model Forecast Dataset for the Conterminous United States

Watershed‐scale hydrologic simulation models generally require climate data inputs including precipitation and temperature. These climate inputs can be derived from downscaled global climate simulations which have the potential to drive runoff forecasts at the scale of local watersheds. While a simulation designed to drive a local watershed model would ideally be constructed at an appropriate scale, global climate simulations are, by definition, arbitrarily determined large rectangular spatial grids. This paper addresses the technical challenge of making climate simulation model results readily available in the form of downscaled datasets that can be used for watershed scale models. Specifically, we present the development and deployment of a new Coupled Model Intercomparison Project phase 5 (CMIP5) based database which has been prepared through a scaling and weighted averaging process for use at the level of U.S. Geological Survey (USGS) Hydrologic Unit Code (HUC)‐8 watersheds. The resulting dataset includes 2,106 virtual observation sites (watershed centroids) each with 698 associated time series datasets representing average monthly temperature and precipitation between 1950 and 2099 based on 234 unique climate model simulations. The new dataset is deployed on a HydroServer and distributed using WaterOneFlow web services in the WaterML format. These methods can be adapted for downscaled General Circulation Model (GCM) results for specific drainage areas smaller than HUC‐8. Two example use cases for the dataset also are presented.     

Sensitivity of Landsat‐Scale Energy Balance to Aerodynamic Variability in Mountains and Complex Terrain

The speed and direction of air flow through complex terrain are difficult to define. Both impact sensible and latent heat flux exchanges at the surface. Evapotranspiration (ET) models such as Mapping EvapoTranspiration at high Resolution with Internalized Calibration (METRIC?) estimate ET as a residual of the surface energy process and are thus sensitive to aerodynamics, including terrain‐induced impacts on roughness governing convective heat transfer (H). There is a need to explore the sensitivities of H estimation and thereby ET estimation to wind speed and terrain roughness in mountainous areas and to determine the merit of operating complex mesoscale wind field models in conjunction with the energy balance process. A sensitivity analysis is explored in METRIC where we increased wind speed in proportion to a relative elevation parameter and we increased aerodynamic roughness to assimilate impacts of relative terrain roughness, estimated in proportion to standard deviation of elevation within a 3 km locality. These aerodynamic modifications increased convective heat transfer in complex terrain and reduced estimated ET. In other sensitivity runs, we reduced estimated wind speed on estimated leeward slopes. Estimated ET with and without these sensitivity adjustments is shown for mountainous areas of Montana and Nevada. Changes in ET ranged from little change (<5%) for lower slopes to about 30% reductions on windward slopes and 25% increases on leeward slopes for some mid to high elevations in the Montana application.     

Effect of Snow Cover Conditions on the Hydrologic Regime: Case Study in a Pluvial‐Nival Watershed,Japan1

Abstract: Hydrologic monitoring in a small forested and mountainous headwater basin in Niigata Prefecture has been undertaken since 2000. An important characteristic of the basin is that the hydrologic regime contains pluvial elements year‐round, including rain‐on‐snow, in addition to spring snowmelt. We evaluated the effect of different snow cover conditions on the hydrologic regime by analyzing observed data in conjunction with model simulations of the snowpack. A degree‐day snow model is presented and applied to the study basin to enable estimation of the basin average snow water equivalent using air temperature at three representative elevations. Analysis of hydrological time series data and master recession curves showed that flow during the snowmelt season was generated by a combination of ground water flow having a recession constant of 0.018/day and diurnal melt water flow having a recession constant of 0.015/hour. Daily flows during the winter/snowmelt season showed greater persistence than daily flows during the warm season. The seasonal water balance indicated that the ratio of runoff to precipitation during the cold season (December to May) was about 90% every year. Seasonal snowpack plays an important role in defining the hydrologic regime, with winter precipitation and snowmelt runoff contributing about 65% of the annual runoff. The timing of the snowmelt season, indicated by the date of occurrence of the first significant snowmelt event, was correlated with the occurrence of low flow events. Model simulations showed that basin average snow water equivalent reached a peak around mid‐February to mid‐March, although further validation of the model is required at high elevation sites.     

County‐Scale Rainwater Harvesting Feasibility in the United States: Climate,Collection Area,Density, and Reuse Considerations

Roof rainwater harvesting (RWH) has the potential to augment water supplies for urban and suburban uses throughout the United States (U.S.). Studies of the performance of RWH at the building and city scales in the U.S. are available, but a countrywide overview of the potential performance of RWH at the county scale has not been done before. Three approaches were taken: (1) assess the viability of RWH in terms of the rainfall that could be captured in relation to the water demand in each county (excluding agriculture), (2) evaluate the performance of a “typical” domestic RWH system across all counties with metrics related to its ability to supply the potable and nonpotable demand, and (3) evaluate the effect of adding a 50% rainwater reuse component to the analysis. We find RWH could be a viable supplemental water source in the U.S., particularly in counties of the Pacific Northwest, Central, and Eastern regions (percent demand covered >50%). Low population density counties have the potential to meet their annual water needs with RWH, while high‐density counties could only source a small portion (~20%) of their annual demand with RWH. Typical RWH systems in counties in the Central and Eastern U.S. performed better than in Western counties. Adding a reuse component can be a key factor in making RWH attractive in many areas of the country. This work can inform future water infrastructure investment and planning in the U.S.     

Intercomparison of Satellite Remote Sensing‐Based Flood Inundation Mapping Techniques

The objective of this study was to determine the accuracy of five different digital image processing techniques to map flood inundation extent with Landsat 8–Operational Land Imager satellite imagery. The May 2016 flooding event in the Hempstead region of the Brazos River, Texas is used as a case study for this first comprehensive comparison of classification techniques of its kind. Five flood water classification techniques (i.e., supervised classification, unsupervised classification, delta‐cue change detection, Normalized Difference Water Index [NDWI], modified NDWI [MNDWI]) were implemented to characterize flooded regions. To identify flood water obscured by cloud cover, a digital elevation model (DEM)–based approach was employed. Classified floods were compared using an Advanced Fitness Index to a “reference flood map” created based on manual digitization, as well as other data sources, using the same satellite image. Supervised classification yielded the highest accuracy of 86.4%, while unsupervised, MNDWI, and NDWI closely followed at 79.6%, 77.3%, and 77.1%, respectively. Delta‐cue change detection yielded the lowest accuracy with 70.1%. Thus, supervised classification is recommended for flood water classification and inundation map generation under these settings. The DEM‐based approach used to identify cloud‐obscured flood water pixels was found reliable and easy to apply. It is therefore recommended for regions with relatively flat topography.     

Glacier Variability (1967‐2006) in the Teton Range,Wyoming, United States1

Edmunds, Jake, Glenn Tootle, Greg Kerr, Ramesh Sivanpillai, and Larry Pochop, 2011. Glacier Variability (1967‐2006) in the Teton Range, Wyoming, United States. Journal of the American Water Resources Association (JAWRA) 48(1): 187‐196. DOI: 10.1111/j.1752‐1688.2011.00607.x Abstract: Glacier area and volume changes were quantified through the use of historical aerial photographs in Wyoming’s Teton Range. Glacier area changes in the Teton Range were estimated for three glaciers using unrectified aerial photography from 1967 to 2006. The total surface area of the three glaciers was 0.53 km² in 1967 and 0.40 km² in 2006, a decrease of 25% during the 39‐year period. The smallest glacier, Teepe, experienced the greatest area loss (60 ± 3%), whereas the largest glacier, Teton Glacier, lost 17 ± 3% of the 1967 area. For the current research, aerial photography from 1967 to 2002 was used to estimate glacier volume loss using stereoscopy techniques. The aerial photographs provide a finer resolution when compared with other datasets including satellite imagery (e.g., Landsat). Volume loss for the three glaciers was estimated to be 3.20 ± 0.46 million cubic meters over the period of 1967 to 2002. In assessing the primary climatic driver of the glacier ice loss, observed summer (June, July, and August) temperature data showed a statistically significant increase in temperatures when comparing the period of study (1968 to 2006) with historical temperatures from 1911 to 1967. When comparing spring (April 1st Snow Water Equivalent) snowpack for the period of study with historical records beginning in 1931, a significant difference in snowpack was not observed.     

Wind‐Powered Desalination: An Estimate of Saline Groundwater in the United States1

James Androwski, Abraham Springer, Thomas Acker, and Mark Manone, 2011. Wind‐Powered Desalination: An Estimate of Saline Groundwater in the United States. Journal of the American Water Resources Association (JAWRA) 47(1):93‐102. DOI: 10.1111/j.1752‐1688.2010.00493.x Abstract: Increasing scarcity of freshwater resources in many regions of the world is leading water resource managers to consider desalination as a potential alternative to traditional freshwater supplies. Desalination technologies are energy intensive and expensive to implement making desalination using renewable energy resources a potentially attractive option. Unfortunately, saline groundwater resources are not well characterized for many regions hindering consideration of such technologies. In this assessment, we estimate the saline groundwater resources of the principal aquifers of the United States using a geographic information system and correlate these resources to wind resources potentially sufficient to supply the energy demand of desalination equipment. We estimate that 3.1 × 10¹⁴ m³ saline groundwater, total volume, are contained in 28 of the country’s principal aquifers known to contain saline groundwater. Of this volume, 1.4 × 10¹⁴ m³ saline groundwater are co‐located with wind resources sufficient for electrical generation to desalinate groundwater.     

Implications of Upstream Flow Availability for Watershed Surface Water Supply across the Conterminous United States

Although it is well established that the availability of upstream flow (AUF) affects downstream water supply, its significance has not been rigorously categorized and quantified at fine resolutions. This study aims to fill this gap by providing a nationwide inventory of AUF and local water resource, and assessing their roles in securing water supply across the 2,099 8‐digit hydrologic unit code watersheds in the conterminous United States (CONUS). We investigated the effects of river hydraulic connectivity, climate variability, and water withdrawal, and consumption on water availability and water stress (ratio of demand to supply) in the past three decades (i.e., 1981–2010). The results show that 12% of the CONUS land relied on AUF for adequate freshwater supply, while local water alone was sufficient to meet the demand in another 74% of the area. The remaining 14% highly stressed area was mostly found in headwater areas or watersheds that were isolated from other basins, where stress levels were more sensitive to climate variability. Although the constantly changing water demand was the primary cause of escalating/diminishing stress, AUF variation could be an important driver in the arid south and southwest. This research contributes to better understanding of the significance of upstream–downstream water nexus in regional water availability, and this becomes more crucial under a changing climate and with intensified human activities.     

Evaluation of Satellite‐Derived Rainfall Data for Multiple Physio‐Climatic Regions in the Santiago River Basin,Mexico

Assessment of water resources requires reliable rainfall data, and rain gauge networks may not provide adequate spatial representation due to limited point measurements. The Tropical Rainfall Measuring Mission (TRMM) provides rainfall data at global scale, and has been used with good results. However, TRMM data are an indirect measurement of rainfall, and therefore must be validated for its proper use. In this work, a validation scheme was designed and implemented to compare the TRMM Version 7 (V7) monthly rainfall product at different time frames with data measured in two hydrologic subregions of the Santiago River Basin (SRB) in Mexico: Río Alto Santiago and Río Bajo Santiago (RBS). Additionally, three physio‐climatic regions provide an assessment of the interplay of topography, distance from coastal regions, and seasonal weather patterns on the correspondence between both datasets. The TRMM V7 rainfall product exhibited good agreement with the rain gauge data particularly for the RBS and for the whole SRB during wettest summer and autumn seasons. However, strong regional dependence was observed due to differences in climate and topography. Overall, in spite of some noted underestimations, the monthly TRMM V7 rainfall product was found to provide useful information that can be used to complement limited monitoring as is the case of RBS. An improved combined rainfall product could be generated and thus gaining the most benefits from both data sources.     

Trends in Pesticide Concentrations in Streamsof the Western United States, 1993‐20051

Johnson, Henry M., Joseph L. Domagalski, and Dina K. Saleh, 2010. Trends in Pesticide Concentrations in Streams of the Western United States, 1993‐2005. Journal of the American Water Resources Association (JAWRA) 00(0):1‐22. DOI: 10.1111/j.1752‐1688.2010.00507.x Abstract: Trends in pesticide concentrations for 15 streams in California, Oregon, Washington, and Idaho were determined for the organophosphate insecticides chlorpyrifos and diazinon and the herbicides atrazine, s‐ethyl diproplythiocarbamate (EPTC), metolachlor, simazine, and trifluralin. A parametric regression model was used to account for flow, seasonality, and antecedent hydrologic conditions and thereby estimate trends in pesticide concentrations in streams arising from changes in use amount and application method in their associated catchments. Decreasing trends most often were observed for diazinon, and reflect the shift to alternative pesticides by farmers, commercial applicators, and homeowners because of use restrictions and product cancelation. Consistent trends were observed for several herbicides, including upward trends in simazine at urban‐influenced sites from 2000 to 2005, and downward trends in atrazine and EPTC at agricultural sites from the mid‐1990s to 2005. The model provided additional information about pesticide occurrence and transport in the modeled streams. Two examples are presented and briefly discussed: (1) timing of peak concentrations for individual compounds varied greatly across this geographic gradient because of different application periods and the effects of local rain patterns, irrigation, and soil drainage and (2) reconstructions of continuous diazinon concentrations at sites in California are used to evaluate compliance with total maximum daily load targets.     

Sensitivity of Thermal Habitat of a Trout Stream to Potential Climate Change,Wisconsin, United States1

Abstract: Airborne thermal remote sensing from four flights on a single day from a single‐engine airplane was used to collect thermal infrared data of a 10.47‐km reach of the upper East Branch Pecatonica River in southwest Wisconsin. The study uses a one‐dimensional stream temperature model calibrated with the longitudinal profiles of stream temperature created from the four thermal imaging flights and validated with three days of continuous stream temperature data from instream data loggers on the days surrounding the thermal remote‐sensing campaign. Model simulations were used to quantify the sensitivity of stream thermal habitat to increases in air and groundwater temperature and changes in base flow. The simulations indicate that stream temperatures may reach critical maximum thresholds for brook trout (Salvelinus fontinalis) and brown trout (Salmo trutta) mortality, particularly if both air temperature increases and base flow declines. The approach demonstrates that thermal infrared data can greatly assist stream temperature model validation due to its high spatial resolution, and that this spatially continuous stream temperature data can be used to pinpoint spatial heterogeneity in groundwater inflow to streams. With this spatially distributed data on thermal heterogeneity and base‐flow accretion, stream temperature models considering various climate change scenarios are able to identify thermal refugia that will be critical for fisheries management under a changing climate.     

Estimation of Evapotranspiration Across the Conterminous United States Using a Regression With Climate and Land‐Cover Data1

Sanford, Ward E. and David L. Selnick, 2012. Estimation of Evapotranspiration Across the Conterminous United States Using a Regression with Climate and Land‐Cover Data. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12010 Abstract: Evapotranspiration (ET) is an important quantity for water resource managers to know because it often represents the largest sink for precipitation (P) arriving at the land surface. In order to estimate actual ET across the conterminous United States (U.S.) in this study, a water‐balance method was combined with a climate and land‐cover regression equation. Precipitation and streamflow records were compiled for 838 watersheds for 1971‐2000 across the U.S. to obtain long‐term estimates of actual ET. A regression equation was developed that related the ratio ET/P to climate and land‐cover variables within those watersheds. Precipitation and temperatures were used from the PRISM climate dataset, and land‐cover data were used from the USGS National Land Cover Dataset. Results indicate that ET can be predicted relatively well at a watershed or county scale with readily available climate variables alone, and that land‐cover data can also improve those predictions. Using the climate and land‐cover data at an 800‐m scale and then averaging to the county scale, maps were produced showing estimates of ET and ET/P for the entire conterminous U.S. Using the regression equation, such maps could also be made for more detailed state coverages, or for other areas of the world where climate and land‐cover data are plentiful.     

Sources and Delivery of Nutrients to the Northwestern Gulf of Mexico from Streams in the South‐Central United States1

Rebich, Richard A., Natalie A. Houston, Scott V. Mize, Daniel K. Pearson, Patricia B. Ging, and C. Evan Hornig, 2011. Sources and Delivery of Nutrients to the Northwestern Gulf of Mexico From Streams in the South‐Central United States. Journal of the American Water Resources Association (JAWRA) 47(5):1061‐1086. DOI: 10.1111/j.1752‐1688.2011.00583.x Abstract: SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were developed to estimate nutrient inputs [total nitrogen (TN) and total phosphorus (TP)] to the northwestern part of the Gulf of Mexico from streams in the South‐Central United States (U.S.). This area included drainages of the Lower Mississippi, Arkansas‐White‐Red, and Texas‐Gulf hydrologic regions. The models were standardized to reflect nutrient sources and stream conditions during 2002. Model predictions of nutrient loads (mass per time) and yields (mass per area per time) generally were greatest in streams in the eastern part of the region and along reaches near the Texas and Louisiana shoreline. The Mississippi River and Atchafalaya River watersheds, which drain nearly two‐thirds of the conterminous U.S., delivered the largest nutrient loads to the Gulf of Mexico, as expected. However, the three largest delivered TN yields were from the Trinity River/Galveston Bay, Calcasieu River, and Aransas River watersheds, while the three largest delivered TP yields were from the Calcasieu River, Mermentau River, and Trinity River/Galveston Bay watersheds. Model output indicated that the three largest sources of nitrogen from the region were atmospheric deposition (42%), commercial fertilizer (20%), and livestock manure (unconfined, 17%). The three largest sources of phosphorus were commercial fertilizer (28%), urban runoff (23%), and livestock manure (confined and unconfined, 23%).