USGS Science⌕ Search

SEARCH · USGS Science

Results for “Soil Systems”

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 1,369 records · Page 76Linked to original sources
USGS1996

Geochemical and isotopic composition of ground water with emphasis on sources of sulfate in the upper Floridan Aquifer in parts of Marion, Sumter, and Citrus counties, Florida

In inland areas of northwest central Florida, sulfate concentrations in the Upper Floridan aquifer are extremely variable and sometimes exceed drinking water standards (250 milligrams per liter). This is unusual because the aquifer is unconfined and near the surface, allowing for active recharge. The sources of sulfate and geochemical processes controlling ground-water composition were evaluated in this area. Water was sampled from thirty-three wells in parts of Marion, Sumter, and Citrus Counties, within the Southwest Florida Water Management District; these included at least a shallow and a deep well at fifteen separate locations. Ground water was analyzed for major ions, selected trace constituents, dissolved organic carbon, and stable isotopes (sulfur-34 of sulfate and sulfide, carbon-13 of inorganic carbon, deuterium, and oxygen-18). Sulfate concentrations ranged from less than 0.2 to 1,400 milligrams per liter, with higher sulfate concentrations usually in water from deeper wells. The samples can be categorized into a low sulfate group (less than 30 milligrams per liter) and a high sulfate group (greater than 30 milligrams per liter). For the high sulfate water, concentrations of calcium and magnesium increased concurrently with sulfate. Chemical and isotopic data and mass-balance modeling indicate that the composition of high sulfate waters is controlled by dedolomitization reactions (dolomite dissolution and calcite precipitation, driven by dissolution of gypsum). Gypsum occurs deeper in the aquifer than open intervals of sampled wells. Upward flow has been documented in deeper parts of the aquifer in the study area, which may be driven by localized discharge areas or rapid flow in shallow parts of the aquifer. Mixing between shallow ground water and sulfate-rich water that dissolved gypsum at the base of the aquifer is probably responsible for the range of concentrations observed in the study area. Other solutes that increased with sulfate apparently originate from the gypsum itself, from other mineral assemblages found deeper in the aquifer in association with gypsum, and from residual seawater from less- flushed, deeper parts of the aquifer. These ions are subsequently transported with sulfate to shallower parts of the aquifer where gypsum is not present. The composition of low sulfate ground water is controlled by differences in the extent of microbially mediated reactions, which produce carbon dioxide. This, in turn, influences the extent of calcite dissolution. Ground waters which underwent limited microbial reactions contained dissolved oxygen and were usually in ridge areas where recharge typically is rapid. Anaerobic waters were in lower lying areas of Sumter County, where soils are poorly drained and aquifer recharge is slow. Anaerobic waters had higher concentrations of calcium, bicarbonate, sulfide, dissolved organic carbon, iron, manganese, and silica, and had lower concentrations of nitrate than aerobic ground waters. For low sulfate waters, sulfate generally originates from meteoric sources (atmospheric precipitation), with variable amounts of oxidation of reduced sulfur and sulfate reduction. Sulfide is sometimes removed from solution, probably by precipitation of a sulfide minerals such as pyrite. In areas where deep ground water has low sulfate concentrations, the shallow flow system is apparently deeper than where high sulfate concentrations occur, and upwelling sulfate-rich water is negligible. The range of sulfate concentrations observed in the study areas and differences in sulfate concentrations with depth indicate a complex interaction between shallow and deep ground-water flow systems.

Water-Resources Investigations Report↗
USGS2017

Integrating landslide and liquefaction hazard and loss estimates with existing USGS real-time earthquake information products

The U.S. Geological Survey (USGS) has made significant progress toward the rapid estimation of shaking and shakingrelated losses through their Did You Feel It? (DYFI), ShakeMap, ShakeCast, and PAGER products. However, quantitative estimates of the extent and severity of secondary hazards (e.g., landsliding, liquefaction) are not currently included in scenarios and real-time post-earthquake products despite their significant contributions to hazard and losses for many events worldwide. We are currently running parallel global statistical models for landslides and liquefaction developed with our collaborators in testing mode, but much work remains in order to operationalize these systems. We are expanding our efforts in this area by not only improving the existing statistical models, but also by (1) exploring more sophisticated, physics-based models where feasible; (2) incorporating uncertainties; and (3) identifying and undertaking research and product development to provide useful landslide and liquefaction estimates and their uncertainties. Although our existing models use standard predictor variables that are accessible globally or regionally, including peak ground motions, topographic slope, and distance to water bodies, we continue to explore readily available proxies for rock and soil strength as well as other susceptibility terms. This work is based on the foundation of an expanding, openly available, case-history database we are compiling along with historical ShakeMaps for each event. The expected outcome of our efforts is a robust set of real-time secondary hazards products that meet the needs of a wide variety of earthquake information users. We describe the available datasets and models, developments currently underway, and anticipated products.

Conference Paper↗
USGS2003

Geologic map of the Rio Puerco quadrangle, Bernalillo and Valencia Counties, New Mexico

The Rio Puerco quadrangle is located southwest of Albuquerque in central New Mexico and covers part of the western part of the Isleta Reservation. The U.S. Geological Survey, the New Mexico Bureau of Geology and Mineral Resources, and the University of New Mexico have conducted geologic mapping on the Isleta Reservation and vicinity as part of the Middle Rio Grande Basin Project. The map area contains surficial deposits, calcic soils, fluvial deposits of the Rio Puerco, deposits of the Santa Fe Group, and three volcanic fields. The area is characterized by predominantly north-trending normal faults with generally down-to-the-east movement. Post-Santa Fe Group deposits are composed of surficial deposits (Pleistocene-Holocene) and fluvial deposits of the Rio Puerco (Pleistocene-Holocene). The surficial deposits are divided into eolian, alluvial, colluvial, and landslide deposits. The fluvial deposits of the Rio Puerco consist of four terrace and present channel deposits. The Santa Fe Group is divided into lower and upper parts. The lower part of the Santa Fe Group is exposed near the southwestern corner of the study area where deposits consist of reddish-brown mudstone and sandstone correlated to the Popotosa Formation (Unit 1) of Lozinsky and Tedford (1991). They interpreted deposition of the unit in a basin-floor playa setting. The Popotosa Formation is in fault contact to the east with deposits of the upper Santa Fe Group. The upper Santa Fe Group is derived from major tributary fluvial systems (ancestral Rio Puerco Puerco and possibly the Rio San Jose drainages) draining the adjacent Colorado Plateau and Sierra Nacimiento and correlated to parts of Kelley's (1977) Ceja Formation of the Santa Fe Group and equivalent to Machette's (1978) Sierra Ladrones Formation, Connell's Arroyo Ojito Formation (Connell and others, 1999, and Maldonado's lithofacies of the Isleta Reservation (Maldonado and Atencio,1998a, b). The group also locally includes a fine- grained unit (lower Pleistocene) referred to here as the sand, silt, and clay of Chavez Grant (Qsc). The Ceja Formation of the Santa Fe Group as defined here is divided into the following units in descending stratigraphic order: (1) upper sand and gravel unit (upper Pliocene), (2) middle silt, sand, and clay unit (upper Pliocene), and (3) lower sand and gravel unit (Pliocene). The three volcanic fields in the map area are: (1) basalt of Cat Hills, dated at 98-110 ka and composed of seven lava flows and four cinder cones; the flows overlie calcic soils that overlie the upper sand and gravel unit of the Ceja Formation; (2) lava flow of Cat Mesa, dated at about 3 Ma and interfingers with the upper part of the Ceja Formation; (3) diabase of Mohinas Mountain, dated at 8.3 Ma (Baldridge and others, 1987) and intrudes the Popotosa Formation. Numerous high-angle faults cut the area but are mostly buried. The faults generally trend north but deviate to the northwest and northeast. The major normal faults are the Cat Mesa and Mohinas Mountain faults.

Miscellaneous Field Studies Map↗
USGS2012

Carbon and sediment accumulation in the Everglades (USA) during the past 4000 years: rates, drivers, and sources of error

Tropical and sub-tropical wetlands are considered to be globally important sources for greenhouse gases but their capacity to store carbon is presumably limited by warm soil temperatures and high rates of decomposition. Unfortunately, these assumptions can be difficult to test across long timescales because the chronology, cumulative mass, and completeness of a sedimentary profile are often difficult to establish. We therefore made a detailed analysis of a core from the principal drainage outlet of the Everglades of South Florida, to assess these problems and determine the factors that could govern carbon accumulation in this large sub-tropical wetland. Accelerator mass spectroscopy dating provided direct evidence for both hard-water and open-system sources of dating errors, whereas cumulative mass varied depending upon the type of method used. Radiocarbon dates of gastropod shells, nevertheless, seemed to provide a reliable chronology for this core once the hard-water error was quantified and subtracted. Long-term accumulation rates were then calculated to be 12.1 g m -2 yr -1 for carbon, which is less than half the average rate reported for northern and tropical peatlands. Moreover, accumulation rates remained slow and relatively steady for both organic and inorganic strata, and the slow rate of sediment accretion ( 0.2 mm yr -1 ) tracked the correspondingly slow rise in sea level (0.35 mm yr -1 ) reported for South Florida over the past 4000 years. These results suggest that sea level and the local geologic setting may impose long-term constraints on rates of sediment and carbon accumulation in the Everglades and other wetlands.

Florida↗
USGS2017

Bacterial, fungal, and plant communities exhibit no biomass or compositional response to two years of simulated nitrogen deposition in a semiarid grassland

Nitrogen (N) deposition affects myriad aspects of terrestrial ecosystem structure and function, and microbial communities may be particularly sensitive to anthropogenic N inputs. However, our understanding of N deposition effects on microbial communities is far from complete, especially for drylands where data are comparatively rare. To address the need for an improved understanding of dryland biological responses to N deposition, we conducted a two-year fertilization experiment in a semiarid grassland on the Colorado Plateau in the southwestern United States. We evaluated effects of varied levels of N inputs on archaeal, bacterial, fungal and chlorophyte community composition within three microhabitats: biological soil crusts (biocrusts), soil below biocrusts, and the plant rhizosphere. Surprisingly, N addition did not affect the community composition or diversity of any of these microbial groups; however, microbial community composition varied significantly among sampling microhabitats. Further, while plant richness, diversity, and cover showed no response to N addition, there were strong linkages between plant properties and microbial community structure. Overall, these findings highlight the potential for some dryland communities to have limited biotic ability to retain augmented N inputs, possibly leading to large N losses to the atmosphere and to aquatic systems.

Environmental Microbiology↗
USGS2011

Effects of baseline conditions on the simulated hydrologic response to projected climate change

Changes in temperature and precipitation projected from five general circulation models, using one late-twentieth-century and three twenty-first-century emission scenarios, were downscaled to three different baseline conditions. Baseline conditions are periods of measured temperature and precipitation data selected to represent twentieth-century climate. The hydrologic effects of the climate projections are evaluated using the Precipitation-Runoff Modeling System (PRMS), which is a watershed hydrology simulation model. The Almanor Catchment in the North Fork of the Feather River basin, California, is used as a case study. Differences and similarities between PRMS simulations of hydrologic components (i.e., snowpack formation and melt, evapotranspiration, and streamflow) are examined, and results indicate that the selection of a specific time period used for baseline conditions has a substantial effect on some, but not all, hydrologic variables. This effect seems to be amplified in hydrologic variables, which accumulate over time, such as soil-moisture content. Results also indicate that uncertainty related to the selection of baseline conditions should be evaluated using a range of different baseline conditions. This is particularly important for studies in basins with highly variable climate, such as the Almanor Catchment.

California↗
USGS2008

Salinity trends in the upper Colorado River basin upstream from the Grand Valley Salinity Control Unit, Colorado, 1986-2003

In 1974, the Colorado River Basin Salinity Control Act was passed into law. This law was enacted to address concerns regarding the salinity content of the Colorado River. The law authorized various construction projects in selected areas or "units" of the Colorado River Basin intended to reduce the salinity load in the Colorado River. One such area was the Grand Valley Salinity Control Unit in western Colorado. The U. S. Geological Survey has done extensive studies and research in the Grand Valley Salinity Control Unit that provide information to aid the U.S. Bureau of Reclamation and the Natural Resources Conservation Service in determining where salinity-control work may provide the best results, and to what extent salinity-control work was effective in reducing salinity concentrations and loads in the Colorado River. Previous studies have indicated that salinity concentrations and loads have been decreasing downstream from the Grand Valley Salinity Control Unit, and that the decreases are likely the result of salinity control work in these areas. Several of these reports; however, also document decreasing salinity loads upstream from the Grand Valley Salinity Control Unit. This finding was important because only a small amount of salinity-control work was being done in areas upstream from the Grand Valley Salinity Control Unit at the time the findings were reported (late 1990’s). As a result of those previous findings, the U.S. Bureau of Reclamation entered into a cooperative agreement with the U.S. Geological Survey to investigate salinity trends in selected areas bracketing the Grand Valley Salinity Control Unit and regions upstream from the Grand Valley Salinity Control Unit. The results of the study indicate that salinity loads were decreasing upstream from the Grand Valley Salinity Control Unit from 1986 through 2003, but the rates of decrease have slowed during the last 10 years. The average rate of decrease in salinity load upstream from the Grand Valley Salinity Control Unit was 10,700 tons/year. This accounts for approximately 27 percent of the decrease observed downstream from the Grand Valley Salinity Control Unit. Salinity loads were decreasing at the fastest rate (6,950 tons/year) in Region 4, which drains an area between the Colorado River at Cameo, Colorado (station CAMEO) and Colorado River above Glenwood Springs, Colorado (station GLEN) streamflow-gaging stations. Trends in salinity concentration and streamflow were tested at station CAMEO to determine if salinity concentration, streamflow, or both are controlling salinity loads upstream from the Grand Valley Salinity Control Unit. Trend tests of individual ion concentrations were included as potential indicators of what sources (based on mineral composition) may be controlling trends in the upper Colorado. No significant trend was detected for streamflow from 1986 to 2003 at station CAMEO; however, a significant downward trend was detected for salinity concentration. The trend slope indicates that salinity concentration is decreasing at a median rate of about 3.54 milligrams per liter per year. Five major ions (calcium, magnesium, sodium, sulfate, and chloride) were tested for trends. The results indicate that processes within source areas with rock and soil types (or other unidentified sources) bearing calcium, sodium, and sulfate had the largest effect on the downward trend in salinity load upstream from station CAMEO. Downward trends in salinity load resulting from ground-water sources and/or land-use change were thought to be possible reasons for the observed decreases in salinity loads; however, the cause or causes of the decreasing salinity loads are not fully understood. A reduction in the amount of ground-water percolation from Region 4 (resulting from work done through Federal irrigation system improvement programs as well as privately funded irrigation system improvements) has helped reduce annual salinity load from Region 4 by approximately 7,400 tons. This amount is equal to about 5.9 percent of the total decrease (125,000 tons, or about 6,950 tons/year) estimated to have occurred in Region 4 during water years 1986 through 2003. A geographic information system was used to quantify the change in the amount of irrigated land upstream from the Grand Valley Salinity Control Unit from 1993 through 2000. These data indicated that the amount of irrigated land did not change substantially, thus indicating that the downward trends in salinity load did not result from land-use change.

Colorado↗
USGS2006

An introduced predator alters Aleutian Island plant communities by thwarting nutrient subsidies

The ramifying effects of top predators on food webs traditionally have been studied within the framework of trophic cascades. Trophic cascades are compelling because they embody powerful indirect effects of predators on primary production. Although less studied, indirect effects of predators may occur via routes that are not exclusively trophic. We quantified how the introduction of foxes onto the Aleutian Islands transformed plant communities by reducing abundant seabird populations, thereby disrupting nutrient subsidies vectored by seabirds from sea to land. We compared soil and plant fertility, plant biomass and community composition, and stable isotopes of nitrogen in soil, plants, and other organisms on nine fox-infested and nine historically fox-free islands across the Aleutians. Additionally, we experimentally augmented nutrients on a fox-infested island to test whether differences in plant productivity and composition between fox-infested and fox-free islands could have arisen from differences in nutrient inputs between island types. Islands with historical fox infestations had soils low in phosphorus and nitrogen and plants low in tissue nitrogen. Soils, plants, slugs, flies, spiders, and bird droppings on these islands had low d15N values indicating that these organisms obtained nitrogen from internally derived sources. In contrast, soils, plants, and higher trophic level organisms on fox-free islands had elevated d15N signatures indicating that they utilized nutrients derived from the marine environment. Furthermore, soil phosphorus (but not nitrogen) and plant tissue nitrogen were higher on fox-free than fox-infested islands. Nutrient subsidized fox-free islands supported lush, high biomass plant communities dominated by graminoids. Fox-infested islands were less graminoid dominated and had higher cover and biomass of low-lying forbs and dwarf shrubs. While d15N profiles of soils and plants and graminoid biomass varied with island size and distance from shore, after accounting for these effects differences between fox-infested and fox-free islands still existed. Fertilization over four years caused a 24-fold increase in graminoid biomass and a shift toward a more graminoid dominated plant community typical of fox-free islands. These results indicate that apex predators can influence plant productivity and composition through complex interaction web pathways involving both top-down forcing and bottom-up nutrient exchanges across systems. ?? 2006 by the Ecological Society of America.

Ecological Monographs↗
USGS2017

Methods for estimating selected low-flow frequency statistics and harmonic mean flows for streams in Iowa

A statewide study was conducted to develop regression equations for estimating six selected low-flow frequency statistics and harmonic mean flows for ungaged stream sites in Iowa. The estimation equations developed for the six low-flow frequency statistics include: the annual 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years, the annual 30-day mean low flow for a recurrence interval of 5 years, and the seasonal (October 1 through December 31) 1- and 7-day mean low flows for a recurrence interval of 10 years. Estimation equations also were developed for the harmonic-mean-flow statistic. Estimates of these seven selected statistics are provided for 208 U.S. Geological Survey continuous-record streamgages using data through September 30, 2006. The study area comprises streamgages located within Iowa and 50 miles beyond the State's borders. Because trend analyses indicated statistically significant positive trends when considering the entire period of record for the majority of the streamgages, the longest, most recent period of record without a significant trend was determined for each streamgage for use in the study. The median number of years of record used to compute each of these seven selected statistics was 35. Geographic information system software was used to measure 54 selected basin characteristics for each streamgage. Following the removal of two streamgages from the initial data set, data collected for 206 streamgages were compiled to investigate three approaches for regionalization of the seven selected statistics. Regionalization, a process using statistical regression analysis, provides a relation for efficiently transferring information from a group of streamgages in a region to ungaged sites in the region. The three regionalization approaches tested included statewide, regional, and region-of-influence regressions. For the regional regression, the study area was divided into three low-flow regions on the basis of hydrologic characteristics, landform regions, and soil regions. A comparison of root mean square errors and average standard errors of prediction for the statewide, regional, and region-of-influence regressions determined that the regional regression provided the best estimates of the seven selected statistics at ungaged sites in Iowa. Because a significant number of streams in Iowa reach zero flow as their minimum flow during low-flow years, four different types of regression analyses were used: left-censored, logistic, generalized-least-squares, and weighted-least-squares regression. A total of 192 streamgages were included in the development of 27 regression equations for the three low-flow regions. For the northeast and northwest regions, a censoring threshold was used to develop 12 left-censored regression equations to estimate the 6 low-flow frequency statistics for each region. For the southern region a total of 12 regression equations were developed; 6 logistic regression equations were developed to estimate the probability of zero flow for the 6 low-flow frequency statistics and 6 generalized least-squares regression equations were developed to estimate the 6 low-flow frequency statistics, if nonzero flow is estimated first by use of the logistic equations. A weighted-least-squares regression equation was developed for each region to estimate the harmonic-mean-flow statistic. Average standard errors of estimate for the left-censored equations for the northeast region range from 64.7 to 88.1 percent and for the northwest region range from 85.8 to 111.8 percent. Misclassification percentages for the logistic equations for the southern region range from 5.6 to 14.0 percent. Average standard errors of prediction for generalized least-squares equations for the southern region range from 71.7 to 98.9 percent and pseudo coefficients of determination for the generalized-least-squares equations range from 87.7 to 91.8 percent. Average standard errors of prediction for weighted-least-squares equations developed for estimating the harmonic-mean-flow statistic for each of the three regions range from 66.4 to 80.4 percent. The regression equations are applicable only to stream sites in Iowa with low flows not significantly affected by regulation, diversion, or urbanization and with basin characteristics within the range of those used to develop the equations. If the equations are used at ungaged sites on regulated streams, or on streams affected by water-supply and agricultural withdrawals, then the estimates will need to be adjusted by the amount of regulation or withdrawal to estimate the actual flow conditions if that is of interest. Caution is advised when applying the equations for basins with characteristics near the applicable limits of the equations and for basins located in karst topography. A test of two drainage-area ratio methods using 31 pairs of streamgages, for the annual 7-day mean low-flow statistic for a recurrence interval of 10 years, indicates a weighted drainage-area ratio method provides better estimates than regional regression equations for an ungaged site on a gaged stream in Iowa when the drainage-area ratio is between 0.5 and 1.4. These regression equations will be implemented within the U.S. Geological Survey StreamStats web-based geographic-information-system tool. StreamStats allows users to click on any ungaged site on a river and compute estimates of the seven selected statistics; in addition, 90-percent prediction intervals and the measured basin characteristics for the ungaged sites also are provided. StreamStats also allows users to click on any streamgage in Iowa and estimates computed for these seven selected statistics are provided for the streamgage.

Iowa↗
USGS2003

Environmental water-quality zones for streams: A regional classification scheme

Various approaches have been used to classify large geographical areas into smaller regions of similar water quality or extrapolate water-quality data from a few streams to other unmonitored streams. A combination of some of the strengths of existing techniques is used to develop a new approach for these purposes. In this new approach, referred to here as SPARTA (SPAtial Regression-Tree Analysis), environmental characteristics for each monitored stream are first quantified using a Geographic Information System (GIS) and then regression-tree analysis is used to determine which characteristics are most statistically important in describing the distribution of a specific water-quality constituent. GIS coverages of only the most statistically significant environmental characteristics are then used to subdivide the area of interest into relatively homogeneous environmental water-quality zones. Results from the regression-tree analysis not only define the most important environmental characteristics, but also describe how to subdivide the coverage of the specific characteristic (for example, areas with <26% or ???26% soil clay content). The resulting regionalization scheme is customized for each water-quality constituent based on the environmental characteristics most statistically related to that constituent. SPARTA was used to delineate areas of similar phosphorus, nitrogen, and sediment concentrations (by including land-use characteristics) and areas of similar potential water quality (by excluding land-use characteristics). The SPARTA approach reduced the variability in water-quality concentrations (phosphorus, total nitrogen, Kjeldahl nitrogen, and suspended sediment) within similarly classified zones from that obtained using the US Environmental Protection Agency's nutrient ecoregions.

Environmental Management↗
USGS2008

Tile drainage as karst: Conduit flow and diffuse flow in a tile-drained watershed

The similarity of tiled-drained watersheds to karst drainage basins can be used to improve understanding of watershed-scale nutrient losses from subsurface tile drainage networks. In this study, short-term variations in discharge and chemistry were examined from a tile outlet collecting subsurface tile flow from a 963 ha agricultural watershed. Study objectives were to apply analytical techniques from karst springs to tile discharge to evaluate water sources and estimate the loads of agricultural pollutants discharged from the tile with conduit, intermediate and diffuse flow regimes. A two-member mixing model using nitrate, chloride and specific conductance was used to distinguish rainwater versus groundwater inputs. Results indicated that groundwater comprised 75% of the discharge for a three-day storm period and rainwater was primarily concentrated during the hydrograph peak. A contrasting pattern of solute concentrations and export loads was observed in tile flow. During base flow periods, tile flow consisted of diffuse flow from groundwater sources and contained elevated levels of nitrate, chloride and specific conductance. During storm events, suspended solids and pollutants adhered to soil surfaces (phosphorus, ammonium and organic nitrogen) were concentrated and discharged during the rapid, conduit flow portion of the hydrograph. During a three-day period, conduit flow occurred for 5.6% of the time but accounted for 16.5% of the total flow. Nitrate and chloride were delivered primarily with diffuse flow (more than 70%), whereas 80-94% of total suspended sediment, phosphorus and ammonium were exported with conduit and intermediate flow regimes. Understanding the water sources contributing to tile drainage and the manner by which pollutant discharge occurs from these systems (conduit, intermediate or diffuse flow) may be useful for designing, implementing and evaluating non-point source reduction strategies in tile-drained landscapes. ?? 2007 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2018

Metamodeling and mapping of nitrate flux in the unsaturated zone and groundwater, Wisconsin, USA

Nitrate contamination of groundwater in agricultural areas poses a major challenge to the sustainability of water resources. Aquifer vulnerability models are useful tools that can help resource managers identify areas of concern, but quantifying nitrogen (N) inputs in such models is challenging, especially at large spatial scales. We sought to improve regional nitrate (NO 3 − ) input functions by characterizing unsaturated zone NO 3 − transport to groundwater through use of surrogate, machine-learning metamodels of a process-based N flux model. The metamodels used boosted regression trees (BRTs) to relate mappable landscape variables to parameters and outputs of a previous “vertical flux method” (VFM) applied at sampled wells in the Fox, Wolf, and Peshtigo (FWP) river basins in northeastern Wisconsin. In this context, the metamodels upscaled the VFM results throughout the region, and the VFM parameters and outputs are the metamodel response variables. The study area encompassed the domain of a detailed numerical model that provided additional predictor variables, including groundwater recharge, to the metamodels. We used a statistical learning framework to test a range of model complexities to identify suitable hyperparameters of the six BRT metamodels corresponding to each response variable of interest: NO 3 − source concentration factor (which determines the local NO 3 − input concentration); unsaturated zone travel time; NO 3 − concentration at the water table in 1980, 2000, and 2020 (three separate metamodels); and NO 3 − “extinction depth”, the eventual steady state depth of the NO 3 − front. The final metamodels were trained to 129 wells within the active numerical flow model area, and considered 58 mappable predictor variables compiled in a geographic information system (GIS). These metamodels had training and cross-validation testing R 2 values of 0.52 – 0.86 and 0.22 – 0.38, respectively, and predictions were compiled as maps of the above response variables. Testing performance was reasonable, considering that we limited the metamodel predictor variables to mappable factors as opposed to using all available VFM input variables. Relationships between metamodel predictor variables and mapped outputs were generally consistent with expectations, e.g. with greater source concentrations and NO 3 − at the groundwater table in areas of intensive crop use and well drained soils. Shorter unsaturated zone travel times in poorly drained areas likely indicated preferential flow through clay soils, and a tendency for fine grained deposits to collocate with areas of shallower water table. Numerical estimates of groundwater recharge were important in the metamodels and may have been a proxy for N input and redox conditions in the northern FWP, which had shallow predicted NO 3 − extinction depth. The metamodel results provide proof-of-concept for regional characterization of unsaturated zone NO 3 − transport processes in a statistical framework based on readily mappable GIS input variables.

Wisconsin↗
USGS2009

Quantifying the spatial details of carbon sequestration potential and performance

Upscaling the spatial and temporal changes of carbon stocks and fluxes from sites to regions is challenging owing to the spatial and temporal variances and covariance of driving variables and the uncertainties in both the model and the input data. Although various modeling approaches have been developed to facilitate the upscaling process, few deal with error transfer from model input to output, and error propagation in time and space. The author has developed the General Ensemble Biogeochemical Modelling System (GEMS) for upscaling carbon stocks and fluxes from sites to regions with measures of uncertainty. This chapter describes the GEMS model, its application to regional- and larger-scale areas, and the new results that demonstrate the challenges of upscaling. GEMS relies on site-scale biogeochemical models to simulate carbon dynamics at the site scale. The spatial deployment of the site-scale model in GEMS is based on the spatial and temporal joint frequency distribution of major driving variables (e.g., land cover and land use change, climate, soils, disturbances, and management). At the site scale, GEMS uses stochastic ensemble simulations to incorporate input uncertainty to quantify uncertainty transfer from input to output, and to identify trends in both input data and simulation results. It permits one to simulate the range of possible permutations of input values and identify the trends and variance in both the input data and results. Using data assimilation techniques, GEMS simulations can be constrained by field and satellite observations, including estimates of net primary production (NPP) from the Moderate Resolution Imaging Spectroradiometer (MODIS), grain yield and cropping practices, and forest inventories. The modeling philosophy embedded in GEMS makes it ideal for assimilating information with various uncertainties to support estimating the spatial details of carbon sequestration potential as well as dynamic monitoring of the performance of carbon sequestration activities over large areas. As a case study, GEMS is applied to simulate the spatial and temporal details of carbon sources, sinks, and uncertainty in the Ridge and Valley ecoregion in the eastern United States

Book chapter↗
USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗
USGS2005

Water-quality assessment of part of the Upper Mississippi River Basin, Minnesota and Wisconsin — Ground-water quality along a flow system in the Twin Cities metropolitan area, Minnesota, 1997-98

As part of a national analysis of the effects of land use on ground-water quality for the National Water-Quality Assessment Program, the U.S. Geological Survey sampled wells along a flow system in surficial glacial aquifers in the northwestern part of the Twin Cities metropolitan area during 1997 and 1998. In addition, a reconnaissance steady-state ground-water model was developed to estimate flowpaths and dates of ground-water recharge using a particle-tracking routine. Sediment samples collected during drilling had high horizontal hydraulic conductivities (ranging from about 3.1 to about 190 feet per day, based on sediment-size analysis of well cuttings), small organic carbon concentrations (ranging from less than 0.2 to 160 grams per kilogram), and pH values that were mostly alkaline (ranging from 4.9 to 8.2). Water samples were analyzed for physical properties, major ions, iron, manganese, nutrients, organic carbon, radon, pesticides, volatile organic compounds (VOCs), chlorofluorocarbons, tritium, and isotopes of nitrogen, hydrogen, and oxygen. Most of the water samples had small dissolved-oxygen concentrations (less than 1 milligram per liter). Calcium, magnesium, sodium, bicarbonate, sulfate, and chloride were the primary dissolved constituents in water samples. Nitrite plus nitrate as nitrogen (nitrate) concentrations were less than the U.S. Environmental Protection Agency (USEPA) Maximum Contaminant Level of 10 mg/L. Nitrogen isotope ratios indicated that the sources of nitrate primarily were soils, animal waste, or denitrification that increased nitrogen isotope ratios in nitrate. Small concentrations of pesticides were detected in the shallow parts of the aquifer. The herbicide prometon was the most frequently detected pesticide. Herbicides applied to control grasses and weeds in corn (atrazine, simazine, and metolachlor) also were frequently detected in water samples. All pesticide and VOCs detected were below USEPA Maximum Contaminant Levels or Health Advisory Limits. Chlorofluorocarbon compounds and tritium concentrations were used to estimate dates of recharge of ground-water samples. In general, shallower ground-water samples were more recently recharged although most water sampled from the aquifer was recharged after 1955. Although land use had substantial effects on ground-water quality, the distribution of contaminants in the aquifer also is affected by complex combinations of factors and processes that include sources of natural and anthropogenic contaminants, three-dimensional advective flow, physical and hydrologic settings, age and evolution of ground water, and transformation of chemical compounds along the flow system. Compounds such as nitrate and dissolved oxygen were greatest in water samples from the upgradient end of the flow system and near the water table. Specific conductance and dissolved solids increased along the flow system and with depth due to increase in residence time in the flow system and dissolution of aquifer materials.

Minnesota↗
USGS1997

Vertical accretion and shallow subsidence in a mangrove forest of southwestern Florida, U.S.A

Simultaneous measurements of vertical accretion from artificial soil marker horizons and soil elevation change from sedimentation-erosion table (SET) plots were used to evaluate the processes related to soil building in range, basin, and overwash mangrove forests located in a low-energy lagoon which recieves minor inputs of terregenous sediments. Vertical accretion measures reflect the contribution of surficial sedimentation (sediment deposition and surface root growth). Measures of elevation change reflect not only the contributions of vertical accretion but also those of subsurface processes such as compaction, decomposition and shrink-swell. The two measures were used to calculate amounts of shallow subsidence (accretion minus elevation change) in each mangrove forest. The three forest types represent different accretionary envrionments. The basin forest was located behind a natural berm. Hydroperiod here was controlled primarily by rainfall rather than tidal exchange, although the basin flooded during extreme tidal events. Soil accretion here occurred primarily by autochthonous organic matter inputs, and elevation was controlled by accretion and shrink-swell of the substrate apparently related to cycles of flooding-drying and/or root growth-decomposition. This hydrologically-restricted forest did not experience an accretion or elevation deficit relative to sea-level rise. The tidally dominated fringe and overwash island forests accreted through mineral sediment inputs bound in place by plant roots. Filamentous turf algae played an important role in stabilizing loose muds in the fringe forest where erosion was prevalent. Elevation in these high-energy environments was controlled not only by accretion but also by erosion and/or shallow subsidence. The rate of shallow subsidence was consistently 3-4 mm y-1 in the fringe and overwash island forests but was negligible in the basin forest. Hence, the vertical development of mangrove soils was influenced by both surface and subsurface processes and the procces controlling soil elevation differed among forest types. The mangrove ecosystem at Rookery Bay has remained stable as sea level has risen during the past 70 years. Yet, lead-210 accretion data suggest a substantial accretion deficit has occurred in the past century (accretion was 10-20 cm < sea-level rise from 1930 to 1990) in the fringe and island forests at Rookery Bay. In contrast, our measures of elevation change mostly equalled the estimates of sea-level rise and shallow subsidence. These data suggest that (1) vertical accretion in this system is driven by local sea-level rise and shallow subsidence, and (2) the mangrove forests are mostly keeping pace with sea-level rise. Thus, the vulnerability of this mangrove ecosystem to sea-level rise is best described in terms of an elevation deficit (elevation change minus sea-level rise) based on annual measures rather than an accretion deficit (accretion minus sea-level rise) based on decadal measures.

Mangroves and Salt Marshes↗
USGS2023

Global methane emissions from rivers and streams

Methane (CH 4 ) is a potent greenhouse gas and its concentrations have tripled in the atmosphere since the industrial revolution. There is evidence that global warming has increased CH 4 emissions from freshwater ecosystems 1 , 2 , providing positive feedback to the global climate. Yet for rivers and streams, the controls and the magnitude of CH 4 emissions remain highly uncertain 3 , 4 . Here we report a spatially explicit global estimate of CH 4 emissions from running waters, accounting for 27.9 (16.7–39.7) Tg CH 4 per year and roughly equal in magnitude to those of other freshwater systems 5 , 6 . Riverine CH 4 emissions are not strongly temperature dependent, with low average activation energy ( E M = 0.14 eV) compared with that of lakes and wetlands ( E M = 0.96 eV) 1 . By contrast, global patterns of emissions are characterized by large fluxes in high- and low-latitude settings as well as in human-dominated environments. These patterns are explained by edaphic and climate features that are linked to anoxia in and near fluvial habitats, including a high supply of organic matter and water saturation in hydrologically connected soils. Our results highlight the importance of land–water connections in regulating CH 4 supply to running waters, which is vulnerable not only to direct human modifications but also to several climate change responses on land.

Nature↗
USGS1972

Falling-stream turbidimeter as a means of measuring sediment concentrations in streams

An optical (photocell) sensing device was used to measure the relative transparency characteristics of sediment suspensions with a view toward improving our ability to measure the temporal variations of suspended-sediment concentration moving in streams. The instrument used was a commercial "falling stream" turbidimeter that measured the relative transparency of a gravity induced smooth flat sample stream formed as the flow drops vertically from an orifice. Output results were recorded on a strip-chart recorder as a function of concentration, up to about 100,000 mg/2, for a range of particle sizes and for three kinds of sediment--a fine sand, a silty alluvial soil, and a clayey residual soil. The sensitivity range was improved by use of both the 1/8- and 1/16-inch thick orifices. The use of the two orifice or nozzle sizes provided data useful in evaluation of the effect of particle size and type of sediment at different concentrations. The unique relative transparency-concentration relationships were more linear for coarse than for fine sediments of a given kind. It was concluded that curvature in this relationship resulted from light reflections from the surface of the particles and (or) the fact that some particles moved through the system while in the "shadow" of other particles closer to the light source. The ratio of suspended-sediment concentration to the square root of the median particle size of sand-sized sediments was found to be a useful parameter for elimination of the effect of particle size in the relative transparency-concentration relationships. Thus it was possible to evaluate the effect of the different type of sediment on the relative transparency independent of particle size. The use of this parameter to eliminate the effect of particle size was unsuccessful for finer sediments. -

Open-File Report↗