USGS ScienceSearch

SEARCH · USGS Science

Results for “Carbon Balance and Management”

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.

72 records · Page 4Linked to original sources

From pools to flow: The PROMISE framework for new insights on soil carbon cycling in a changing world

Soils represent the largest terrestrial reservoir of organic carbon, and the balance between soil organic carbon (SOC) formation and loss will drive powerful carbon‐climate feedbacks over the coming century. To date, efforts to predict SOC dynamics have rested on pool‐based models, which assume classes of SOC with internally homogenous physicochemical properties. However, emerging evidence suggests that soil carbon turnover is not dominantly controlled by the chemistry of carbon inputs, but rather by restrictions on microbial access to organic matter in the spatially heterogeneous soil environment. The dynamic processes that control the physicochemical protection of carbon translate poorly to pool‐based SOC models; as a result, we are challenged to mechanistically predict how environmental change will impact movement of carbon between soils and the atmosphere. Here, we propose a novel conceptual framework to explore controls on belowground carbon cycling: P robabilistic R epresentation of O rganic M atter I nteractions within the S oil E nvironment (PROMISE). In contrast to traditional model frameworks, PROMISE does not attempt to define carbon pools united by common thermodynamic or functional attributes. Rather, the PROMISE concept considers how SOC cycling rates are governed by the stochastic processes that influence the proximity between microbial decomposers and organic matter, with emphasis on their physical location in the soil matrix. We illustrate the applications of this framework with a new biogeochemical simulation model that traces the fate of individual carbon atoms as they interact with their environment, undergoing biochemical transformations and moving through the soil pore space. We also discuss how the PROMISE framework reshapes dialogue around issues related to SOC management in a changing world. We intend the PROMISE framework to spur the development of new hypotheses, analytical tools, and model structures across disciplines that will illuminate mechanistic controls on the flow of carbon between plant, soil, and atmospheric pools.

Global Change Biology

The effects of cheatgrass invasion on US Great Basin carbon storage depend on interactions between plant community composition, precipitation seasonality, and soil climate regime

Annual-grass invasions are transforming desert ecosystems in ways that affect ecosystem carbon (C) balance, but previous studies do not agree on the pattern, magnitude and direction of changes. A recent meta-analysis of 41 articles and 386 sites concludes that invasion by annual grasses such as cheatgrass ( Bromus tectorum L) reduces C in biomass across the Great Basin (Nagy et al., 2021). Reanalysis reveals that whether cheatgrass affects biomass C stocks is not generalizable, but rather depends on the considerable variation in climate across the subject sites. Our analysis suggests that accurate Great Basin-scale estimates of cheatgrass effects on C balance are not yet possible. Addition of climate variables to the meta-analysis reveals that cheatgrass invasion (a) reduced C in above-ground biomass in relatively summer-wet sites but not in summer-dry sites, (b) increased surface soil C in sites with intermediate resistance and resilience classifications (R&R) but not in low R&R sites—that is, mesic/aridic soil climates and (c) did not affect deep soil C. Considering that cheatgrass has expanded most in relatively summer-dry sites and mesic/aridic sites, omission of climate factors leads to model overestimates of cheatgrass effects on C when extrapolating to larger areas. Estimates of cheatgrass effects on C would also be improved if the analysis considered that (a) perennial grasslands are a common community state in the Great Basin that have intermediary C relative to annual grasslands and sagebrush stands, that is the omission of perennial grasslands from analysis inflates the baseline C storage of uninvaded Great Basin ecosystems, and( b) cheatgrass does not often exist in stable monocultures and soil carbon can reflect current or recent presence of other species. Synthesis and applications . Invasions often reveal heterogeneity in ecosystem structure and function that is not otherwise evident, and the heterogeneity can influence estimation of the net impacts of the invaders. For cheatgrass and other invaders, we propose that formally accounting for the spatial variability of invasion on ecosystem functions will improve the estimation of their net effect on ecosystem C, and thus improve prospects for adjusting management practices to optimize C sequestration.

western United States

Dynamics of land-use change and conservation in the Prairie Pothole Region of the United States: environmental and economic implications with linkages to rural community well-being

Rural America has changed dramatically over the last century, from having over half the population living in rural settings to only 20 percent residing in a rural area today, and outmigration of younger populations from rural communities remains a constant issue for local governing officials. A declining tax base and concurrent rising costs for maintenance and repair of aging infrastructure add further challenges to policy decisions. Reduced enrollment has caused school closures or mergers. Farm consolidation and technical advances reduced the demand for local labor. On the positive side, however, record-high commodity prices have amplified farm income to new heights. The increased revenues can lead to farmers spending additional money within the local region, while at the same time increased transportation of products has impacted local infrastructure such as roads and bridges. Such dynamics present challenges for municipal leaders charged with promoting economic development and balanced spending, while at the same time maintaining the way of life and rural character that are so important to area residents. The Prairie Pothole Region (PPR) of the United States covers much of the Northern Great Plains, including parts of North and South Dakota, Minnesota, Iowa, and a small part of Montana, and extends across a broad swath of Alberta and Saskatchewan. The region is defined largely by its rural character but has experienced extensive land conversion over the last century, with agricultural areas replacing native prairie habitat. Additional pressures arise from oil and gas development, global markets for agricultural production, and increased demands for biofuel feedstocks. Record-high commodity prices increase pressure on the native prairie as farmers look for new cropland acres. The volatility of commodity prices has raised fears over the intensity of land conversion to row-crop agriculture, the economic health and resiliency of rural communities, and ultimately, population dynamics and outmigration of younger generations. Land-use pressures are increased by the exponential growth of oil and gas production in the region, where some 8,200 wells are now in production within the Williston Basin of North Dakota, accompanied by increased population pressures on housing and municipal services. The U.S. Department of Agriculture’s Conservation Reserve Program (CRP)--a cropland retirement program with close to 4.8 million acres enrolled in the PPR--faces uncertainty in upcoming legislative actions, with a large majority of property enrollments scheduled to expire by 2017. The CRP historically has provided improved habitat conditions, reductions of soil damage through erosion and loss of nutrients, and sequestration of millions of tons of atmospheric carbon. In turn, wildlife-related recreation levels have increased in many parts of the PPR, with money spent in local communities. Contemporary resource-management and rural-development planning increasingly emphasize the need for diversification and integration of resource-extractive industries with nonmarket-based recreational and amenity values that tie into quality of life. Ultimately, each community is unique in its environmental, social, economic, and fiscal endowments. One rural-development policy may work better in one community than another. In addition, rural-development issues such as migration, job growth, and taxes are diverse in themselves. The goal of this report is to qualitatively and quantitatively discuss the economic impacts of land-use decisions in rural areas, particularly in the PPR.

Professional Paper

Water-budgets and recharge-area simulations for the Spring Creek and Nittany Creek Basins and parts of the Spruce Creek Basin, Centre and Huntingdon Counties, Pennsylvania, Water Years 2000–06

This report describes the results of a study by the U.S. Geological Survey in cooperation with ClearWater Conservancy and the Pennsylvania Department of Environmental Protection to develop a hydrologic model to simulate a water budget and identify areas of greater than average recharge for the Spring Creek Basin in central Pennsylvania. The model was developed to help policy makers, natural resource managers, and the public better understand and manage the water resources in the region. The Groundwater and Surface-water FLOW model (GSFLOW), which is an integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Groundwater Flow Model (MODFLOW-NWT), was used to simulate surface water and groundwater in the Spring Creek Basin for water years 2000–06. Because the groundwater and surface-water divides for the Spring Creek Basin do not coincide, the study area includes the Nittany Creek Basin and headwaters of the Spruce Creek Basin. The hydrologic model was developed by the use of a stepwise process: (1) develop and calibrate a PRMS model and steady-state MODFLOW-NWT model; (2) re-calibrate the steady-state MODFLOW-NWT model using potential recharge estimates simulated from the PRMS model, and (3) integrate the PRMS and MODFLOW-NWT models into GSFLOW. The individually calibrated PRMS and MODFLOW-NWT models were used as a starting point for the calibration of the fully coupled GSFLOW model. The GSFLOW model calibration was done by comparing observations and corresponding simulated values of streamflow from 11 streamgages and groundwater levels from 16 wells. The cumulative water budget and individual water budgets for water years 2000–06 were simulated by using GSFLOW. The largest source and sink terms are represented by precipitation and evapotranspiration, respectively. For the period simulated, a net surplus in the water budget was computed where inflows exceeded outflows by about 1.7 billion cubic feet (0.47 inches per year over the basin area); storage increased by about the same amount to balance the budget. The rate and distribution of recharge throughout the Spring Creek, Nittany Creek, and Spruce Creek Basins is variable as a result of the high degree of hydrogeologic heterogeneity and karst features. The greatest amount of recharge was simulated in the carbonate-bedrock valley, near the toe slopes of Nittany and Tussey Mountains, in the Scotia Barrens, and along the area coinciding with the Gatesburg Formation. Runoff extremes were observed for water years 2001 (dry year) and 2004 (wet year). Simulated average recharge rates (water reaching the saturated zone as defined in GSFLOW) for 2001 and 2004 were 5.4 in/yr and 22.0 in/yr, respectively. Areas where simulations show large variations in annual recharge between wet and dry years are the same areas where simulated recharge was large. Those areas where rates of groundwater recharge are much higher than average, and are capable of accepting substantially greater quantities of recharge during wet years, might be considered critical for maintaining the flow of springs, stream base flow, or the source of water to supply wells. The slopes of the Bald Eagle, Tussey, and Nittany Mountains are relatively insensitive to variations in recharge, primarily because of reduced infiltration rates and steep slopes.

Pennsylvania

Precipitation, Ground-water Hydrology, and Recharge Along the Eastern Slopes of the Sandia Mountains, Bernalillo County, New Mexico

The spatial and temporal distribution of recharge to carbonate and clastic aquifers along the eastern slopes of the Sandia Mountains was investigated by using precipitation, water-level, dissolved chloride, and specific-conductance data. The U.S. Geological Survey (USGS), in cooperation with the Bernalillo County Public Works Division, conducted a study to assess ground-water conditions and provide technical data that could be used as a basis for management and future planning of eastern Bernalillo County water resources. The intent of the investigation was to improve the current understanding of subsurface mechanisms controlling recharge dynamics in a geologically complex aquifer system. In the Sandia Mountains, precipitation events are generally limited to snowfalls in winter months and monsoon rainfall in late summer. Monthly meteorological data from weather stations in the study area indicate that monsoon rainfall during July and August constitutes close to one-third of annual precipitation totals. Following precipitation and snowmelt events, daily ground-water level data show low-amplitude, long-duration peaks in hydrographs of wells north and west of the Tijeras Fault. Hydrographs of monthly and biannual water-level data from across the study area show seasonal variation and water-level fluctuations in excess of 30 ft during a period of below-average precipitation. Water level observations in 67 percent of wells showing drought-induced water-level declines rebounded to at or near predrought conditions within 6 months of return to normal climate conditions. Cross-correlation of annual hydrologic data shows aquifer response to periods of monsoon recharge to persist from 1 to 6 months following events. The lag time between precipitation input and response of water levels or solute concentrations was largest near the Tijeras and Gutierrez Faults. These results indicate regional faults hydrologically isolate the Tijeras Graben from groundwater recharge originating at high elevations along the eastern slopes of the Sandia Mountains. Recharge rates calculated by using the chloride-mass-balance method for five springs located at the base of the Sandia Mountains ranged from 1 to 23 percent of annual precipitation.

Scientific Investigations Report

Simulation of flow and eutrophication in the central Salem River, New Jersey

The central Salem River in New Jersey is subject to periods of water-quality impairment, marked by elevated concentrations of phosphorus and chlorophyll- a , and low concentrations of and large diurnal swings in concentrations of dissolved oxygen. These seasonal eutrophic conditions are controlling factors for water quality in lower reaches, where the river is more lacustrine than in upper reaches, as a result of downstream damming. This biological productivity is supported by nutrient wash-off from agricultural areas in the surrounding watershed. To investigate this impairment, flow measurement and water-quality sampling were conducted during 2007–08 in support of development of a one-dimensional surface-water-quality model that simulates nutrient cycling and transformation processes. The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection, used the U.S. Environmental Protection Agency Water Quality Analysis Simulation Program (WASP) to develop a receiving-water-quality model of the central Salem River between Woodstown and Deepwater, New Jersey, from April 2007 to October 2008. The main-stem river and largest tributary were simulated. In the flow model, kinematic wave flow is used to simulate flow in upper reaches and ponded weir flow is used to simulate flow in lower reaches. The water-quality model makes use of a mass-balance equation to simulate the fate and transport of nutrients, phytoplankton chlorophyll- a , dissolved oxygen, and oxygen demands (an indicator rather than a substance) in the river. Model input included channel characteristics, boundary conditions for flow and water quality, environmental parameters, vertical dispersion coefficients, settling rates, and kinetic constants. Inputs were estimated where field data were lacking, notably for tributary flows and nutrient loads. The model was calibrated to observed flow variables and concentrations of dissolved oxygen, chlorophyll- a , and nutrients at sampling locations, with emphasis on growing-season conditions. Calibration was achieved through graphical and statistical comparison of simulated results to observed data. Sensitivity analyses were performed, and model limitations and applicability were evaluated. Simulated results closely matched observed data in most cases, although some were overpredicted slightly. The most important causes of overprediction were estimated tributary flows for the flow model and estimated tributary watershed loads for the water-quality model. Calibration of dissolved-oxygen concentrations was closer, and predicted diurnal variations were consistent with high algal photosynthesis/respiration, although lack of continuous dissolved-oxygen data precluded verifying these predictions. A similar caveat applies to predicted diurnal variations in chlorophyll- a . Simulated limitations on algal growth were consistent with those based on observed data and indicated phosphorus was the main limiting nutrient, except during certain periods when nitrogen was limiting. Two water-quality management scenarios were simulated with the model to assess the effect of point- and nonpoint-source nutrient reductions on water-quality conditions in the river. Scenarios involved (1) a return of watershed land use to predevelopment natural conditions and (2) an extreme reduction in nutrient input. Although the extreme-nutrient-reduction scenario yielded improvements in water quality, the natural-conditions scenario yielded the largest improvements as indicated by minimal violations of surface-water-quality standards or thresholds. However, years may be needed to attain the full benefit of these management scenarios as a result of accumulation of phosphorus and organic carbon in riverbed sediments in lacustrine reaches. The results of this study indicate that the quality of water in the central Salem River will improve if management policies that mitigate the effects of nutrient-loading practices in the watershed, particularly those related to agriculture, are implemented.

New Jersey

U.S. Geological Survey program of offshore resource and geoenvironmental studies, Atlantic-Gulf of Mexico region, from September 1, 1976, to December 31, 1978

Mineral and energy resources of the continental margins of the United States arc important to the Nation's commodity independence and to its balance of payments. These resources are being studied along the continental margins of the Atlantic Ocean and the Gulf of Mexico in keeping with the mission of the U.S. Geological Survey to survey the geologic structures, mineral resources, and products of the national domain.'(Organic Act of 1879). An essential corollary to these resource studies is the study of potential geologic hazards that may be associated with offshore resource exploration and exploitation. In cooperation with the U.S. Bureau of Land Management, the Geological Survey, through its Atlantic-Gulf of Mexico Marine Geology Program, carries out extensive research to evaluate hazards from sediment mobility, shallow gas, and slumping and to acquire information on the distribution and concentration of trace metals and biogenic and petroleum-derived hydrocarbons in sea-floor sediments. All these studies arc providing needed background information, including information on pollutant dispersal, on the nearshore, estuarine, and lacustrine areas that may be near pipeline and nuclear powerplant sites. Users of these data include the Congress, many Federal agencies, the coastal States, private industry, academia, and the concerned public. The results of the regional structural, stratigraphic, and resource studies carried out under the Atlantic-Gulf of Mexico Marine Geology Program have been used by the Geological Survey and the Bureau of Land Management to select areas for future leasing and to aid in the evaluation of tracts nominated for leasing. Resource studies have concentrated mostly on the Atlantic Outer Continental Shelf frontier areas. Geologic detailing of five major basins along the U.S. Atlantic margin, where sediments are as much as 14 km thick, have been revealed by 25,000 km of 24-and 48-channel common-depth-point seismic data, 187,000 km of acromagnetic data, and 39,000 km of gravity data, plus 10,000 samples and logs obtained from U.S. Geological Survey and industry drilling (for example, coreholes of the Atlantic Slope Program, Joint Oceanographic Institutions Deep Earth Sampling, Continental Offshore Stratigraphic Tests, and the Atlantic Margin Coring Program). A sedimentary section of Jurassic and Cretaceous age grades from terrigenous clastic rocks nearshore to carbonate rocks offshore; this section is part of an extensive buried bank-platform complex that could contain large reserves of natural gas and oil. The volume of sediment deposited offshore far exceeds the volume deposited onshore where extensive accumulations of oil, gas, and minerals have been found. Commercial exploratory drilling offshore thus far has been limited to the Baltimore Canyon Trough area off New Jersey, where at least two holes have found gas; leasing has taken place in the Southeast Georgia Embayment, where drilling was scheduled to begin in 1979, and is imminent in the Georges Bank area off New England. In addition, hydrogeologic and hydrochemical data obtained from the drilling studies have delineated freshwater-bearing submarine extensions of land aquifers that are important coastal ground-water resources. Hazards in the Georges Bank area include sand mobility associated with strong currents and storm-driven waves; high concentrations of suspended sediment in the water column that, when mixed with spilled oil, may sink to the bottom; and slumping along the upper slope. In the Baltimore Canyon, high sediment mobility accompanies major winter storms, and slumped material may cover as much as 20 percent of the upper slope. Potentially unstable slope areas are being studied in great detail to provide data on timing, triggering mechanisms, and rates of sediment movement. In the Southeast Georgia Embayment and Blake Plateau Basin, strong Gulf Stream flow poses a major problem to all offshore operations. In the Gulf o

Circular

Long-term dynamics of production, respiration, and net CO 2 exchange in two sagebrush-steppe ecosystems

We present a synthesis of long-term measurements of CO 2 exchange in 2 US Intermountain West sagebrush-steppe ecosystems. The locations near Burns, Oregon (1995–2001), and Dubois, Idaho (1996–2001), are part of the AgriFlux Network of the Agricultural Research Service, United States Department of Agriculture. Measurements of net ecosystem CO 2 exchange ( F c ) during the growing season were continuously recorded at flux towers using the Bowen ratio-energy balance technique. Data were partitioned into gross primary productivity ( P g ) and ecosystem respiration ( R e ) using the light-response function method. Wintertime fluxes were measured during 1999/2000 and 2000/2001 and used to model fluxes in other winters. Comparison of daytime respiration derived from light-response analysis with nighttime tower measurements showed close correlation, with daytime respiration being on the average higher than nighttime respiration. Maxima of P g and R e at Burns were both 20 g CO 2 ·m −2 ·d −1 in 1998. Maxima of P g and R e at Dubois were 37 and 35 g CO 2 ·m −2 ·d −1 , respectively, in 1997. Mean annual gross primary production at Burns was 1 111 (range 475–1 715) g CO 2 ·m −2 ·y −1 or about 30% lower than that at Dubois (1 602, range 963–2 162 g CO 2 ·m −2 ·y −1 ). Across the years, both ecosystems were net sinks for atmospheric CO 2 with a mean net ecosystem CO 2 exchange of 82 g CO 2 ·m −2 ·y −1 at Burns and 253 g CO 2 ·m −2 ·y −1 at Dubois, but on a yearly basis either site could be a C sink or source, mostly depending on precipitation timing and amount. Total annual precipitation is not a good predictor of carbon sequestration across sites. Our results suggest that F c should be partitioned into P g and R e components to allow prediction of seasonal and yearly dynamics of CO 2 fluxes.

Rangeland Ecology and Management

Above- and belowground biomass carbon stock and net primary productivity maps for tidal herbaceous marshes of the United States

Accurate assessments of greenhouse gas emissions and carbon sequestration in natural ecosystems are necessary to develop climate mitigation strategies. Regional and national-level assessments of carbon sequestration require high-resolution data to be available for large areas, increasing the need for remote sensing products that quantify carbon stocks and fluxes. The Intergovernmental Panel on Climate Change (IPCC) provides guidelines on how to quantify carbon flux using land cover land change and biomass carbon stock information. Net primary productivity (NPP), carbon uptake, and storage in vegetation, can also be used to model net carbon sequestration and net carbon export from an ecosystem (net ecosystem carbon balance). While biomass and NPP map products for terrestrial ecosystems are available, there are currently no conterminous United States (CONUS) biomass carbon stock or NPP maps for tidal herbaceous marshes. In this study, we used peak soil adjusted vegetation index (SAVI) values, derived from Landsat 8 composites, and five other vegetation indices, plus a categorical variable for the CONUS region (Pacific Northwest, California, Northeast, Mid-Atlantic, South Atlantic-Gulf, or Everglades), to model spatially explicit aboveground peak biomass stocks in tidal marshes (i.e., tidal palustrine and estuarine herbaceous marshes) for the first time. Tidal marsh carbon conversion factors, root-to-shoot ratios, and vegetation turnover rates, were compiled from the literature and used to convert peak aboveground biomass to peak total (above- and belowground) biomass and NPP. An extensive literature search for aboveground turnover rates produced sparse and variable values; therefore, we used an informed assumption of a turnover rate of one crop per year for all CONUS tidal marshes. Due to the lack of turnover rate data, the NPP map is identical to the peak biomass carbon stock map. In reality, it is probable that turnover rate varies by region, given seasonal length differences; however, the NPP map provides the best available information on spatially explicit CONUS tidal marsh NPP. This study identifies gaps in the scientific knowledge, to support future studies in addressing this lack of turnover data. Across CONUS, average total peak biomass carbon stock in tidal marshes was 848 g C m −2 (871 g C m −2 in palustrine and 838 g C m −2 in estuarine marshes), and based on a median biomass turnover rate of 1, it is expected that the mean NPP annual flux for tidal marshes is similar (e.g., 848 g C m −2 y −1 ). Peak biomass carbon stocks in tidal marshes were lowest in the Florida Everglades region and highest in the California regions. These are the first fine-scale national maps of biomass carbon and NPP for tidal wetlands, spanning all of CONUS. These estimates of CONUS total peak biomass carbon stocks and NPP rates for tidal marshes can support regional- and national-scale assessments of greenhouse gas emissions, as well as natural resource management of coastal wetlands, as part of nature-based climate solution efforts.

Remote Sensing

Carbon balance of the terrestrial biosphere in the twentieth century: analyses of CO 2 , climate and land use effects with four process-based ecosystem models

The concurrent effects of increasing atmospheric CO 2 concentration, climate variability, and cropland establishment and abandonment on terrestrial carbon storage between 1920 and 1992 were assessed using a standard simulation protocol with four process-based terrestrial biosphere models. Over the long-term(1920–1992), the simulations yielded a time history of terrestrial uptake that is consistent (within the uncertainty) with a long-term analysis based on ice core and atmospheric CO 2 data. Up to 1958, three of four analyses indicated a net release of carbon from terrestrial ecosystems to the atmosphere caused by cropland establishment. After 1958, all analyses indicate a net uptake of carbon by terrestrial ecosystems, primarily because of the physiological effects of rapidly rising atmospheric CO 2 . During the 1980s the simulations indicate that terrestrial ecosystems stored between 0.3 and 1.5 Pg C yr−1, which is within the uncertainty of analysis based on CO 2 and O 2 budgets. Three of the four models indicated (in accordance with O 2 evidence) that the tropics were approximately neutral while a net sink existed in ecosystems north of the tropics. Although all of the models agree that the long-term effect of climate on carbon storage has been small relative to the effects of increasing atmospheric CO 2 and land use, the models disagree as to whether climate variability and change in the twentieth century has promoted carbon storage or release. Simulated interannual variability from 1958 generally reproduced the El Niño/Southern Oscillation (ENSO)-scale variability in the atmospheric CO 2 increase, but there were substantial differences in the magnitude of interannual variability simulated by the models. The analysis of the ability of the models to simulate the changing amplitude of the seasonal cycle of atmospheric CO 2 suggested that the observed trend may be a consequence of CO 2 effects, climate variability, land use changes, or a combination of these effects. The next steps for improving the process-based simulation of historical terrestrial carbon include (1) the transfer of insight gained from stand-level process studies to improve the sensitivity of simulated carbon storage responses to changes in CO 2 and climate, (2) improvements in the data sets used to drive the models so that they incorporate the timing, extent, and types of major disturbances, (3) the enhancement of the models so that they consider major crop types and management schemes, (4) development of data sets that identify the spatial extent of major crop types and management schemes through time, and (5) the consideration of the effects of anthropogenic nitrogen deposition. The evaluation of the performance of the models in the context of a more complete consideration of the factors influencing historical terrestrial carbon dynamics is important for reducing uncertainties in representing the role of terrestrial ecosystems in future projections of the Earth system.

Global Biogeochemical Cycles

Productivity and CO 2 exchange of Great Plains ecoregions. I. Shortgrass steppe: Flux tower estimates

The shortgrass steppe (SGS) occupies the southwestern part of the Great Plains. Half of the land is cultivated, but significant areas remain under natural vegetation. Despite previous studies of the SGS carbon cycle, not all aspects have been completely addressed, including gross productivity, ecosystem respiration, and ecophysiological parameters. Our analysis of 1998 − 2007 flux tower measurements at five Bowen ratio–energy balance (BREB) and three eddy covariance (EC) sites characterized seasonal and interannual variability of gross photosynthesis and ecosystem respiration. Identification of the nonrectangular hyperbolic equation for the diurnal CO 2 exchange, with vapor pressure deficit (VPD) limitation and exponential temperature response, quantified quantum yield α , photosynthetic capacity A max , and respiration rate r d with variation ranges (19 \< α \ < 51 mmol mol − 1 , 0.48 \< A max \ < 2.1 mg CO 2 m − 2 s − 1 , 0.15 \< r d \ < 0.49 mg CO 2 m − 2 s − 1 ). Gross photosynthesis varied from 1 100 to 2 700 g CO 2 m − 2 yr − 1 , respiration from 900 to 3,000 g CO 2 m − 2 yr − 1 , and net ecosystem production from − 900 to + 700 g CO 2 m − 2 yr − 1 , indicating that SGS may switch from a sink to a source depending on weather. Comparison of the 2004 − 2006 measurements at two BREB and two parallel EC flux towers located at comparable SGS sites showed moderately higher photosynthesis, lower respiration, and higher net production at the BREB than EC sites. However, the difference was not related only to methodologies, as the normalized difference vegetation index at the BREB sites was higher than at the EC sites. Overall magnitudes and seasonal patterns at the BREB and the EC sites during the 3-yr period were similar, with trajectories within the ± 1.5 standard deviation around the mean of the four sites and mostly reflecting the effects of meteorology.

Colorado, Wyoming

Geohydrology of Big Bear Valley, California: phase 1--geologic framework, recharge, and preliminary assessment of the source and age of groundwater

The Big Bear Valley, located in the San Bernardino Mountains of southern California, has increased in population in recent years. Most of the water supply for the area is pumped from the alluvial deposits that form the Big Bear Valley groundwater basin. This study was conducted to better understand the thickness and structure of the groundwater basin in order to estimate the quantity and distribution of natural recharge to Big Bear Valley. A gravity survey was used to estimate the thickness of the alluvial deposits that form the Big Bear Valley groundwater basin. This determined that the alluvial deposits reach a maximum thickness of 1,500 to 2,000 feet beneath the center of Big Bear Lake and the area between Big Bear and Baldwin Lakes, and decrease to less than 500 feet thick beneath the eastern end of Big Bear Lake. Interferometric Synthetic Aperture Radar (InSAR) was used to measure pumping-induced land subsidence and to locate structures, such as faults, that could affect groundwater movement. The measurements indicated small amounts of land deformation (uplift and subsidence) in the area between Big Bear Lake and Baldwin Lake, the area near the city of Big Bear Lake, and the area near Sugarloaf, California. Both the gravity and InSAR measurements indicated the possible presence of subsurface faults in subbasins between Big Bear and Baldwin Lakes, but additional data are required for confirmation. The distribution and quantity of groundwater recharge in the area were evaluated by using a regional water-balance model (Basin Characterization Model, or BCM) and a daily rainfall-runoff model (INFILv3). The BCM calculated spatially distributed potential recharge in the study area of approximately 12,700 acre-feet per year (acre-ft/yr) of potential in-place recharge and 30,800 acre-ft/yr of potential runoff. Using the assumption that only 10 percent of the runoff becomes recharge, this approach indicated there is approximately 15,800 acre-ft/yr of total recharge in Big Bear Valley. The INFILv3 model was modified for this study to include a perched zone beneath the root zone to better simulate lateral seepage and recharge in the shallow subsurface in mountainous terrain. The climate input used in the INFILv3 model was developed by using daily climate data from 84 National Climatic Data Center stations and published Parameter Regression on Independent Slopes Model (PRISM) average monthly precipitation maps to match the drier average monthly precipitation measured in the Baldwin Lake drainage basin. This model resulted in a good representation of localized rain-shadow effects and calibrated well to measured lake volumes at Big Bear and Baldwin Lakes. The simulated average annual recharge was about 5,480 acre-ft/yr in the Big Bear study area, with about 2,800 acre-ft/yr in the Big Bear Lake surface-water drainage basin and about 2,680 acre-ft/yr in the Baldwin Lake surface-water drainage basin. One spring and eight wells were sampled and analyzed for chemical and isotopic data in 2005 and 2006 to determine if isotopic techniques could be used to assess the sources and ages of groundwater in the Big Bear Valley. This approach showed that the predominant source of recharge to the Big Bear Valley is winter precipitation falling on the surrounding mountains. The tritium and uncorrected carbon-14 ages of samples collected from wells for this study indicated that the groundwater basin contains water of different ages, ranging from modern to about 17,200-years old.The results of these investigations provide an understanding of the lateral and vertical extent of the groundwater basin, the spatial distribution of groundwater recharge, the processes responsible for the recharge, and the source and age of groundwater in the groundwater basin. Although the studies do not provide an understanding of the detailed water-bearing properties necessary to determine the groundwater availability of the basin, they do provide a framework for the future development of a groundwater model that would help to improve the understanding of the potential hydrologic effects of water-management alternatives in Big Bear Valley.

California

Mass loading of selected major and trace elements in Lake Fork Creek near Leadville, Colorado, September-October 2001

A mass-loading study of Lake Fork Creek of the Arkansas River between Sugarloaf Dam and the mouth was completed in September-October 2001 to help ascertain the following: (1) variation of pH and aqueous constituent concentrations (calcium, sulfate, alkalinity, aluminum, cadmium, copper, iron, manganese, lead, and zinc) and their relation to toxicity standards along the study reach; (2) location and magnitude of sources of metal loading to Lake Fork Creek; (3) amount and locations of metal attenuation; (4) the effect of streamside wetlands on metal transport from contributing mine tunnels; and (5) the effect of organic-rich inflow from the Leadville National Fish Hatchery on water quality in Lake Fork Creek. The study was done in cooperation with the Bureau of Land Management, U.S. Department of Agriculture Forest Service, and U.S. Fish and Wildlife Service. Constituent concentrations and pH showed variable patterns over the study reach. Hardness-based acute and chronic toxicity standards were exceeded for some inflows and some constituents. However, stream concentrations did not exceed standards except for zinc starting in the upper parts of the study reach and extending to just downstream from the inflow from the Leadville National Fish Hatchery. Dilution from that inflow lowered stream zinc concentrations to less than acute and chronic toxicity standards. The uppermost 800 meters of the study reach that contained inflow from the Bartlett, Dinero, and Nelson mine tunnels and the Dinero wetland was the greatest source of loading for manganese and zinc. A middle section of the study reach that extended approximately 2 kilometers upstream from the National Fish Hatchery inflow to just downstream from that inflow was the largest source of aluminum, copper, iron, and lead loading. The loading was partially from the National Fish Hatchery inflow but also from unknown sources upstream from that inflow, possibly ground water. The largest sources for calcium and sulfate load to the stream were the parts of the study reach containing inflow from the tribu-taries Halfmoon Creek (calcium) and Willow Creek (sulfate). The Arkansas River and its tributaries upstream from Lake Fork Creek were the source of most of the calcium (70 percent), sulfate (82 percent), manganese (77 percent), lead (78 percent), and zinc (95 percent) loads in the Arkansas River downstream from the Lake Fork confluence. In contrast, Lake Fork Creek was the major source of aluminum (68 percent), copper (65 percent), and iron (87 percent) loads to the Arkansas River downstream from the confluence. Attenuation was not important for calcium, sulfate, or iron. However, other metals loads were reduced up to 81 percent over the study reach (aluminum, 25 percent; copper, 20 percent; manganese, 81 percent; lead, 30 percent; zinc, 72 percent). Metal attenuation in the stream occurred primarily in three locations (1) the irrigation diversion ditch; (2) the beaver pond complex extending from upstream from the Colorado Gulch inflow to just downstream from that inflow; and (3) the stream reach that included the inflow from Willow Creek. The most likely attenuation mechanism is precipitation of metal oxides and hydroxides (primarily manganese), and sorption or coprecipitation of trace elements with the precipitating phase. A mass-balance calculation indicated that the wetland between the Dinero Tunnel and Lake Fork Creek removed iron, had little effect on zinc mass load, and was a source for, or was releasing, aluminum and manganese. In contrast, the wetland that occurred between the Siwatch Tunnel and Lake Fork Creek removed aluminum, iron, manganese, and zinc from the tunnel drainage before it entered the creek. Inflow from the National Fish Hatchery increased dissolved organic carbon concentrations in Lake Fork Creek and slightly changed the composition of the dissolved organic carbon. However, dissolved organic carbon loads increased in the stream reach downs

Colorado

Applied Geochemistry Special Issue on Environmental geochemistry of modern mining

Environmental geochemistry is an integral part of the mine-life cycle, particularly for modern mining. The critical importance of environmental geochemistry begins with pre-mining baseline characterization and the assessment of environmental risks related to mining, continues through active mining especially in water and waste management practices, and culminates in mine closure. The enhanced significance of environmental geochemistry to modern mining has arisen from an increased knowledge of the impacts that historical and active mining can have on the environment, and from new regulations meant to guard against these impacts. New regulations are commonly motivated by advances in the scientific understanding of the environmental impacts of past mining. The impacts can be physical, chemical, and biological in nature. The physical challenges typically fall within the purview of engineers, whereas the chemical and biological challenges typically require a multidisciplinary array of expertise including geologists, geochemists, hydrologists, microbiologists, and biologists. The modern mine-permitting process throughout most of the world now requires that potential risks be assessed prior to the start of mining. The strategies for this risk assessment include a thorough characterization of pre-mining baseline conditions and the identification of risks specifically related to the manner in which the ore will be mined and processed, how water and waste products will be managed, and what the final configuration of the post-mining landscape will be. In the Fall 2010, the Society of Economic Geologists held a short course in conjunction with the annual meeting of the Geological Society of America in Denver, Colorado (USA) to examine the environmental geochemistry of modern mining. The intent was to focus on issues that are pertinent to current and future mines, as opposed to abandoned mines, which have been the focus of numerous previous short courses. The geochemical challenges of current and future mines share similarities with abandoned mines, but differences also exist. Mining and ore processing techniques have changed; the environmental footprint of waste materials has changed; environmental protection has become a more integral part of the mine planning process; and most historical mining was done with limited regard for the environment. The 17 papers in this special issue evolved from the Society of Economic Geologists’ short course. The relevant geochemical processes encompass the source, transport, and fate of contaminants related to the life cycle of a mine. Contaminants include metals and other inorganic species derived from geologic sources such as ore and solid mine waste, and substances brought to the site for ore processing, such as cyanide to leach gold. Factors, such as mine-waste mineralogy, hydrologic setting, mine-drainage chemistry, and microbial activity, that affect the hydrochemical risks from mining are reviewed by Nordstrom et al. In another paper, Nordstrom discusses baseline characterization at mine sites in a regulatory framework, and emphasizes the influence of mineral deposits in producing naturally elevated concentrations of many trace elements in surface water and groundwater. Surface water quality in mineralized watersheds is influenced by a number of processes that act on daily (diel) cycles and can produce dramatic variations in trace element concentrations as described by Gammons et al. Pre-mining baseline characterization studies should strive to capture the magnitude of these diel variations. Desbarats et al., using a case study of mine drainage from a gold mine, illustrate how elements that commonly occur as negatively charged species (anions) in solution, such as arsenic as arsenate, behave in an opposite fashion than most metals, which occur as positively charged species (cations). Significant improvement in the understanding of factors that influence the toxicity of metals to aquatic organisms in surface water has highlighted the importance of aqueous chemistry, particularly dissolved organic carbon, as described by Smith et al. Stream sediment contamination is another important pathway for affecting aquatic organisms, as reviewed by Besser et al. Understanding and predicting environmental consequences from mining begins with knowing the mineralogy and mineral reactivity of the ore, the wastes, and of secondary minerals formed later. Jamieson et al. review the importance of mineralogical studies in mine planning and remediation. A number of types of site-specific studies are needed to identify environmental risks related to individual mines. Lapakko reviews the general framework of mine waste characterization studies that are integral to the mine planning process. Hageman et al. present a comparative study of several static tests commonly used to characterize mine waste. The mining and ore processing practices employed at a specific mine site will vary on the basis of the commodities being targeted, the geology of the deposit, the geometry of the deposit, and the mining and ore processing methods used. Thus, these factors, in addition to the waste management practices used, can result in a variety of end-member mine waste features, each of which has its own set of challenges. Open pit mines and underground mines require waste rock to be removed to access ore. Waste rock presents unique problems because the rock is commonly mineralized at sub-economic grades and has not been processed to remove potentially problematic minerals, such as pyrite. Amos et al. examine the salient aspects of the geochemistry of waste rock. Mill tailings – the waste material after ore minerals have been removed – are a volumetrically important solid waste at many mine sites. Their fine grain size and the options for their management make their behavior in the environment distinct from that of waste rock. Lindsay et al. describe some of these differences through three case-study examples. Subaqueous disposal of tailings is another option described by Moncur et al. Cyanide leaching for gold extraction is a common method throughout the world. Johnson describes environmental aspects of cyanidation. Uranium mining presents unique environmental challenges, particularly since in-situ recovery has seen widespread use. Campbell et al. review the environmental geochemistry of uranium mining and current research on bioremediation. Ore concentrates from many types of metal mining undergo a pyrometallurgical technique known as smelting to extract the metal. Slag is the result of smelting, and it may be an environmental liability or a valuable byproduct, as described by Piatak et al. Finally, the open pits that result from surface mining commonly reach below the water table. At the end of mining, these pits may fill to form lakes that become part of the legacy of the mine. Castendyk et al., in two papers, review theoretical aspects of the environmental limnology of pit lakes. They also describe approaches that have been used to model pit lake water balance, wall-rock contributions to pit lake chemistry, pit lake water quality, and limnological processes, such as vertical mixing, through the use of three case studies.

Applied Geochemistry

Hydrology and geochemistry of a slag-affected aquifer and chemical characteristics of slag-affected ground water, northwestern Indiana and northeastern Illinois

Slag is a by-product of steel manufacturing and a ubiquitous fill material in northwestern Indiana. Ground water associated with slag deposits generally is characterized by high pH and elevated concentrations of many inorganic water-quality constituents. The U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, conducted a study in northwestern Indiana from June 1995 to September 1996 to improve understanding of the effects of slag deposits on the water quality of a glacial-outwash aquifer. The Bairstow Landfill, a slag-fill deposit overlying the Calumet aquifer near Hammond, Indiana, was studied to represent conditions in slag-deposit settings that are common in northwestern Indiana. Ground water from 10 observation wells, located in four nests at the site, and surface water from the adjacent Lake George were analyzed for values of field-measured parameters and concentrations of major ions, nutrients, trace elements, and bulk properties. Solid-phase samples of slag and aquifer sediment collected during drilling were examined with X-ray diffraction and geochemical digestion and analysis. Concentrations of calcium, potassium, sodium, and sulfate were highest in wells screened partly or fully in slag. Potassium concentrations in ground water ranged from 2.9 to 120 milligrams per liter (mg/L), were highest in water from slag deposits, and decreased with depth. The highest concentrations for aluminum, barium, molybdenum, nickel, and selenium were in water from the slag. Silica concentrations were highest in wells screened directly beneath the slag-aquifer interface, and magnesium concentrations were highest in intermediate and deep aquifer wells. Silica concentrations in shallow and intermediate aquifer wells ranged from 27 to 41 mg/L and were about 10 times greater than those in water from slag deposits. The highest concentrations for chromium, lead, and zinc were in ground water from immediately below the slag-aquifer interface. The solid-phase analyses indicated that calcite, dolomite, and quartz generally were present throughout the slag-aquifer system; barian celestite, cristobalite, manganese-bearing calcite, and minrecordite were present in fewer samples. Trace elements that are liberated from the slag may be incorporated as impurities during precipitation of major minerals, sorbed onto clays and other grainsize fractions not analyzed as part of this study, or present in low-abundance minerals that were not detected by the X-ray analysis. Mass-balance and speciation programs were used to identify geochemical processes that may be occurring as water infiltrates through the slag, flows into the aquifer, and discharges into Lake George. The geochemical models indicate that precipitation of calcite may be occurring where slag-affected water enters the aquifer. Models also indicate that dolomite precipitation and clay-mineral dissolution may be occurring at the slag-aquifer interface; however, dolomite precipitation is generally believed to require geologically long time periods. Silica may be dissolving where slag-affected ground water enters the aquifer and may be precipitating where slag-affected ground water discharges to the lakebed of Lake George. In addition to the site-specific study, a statistical analysis of regional water quality was done to compare ground water in wells affected and unaffected by slag. When com-pared to wells in background locations in the Calumet aquifer, wells screened in slag across northwestern Indiana and northeastern Illinois generally had relatively higher pH and specific-conductance values and relatively higher concentrations of alkalinity, dissolved solids, suspended solids, total organic carbon, calcium, potassium, sodium, chloride, aluminum, barium, and possibly magnesium, sulfate, chromium, cobalt, copper, cyanide, manganese, mercury, nickel, and vanadium. When compared to wells in slag and wells in background locations, ground water from immediately beneath or immediately downgradient from slag had relatively high concentrations of arsenic and silica. Water-quality characteristics in ground water at the Bairstow Landfill were similar to water-quality characteristics in slag-contact and slag-affected wells throughout northwestern Indiana.

Illinois, Indiana

Spatial and temporal variation of stream chemistry associated with contrasting geology and land-use patterns in the Chesapeake Bay watershed—Summary of results from Smith Creek, Virginia; Upper Chester River, Maryland; Conewago Creek, Pennsylvania; and Difficult Run, Virginia, 2010–2013

Despite widespread and ongoing implementation of conservation practices throughout the Chesapeake Bay watershed, water quality continues to be degraded by excess sediment and nutrient inputs. While the Chesapeake Bay Program has developed and maintains a large-scale and long-term monitoring network to detect improvements in water quality throughout the watershed, fewer resources have been allocated for monitoring smaller watersheds, even though water-quality improvements that may result from the implementation of conservation practices are likely to be first detected at smaller watershed scales. In 2010, the U.S. Geological Survey partnered with the U.S. Environmental Protection Agency and the U.S. Department of Agriculture to initiate water-quality monitoring in four selected small watersheds that were targeted for increased implementation of conservation practices. Smith Creek watershed is an agricultural watershed in the Shenandoah Valley of Virginia that is dominated by cattle and poultry production, and the Upper Chester River watershed is an agricultural watershed on the Eastern Shore of Maryland that is dominated by row-cropping activities. The Conewago Creek watershed is an agricultural watershed in southeastern Pennsylvania that is characterized by mixed agricultural activities. The fourth watershed, Difficult Run, is a suburban watershed in northern Virginia that is dominated by medium density residential development. The objective of this study was to investigate spatial and temporal variations in water chemistry and suspended sediment in these four relatively small watersheds that represent a range of land-use patterns and underlying geology to (1) characterize current water-quality conditions in these watersheds, and (2) identify the dominant sources, sinks, and transport processes in each watershed. The general study design involved two components. The first included intensive routine water-quality monitoring at an existing streamgage within each study area (including continuous water-quality monitoring as well as discrete water-quality sampling) to develop a detailed understanding of the temporal and hydrologic variability in stream chemistry and sediment transport in each watershed. The second component involved extensive water-quality monitoring at various sites throughout each watershed to develop a detailed understanding of spatial patterns. Both components were used to improve understanding of sources and transport processes affecting stream chemistry, including nutrients and suspended sediments, and their implications for detecting long-term trends related to best management practices. This report summarizes the results of monitoring that was performed from April 2010 through September 2013. Individual Small Watershed Summaries Summaries for each of the four small watersheds are presented below. Each watershed has a more descriptive and detailed section in the report, but these summaries may be particularly useful for some watershed managers and stakeholders desiring slightly less technical detail. Smith Creek Smith Creek is a 105.39-mi 2 watershed within the Shenandoah Valley that drains to the North Fork Shenandoah River. The long-term Smith Creek base-flow index is 72.3 percent, indicating that on average, approximately 72 percent of Smith Creek flow was base flow, which suggests that Smith Creek streamflow is dominated by groundwater discharge rather than stormwater runoff. A series of cluster and principal components analyses demonstrated that the majority of the variability in Smith Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed for turbidity, suspended sediments, total nitrogen, ammonium, orthophosphate, iron, total phosphorus, and the ratio of calcium to magnesium. Statistically significant inverse correlations with flow were observed for specific conductance, magnesium, δ 15 N of nitrate, pH, bicarbonate, calcium, and δ 18 O of nitrate. Of particular note, flow and nitrate were not statistically significantly correlated, likely because of the relatively complex concentration-discharge relationship observed in continuous and discrete datasets. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, turbidity, orthophosphate, total phosphorus, suspended-sediment concentration, and silica were higher during the warm season, but pH, dissolved oxygen, and sulfate were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loads in Smith Creek using the continuous water-quality monitors. The mean Smith Creek in-stream sediment load was approximately 6,900 tons per year, with nearly 90 percent of the sediment load over the 3-year study period contributed during the eight largest storm events during that period. The Smith Creek total phosphorus load was approximately 21,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods, although a substantial phosphorus load still occurs during base-flow conditions. The Smith Creek total nitrogen load was approximately 400,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions (as was the case for sediment and total phosphorus) and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Smith Creek watershed revealed how the complex geology and hydrology interacted to result in variable water chemistry. During relatively dry and low base-flow periods, much of the discharge in Smith Creek was contributed by a single dominant spring—Lacey Spring. During wetter base-flow periods, the flows in Smith Creek were largely generated by a mixture of headwater springs and forested mountain tributaries with very different geochemical composition. The headwater springs generally issued from limestone bedrock and were characterized as having relatively high nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The undeveloped, high-gradient, forested mountain sites were generally characterized by low ionic strength waters with low nutrient concentrations. Nitrate isotope data from the limestone springs generally were consistent with manure-derived nitrogen sources (such as cattle and poultry), although the possibility of other mixed sources cannot be excluded. Nitrate isotope data from the undeveloped, high-gradient forested mountain sites were more consistent with nitrogen from undisturbed soils, atmospheric deposition, or nitrogen fixation. Regardless of the nitrogen source, oxygen isotope data indicate that the nitrate was largely a result of nitrification. Land-use data indicate that manure sources of nitrogen dominated watershed nitrogen inputs. Phosphorus sources were less well studied. The presence of a single point-source discharge near the town of New Market contributed the majority of the phosphorus to Smith Creek under base-flow conditions, but nonpoint sources of phosphorus dominated the loading to Smith Creek during stormflow periods. Implementation of conservation practices increased in the Smith Creek watershed during the study period, and even though a broad range of practice types was implemented, the most common practices included stream fencing (for cattle exclusion), the development of nutrient management plans, conservation crop rotation, and the planting of cover crops. While the implementation of these conservation practices is encouraging, results indicate small increases in nitrate concentrations at the streamgage over the last 29 years, concurrent with small decreases in nitrate fluxes. It will likely be years before the cumulative effect of these practices can be detected in the Smith Creek water quality, and the magnitude of the effect of these conservation practices detected in Smith Creek will depend largely on whether nutrient loading (of manure and commercial fertilizer) is reduced over time. Upper Chester River The Upper Chester River watershed includes the 36-square-mile (mi 2 ) watershed area around several nontidal tributaries that drain into the tidal Chester River. The streamgage is on Chesterville Branch, the largest nontidal tributary (approximately 6.12 mi 2 ) and is the site for continuous water-quality monitoring for this project. The base-flow index at Chesterville Branch is about 72 percent and indicates that, as in most of the Coastal Plain, groundwater is the greatest contributor to streamflow. As such, more than 90 percent of the nitrogen in the stream is in the form of nitrate from groundwater. Continuous and discrete data collected at Chesterville Branch show the effects of streamflow and season on water quality. Significantly positive correlations with flow were observed for ammonium, dissolved and total phosphorus, sediment, and turbidity as runoff carried these constituents from the land surface into Chesterville Branch. Other constituents that increased significantly with flow include potassium, sulfate, iron, and manganese, which are likely contributed from near-stream areas and ponds with high organic-matter content. Total nitrogen, pH, and specific conductance, along with chemical constituents associated with groundwater inputs including nitrate, calcium, ratio of calcium to magnesium, silica, bicarbonate, and sodium, were negatively correlated with flow because concentrations of these constituents were diluted by runoff. Seasonal differences in water chemistry, which are most likely related to increased biologic effects on the uptake and release of chemicals in the stream and near-stream areas, also were observed. Water temperature, orthophosphate, δ 15 N of nitrate, bicarbonate, sodium, and the ratio of sodium to chloride were higher during the warm season, and dissolved oxygen, total nitrogen, nitrate, magnesium, sulfate, and manganese were higher during the cool season. Surrogate-regression models developed by using continuous water-quality data showed that the annual sediment load for the 2013 water year was about 2,600 tons, with more than 90 percent of this sediment contributed during two storms. The total phosphorus load in 2013 was about 13,000 pounds with more than 90 percent contributed during the same two storms as sediment. The load of total nitrogen, 140,000 pounds, accumulated steadily throughout the 2013 water year as nitrate in groundwater continuously discharged into the stream. The same two large storms that contributed 90 percent of the suspended-sediment and total phosphorus load only contributed about 20 percent of the annual total nitrogen load. Extensive water-quality monitoring of stream base flow throughout the Upper Chester River watershed identified how differences in land use and hydrogeology affected water chemistry. In parts of the watershed with well-drained soil and thick sandy aquifer sediments, concentrations of nitrate and other chemicals associated with fertilizer and lime application increased in streams as agricultural land use increased. More than 90 percent of the nitrogen in streams from these areas was in the form of nitrate, and concentrations ranged from about 5 milligrams per liter (mg/L) to 8 mg/L as nitrogen in the two largest tributaries. Stream nitrate concentrations were about 1 mg/L as nitrogen where soils were more poorly drained, the surficial aquifer sediments were thinner, and forests and wetlands were more widespread than agriculture. Nitrate isotope data were consistent with inorganic fertilizers ± atmospheric deposition and N 2 fixation as sources of nitrogen, and with nitrification as the dominant nitrate-forming process. Nitrate reduction was indicated by elevated δ 15 N and δ 18 O values in some samples from streams draining watersheds with poorly drained soils. An analysis of land-use data and SPARROW modeling input data attributed almost 90 percent of the nitrogen sources in the Upper Chester River watershed to inorganic fertilizer and fixation of atmospheric nitrogen by legumes, which is in agreement with the isotopic characteristics of nitrate in this watershed. Local sources of manure are limited in this area. Total phosphorus concentrations during base flow ranged from below detection to about 0.2 mg/L. Stream phosphorus concentrations during base flow were generally lower than those measured during storms because most phosphorus transport likely occurs as phosphorus attached to sediment particles during runoff. Because manure is not widely used in this area, the major source of phosphorus is likely fertilizer. The implementation of conservation practices in the Upper Chester River watershed increased substantially during the study period, with a total implementation of 1,194 U.S. Department of Agriculture-compliant practices. The most frequently used practices were oriented towards nutrient and sediment control, including cover crops, nutrient management planning, conservation crop rotation, conservation tillage, and irrigation management. The current Chesapeake Bay model for this area predicts that implementation of best management practices should result in a 13-percent decrease in overall delivery of nitrogen to the Upper Chester River. Because most nitrogen travels through the groundwater system for years to decades before being discharged to streams, the time period of monitoring was not sufficient to see the effects of these practices on water quality. The magnitude of the effect that may eventually be detected will depend on the degree to which nitrate leaching into the groundwater system is reduced over time. Loadings of phosphorus and sediment are primarily transported during large runoff events and are difficult to control and analyze for trends because of their timing and episodic nature. Conewago Creek Conewago Creek has two primary monitoring locations—one near the middle of the 47-mi 2 watershed and the other near the outlet just upstream of the Susquehanna River. The base-flow index was 47.3 percent for 2012–2013, indicating that on average, approximately 53 percent of the streamflow in Conewago Creek exited the watershed as surface flow, which suggests that the stormwater runoff was somewhat greater than groundwater discharge (base flow). A series of cluster and principal components analyses demonstrated that the majority of the variability in the Conewago Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed at both monitoring sites for ammonium, total phosphorus, orthophosphate, iron, and manganese; additionally, at the upstream monitoring station, total nitrogen demonstrated a statistically significant positive correlation with flow. Statistically significant inverse correlations with flow were observed at both sites for water temperature, specific conductance (at the downstream site only), sulfate, chloride, calcium, and magnesium. Statistically significant seasonal patterns were observed for several water-quality constituents. Water temperature, phosphorus (upstream site only), and orthophosphate were higher during the warm season, and nitrate and total nitrogen (upstream site only) were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient load in Conewago Creek by using the continuous water-quality monitors and water-quality samples. Conewago Creek sediment load was approximately 9,900 tons in 2012 and approximately 18,900 tons in 2013, with nearly 80 percent of the sediment load in 2013 contributed by the three largest storm events. Annual total nitrogen loads could not be estimated due to poor model performance. The addition of continued monitoring or a continuously recording nitrate sensor could improve estimates of total nitrogen loads. During 2012 and 2013, phosphorus loads in Conewago Creek were approximately 50,000 pounds in each year. Combining data from one high-flow synoptic sampling with the data from routine sampling revealed how the geology and hydrology interact to result in variable water chemistry throughout the Conewago Creek watershed. The areas above the upstream gage in the headwaters are generally underlain by forested non-carbonate bedrock and are characterized by relatively low nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The more developed, agricultural areas below the upstream site were generally characterized by higher ionic strength waters with higher nutrient and metal concentrations. An analysis of land-use data and SPAtially Referenced Regressions On Watershed (SPARROW) modeling data indicates that manure sources of nitrogen dominate the input of nitrogen to the watershed. Implementation of conservation practices increased in the Conewago Creek watershed during the study period, and while a broad range of practice types were implemented, the most common practices included residue and tillage management, cover crops, nutrient management, terracing, and stream fencing (for animal exclusion or bank restoration). While the implementation of these conservation practices is encouraging, the cumulative effects of these practices probably will not be detected in Conewago Creek water quality for several years. The magnitude of the effects of these conservation practices on water quality in Conewago Creek will depend largely on the extent to which nutrient loading (septic, manure, and commercial fertilizer) and sediment-producing activities are reduced over time. Difficult Run The Difficult Run watershed is a 57.82-mi 2 watershed that drains to the Potomac River. The long-term Difficult Run base-flow index (from 1936 to 2010) was 57.9, indicating that approximately 58 percent of streamflow exited the watershed as base flow and 42 percent as stormflow; however, with continued development and urbanization of the watershed, the base-flow index has decreased to 50 percent during the last 20 years. This base-flow index was less than those of the other watersheds evaluated in this study, likely because the Difficult Run watershed largely is underlain by crystalline piedmont metamorphic rocks and has a greater proportion of impervious urban land cover. A series of cluster and principal components analyses indicated that most of the variability in Difficult Run water quality could be attributed to hydrologic variability and seasonality. Statistically significant positive correlations with flow were observed for turbidity, dissolved oxygen, suspended sediments, ammonium, orthophosphate, iron, and total phosphorus. Statistically significant inverse correlations with flow were observed for water temperature, pH, specific conductance, bicarbonate, calcium, magnesium, nitrate, δ 15 N of nitrate, and silica. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, ammonium, orthophosphate, and δ 15 N of nitrate were higher during the warm season, and dissolved oxygen, nitrate, and manganese were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loading rates. The Difficult Run sediment load was approximately 8,000 tons per year, with greater than 95 percent of the sediment load in the 2013 water year contributed by the seven largest storm events. The total phosphorus load in Difficult Run was approximately 14,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods. The total nitrogen load in Difficult Run is estimated to have been approximately 140,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions than that of phosphorus and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Difficult Run watershed revealed relatively uniform generation of flow per unit of watershed area, as well as spatial variation in water quality that is strongly related to land-use activities. Elevated nitrate concentrations were observed in a subset of monitoring sites that are inversely correlated with population density and positively correlated to the septic system density within each subwatershed. The majority of the elevated nitrate concentrations for these sites are hypothesized to be caused by nitrate leaching from septic systems, more so than homeowner fertilizer usage among these subwatersheds that have lower population densities than other parts of the watershed. Nitrate isotope data, temporal patterns in the water-quality data, mass-balance computations, and a separate land-use analysis all generally indicate that leachate from septic systems was the likely source of the elevated nitrate. Another group of water-quality sites have relatively low nitrogen concentrations, are located in areas that are served by city sewer lines, and have experienced stream restoration activities. A final group of sites drained the areas with the highest imperviousness and had strongly elevated specific conductance, chloride, and sodium, which were likely caused by a combination of road salting and other anthropogenic sources draining these urbanized areas in the watershed. A fourth group of sites represents a mixture of water sources and had water quality similar to that at the Difficult Run streamgage. Analysis of the nitrate isotope data generally indicates a broad range of composition indicative of mixed natural and anthropogenic nitrogen sources. Implementation of conservation practices increased in the Difficult Run watershed during the study period, and while a broad range of practice types was implemented, the most common practices included stream restoration. While the implementation of these conservation practices is encouraging, the cumulative effect of these practices probably will not be detected in Difficult Run water quality for several years.

Maryland, Pennsylvania, Virginia

A spatially referenced regression model (SPARROW) for suspended sediment in streams of the Conterminous U.S.

Suspended sediment has long been recognized as an important contaminant affecting water resources. Besides its direct role in determining water clarity, bridge scour and reservoir storage, sediment serves as a vehicle for the transport of many binding contaminants, including nutrients, trace metals, semi-volatile organic compounds, a nd numerous pesticides (U.S. Environmental Protection Agency, 2000a). Recent efforts to addr ess water-quality concerns through the Total Maximum Daily Load (TMDL) process have iden tified sediment as the single most prevalent cause of impairment in the Nation’s streams a nd rivers (U.S. Environmental Protection Agency, 2000b). Moreover, sediment has been identified as a medium for the tran sport and sequestration of organic carbon, playing a potentia lly important role in understa nding sources and sinks in the global carbon budget (Stallard, 1998). A comprehensive understanding of sediment fate a nd transport is considered essential to the design and implementation of effective plans for sediment management (Osterkamp and others, 1998, U.S. General Accounting Office, 1990). An exte nsive literature addr essing the problem of quantifying sediment transport has produced a nu mber of methods for estimating its flux (see Cohn, 1995, and Robertson and Roerish, 1999, for us eful surveys). The accuracy of these methods is compromised by uncertainty in the concentration measurements and by the highly episodic nature of sediment movement, particul arly when the methods are applied to smaller basins. However, for annual or decadal flux es timates, the methods are generally reliable if calibrated with extended periods of data (Robertson and Roerish, 1999). A substantial literature also supports the Universal Soil Loss Equation (U SLE) (Soil Conservation Service, 1983), an engineering method for estimating sheet and rill erosion, although the empirical credentials of the USLE have recently been questioned (Tri mble and Crosson, 2000). Conversely, relatively little direct evidence is available concerning the fate of sediment. The common practice of quantifying sediment fate with a sediment deliv ery ratio, estimated from a simple empirical relation with upstream basin area, does not artic ulate the relative importance of individual storage sites within a basin (Wolman, 1977). Rates of sediment deposition in reservoirs and flood plains can be determined from empirical measurement s , but only a limited number of sites have been monitored, and net rates of deposition or loss from other potential sinks and sources is largely unknown (Stallard, 1998). In particular, little is known about how much sediment loss from fields ultimately makes its way to stream channels, and how much sediment is subsequently stored in or lost from th e streambed (Meade and Parker, 1985, Trimble and Crosson, 2000). This paper reports on recent progress made to a ddress empirically the question of sediment fate and transport on a national scale. The model pres ented here is based on the SPAtially Referenced Regression On Watershed attr ibutes (SPARROW) methodology, fi rst used to estimate the distribution of nutrients in str eams and rivers of the United Stat es, and subsequently shown to describe land and stream processes affecting the delivery of nutrients (Smith and others, 1997, Alexander and others, 2000, Preston and Brakeb ill, 1999). The model makes use of numerous spatial datasets, available at the national level, to explain long-term sediment water-quality conditions in major streams and rivers throughou t the United States. Sediment sources are identified using sediment erosion rates from the National Resources I nventory (NRI) (Natural Resources Conservation Service, 2000) and apportioned over the landscape according to 30- meter resolution land-use information from th e National Land Cover Data set (NLCD) (U.S. Geological Survey, 2000a). More than 76,000 reservoirs from the National Inventory of Dams (NID) (U.S. Army Corps of Engin eers, 1996) are identified as pot ential sediment sinks. Other, non-anthropogenic sources and sinks are identified using soil in formation from the State Soil Survey Geographic (STATSGO) data base (Schwarz and Alexander, 1995) and spatial coverages representing surficial rock t ype and vegetative cover. The SPA RROW model empirically relates these diverse spatial datasets to estimates of long-term, mean annual sediment flux computed from concentration and flow measurements co llected over the period 1985 -95 from more than 400 monitoring stations maintained by the Na tional Stream Quality Accounting Network (Alexander and others, 1998), the National Wa ter Quality Assessment Program, and U.S. Geological Survey District offices (Turcios and Gray, in press). Th e calibrated model is used to estimate sediment flux for over 60,000 stream segments included in the River Reach File 1 (RF1) stream network (Alexander and others, 1999). SPARROW uses statis tical methods to calibrate a simple, structural model of riverine water quality, one that imposes mass ba lance in accounting for changes in contaminant flux. As applied here, the mass-balance approach facilitates the interpretation of model results in terms of physical processes affecting sediment transport, and makes possible the estimation of various rates of sediment generation and loss associated with stream channels and features of the landscape. The statistical approach provides a basi s for assessing the error of these inferred rates and of the error in extrapolated estimates of sediment flux made for streams in the RF1 network. An important implication of the holistic modeling approach adopted in this analysis is that estimates of sediment production and loss ar e based on, and therefore consistent with, measurements of in-stream flux. Other ancillary information, such as direct measurements of long-term sediment storage and release from rese rvoirs (Steffen, 1996), is incorporated into the analysis by specifying additional equations expl aining these ancillary variables. The imposition of cross-equation constraints affords this info rmation a statistically consistent weight in explaining in-stream sediment flux. Thus, the me thodology described here represents a general framework for synthesizing a wide spectrum of available information relevant to the understanding of sediment fate and transport.

Conterminous United States

Channel change and sediment transport in the Puyallup River watershed through 2022

The Puyallup River drains a 990 square mile watershed in western Washington, with headwaters on the glacier-covered flanks of Mount Rainier. Major tributaries include the White, Carbon, and Mowich Rivers. In the levee-confined reaches of the lower watershed, loss of flood conveyance due to sand and gravel deposition has been a chronic issue. Over much of the 20th century, flood conveyance was maintained through sediment removal, but this practice ended in the late 1990s. Flood hazard management activities since the 1990s have primarily involved levee removal or setback projects. Assessments of 1984-2009 repeat cross sections suggested that sediment deposition rates were particularly high in reaches with recent levee setbacks. However, there have been no assessments of recent deposition rates since the 2009 surveys. There are also concerns that intensifying flood hydrology or increased sediment delivery from Mount Rainier may exacerbate deposition. However, assessment of those risks has been hindered by limited understanding of watershed-scale sediment delivery and routing, particularly for coarse sand and gravel. The U.S. Geological Survey, in cooperation with Pierce County, initiated this study to improve understanding of sediment deposition in the lower Puyallup River watershed. This work is primarily based on differencing of multiple aerial lidar datasets collected during 2002–2022, supplemented by early 1990 photogrammetric elevation datasets, geomorphic assessments of streamgage data, historical topographic surveys from 1907, and previously collected sediment transport measurements. Analyses cover the Puyallup, Carbon, and Mowich Rivers, but do not include the White River. During 2004–2020, repeat aerial lidar indicates that 1.3 ± 0.3 million yd3 of sediment accumulated in the lower 20 valley miles (VMs) of the Puyallup River, averaging 80,000 ± 20,000 cubic yards per year (yd3/yr). Deposition was observed during both 2004–11 and 2011–20 lidar differencing intervals. This continued a long-term depositional trend that extends back to at least 1977. From 2004 to 2011, deposition rates along the Soldiers Home levee setback reach, the only setback project downstream of VM 20 completed prior to 2011, were approximately four times higher than in adjacent unmodified reaches. From 2011 to 2020, two additional setback projects were completed; volumetric deposition rates over all three setback reaches were similar to adjacent unmodified reaches, suggesting elevated setback deposition in the 2004–11 interval may have been influenced by an extreme flood in November 2006. These levee setback projects increased the local cross-sectional area of the floodway, used as a rough proxy for relative flood conveyance, by 50 to 200 percent above 2004 conditions. If deposition continued at recent rates, cross-sectional area over the levee setback reaches would be reduced back to 2004 values by 2050-90. Deposition also occurred over the lower six VMs of the Carbon River during 2004–20, though volumes (0.15 ± 0.09 million yd3) were an order of magnitude lower than along the Puyallup River. Relatively lower deposition rates in the Carbon River are most likely the combined result of modestly lower incoming sediment loads, modestly steeper channel slope, and the additional sediment transport capacity provided by two large non-glacial tributaries that enter the Carbon River near VM 5. Upstream of the depositional reaches described above, 2002–22 sediment storage trends along the Puyallup, Carbon, and Mowich Rivers were predominately negative (net erosion) up to the Mount Rainier National Park boundary. Net erosion was the result of bank and bluff erosion exceeding deposition across wetted channel and bare gravel areas, as opposed to uniform vertical downcutting. Net erosion along these river valleys delivered 3.4 ± 0.6 million yd3 to the river system, equivalent to 190,000 ± 35,000 yd3/yr. Most of that volume was supplied by erosion of relatively low (4–10 ft) surfaces along the Puyallup and Mowich Rivers and tall (300 ft) glacial bluffs along the lower Carbon River. Substantial aggradation from 1984 to 2009 reported by Czuba and others (2010) along reaches of the Puyallup River (VM 19–22) where levee confinement has recently been removed was most likely an artifact of methodologic bias. The Puyallup, Mowich, and Carbon Rivers drain five distinct glaciated watersheds on the flanks of Mount Rainier, four of which were assessed in this study. All four watersheds were impacted by an extreme November 2006 rainstorm. Between 2002 and 2008, debris flows occurred in all four headwater areas, collectively eroding at least 2.1 million yd3 of sediment. These debris flows formed distinct deposits one to two miles downstream of source areas, depositing 30-50 percent of the material eroded upstream. From 2008 to 2022, no headwater debris flows were observed and overall rates of geomorphic change in the headwaters were low. Rivers eroded into debris flow deposits emplaced over the 2002–08 interval, but re-deposited equivalent volumes of material within a half mile downstream. Stage-discharge relations at five streamgages on upland rivers draining Mount Rainier show either net channel incision or dynamic variability with no long-term trend over the past 60–100 years. Observations of pervasive river valley erosion and stable or incising trends at long-term streamgages in the upper watershed do not support prior claims of widespread and accelerating aggradation of upland rivers draining Mount Rainier. Erosion and deposition volumes estimated in this report were combined with sediment transport estimates from limited suspended sediment and bedload measurements, estimates of sub-glacial erosion rates, and sediment delivery from non-glacial tributaries to construct watershed-scale sediment budgets for the Puyallup River watershed. During 2004–20, the estimated sediment load entering the depositional lowlands was well balanced by estimated inputs from, in order of relative magnitude, subglacial erosion (33–60 percent of total sediment load), erosion along the major river valleys (25–45 percent), erosion in recently deglaciated headwater areas (7–17 percent) and non-glacial tributaries (3–9 percent). These results are specific to the study period and represent total sediment loads, most of which is fine material carried in suspension. The relative sourcing of sand and gravel may be different than implied by this sediment budget. Downstream of VM 12, comparison of 1907 and 2009 channel surveys show net lowering of the channel thalweg of 4–12 ft. A long-term gage near VM 22 shows lowering of 4–5 ft through the 1960s. Lowering at both locations was inferred to be a channel response to the substantial straightening, and so steepening, of the river during major phases of levee construction through the early and mid-20th century. Application of a simple empirical bedload-discharge power-law relation to an ensemble of model-estimated daily mean discharge records in the lower Puyallup River between 1977 and 2100 projects that annual bedload transport capacity in the lower Puyallup River will increase by 20–60 percent by the middle of the 21st century. Actual changes in bedload transport and deposition rates will depend on concurrent changes in sediment supply and local hydraulics governing deposition. This report presents several key conclusions. First, the persistence and spatial patterns of sand and gravel deposition along the lower Puyallup River support prior claims that deposition is fundamentally caused by decreases in channel slope moving downstream. Given this underlying cause and the abundance of sand and gravel available to be transported downstream, deposition is likely to continue for the foreseeable future. Second, despite continued sediment deposition, recent levee setback projects in the lower Puyallup River will likely provide several decades of flood conveyance benefits relative to a no-action alternative. Third, while the rivers linking Mount Rainier to the Puget Sound lowlands have often been discussed as conduits that either pass or accumulate sediment from Mount Rainier, observations from 2002–22 show these river valleys acting as substantial sediment sources, delivering three times more sediment than recently deglaciated headwater areas on Mount Rainier. While the persistence and underlying cause of recent river valley erosion remain unknown, sediment storage dynamics along these river valleys are likely to be a major control on sand and gravel delivery to the lower watershed.

EarthArXiv