USGS Science⌕ Search

SEARCH · USGS Science

Results for “Input”

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 595 records · Page 33Linked to original sources

Precipitation-runoff and streamflow-routing models for the Willamette River basin, Oregon

Precipitation-runoff and streamflow-routing models were constructed and assessed as part of a water-quality study of the Willamette River Basin. The study was a cooperative effort between the U.S. Geological Survey (USGS) and the Oregon Department of Environmental Quality (ODEQ) and was coordinated with the USGS National Water-Quality Assessment (NAWQA) study of the Willamette River. Routing models are needed to estimate streamflow so that water-quality constituent loads can be calculated from measured concentrations and so that sources, sinks, and downstream changes in those loads can be identified. Runoff models are needed to estimate ungaged-tributary inflows for routing models and to identify flow contributions from different parts of the basin. The runoff and routing models can be run either separately or together to simulate streamflow at various locations and to examine streamflow contributions from overland flow, shallow-subsurface flow, and ground-water flow. The 11,500-square-mile Willamette River Basin was partitioned into 21 major basins and 253 subbasins. For each subbasin, digital data layers of land use, soils, geology, and topography were combined in a geographic information system (GIS) to define hydrologic response units (HRU's), the basic computational unit for the Precipitation-Runoff Modeling System (PRMS). Spatial data layers were also used to calculate noncalibrated PRMS parameter values. Other PRMS parameter values were obtained from 10 nearby calibrated subbasins of representative location and character. About 760 miles of the Willamette River system were partitioned into 4 main-stem networks and 17 major tributary networks for streamflow routing. Data from time-of-travel studies, discharge measurements, and flood analyses were used to develop equations that related stream cross-sectional area to discharge and stream width to discharge. These relations were derived for all 21 stream networks at approximately 3-mile intervals and used in the Diffusion Analogy Flow model (DAFLOW) in streamflow routing. Ten representative runoff models and 11 network-routing models were calibrated for water years 1972-75 and verified for water years 1976-78. These were the periods with the most complete and widespread streamflow record for the Willamette River Basin. Observed and estimated daily precipitation and daily minimum and maximum air temperature were used as input to the runoff models. The resulting coefficient of determination (R2) for the representative runoff models ranged from 0.69 to 0.93 for the calibration period and from 0.63 to 0.92 for the verification period; absolute errors ranged from 18 to 39 percent and from 27 to 51 percent, respectively. Bias error for the runoff modeling ranged from + 13 to -32 percent. Observed daily streamflow data were used as input to the network-routing models where available, and simulated streamflows from runoff model results were used for ungaged areas. Absolute error for the network-routing models ranged from about 21 percent for the Molalla River model, for which 70 percent of the subbasin was ungaged, to about 4 percent for the Willamette main-stem model (Albany to Salem), for which only 9 percent of the subbasin was ungaged. With an input of current streamflow, precipitation, and air temperature data the combined runoff and routing models can provide current estimates of streamflow at almost 500 locations on the main stem and major tributaries of the Willamette River with a high degree of accuracy. Relative contributions of surface runoff, subsurface flow, and ground-water flow can be assessed for 1 to 10 HRU classes in each of 253 subbasins identified for precipitation-runoff modeling. Model outputs were used with a water-quality model to simulate the movement of dye in the Pudding River as an example

Oregon↗

Phosphorus dynamics in Delavan Lake Inlet, southeastern Wisconsin, 1994

Removal of rough fish from Delavan Lake Inlet in southeastern Wisconsin transformed the inlet from a shallow, turbid-water system with few macrophytes to a clear-water system with an abundance of macrophytes and increased phosphorus concentrations. To understand the changes in the phosphorus dynamics, the U.S. Geological Survey and the U.S. Army Corps of Engineers, in cooperation with the Town of Delavan, constructed a detailed phosphorus budget for the inlet for April through September 1994. The budget included inputs from the drainage basin, atmosphere, sediments (as a function of pH and the fraction of the day under anoxic conditions), and ground water. Laboratory studies estimated phosphorus flux rates under various conditions. Field studies documented ambient conditions in the inlet and fluxes of phosphorus into and out of the inlet. Phosphorus released from the sediments was estimated by coupling the flux rates, estimated in the laboratory, with ambient conditions in the inlet. The detailed phosphorus budget indicated that the increase in phosphorus concentrations was caused primarily by elevated pH resulting from increased photosynthetic activity of the macrophytes and a high release of phosphorus from the sediments. The release of phosphorus from the sediments was the largest source of phosphorus to the inlet in the spring and summer of 1994 and in other years of low to near normal runoff; however, in years of high runoff, phosphorus input from the inlet's drainage basin was the largest source of phosphorus. A less-detailed phosphorus budget constructed for the period from February 1993 to September 1994 demonstrated that, over the entire year, runoff from the drainage basin was the dominant source in the phosphorus budget. During April-September 1994, the input of phosphorus from the inlet may especially affect the summer productivity in Delavan Lake because almost 80 percent of the phosphorus loading during this time was in the form of dissolved orthophosphate.

Wisconsin↗

Geohydrology of the unsaturated zone and simulated time of arrival of landfill leachate at the water table, Municipal Solid Waste Landfill Facility, U.S. Army Air Defense Artillery Center and Fort Bliss, El Paso County, Texas

The U.S. Air Defense Artillery Center and Fort Bliss Municipal Solid Waste Landfill Facility (MSWLF) is located about 10 miles northeast of downtown El Paso, Texas. The landfill is built on the Hueco Bolson, a deposit that yields water to five public-supply wells within 1.1 miles of the landfill boundary on all sides. The bolson deposits consist of lenses and mixtures of sand, clay, silt, gravel, and caliche. The unsaturated zone at the landfill is about 300 feet thick. The Hydrologic Evaluation of Landfill Performance (HELP) and the Multimedia Exposure Assessment Model for Evaluating the Land Disposal of Wastes (MULTIMED) computer models were used to simulate the time of first arrival of landfill leachate at the water table. Site-specific data were collected for model input. At five sites on the landfill cover, hydraulic conductivity was measured by an in situ method; in addition, laboratory values were obtained for porosity, moisture content at field capacity, and moisture content at wilting point. Twenty-seven sediment samples were collected from two adjacent boreholes drilled near the southwest corner of the landfill. Of these, 23 samples were assumed to represent the unsaturated zone beneath the landfill. The core samples were analyzed in the laboratory for various characteristics required for the HELP and MULTIMED models: initial moisture content, dry bulk density, porosity, saturated hydraulic conductivity, moisture retention percentages at various suction values, total organic carbon, and pH. Parameters were calculated for the van Genuchten and Brooks-Corey equations that relate hydraulic conductivity to saturation. A reported recharge value of 0.008 inch per year was estimated on the basis of soil- water chloride concentration. The HELP model was implemented using input values that were based mostly on site-specific data or assumed in a conservative manner. Exceptions were the default values used for waste characteristics. Flow through the landfill was assumed to be at steady state. The HELP-estimated landfill leakage rate was 101.6 millimeters per year, approximately 500 times the estimated recharge rate for the area near the landfill. The MULTIMED model was implemented using input values that were based mainly on site-specific data and some conservatively assumed values. Landfill leakage was assumed to begin when the landfill was established and to continue at a steady-state rate of 101.6 millimeters per year as estimated by the HELP model. By using an assumed solute concentration in the leachate of 1 milligram per liter and assuming no delay or decay of solute, the solute serves as a tracer to indicate the first arrival of landfill leachate. The simulated first arrival of leachate at the water table was 204 to 210 years after the establishment of the landfill.

Texas↗

Characterization of water quality and simulation of temperature, nutrients, biochemical oxygen demand, and dissolved oxygen in the Wateree River, South Carolina, 1996-98

In May 1996, the U.S. Geological Survey entered into a cooperative agreement with the Kershaw County Water and Sewer Authority to characterize and simulate the water quality in the Wateree River, South Carolina. Longitudinal profiling of dissolved-oxygen concentrations during the spring and summer of 1996 revealed dissolved-oxygen minimums occurring upstream from the point-source discharges. The mean dissolved-oxygen decrease upstream from the effluent discharges was 2.0 milligrams per liter, and the decrease downstream from the effluent discharges was 0.2 milligram per liter. Several theories were investigated to obtain an improved understanding of the dissolved-oxygen dynamics in the upper Wateree River. Data suggest that the dissolved-oxygen concentration decrease is associated with elevated levels of oxygen-consuming nutrients and metals that are flowing into the Wateree River from Lake Wateree. Analysis of long-term streamflow and water-quality data collected at two U.S. Geological Survey gaging stations suggests that no strong correlation exists between streamflow and dissolved-oxygen concentrations in the Wateree River. However, a strong negative correlation does exist between dissolved-oxygen concentrations and water temperature. Analysis of data from six South Carolina Department of Health and Environmental Control monitoring stations for 1980.95 revealed decreasing trends in ammonia nitrogen at all stations where data were available and decreasing trends in 5-day biochemical oxygen demand at three river stations. The influence of various hydrologic and point-source loading conditions on dissolved-oxygen concentrations in the Wateree River were determined by using results from water-quality simulations by the Branched Lagrangian Transport Model. The effects of five tributaries and four point-source discharges were included in the model. Data collected during two synoptic water-quality samplings on June 23.25 and August 11.13, 1997, were used to calibrate and validate the Branched Lagrangian Transport Model. The data include dye-tracer concentrations collected at six locations, stream-reaeration data collected at four locations, and water-quality and water-temperature data collected at nine locations. Hydraulic data for the Branched Lagrangian Transport Model were simulated by using the U.S. Geological Survey BRANCH one-dimensional, unsteady-flow model. Data that were used to calibrate and validate the BRANCH model included time-series of water-level and streamflow data at three locations. The domain of the hydraulic model and the transport model was a 57.3- and 43.5-mile reach of the river, respectively. A sensitivity analysis of the simulated dissolved-oxygen concentrations to model coefficients and data inputs indicated that the simulated dissolved-oxygen concentrations were most sensitive to changes in the boundary concentration inputs of water temperature and dissolved oxygen followed by sensitivity to the change in streamflow. A 35-percent increase in streamflow resulted in a negative normalized sensitivity index, indicating a decrease in dissolved-oxygen concentrations. The simulated dissolved-oxygen concentrations showed no significant sensitivity to changes in model input rate kinetics. To demonstrate the utility of the Branched Lagrangian Transport Model of the Wateree River, the model was used to simulate several hydrologic and water-quality scenarios to evaluate the effects on simulated dissolved-oxygen concentrations. The first scenario compared the 24-hour mean dissolved-oxygen concentrations for August 13, 1997, as simulated during the model validation, with simulations using two different streamflow patterns. The mean streamflow for August 13, 1997, was 2,000 cubic feet per second. Simulations were run using mean streamflows of 1,000 and 1,400 cubic feet per second while keeping the water-quality boundary conditions the same as were used during the validation simulations.

South Carolina↗

The effects of tropical cyclone-generated deposition on the sustainability of the Pearl River marsh, Louisiana: The importance of the geologic framework

Shoreline retreat is a tremendously important issue along the coast of the northern Gulf of Mexico, especially in Louisiana. Although this marine transgression results from a variety of causes, the crucial factor is the difference between marsh surface elevation and rising sea levels. In most cases, the primary cause of a marsh's inability to keep up with sea level is the lack of input of inorganic material. Although tropical cyclones provide an important source of such sediment, little effort has been made to determine the point of origin of the deposited material. In this study we use sedimentary, geochemical and biogeochemical data to identify the bed of the Pearl River and/or Lake Borgne as the source of a ~5 cm thick clastic layer deposited on the surface of the Pearl River marsh on the Louisiana/Mississippi border. Radiochemical chronologies and sedimentary evidence indicate that this layer was associated with the passage of Hurricane Katrina in 2005. As this material would otherwise have been lost to the system, this deposition indicates a net gain to marsh surface elevation. Accretion rates, determined from 137Cs and 14C profiles and the use of the Katrina layer as a stratigraphic marker, indicate that short-term (~50 years) rates are as much as an order of magnitude higher than the long- term (1000s of years) rates. We suggest that the marsh's geologic setting in an incised river valley with steep vertical constraints and a large fluvial discharge, promotes rapid accretion rates, with rates accelerating as the sea moves inland, due to extended hydroperiods and the input of clastic material from both the marine and terrestrial sides. These rates are especially large when compared to accretion occurring in the more common open marshes fringing the Gulf that lack fluvial input. The difference is particularly large when related to marsh recovery/regrowth following the deposition of thick hurricane-generated clastic layers. Given the number of similar incised river valleys along the Gulf Coast, we believe that understanding the processes controlling marsh accretion in such environments is essential in evaluating marsh sustainability on a regional basis.

Louisiana↗

Recent nitrogen storage and accumulation rates in mangrove soils exceed historic rates in the urbanized San Juan Bay Estuary (Puerto Rico, United States)

Tropical mangrove forests have been described as “coastal kidneys,” promoting sediment deposition and filtering contaminants, including excess nutrients. Coastal areas throughout the world are experiencing increased human activities, resulting in altered geomorphology, hydrology, and nutrient inputs. To effectively manage and sustain coastal mangroves, it is important to understand nitrogen (N) storage and accumulation in systems where human activities are causing rapid changes in N inputs and cycling. We examined N storage and accumulation rates in recent (1970 – 2016) and historic (1930 – 1970) decades in the context of urbanization in the San Juan Bay Estuary (SJBE, Puerto Rico), using mangrove soil cores that were radiometrically dated. Local anthropogenic stressors can alter N storage rates in peri-urban mangrove systems either directly by increasing N soil fertility or indirectly by altering hydrology (e.g., dredging, filling, and canalization). Nitrogen accumulation rates were greater in recent decades than historic decades at Piñones Forest and Martin Peña East. Martin Peña East was characterized by high urbanization, and Piñones, by the least urbanization in the SJBE. The mangrove forest at Martin Peña East fringed a poorly drained canal and often received raw sewage inputs, with N accumulation rates ranging from 17.7 to 37.9 g m –2 y –1 in recent decades. The Piñones Forest was isolated and had low flushing, possibly exacerbated by river damming, with N accumulation rates ranging from 18.6 to 24.2 g m –2 y –1 in recent decades. Nearly all (96.3%) of the estuary-wide mangrove N (9.4 Mg ha –1 ) was stored in the soils with 7.1 Mg ha –1 sequestered during 1970–2017 (0–18 cm) and 2.3 Mg ha –1 during 1930–1970 (19–28 cm). Estuary-wide mangrove soil N accumulation rates were over twice as great in recent decades (0.18 ± 0.002 Mg ha –1 y –1 ) than historically (0.08 ± 0.001 Mg ha –1 y –1 ). Nitrogen accumulation rates in SJBE mangrove soils in recent times were twofold larger than the rate of human-consumed food N that is exported as wastewater (0.08 Mg ha –1 y –1 ), suggesting the potential for mangroves to sequester human-derived N. Conservation and effective management of mangrove forests and their surrounding watersheds in the Anthropocene are important for maintaining water quality in coastal communities throughout tropical regions.

Puerto Rico, San Juan Bay Estuary↗

A 10-year continuous daily simulation of chloride flux from a suburban watershed in Fairfax County, Virginia, USA

Increasing levels of chloride in surface water are associated with detrimental effects on water quality, aquatic ecosystems, infrastructure, and human health. Numerous mass-balance studies have inferred watershed transport processes by interpreting chloride inputs and outputs, but few represent internal dynamics explicitly. We constructed a coupled water/chloride mass balance model to gain insights into storage, residence time, and transport processes in a 10-km 2 urban watershed. The model, which operates over a 10-year period at a daily time scale, represents storage in a dynamic soil-moisture reservoir, quick-flow runoff from storm events, and slow-flow runoff that sustains streamflow in dry weather. The calibrated model accurately represented (a)the observed transition from a streamflow enrichment regime in cold months to a dilution regime in warmer months, (b) the observed tendency for late-summer concentrations to be higher after winters with heavy snowfall, and (c) a period-of-record downward trend in chloride concentration likely associated with a downward trend in annual snowfall. Estimated chloride inputs averaged 195 metric tons per year, while the average output was 270 metric tons per year. In contrast, estimated storage was only 107 metric tons. The estimated mean residence time in groundwater was 1.27 years. This short residence time indicates that efforts to reduce inputs will manifest as decreased concentrations in streamflow on a management-relevant time scale of several years. The coupled mass balance model yielded insights into internal watershed dynamics that would not be possible from simple input/output analysis; such models can be useful tools for gaining insight into small watershed hydrology and pollutant transport.

Virginia↗

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal↗

Integration of vegetation classification with land cover mapping: Lessons from regional mapping efforts in the Americas

Aims : Natural resource management and biodiversity conservation rely on inventories of vegetation that span multiple management or political jurisdictions. However, while remote sensing data and analytical tools have enabled production of maps at increasing spatial resolution and reliability, there are limited examples where national or continental-scaled maps are produced to represent vegetation at high thematic detail. We illustrate two examples that have bridged the gap between traditional land cover mapping and modern vegetation classification. Study area : Our two case studies include national ( USA ) and continental (North and South America) vegetation and land cover mapping. These studies span conditions from subpolar to tropical latitudes of the Americas. Methods : Both case studies used a supervised modeling approach with the International Vegetation Classification ( IVC ) to produce maps that provide for greater thematic detail. Georeferenced locations for these vegetation types are used by machine learning algorithms to train a predictive model and generate a distribution map. Results : The USA LANDFIRE (Landscape Fire and Resource Management Planning Tools Project) case study illustrates how a history of vegetation-based classification and availability of key inputs can come together to generate standard map products covering more than 9.8 million km 2 that are unsurpassed anywhere in the world in terms of spatial and thematic resolution. That being said, it also remains clear that mapping at the thematic resolution of the IVC Group and finer resolution require very large and spatially balanced inputs of georeferenced samples. Even with extensive prior data collection efforts, these remain a key limitation. The NatureServe effort for the Americas - encompassing 22% of the global land surface - demonstrates methods and outputs suitable for worldwide application at continental scales. Conclusions : Continued collection of input data used in the case studies could enable mapping at these spatial and thematic resolutions around the globe.

Vegetation Classification and Survey↗

Metal uptake by phytoplankton during a bloom in South San Francisco Bay: Implications for metal cycling in estuaries

The 1994 spring phytoplankton bloom in South San Francisco Bay caused substantial reductions in concentrations of dissolved Cd, Ni, and Zn, but not Cu. We estimate that the equivalent of ~60% of the total annual input of Cd, Ni, and Zn from local waste‐water treatment plants is cycled through the phytoplankton in South Bay. The results suggest that processes that affect phytoplankton bloom frequency or intensity in estuaries (e.g. nutrient enrichment) may also affect metal trapping. The bloom was characterized by hydrographic surveys conducted at weekly intervals for 9 weeks. Metal samples were collected from the water column on three occasions, timed to bracket the period when the bloom was predicted. Factors that might confound observations of biological influences, such as freshwater inputs, were relatively constant during the study. Before the bloom, concentrations of dissolved Cd were 0.81 ± 0.02 nmol kg −1 , Zn concentrations were 19.8 ± 1.5 nmol kg −1 , Ni were 42 ± 1.4 nmol kg −1 , and Cu were 37 ± 1.4 nmol kg −1 . These values are elevated relative to riverine and coastal end‐members, reflecting inputs from wastewater and(or) sediments. At the height of the bloom, dissolved Zn, Cd, and Ni were reduced to 19, 50, and 75% of their prebloom concentrations, respectively. Dissolved Cu concentrations increased 20%. The mass of Cd taken up by phytoplankton was similar to the mass of Cd removed from solution if particle settling was considered, and Cd concentrations estimated in phytoplankton were higher than concentrations in suspended particulate material (SPM). Particulate concentrations of Zn and Ni during the bloom appeared to be dominated by the influence of changes in resuspension of Zn and Ni‐rich sediments.

California↗

GSFLOW-GRASS v1.0.0: GIS-enabled hydrologic modeling of coupled groundwater–surface-water systems

The importance of water moving between the atmosphere and aquifers has led to efforts to develop and maintain coupled models of surface water and groundwater. However, developing inputs to these models is usually time-consuming and requires extensive knowledge of software engineering, often prohibiting their use by many researchers and water managers, thus reducing these models' potential to promote science-driven decision-making in an era of global change and increasing water resource stress. In response to this need, we have developed GSFLOW–GRASS, a bundled set of open-source tools that develops inputs for, executes, and graphically displays the results of GSFLOW, the U.S. Geological Survey's coupled groundwater and surface-water flow model. In order to create a robust tool that can be widely implemented over diverse hydro(geo)logic settings, we built a series of GRASS GIS extensions that automatically discretizes a topological surface-water flow network that is linked with an underlying gridded groundwater domain. As inputs, GSFLOW–GRASS requires at a minimum a digital elevation model, a precipitation and temperature record, and estimates of channel parameters and hydraulic conductivity. We demonstrate the broad applicability of the toolbox by successfully testing it in environments with varying degrees of drainage integration, landscape relief, and grid resolution, as well as the presence of irregular coastal boundaries. These examples also show how GSFLOW–GRASS can be implemented to examine the role of groundwater–surface-water interactions in a diverse range of water resource and land management applications.

Geoscientific Model Development↗

Improving predictive power of physically based rainfall-induced shallow landslide models: a probablistic approach

Distributed models to forecast the spatial and temporal occurrence of rainfall-induced shallow landslides are deterministic. These models extend spatially the static stability models adopted in geotechnical engineering and adopt an infinite-slope geometry to balance the resisting and the driving forces acting on the sliding mass. An infiltration model is used to determine how rainfall changes pore-water conditions, modulating the local stability/instability conditions. A problem with the existing models is the difficulty in obtaining accurate values for the several variables that describe the material properties of the slopes. The problem is particularly severe when the models are applied over large areas, for which sufficient information on the geotechnical and hydrological conditions of the slopes is not generally available. To help solve the problem, we propose a probabilistic Monte Carlo approach to the distributed modeling of shallow rainfall-induced landslides. For the purpose, we have modified the Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Analysis (TRIGRS) code. The new code (TRIGRS-P) adopts a stochastic approach to compute, on a cell-by-cell basis, transient pore-pressure changes and related changes in the factor of safety due to rainfall infiltration. Infiltration is modeled using analytical solutions of partial differential equations describing one-dimensional vertical flow in isotropic, homogeneous materials. Both saturated and unsaturated soil conditions can be considered. TRIGRS-P copes with the natural variability inherent to the mechanical and hydrological properties of the slope materials by allowing values of the TRIGRS model input parameters to be sampled randomly from a given probability distribution. The range of variation and the mean value of the parameters can be determined by the usual methods used for preparing the TRIGRS input parameters. The outputs of several model runs obtained varying the input parameters are analyzed statistically, and compared to the original (deterministic) model output. The comparison suggests an improvement of the predictive power of the model of about 10% and 16% in two small test areas, i.e. the Frontignano (Italy) and the Mukilteo (USA) areas, respectively. We discuss the computational requirements of TRIGRS-P to determine the potential use of the numerical model to forecast the spatial and temporal occurrence of rainfall-induced shallow landslides in very large areas, extending for several hundreds or thousands of square kilometers. Parallel execution of the code using a simple process distribution and the Message Passing Interface (MPI) on multi-processor machines was successful, opening the possibly of testing the use of TRIGRS-P for the operational forecasting of rainfall-induced shallow landslides over large regions.

Geoscientific Model Development and Discussions↗

Geochemical mass balances of major elements in Lake Baikal

Major element mass balances for Lake Baikal are calculated with mostly previously published data for soluble fluxes and new, unpublished data for riverine suspended particulate matter chemistry. Physical transport seems to be the most important riverine process. The elements Ca, Mg, and Na seem to be very mobile in the weathering mantle and K and Si seem to be relatively mobile. A comparison of elemental input-output budgets and mass accumulation rates (MAR) in bottom sediments shows that most major elements, except Ca, Si, and Mn, have comparable riverine particulate matter fluxes and MARs. The addition of wet atmospheric deposition fluxes results in an excess of Ca, Mg, and Na entering the lake. The additive effect of these excess inputs during a 40-year period amounts to undetectable concentration increases in the water column. If erosion of weathered bedrock is the source of most dissolved and all particulate matter transported to the lake, theoretical elemental fluxes can be calculated with Al as the conservative element. Flux ratios (observed/theoretical) range from 0.7 to 2.2, but most fall within the acceptable range of 0.7-1.5. Major rock-forming elements are carried by rivers as weathering products and there are minimal biogeochemical processes that modify these inputs as suspended particulate matter accumulates in the bottom sediments of the lake.

Limnology and Oceanography↗

Hydrologic influence on methane and carbon dioxide dynamics at two north-central Minnesota lakes

Annual emissions of (CH4 + CO2) to the atmosphere were proportional to net hydrologic inputs of C, mostly by groundwater, at two lakes in the Shingobee River watershed in north-central Minnesota. Williams Lake (WL), a closed basin lake near the top of the watershed, had a hydraulic residence time of 2-4 yr and groundwater exchange of about +2 mol dissolved inorganic carbon (DIC) and -0.1 mol dissolved organic carbon (DOC) m-2 lake area yr-1. The Shingobee River flows through Shingobee Lake (SL) that had a hydraulic residence of 0.3-0.4 yr and received net groundwater plus surface-water inputs of +5.3 to +7.3 mol DIC and fewer than +1.3 mol (DOC + particulate organic carbon) m-2 yr-1. Approximately 60-80% of net annual C input to SL was from groundwater. Lake storage of CH4 and CO2 was greatest in late winter, with maximum emissions to the atmosphere immediately following ice melt. The lakes emitted CH4 continuously during open water, having annual losses of -1.6 mol CH4 m-2 yr-1 at WL and -1.9 mol CH4 m-2 yr-1 at SL. Although the WL epilimnion was CO2 depleted throughout summer, net annual CO2 exchange with the atmosphere was near zero because springtime emission offset summertime uptake. CO2 supersaturation resulted in emission of -8.0 mol CO2 m-2 yr-1 at SL.

Limnology and Oceanography↗

Groundwater geochemistry of Isla de Mona, Puerto Rico

In this study, we explore the differences between the hydrogeochemical processes observed in a setting that is open to input from the land surface and in a setting that is closed with respect to input from the land surface. The closed setting was a water-filled passage in a cave. Samples of groundwater and of a solid that appeared to be suspended in the relatively fresh region of saline-freshwater mixing zone were collected. The solid was determined to be aragonite. Based on the analyses of the composition and saturation state of the groundwater, the mixing of fresh and saline water and precipitation of aragonite are the controlling geochemical processes in this mixing zone. We found no evidence of sulfate reduction. Thus, this mixing zone is similar to that observed in Caleta Xel Ha, Quintana Roo, also a system that is closed with respect to input from the land surface. The open setting was an unconfined aquifer underlying the coastal plain along which four hand-dug wells are located. Two wells are at the downgradient ends of inferred flowpaths and one is along a flowpath. The composition of the groundwater in the downgradient wells is sulfide-rich and brackish. In contrast, at the well located along a flow line, the groundwater is oxygenated and brackish. All groundwater is oversaturated with respect to calcite, aragonite, and dolomite. The composition is attributed to mixing of fresh and saline groundwater, CO2 outgassing, and sulfate reduction. This mixing zone is geochemically similar to that observed in blue holes and cenotes.

Journal of Cave and Karst Studies↗

Geochemistry of Florida Bay sediments: Nutrient history at five sites in eastern and central Florida Bay

Recent seagrass dieoff and massive microalgal blooms have focused attention on the health of the Florida Bay ecosystem. Changes in nutrient input and the nutrient dynamics of Florida Bay are hypothesized to be linked to these problems, but crucial baseline information is still lacking. Efforts to restore Florida Bay to its natural condition will require information on the nutrient history of the bay. The purpose of this study was to examine distributions of organic C, total N, and total P in carbonate sediments from sites of continuous and known sedimentation rate (210Pb and 137Cs dated), in eastern and central Florida Bay. These sediments provide a record of historical changes in the C, N, and P load to the eastern and central bay. Analyses were conducted on sediments from cores collected at five sites, and on buried seagrass fragments at two sites. At three of the sites, sediments from seagrass-covered and adjacent barren areas were examined to determine differences in sedimentary geochemistry. Stable isotope analyses (??13C and ??15N) of sedimentary organic C and total N and of buried seagrass fragments were also carried out at two sites to examine possible changes in nutrient sources to the estuary. Results were consistent with recent increases in N and P in eastern Florida Bay, beginning in the early to mid 1980's. The timing of the increase in nutrient load observed in the sediment data directly preceded the first observations of massive microalgal blooms and seagrass dieoff in Florida Bay in 1987. The observed nutrification was greater for P than N, and was most pronounced at the most northeasterly site sampled (Pass Key). Isotope data (??15N) suggested that an increase in algal production accompanied the increase in N load at the Pass Key site. Along record of organic C, total N, and total P distributions from Whipray Basin in central Florida Bay showed historical peaks (mid 1700's and late 1800's) in organic C and total N, but not total P; these enrichments were nearly equivalent to recent inputs to the estuary. Barren areas were observed to have generally lower concentrations of organic C, total N, and total P in near surface sediments compared to seagrass-covered areas, but had generally similar concentrations in deeper sediments. This suggested that barren areas adjacent to seagrass-covered sites were places where relict sediment was physically transported and covered seagrass beds. This dataset provides an historical view of changes in nutrient inputs to Florida Bay, and baseline information needed for nutrient modeling of the bay.

Journal of Coastal Research↗

Operational modeling system with dynamic-wave routing

A near real-time streamflow-simulation system utilizing continuous-simulation rainfall-runoff generation with dynamic-wave routing is being developed by the U.S. Geological Survey in cooperation with the Du Page County Department of Environmental Concerns for a 24-kilometer reach of Salt Creek in Du Page County, Illinois. This system is needed in order to more effectively manage the Elmhurst Quarry Flood Control Facility, an off-line stormwater diversion reservoir located along Salt Creek. Near real time simulation capabilities will enable the testing and evaluation of potential rainfall, diversion, and return-flow scenarios on water-surface elevations along Salt Creek before implementing diversions or return-flows. The climatological inputs for the continuous-simulation rainfall-runoff model, Hydrologic Simulation Program - FORTRAN (HSPF) are obtained by Internet access and from a network of radio-telemetered precipitation gages reporting to a base-station computer. The unit area runoff time series generated from HSPF are the input for the dynamic-wave routing model. Full Equations (FEQ). The Generation and Analysis of Model Simulation Scenarios (GENSCN) interface is used as a pre- and post-processor for managing input data and displaying and managing simulation results. The GENSCN interface includes a variety of graphical and analytical tools for evaluation and quick visualization of the results of operational scenario simulations and thereby makes it possible to obtain the full benefit of the fully distributed dynamic routing results.

Conference Paper↗

Nitrogen mineralization in a mature boreal forest, Isle Royale, Michigan

The 115-ha boreal Wallace Lake watershed, located on Isle Royale in the Lake Superior Basin, receives moderate anthropogenic atmospheric N inputs (3 kg ha-1 yr-1). Ecosystem response to atmospheric N inputs is, in part, determined by internal production and consumption of this limiting nutrient. The objectives of this 2-yr study on plots dominated by birch, spruce, or alder were to quantify forest floor and surface soil net and gross N mineralization rates, and examine potential effects of changes in temperature and moisture on these processes. Gross N mineralization rates were 23 times net mineralization rates, and increased with forest floor/soil temperature and moisture. Substrate quality was a likely factor in higher gross mineralization rates beneath birch and spruce. Ammonium immobilization increased with forest floor/soil temperature and moisture. Higher net N mineralization rates beneath alder resulted from lower microbial immobilization rather than greater gross N mineralization. The greatest differences between gross N mineralization and immobilization occurred in early summer. Ammonium immobilization averaged 62% of gross N mineralization. Net and gross nitrification rates differed by vegetation type, were highest in spring and fall, and increased with moisture. Gross nitrification was 19 times net nitrification rates. Nitrate immobilization increased with soil moisture, and equalled or exceeded gross nitrification. Net and gross N mineralization or nitrification rates were not correlated. Seasonal variation in forest floor and surface soil N cycling coupled with high ecosystem retention of precipitation N inputs suggest streamwater N concentrations and output reflect trends in soil processes.

Conference Paper↗