USGS ScienceSearch

SEARCH · USGS Science

Results for “Research Square”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11Linked to original sources

Moon meteoritic seismic hum: Steady state prediction

We use three different statistical models describing the frequency of meteoroid impacts on Earth to estimate the seismic background noise due to impacts on the lunar surface. Because of diffraction, seismic events on the Moon are typically characterized by long codas, lasting 1 h or more. We find that the small but frequent impacts generate seismic signals whose codas overlap in time, resulting in a permanent seismic noise that we term the "lunar hum" by analogy with the Earth's continuous seismic background seismic hum. We find that the Apollo era impact detection rates and amplitudes are well explained by a model that parameterizes (1) the net seismic impulse due to the impactor and resulting ejecta and (2) the effects of diffraction and attenuation. The formulation permits the calculation of a composite waveform at any point on the Moon due to simulated impacts at any epicentral distance. The root-mean-square amplitude of this waveform yields a background noise level that is about 100 times lower than the resolution of the Apollo long-period seismometers. At 2 s periods, this noise level is more than 1000 times lower than the low noise model prediction for Earth's microseismic noise. Sufficiently sensitive seismometers will allow the future detection of several impacts per day at body wave frequencies. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Evaluation of procedures for prediction of unconventional gas in the presence of geologic trends

This study extends the application of local spatial nonparametric prediction models to the estimation of recoverable gas volumes in continuous-type gas plays to regimes where there is a single geologic trend. A transformation is presented, originally proposed by Tomczak, that offsets the distortions caused by the trend. This article reports on numerical experiments that compare predictive and classification performance of the local nonparametric prediction models based on the transformation with models based on Euclidean distance. The transformation offers improvement in average root mean square error when the trend is not severely misspecified. Because of the local nature of the models, even those based on Euclidean distance in the presence of trends are reasonably robust. The tests based on other model performance metrics such as prediction error associated with the high-grade tracts and the ability of the models to identify sites with the largest gas volumes also demonstrate the robustness of both local modeling approaches. ?? International Association for Mathematical Geology 2009.

Natural Resources Research

Determination of three‐dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model

Geiger's method of locating local earthquakes has been extended to include the effect of P velocity variation along the ray paths in three dimensions. The crustal structure was modeled by rectangular blocks, and a parameter was assigned to each block describing the perturbation of P wave slowness in the block. On the basis of an initial model, a set of linear equations for the observed first P arrival times was formulated in terms of the source and medium parameters. The source parameters for all the earthquakes in the data set and medium parameters for all the blocks penetrated by the seismic rays were then determined simultaneously by the damped least squares method. A computer program has been written for the case of a homogeneous initial medium model with constant P velocity. It has been tested with artificial data and applied to a set of local earthquakes recorded by a dense seismic array in Bear Valley, California. The resultant velocity distribution in the top 5 km shows a narrow low‐velocity zone of about 5 km/s in the San Andreas fault zone sandwiched between high‐velocity regions of about 6 km/s.

California

An analysis of autocorrelation and bias in home range estimation

Home range estimation is routine practice in ecological research. While advances in animal tracking technology have increased our capacity to collect data to support home range analysis, these same advances have also resulted in increasingly autocorrelated data. Consequently, the question of which home range estimator to use on modern, highly autocorrelated tracking data remains open. This question is particularly relevant given that most estimators assume independently sampled data. Here, we provide a comprehensive evaluation of the effects of autocorrelation on home range estimation. We base our study on an extensive data set of GPS locations from 369 individuals representing 27 species distributed across five continents. We first assemble a broad array of home range estimators, including Kernel Density Estimation (KDE) with four bandwidth optimizers (Gaussian reference function, autocorrelated-Gaussian reference function AKDE, Silvermans rule of thumb, and least squares cross-validation), Minimum Convex Polygon, and Local Convex Hull methods. Notably, all of these estimators except AKDE assume independent and identically distributed (IID) data. We then employ half-sample cross-validation to objectively quantify estimator performance, and the recently introduced effective sample size for home range area estimation ( N̂ area ) to quantify the information content of each data set. We found that AKDE 95% area estimates were larger than conventional IID-based estimates by a mean factor of 2. The median number of cross-validated locations included in the hold-out sets by AKDE 95% (or 50%) estimates was 95.3% (or 50.1%), confirming the larger AKDE ranges were appropriately selective at the specified quantile. Conversely, conventional estimates exhibited negative bias that increased with decreasing N̂ area . To contextualize our empirical results, we performed a detailed simulation study to tease apart how sampling frequency, sampling duration, and the focal animals movement conspire to affect range estimates. Paralleling our empirical results, the simulation study demonstrated that AKDE was generally more accurate than conventional methods, particularly for small N̂ area . While 72% of the 369 empirical data sets had >1,000 total observations, only 4% had an N̂ area >1,000, where 30% had an N̂ area <30. In this frequently encountered scenario of small N̂ area , AKDE was the only estimator capable of producing an accurate home range estimate on autocorrelated data.

Ecological Monographs

Geomagnetic local and regional harmonic analyses

Procedures are developed for using rectangular and cylindrical harmonic analyses in local and regional areas. Both the linear least squares analysis, applicable when component data are available, and the nonlinear least squares analysis, applicable when only total field data are available, are treated. When component data are available, it is advantageous to work with residual fields obtained by subtracting components derived from a harmonic potential from the observed components. When only total field intensity data are available, they must be used directly. Residual values cannot be used. Cylindrical harmonic analyses are indicated when fields tend toward cylindrical symmetry; otherwise, rectangular harmonic analyses will be more advantageous. Examples illustrating each type of analysis are given.

Journal of Geophysical Research Solid Earth

Maps of active layer thickness in northern Alaska by upscaling P-band polarimetric synthetic aperture radar retrievals

Extensive, detailed information on the spatial distribution of active layer thickness (ALT) in northern Alaska and how it evolves over time could greatly aid efforts to assess the effects of climate change on the region and also help to quantify greenhouse gas emissions generated due to permafrost thaw. For this reason, we have been developing high-resolution maps of ALT throughout northern Alaska. The maps are produced by upscaling from high-resolution swaths of estimated ALT retrieved from airborne P-band synthetic aperture radar (SAR) images collected for three different years. The upscaling was accomplished by using hundreds of thousands of randomly selected samples from the SAR-derived swaths of ALT to train a machine learning regression algorithm supported by numerous spatial data layers. In order to validate the maps, thousands of randomly selected samples of SAR-derived ALT were excluded from the training in order to serve as validation pixels; error performance calculations relative to these samples yielded root-mean-square errors (RMSEs) of 7.5–9.1 cm, with bias errors of magnitude under 0.1 cm. The maps were also compared to ALT measurements collected at a number of in situ test sites; error performance relative to the site measurements yielded RMSEs of approximately 11–12 cm and bias of 2.7–6.5 cm. These data are being used to investigate regional patterns and underlying physical controls affecting permafrost degradation in the tundra biome.

Environmental Research Letters

Seasonal cultivated and fallow cropland mapping using MODIS-based automated cropland classification algorithm

Increasing drought occurrences and growing populations demand accurate, routine, and consistent cultivated and fallow cropland products to enable water and food security analysis. The overarching goal of this research was to develop and test automated cropland classification algorithm (ACCA) that provide accurate, consistent, and repeatable information on seasonal cultivated as well as seasonal fallow cropland extents and areas based on the Moderate Resolution Imaging Spectroradiometer remote sensing data. Seasonal ACCA development process involves writing series of iterative decision tree codes to separate cultivated and fallow croplands from noncroplands, aiming to accurately mirror reliable reference data sources. A pixel-by-pixel accuracy assessment when compared with the U.S. Department of Agriculture (USDA) cropland data showed, on average, a producer’s accuracy of 93% and a user’s accuracy of 85% across all months. Further, ACCA-derived cropland maps agreed well with the USDA Farm Service Agency crop acreage-reported data for both cultivated and fallow croplands with R-square values over 0.7 and field surveys with an accuracy of ≥95% for cultivated croplands and ≥76% for fallow croplands. Our results demonstrated the ability of ACCA to generate cropland products, such as cultivated and fallow cropland extents and areas, accurately, automatically, and repeatedly throughout the growing season.

California

Integration of remote sensing and field observations in evaluating DSSAT model for estimating maize and soybean growth and yield in Maryland, USA

Crop models are useful for evaluating crop growth and yield at the field and regional scales, but their applications and accuracies are restricted by input data availability and quality. To overcome difficulties inherent to crop modeling, input data can be enhanced by the incorporation of remotely sensed and field observations into crop growth models. This approach has been recognized to be an important way to monitor crop growth conditions and to predict yield at the field and regional scale. In recent years, satellite remote sensing has provided high-temporal and high-spatial-resolution data that allow for generating continuous time series of biophysical parameters such as vegetation indices, leaf area index, and phenology. The objectives of this study were to use remote sensing along with field observations as inputs to the Decision Support System for Agro-Technology (DSSAT) model to estimate soybean and maize growth and yield. The study used phenology and leaf area index (LAI) data derived from Planet Fusion (daily, 3 m) satellite imagery along with field observation data on crop growth stage, LAI and yield collected at the United State Department of Agriculture, Agricultural Research Service, Beltsville Agricultural Research Center (BARC), Beltsville, Maryland. For maize, a total of 17 treatments (site years) were used (ten treatments for model calibration and seven treatments for validation), while for soybean (maturity groups three and four), a total of 18 treatments were used (nine for calibration and nine for validation). The calibrated model was tested against an independent, multi-location and multi-year set of phenology and yield data (2017–2020) from BARC fields. The model accurately simulated maize and soybean days to flowering and maturity and produced reasonable yield estimates for most fields and years. Model run for independent locations and years produced good results for phenology and yields for both maize and soybean, as indicated by index of agreement (d) values ranging from 0.65 to 0.93 and normalized root-mean-squared error values ranging from 1 to 20%, except for soybean maturity group four. Overall, model performances with respect to phenology and grain yield for maize and soybean were good and consistent with other DSSAT evaluation studies. The inclusion of remote sensing along with field observations in crop-growth model inputs can provide an effective approach for assessing crop conditions, even in regions lacking ground data.

Maryland

Carbon dioxide emission rate of Kīlauea Volcano: Implications for primary magma and the summit reservoir

We report a CO 2 emission rate of 8500 metric tons per day (t d &minus;1 ) for the summit of Kīlauea Volcano, several times larger than previous estimates. It is based on three sets of measurements over 4 years of synchronous SO 2 emission rates and volcanic CO 2 /SO 2 concentration ratios for the summit correlation spectrometer (COSPEC) traverse. Volcanic CO 2 /SO 2 for the traverse is representative of the global ratio for summit emissions. The summit CO 2 emission rate is nearly constant, despite large temporal variations in summit CO 2 /SO 2 and SO 2 emission rates. Summit CO 2 emissions comprise most of Kīlauea's total CO 2 output (&sim;9000 t d &minus;1 ). The bulk CO 2 content of primary magma determined from CO 2 emission and magma supply rate data is &sim;0.70 wt %. Most of the CO 2 is present as exsolved vapor at summit reservoir depths, making the primary magma strongly buoyant. Turbulent mixing with resident reservoir magma, however, prevents frequent eruptions of buoyant primary magma in the summit region. CO 2 emissions confirm that the magma supply enters the edifice through the summit reservoir. A persistent several hundred parts per million CO 2 anomaly arises from the entry of magma into the summit reservoir beneath a square kilometer area east of Halemaumau pit crater. Since most of the CO 2 in primary magma is degassed in the summit, the summit CO 2 emission rate is an effective proxy for the magma supply rate. Both scrubbing of SO 2 and solubility controls on CO 2 and S in basaltic melt cause high CO 2 /SO 2 in summit emissions and spatially uncorrelated distributions of CO 2 and SO 2 in the summit plume.

Journal of Geophysical Research B: Solid Earth

National Water-Quality Assessment Program, western Lake Michigan drainages: Summaries of liaison committee meeting, Green Bay, Wisconsin, March 28-29, 1995

The Western Lake Michigan Drainages (WMIC) study unit, under investigation since 1991, drains 20,000 square miles (mi2) in eastern Wisconsin and Upper Michigan (fig. 1). The major water-quality issues in the WMIC study unit are: (1) nonpoint-source contamination of surface and ground water by agricultural chemicals, (2) contamination in bottom sediments of rivers and harbors by toxic substances, including polychlorinated biphenyls (PCB's), other synthetic organic compounds, and trace elements, (3) nutrient enrichment of rivers and lakes resulting from nonpoint- and point-source discharges, and (4) acidification and mercury contamination of lakes in poorly buffered watersheds in the northwestern part of the study unit. A study-unit liaison committee, which includes representatives of Federal, State, university, and private and citizen organizations, has met annually since 1991 to review plans and results and guide the investigators toward policy-relevant efforts. The results of research conducted in the WMIC study unit by U.S. Geological Survey (USGS) and non-USGS researchers were presented at the liaison committee meeting held in Green Bay, Wis., on March 28-29, 1995. This report contains summaries of the oral presentations given at the WMIC 1995 liaison committee meeting.

Michigan, Wisconsin

Subsidence and carbon fluxes in the Sacramento/San Joaquin Delta, California

The Sacramento/San Joaquin Delta was once tidal marshland and was formed by the confluence of the Sacramento and the San Joaquin Rivers. The delta covers more than 1,000 square miles and is bordered by the major population centers of the San Francisco Bay area, Sacramento, and Stockton (fig. 1). The delta took its current form by the 1930's when drainage of 100 islands and tracts and construction of 2,250 miles of levees were completed. The land surface on the delta islands has progressively subsided below sea level since the early 20th century. The integrity of the levee network is threatened because subsidence increases the differences in water levels maintained by the levees, which can cause levee failure and consequent flooding of the delta islands. Levee failure affects the quality and beneficial use of water in the delta, which is the major source of California's water supply to millions of residents and acres of agricultural land. This Fact Sheet briefly summarizes the research done by the U.S. Geological Survey, in cooperation with the California Department of Water Resources to assess the causes of subsidence, processes that affect subsidence, and carbon fluxes in the delta.

California

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface

Direct geoelectrical evidence of mass transfer at the laboratory scale

Previous field-scale experimental data and numerical modeling suggest that the dual-domain mass transfer (DDMT) of electrolytic tracers has an observable geoelectrical signature. Here we present controlled laboratory experiments confirming the electrical signature of DDMT and demonstrate the use of time-lapse electrical measurements in conjunction with concentration measurements to estimate the parameters controlling DDMT, i.e., the mobile and immobile porosity and rate at which solute exchanges between mobile and immobile domains. We conducted column tracer tests on unconsolidated quartz sand and a material with a high secondary porosity: the zeolite clinoptilolite. During NaCl tracer tests we collected nearly colocated bulk direct-current electrical conductivity (σ b ) and fluid conductivity (σ f ) measurements. Our results for the zeolite show (1) extensive tailing and (2) a hysteretic relation between σ f and σ b , thus providing evidence of mass transfer not observed within the quartz sand. To identify best-fit parameters and evaluate parameter sensitivity, we performed over 2700 simulations of σ f , varying the immobile and mobile domain and mass transfer rate. We emphasized the fit to late-time tailing by minimizing the Box-Cox power transformed root-mean square error between the observed and simulated σ f . Low-field proton nuclear magnetic resonance (NMR) measurements provide an independent quantification of the volumes of the mobile and immobile domains. The best-fit parameters based on σ f match the NMR measurements of the immobile and mobile domain porosities and provide the first direct electrical evidence for DDMT. Our results underscore the potential of using electrical measurements for DDMT parameter inference.

Water Resources Research

Morphologic dating of fault scarps using airborne laser swath mapping (ALSM) data

Models of fault scarp morphology have been previously used to infer the relative age of different fault scarps in a fault zone using labor-intensive ground surveying. We present a method for automatically extracting scarp morphologic ages within high-resolution digital topography. Scarp degradation is modeled as a diffusive mass transport process in the across-scarp direction. The second derivative of the modeled degraded fault scarp was normalized to yield the best-fitting (in a least-squared sense) scarp height at each point, and the signal-to-noise ratio identified those areas containing scarp-like topography. We applied this method to three areas along the San Andreas Fault and found correspondence between the mapped geometry of the fault and that extracted by our analysis. This suggests that the spatial distribution of scarp ages may be revealed by such an analysis, allowing the recent temporal development of a fault zone to be imaged along its length.

Geophysical Research Letters

Stream temperature prediction in a shifting environment: The influence of deep learning architecture

Stream temperature is a fundamental control on ecosystem health. Recent efforts incorporating process guidance into deep learning models for predicting stream temperature have been shown to outperform existing statistical and physical models. This performance is in part because deep learning architectures can actively learn spatiotemporal relationships that govern how water and energy propagate through a river network. However, exploration of how spatiotemporal awareness and process guidance influence a model's generalizability under shifting environmental conditions such as climate change is limited. Here, we use Explainable Artificial Intelligence (XAI) to interrogate how differing deep learning architectures affect a model's learned spatial and temporal dependencies, and how those learned dependencies affect a model's ability to maintain high accuracy when applied to unseen environmental conditions. Using the Delaware River Basin in the northeastern United States as a test case, we compare two spatiotemporally aware process-guided deep learning models for predicting stream temperature (a recurrent graph convolution network—RGCN, and a temporal convolution graph model—Graph WaveNet). Both models achieve equally high predictive performance when testing data are well represented in the training data (test root mean squared errors of 1.64°C and 1.65°C); however, Graph WaveNet significantly outperforms RGCN in 4 out of 5 experiments where test partitions represent different types of unseen environmental conditions. XAI results show that the architecture of Graph WaveNet leads to learned spatial relationships with greater fidelity to physical processes, and that this fidelity improves the generalizability of the model when applied to shifting and/or unseen environmental conditions.

Delaware River Basin

Freshwater DOM quantity and quality from a two-component model of UV absorbance

We present a model that considers UV-absorbing dissolved organic matter (DOM) to consist of two components (A and B), each with a distinct and constant spectrum. Component A absorbs UV light strongly, and is therefore presumed to possess aromatic chromophores and hydrophobic character, whereas B absorbs weakly and can be assumed hydrophilic. We parameterised the model with dissolved organic carbon concentrations [DOC] and corresponding UV spectra for c. 1700 filtered surface water samples from North America and the United Kingdom, by optimising extinction coefficients for A and B, together with a small constant concentration of non-absorbing DOM (0.80 mg DOC L -1 ). Good unbiased predictions of [DOC] from absorbance data at 270 and 350 nm were obtained ( r 2 = 0.98), the sum of squared residuals in [DOC] being reduced by 66% compared to a regression model fitted to absorbance at 270 nm alone. The parameterised model can use measured optical absorbance values at any pair of suitable wavelengths to calculate both [DOC] and the relative amounts of A and B in a water sample, i.e. measures of quantity and quality. Blind prediction of [DOC] was satisfactory for 9 of 11 independent data sets (181 of 213 individual samples).

Water Research

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth

Nutrient concentrations and yields in undeveloped stream basins of the United States

Data from 85 sites across the United States were used to estimate concentrations and yields of selected nutrients in streams draining relatively undeveloped basins. Flow-weighted concentrations during 1990-1995 were generally low with median basin concentrations of 0.020, 0.087, 0.26, 0.010, and 0.022 milligrams per liter (mg/L) for ammonia as N, nitrate as N, total nitrogen, orthophosphate as P, and total phosphorus, respectively. The flow-weighted concentration of nitrate exceeded 0.6 mg/L in only three basins. Total nitrogen exceeded 1 mg/L in only four basins, and total phosphorus exceeded 0.1 mg/L in only four basins. The median annual basin yield of ammonia as N, nitrate as N, total nitrogen, orthophosphate as P, and total phosphorus was 8.1, 26, 86, 2.8, and 8.5 kilograms per square kilometer, respectively. Concentrations and yields of nitrate tended to be highest in northeastern and mid-Atlantic coastal states and correlated well with areas of high atmospheric nitrogen deposition. Concentrations and yields of total nitrogen were highest in the southeastern part of the nation and in parts of the upper Midwest. In the northeast, nitrate was generally the predominant form of nitrogen, and in the southeast and parts of the upper Midwest, organic nitrogen was the dominant form. Concentrations of total phosphorus were generally highest in the Rocky Mountain and Central Plain states.Data from 85 sites across the United States were used to estimate concentrations and yields of selected nutrients in streams draining relatively undeveloped basins. Flow-weighted concentrations during 1990-1995 were generally low with median basin concentrations of 0.020, 0.087, 0.26, 0.010, and 0.022 milligrams per liter (mg/L) for ammonia as N, nitrate as N, total nitrogen, orthophosphate as P, and total phosphorus, respectively. The flow-weighted concentration of nitrate exceeded 0.6 mg/L in only three basins, Total nitrogen exceeded 1 mg/L in only four basins, and total phosphorus exceeded 0.1 mg/L in only four basins. The median annual basin yield of ammonia as N, nitrate as N, total nitrogen, orthophosphate as P, and total phosphorus was 8.1, 26, 86, 2.8, and 8.5 kilograms per square kilometer, respectively. Concentrations and yields of nitrate tended to be highest in northeastern and mid-Atlantic coastal states and correlated well with areas of high atmospheric nitrogen deposition. Concentrations and yields of total nitrogen were highest in the southeastern part of the nation and in parts of the upper Midwest. In the northeast, nitrate was generally the predominant form of nitrogen, and in the southeast and parts of the upper Midwest, organic nitrogen was the dominant form. Concentrations of total phosphorus were generally highest in the Rocky Mountain and Central Plain states.Data collected across the US from 85 streams draining relatively undeveloped basins were used to identify broad regional and national patterns in nutrient concentrations and yields. The basins of interest were selected from three USGS programs: the Hydrologic Benchmark Network, the National Water-Quality Assessment, and the Research Program. Water samples from most basins were collected on a weekly to bimonthly schedule. While the flow-weighted concentrations of nutrients varied, concentrations were low in most basins. Median flow-weighted concentrations of ammonia, nitrate, total nitrogen, orthophosphate, and total phosphorus were 0.020, 0.087, 0.26, 0.010, and 0.022 mg/l, respectively. Nitrate concentrations tended to be highest in the northeastern US, while southeastern and north-central basins had some of the highest NH3 concentrations. Flow-weighted concentrations of total P were generally highest in the Rocky Mountain and Central Plain states and in parts of the southeast.

Journal of the American Water Resources Associatio