USGS ScienceSearch

SEARCH · USGS Science

Results for “Interfaces”

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 361 records · Page 20Linked to original sources

Pyroclast/snow interactions and thermally driven slurry formation. Part 1: Theory for monodisperse grain beds

Lahars are often produced as pyroclastic flows move over snow. This phenomenon involves a complicated interplay of mechanical and thermal processes that need to be separated to get at the fundamental physics. The thermal physics of pyroclast/snow interactions form the focus of this paper. A theoretical model is developed of heat- and mass transfer at the interface between a layer of uniformly sized pyroclasts and an underlying bed of snow, for the case in which there is no relative shear motion between pyroclasts and snow. A microscale view of the interface is required to properly specify boundary conditions. The physical model leads to the prediction that the upward flux of water vapor - which depends upon emplacement temperature, pyroclast grain size, pyroclast-layer thickness, and snow permeability - is sometimes sufficient to fluidize the pyroclasts. Uniform fluidization is usually unstable to bubble formation, which leads to vigorous convection of the pyroclasts themselves. Thus, predicted threshold conditions for fluidization are tantamount to predicted thresholds for particle convection. Such predictions are quantitatively in good agreement with results of experiments described in part 2 of this paper. Because particle convection commonly causes scour of the snow bed and transformation of the pyroclast layer to a slurry, there exists a 'thermal scour' process for generating lahars from pyroclastic flows moving over snow regardless of the possible role of mechanical scour.

Bulletin of Volcanology

Three-dimensional variable-density flow simulation of a coastal aquifer in southern Oahu, Hawaii, USA

Three-dimensional modeling of groundwater flow and solute transport in the Pearl Harbor aquifer, southern Oahu, Hawaii, shows that the readjustment of the freshwater–saltwater transition zone takes a long time following changes in pumping, irrigation, or recharge in the aquifer system. It takes about 50 years for the transition zone to move 90% of the distance to its new steady position. Further, the Ghyben–Herzberg estimate of the freshwater/saltwater interface depth occurred between the 10 and 50% simulated seawater concentration contours in a complex manner during 100 years of the pumping history of the aquifer. Thus, it is not a good predictor of the depth of potable water. Pre-development recharge was used to simulate the 1880 freshwater-lens configuration. Historical pumpage and recharge distributions were used and the resulting freshwater-lens size and position were simulated through 1980. Simulations show that the transition zone moved upward and landward during the period simulated. Previous groundwater flow models for Oahu have been limited to areal models that simulate a sharp interface between freshwater and saltwater or solute-transport models that simulate a vertical aquifer section. The present model is based on the US Geological Survey’s three-dimensional solute transport (3D SUTRA) computer code. Using several new tools for pre- and post-processing of model input and results have allowed easy model construction and unprecedented visualization of the freshwater lens and underlying transition zone in Hawaii’s most developed aquifer.

Hawaii

In situ allelopathic expression by the invasive amphibious plant, Ludwigia hexapetala (water primrose) across habitat types, seasons, and salinities

Broad infestations of invasive, non-native vegetation have transformed wetlands around the world. Ludwigia hexapetala is a widespread, amphibious invasive plant with a creeping growth habit in open water and an erect growth habit in terrestrial habitats. In the upper San Francisco Estuary of California, L. hexapetala is increasingly terrestrializing into marshes and this expansion may be facilitated by allelopathy. We conducted the first field-based study on L. hexapetala allelopathy to determine whether (1) three allelochemicals known to be exuded by L. hexapetala are expressed in situ, (2) the allelochemicals are detectable in leaves, soil, and water, and (3) allelopathic expression varies by season, salinity, and growth habit (open water “patch” vs. terrestrial marsh “interface” locations). Water, soil, and L. hexapetala leaves were collected in two freshwater sites and two oligohaline sites in the upper San Francisco Estuary in summer 2021, fall 2021, and spring 2022. Myricitrin and quercitrin, known allelochemicals, and salipurposid, a newly identified polyphenol, were detected in water, soil, and leaves. There were significant differences in allelochemical concentrations under fresh versus oligohaline conditions in water and soil, but not leaves. All three allelochemicals generally had higher concentrations in patch versus interface locations, suggesting that L. hexapetala allelopathy plays a greater competitive role in open water than terrestrial habitats. Leaf concentrations of each allelochemical varied seasonally; however, both myricitrin and salipurposid had heightened concentrations in spring. These results suggest that herbicide application in early spring may be most effective in controlling L. hexapetala terrestrialization from open water to marshes.

California

Influence of groundwater recharge and well characteristics on dissolved arsenic concentrations in southeastern Michigan groundwater

Arsenic concentrations exceeding 10 ??g/l, the United States maximum contaminant level and the World Health Organization guideline value, are frequently reported in groundwater from bedrock and unconsolidated aquifers of southeastern Michigan. Although arsenic-bearing minerals (including arsenian pyrite and oxide/hydroxide phases) have been identified in Marshall Sandstone bedrock of the Mississippian aquifer system and in tills of the unconsolidated aquifer system, mechanisms responsible for arsenic mobilization and subsequent transport in groundwater are equivocal. Recent evidence has begun to suggest that groundwater recharge and characteristics of well construction may affect arsenic mobilization and transport. Therefore, we investigated the relationship between dissolved arsenic concentrations, reported groundwater recharge rates, well construction characteristics, and geology in unconsolidated and bedrock aquifers. Results of multiple linear regression analyses indicate that arsenic contamination is more prevalent in bedrock wells that are cased in proximity to the bedrock-unconsolidated interface; no other factors were associated with arsenic contamination in water drawn from bedrock or unconsolidated aquifers. Conditions appropriate for arsenic mobilization may be found along the bedrock-unconsolidated interface, including changes in reduction/oxidation potential and enhanced biogeochemical activity because of differences between geologic strata. These results are valuable for understanding arsenic mobilization and guiding well construction practices in southeastern Michigan, and may also provide insights for other regions faced with groundwater arsenic contamination. ?? Springer-Verlag 2008.

Environmental Geochemistry and Health

Riders on the storm: selective tidal movements facilitate the spawning migration of threatened delta smelt in the San Francisco Estuary

Migration strategies in estuarine fishes typically include behavioral adaptations for reducing energetic costs and mortality during travel to optimize reproductive success. The influence of tidal currents and water turbidity on individual movement behavior were investigated during the spawning migration of the threatened delta smelt, Hypomesus transpacificus , in the northern San Francisco Estuary, California, USA. Water current velocities and turbidity levels were measured concurrently with delta smelt occurrence at sites in the lower Sacramento River and San Joaquin River as turbidity increased due to first-flush winter rainstorms in January and December 2010. The presence/absence of fish at the shoal-channel interface and near the shoreline was quantified hourly over complete tidal cycles. Delta smelt were caught consistently at the shoal-channel interface during flood tides and near the shoreline during ebb tides in the turbid Sacramento River, but were rare in the clearer San Joaquin River. The apparent selective tidal movements by delta smelt would facilitate either maintaining position or moving upriver on flood tides, and minimizing advection down-estuary on ebb tides. These movements also may reflect responses to lateral gradients in water turbidity created by temporal lags in tidal velocities between the near-shore and mid-channel habitats. This migration strategy can minimize the energy spent swimming against strong river and tidal currents, as well as predation risks by remaining in turbid water. Selection pressure on individuals to remain in turbid water may underlie population-level observations suggesting that turbidity is a key habitat feature and cue initiating the delta smelt spawning migration.

California

Rates of dissolution of aluminosilicates in seawater

Dissolution of eight clay minerals, four zeolites, and quartz in seawater has been monitored for 8 1 2 years. For most of the minerals, dissolution can be described as a first-order reaction in which dissolved silica approaches from undersaturation steady concentration values with time. Characteristic reaction rate constants (k 1 ) are of the order of 10 -7 sec -1 . One of the zeolites, clinoptilolite, shows a different dissolution behavior: SiO 2 concentration in solution reaches a high value within one year, followed by a decline to a lower value, suggestive of precipitation of another silicate phase (possibly sepiolite). A mathematical solution is given for a kinetic equation combining the parabolic-rate and first-order rate processes. It is shown that in a wide range of silicate dissolution reactions taking place over long periods of time, the presence of the parabolic-rate dissolution processes cannot be detected, thereby making its inclusion in the kinetic equations unnecessary. The experimental rates of dissolution are comparable to the SiO 2 - dissolution rates in oceanic sediments near the sediment/water interface. But deeper in the sediment, the calculated dissolution rates are significantly lower than the near-interface and experimental values. © 1975.

Earth and Planetary Science Letters

Biogeochemical cycling in an organic-rich coastal marine basin. 7. Sulfur mass balance, oxygen uptake and sulfide retention

Sulfur and oxygen fluxes were quantified in the seasonally varying anoxic marine sedimentary system of Cape Lookout Bight, N.C., U.S.A. Over the three year study period, 1981–1983, the mean annual sulfate reduction rate was determined to be 18.2 ± 1.6 moles · m −2 · y −1 . This value, added to the estimate of the detrital sulfur input of 1.2 ± 4.4 gave a total sulfur input of 19.4 ± 4.7 moles · m −2 · y −1 . The sulfide flux to the sediment-water interface, measured in anaerobic benthic chambers was 4.6 ± 0.5 moles · m −2 · y −1 , and represented 37% of the annual oxygen uptake rate of 25.2 ± 2.8 moles · m −2 · y −1 . The sulfide burial rate, determined to be 15.5 ± 3.1 moles · m −2 · y −1 , was within 5% of the value predicted by summing the fluxes at the sediment-water interface. The C S "> CS weight ratio of sediment below the depth of diagenetic reaction was determined to be 2.75. The sulfide retention rate in these rapidly accumulating sediments (10 cm/yr) was 77 ± 19%. Comparison of this result with previous studies shows that rapid sediment accumulation and the lack of bioturbation control this unusually high degree of sulfide retention.

Geochimica et Cosmochimica Acta

Organic geochemistry and brine composition in Great Salt, Mono, and Walker Lakes

Samples of Recent sediments, representing up to 1000 years of accumulation, were collected from three closed basin lakes (Mono Lake, CA, Walker Lake, NV, and Great Salt Lake, UT) to assess the effects of brine composition on the accumulation of total organic carbon, the concentration of dissolved organic carbon, humic acid structure and diagenesis, and trace metal complexation. The Great Salt Lake water column is a stratified Na-Mg-Cl-SO 4 brine with low alkalinity. Algal debris is entrained in the high density (1.132–1.190 g/cc) bottom brines, and in this region maximum organic matter decomposition occurs by anaerobic processes, with sulfate ion as the terminal electron acceptor. Organic matter, below 5 cm of the sediment-water interface, degrades at a very slow rate in spite of very high pore-fluid sulfate levels. The organic carbon concentration stabilizes at 1.1 wt%. Mono Lake is an alkaline (Na-CO 3 -Cl-SO 4 ) system. The water column is stratified, but the bottom brines are of lower density relative to the Great Salt Lake, and sedimentation of algal debris is rapid. Depletion of pore-fluid sulfate, near l m of core, results in a much higher accumulation of organic carbon, approximately 6 wt%. Walker Lake is also an alkaline system. The water column is not stratified, and decomposition of organic matter occurs by aerobic processes at the sediment-water interface and by anaerobic processes below. Total organic carbon and dissolved organic carbon concentrations in Walker Lake sediments vary with location and depth due to changes in input and pore-fluid sulfate concentrations. Nuclear magnetic resonance studies ( 13 C) of humic substances and dissolved organic carbon provide information on the source of the Recent sedimentary organic carbon (aquatic vs. terrestrial), its relative state of decomposition, and its chemical structure. The spectra suggest an algal origin with little terrestrial signature at all three lakes. This is indicated by the ratio of aliphatic to aromatic carbon and the absence of chemical structures indicative of the lignin of vascular plants. The dissolved organic carbon of the Mono Lake pore fluids is structurally related to humic acid and is also related to carbohydrate metabolism. The alkaline pore fluids, due to high pH, solubilize high molecular weight organic matter from the sediments. This hydrophilic material is a metal complexing agent. Despite very high algal productivities, organic carbon accumulation can be low in stratified lakes if the anoxic bottom waters are hypersaline with high concentrations of sulfate ion. Labile organic matter is recycled to the water column and the sedimentary organic matter is relatively nonsusceptible to bacterial metabolism. As a result, pore-fluid dissolved organic carbon and metal-organic complexation are low.

Geochimica et Cosmochimica Acta

Crystallization of accessory phases in magmas by local saturation adjacent to phenocrysts

Accessory minerals commonly occur attached to or included in the major crystalline phases of felsic and some intermediate igneous rocks. Apatite is particularly common as inclusions, but Fe-Ti oxides, pyrrhotite, zircon, monazite, chevkinite and xenotime are also known from silicic rocks. Accessories may nucleate near the host crystal/ liquid interface as a result of local saturation owing to formation of a differentiated chemical boundary layer in which accessory mineral solubility would be lower than in the surrounding liquid. Differentiation of this boundary layer would be greatest adjacent to ferromagnesian phenocrysts, especially Fe-Ti oxides; it is with oxides that accessories are most commonly associated in rocks. A boundary layer may develop if the crystal grows more rapidly than diffusion can transport incorporated and rejected elements to and from the phenocryst. Diffusion must dominate over convection as a mode of mass transfer near the advancing crystal/liquid interface in order for a boundary layer to exist. Accumulation of essential structural constituent elements of accessory minerals owing to their slow diffusion in evolved silicate melt also may force local saturation, but this is not a process that applies to all cases. Local saturation is an attractive mechanism for enhancing fractionation during crystallization differentiation. If accessory minerals attached to or included in phenocrysts formed because of local saturation, their host phenocrysts must have grown rapidly when accessories nucleated in comparison to lifetimes of magma reservoirs. Some inconsistencies remain in a local saturation origin for accessory phases that cannot be evaluated without additional information.

Geochimica et Cosmochimica Acta

Salt-water-freshwater transient upconing - An implicit boundary-element solution

The boundary-element method is used to solve the set of partial differential equations describing the flow of salt water and fresh water separated by a sharp interface in the vertical plane. In order to improve the accuracy and stability of the numerical solution, a new implicit scheme was developed for calculating the motion of the interface. The performance of this scheme was tested by means of numerical simulation. The numerical results are compared to experimental results for a salt-water upconing under a drain problem.

Journal of Hydrology

Theory and application of an approximate model of saltwater upconing in aquifers

Motion and mixing of salt water and fresh water are vitally important for water-resource development throughout the world. An approximate model of saltwater upconing in aquifers is developed, which results in three non-linear coupled equations for the freshwater zone, the saltwater zone, and the transition zone. The description of the transition zone uses the concept of a boundary layer. This model invokes some assumptions to give a reasonably tractable model, considerably better than the sharp interface approximation but considerably simpler than a fully three-dimensional model with variable density. We assume the validity of the Dupuit-Forchheimer approximation of horizontal flow in each layer. Vertical hydrodynamic dispersion into the base of the transition zone is assumed and concentration of the saltwater zone is assumed constant. Solute in the transition zone is assumed to be moved by advection only. Velocity and concentration are allowed to vary vertically in the transition zone by using shape functions. Several numerical techniques can be used to solve the model equations, and simple analytical solutions can be useful in validating the numerical solution procedures. We find that the model equations can be solved with adequate accuracy using the procedures presented. The approximate model is applied to the Smoky Hill River valley in central Kansas. This model can reproduce earlier sharp interface results as well as evaluate the importance of hydrodynamic dispersion for feeding salt water to the river. We use a wide range of dispersivity values and find that unstable upconing always occurs. Therefore, in this case, hydrodynamic dispersion is not the only mechanism feeding salt water to the river. Calculations imply that unstable upconing and hydrodynamic dispersion could be equally important in transporting salt water. For example, if groundwater flux to the Smoky Hill River were only about 40% of its expected value, stable upconing could exist where hydrodynamic dispersion into a transition zone is the primary mechanism for moving salt water to the river. The current model could be useful in situations involving dense saltwater layers.

Journal of Hydrology

Transverse dispersion of contaminants in fractured permeable formations

Our interest in understanding some of the mechanisms involved in the mineralization of the Great Bend Prairie aquifer of Kansas by salt water originating from Permian bedrock formations, which are fractured sandstones, has prompted this basic study. The fractured permeable formatiodis represented by a simplified conceptual model incorporating two sets of oblique and parallel fractures embedded in permeable blocks. The domain is initially divided into a completely freshwater zone overlying a completely saline water zone. However, the sharp interface originally existing between the saline and fresh water is subject to dispersion because of mixing in fracture intersections and between the fracture flow and the permeable block flow. Simulations based on the use of an appropriate numerical model developed in the present study have helped us to characterize dispersion of the sharp interface and creation of the transition zone. Relationships of transverse and longitudinal dispersion in the domain are also determined.

Journal of Hydrology

Seasonal cycles of dissolved constituents in streamwater in two forested catchments in the mid-Atlantic region of the eastern U.S.A.

Streamwater discharge and chemistry of two small catchments on Catoctin Mountain in north-central Maryland have been monitored since 1982. Repetitive seasonal cycles in stream-water chemistry have been observed each year, along with seasonal cycles in the volume of stream discharge and in groundwater levels. The hypothesis that the observed streamwater chemical cycles are related to seasonal changes in the hydrological flow paths that contribute to streamflow is examined using a combination of data on groundwater levels, shallow and deep groundwater chemistry, streamwater discharge, streamwater chemistry, soil-water chemistry, and estimates of water residence times. The concentrations of constituents derived from rock weathering, particularly bicarbonate and silica, increase in streamwater during the summer when the water table is below the regolith-bedrock interface and stream discharge consists primarily of deep groundwater from the fractured-bedrock aquifer. Conversely, the concentrations in streamwater of atmospherically derived components, particularly sulfate, increase in winter when the water table is above the regolith-bedrock interface and stream discharge consists primarily of shallow groundwater from the regolith. Tritium and chlorofluorocarbon (CFC) measurements suggest that the groundwater in these systems is young, with a residence time of less than several years. The results of this study have implications for the design of large-scale water-quality monitoring programs.

Maryland

Ore microscopy of the Paoli silver-copper deposit, Oklahoma

The Paoli silver-copper deposit is located in south-central Oklahoma, 56 km south-southeast from Norman, Oklahoma. It was mined for high-grade silver-copper near the beginning of this century, and intensive exploratory drilling during the early 1970's delineated unmined portions of the deposit. A collaborative study between the U.S.G.S., the Kansas Geological Survey, and the University of Missouri-Rolla was undertaken to provide new information on the character of red bed copper deposits of the Midcontinent region. The Paoli deposit has been interpreted to occur as a roll-front type of deposit. The silver and copper mineralization occurs within paleochannels in the Permian Wellington Formation. The silver-copper interfaces appear to be controlled by oxidation-reduction interfaces that are marked by grey to red color changes in the host sandstone. Ore microscopic examinations of polished thin sections show that unoxidized ore consists of chalcocite, digenite, chalcopyrite, covellite and pyrite; and oxidized ores are characterized by covellite, bornite, hematite and goethite. In sandstone-hosted ores, chalcocite and digenite replace dolomite and border clastic quartz grains. In siltstone-hosted ores, the copper sulfide grains have varied shapes; most are irregular in shape and 5–25 μm across, others have euhedral shapes suggestive of pyrite crystal replacements, and some are crudely spherical and are 120–200 μm across. Chalcopyrite is the predominant copper sulfide at depth. Covellite and malachite replace chalcocite and digenite near the surface. Silver only occurs as native silver; most as irregularly shaped grains 40–80 μm across, but some as cruciform crystals that are up to 3.5 mm across. The native silver has been deposited after copper sulfides, and locally replaces chalcocite. Surficial nodules of pyrite, malachite and hematite locally are present in outcrops at the oxidation-reduction fronts. Polished sections of the nodules show that malachite forms a cement around quartz sand grains, and brecciated pyrite grains are surrounded by rims of hematite and goethite. Dolomite is the principal sandstone cement. Cathodoluminescence microscopic study of the mineral has shown that it was deposited during seven periods before the copper sulfide mineralization.

Oklahoma

Distribution and effects of shallow gas on bulk estuarine sediment properties

Gas bubble are present in sediments covering approximately 30% of the main stem of Chesapeake Bay, with bubbles occurring at the sediment-water interface in 18% of the main stem sediments. This biogenic gas is found either in the sediments in the lower salinity reaches of the Bay, or confined to sediments which overline infilled palaeodrainage channels formed during the Wisconsinan low sea level stand (approximately 18 ka). Gas associated with the old drainage network does not correlate with present bathymetry or sedimentological patterns. Some differences between the gas-charged and gas-free sediments are: (1) gas-charged sediments have water contents 10-20% higher than comparable gas-free cores; (2) organic matter is better presented with depth in the gas-charged sediments (upwards of 60% more at one depth); (3 monosulphides are dominant sulphide mineral phase within the gas-charged sediments, comprising over 40% of the total sulphur. Within the gas-free sediments monosulphides are significant only near the sediment-water interface and rapidly become negligible with depth, and; (4) cores of gas-charged sediments are highly colour-banded due to preservation of sulphide mineral variations, while gas-free cores are diagenetically altered to pyrite.

Maryland

Use of geoelectrical methods in groundwater pollution surveys in a coastal environment

The pollution of coastal aquifers by old landfills can contaminate valuable and scarce water resources in the freshwater lens utilized seasonably by overcrowded communities. The pollutants will ultimately flow into the sea where they may also cause a coastal water pollution problem. We have detected pollution in the freshwater lens from a sanitary landfill near Provincetown, Cape Cod, using the geoelectrical resistivity method. This survey included Schlumberger geoelectrical depth soundings and a horizontal geoelectrical profile using the Wenner configuration. The geoelectrical survey was conducted at a site along Highway 6 where it passes the coastal town of Provincetown and a sanitary landfill that has been in operation since 1954. The depth soundings suggest the characteristic decrease in resistivity vs. depth from the high resistivity of the unsaturated zone to the low resistivity of the saltwater saturated zone. The freshwater lens is clearly identified by the change in slope of the steeply dipping curve of resistivity versus electrode spacing. Interpretations made using a multilayer program, Geomate, resulted in layer resistivities between 460 and 95 ohm?m for the freshwater lens. A comparison with well water resistivities suggests that a layer resistivity of 230 ohm?m or lower is indicative of pollution in the freshwater lens. The results of the geoelectrical depth soundings were confirmed in the Wenner horizontal profile. Both measurements suggest that the pollutants do not spread evenly as one would expect for a homogeneous and isotropic medium. Instead, a preferred channel for the flow of the pollutants is observed along a path from the landfill toward the shoreline. The depth to the saltwater/freshwater interface or, more specifically, to the low resistivity-high resistivity interface appears to be shallow where the freshwater lens is polluted. This was confirmed by pore water well samples that were highly mineralized. The equilibrium postulated by the Ghyben-Herzberg relation appears to be disturbed in the area of aquifer pollution. This rise in the conductivity boundary is caused by the highly mineralized bottom of the contaminant plume that submerges into the saltwater saturated zone. In the area of high freshwater pollution the groundwater can be subdivided into three layers that show a decrease in resistivity with depth. The formation factor, F, defined as the ratio of bulk aquifer resistivity to pore water resistivity, shows unusually high values between 10 and 12. These high values are unexpected for an unconsolidated sand. Pollution residues are suspected to clog the pores and thus to increase the resistivity. It is possible that iron-oxidizing bacteria and the precipitation of dissolved iron or organic pollutants are the cause of the high values of F. If proven correct, these interesting possibilities could lead to future new applications of the geoelectrical resistivity method in contaminant hydroloy.

Journal of Applied Geophysics

Wetlands under global change

Wetlands are among the ecosystem types most threatened by global change, including both climate change and other anthropogenic factors such as sea level rise, urban development, deforestation, agricultural land use, drainage, levees, tidal flow restrictions, pollution, eutrophication, and fires. Wetlands not only store disproportionate amounts of carbon compared to other terrestrial ecosystems, but they lie at the terrestrial-aquatic interface crucial to understanding landscape and global scale biogeochemical cycles. In this chapter, we focus on the major global change factors affecting wetlands and the responses of different wetland types to those global change factors. Special attention is given to direct responses to increasing atmospheric carbon dioxide levels. Because of their hydrological connections and placement at the terrestrial-aquatic interface, the conservation of wetlands involves accounting for uncertainties related to interacting stressors. While the past decades have seen many important experimental and observational studies of wetland responses to global change factors, large uncertainties remain, especially within tropical regions where even the basic extent of wetland ecosystems is not well documented.

Book chapter

Modeling Carbon Dioxide, pH and Un-Ionized Ammonia Relationships in Serial Reuse Systems

In serial reuse systems, excretion of metabolic carbon dioxide has a significant impact on ambient pH, carbon dioxide, and un-ionized ammonia concentrations. This impact depends strongly on alkalinity, water flow rate, feeding rate, and loss of carbon dioxide to the atmosphere. A reduction in pH from metabolic carbon dioxide can significantly reduce the un-ionized ammonia concentration and increase the carbon dioxide concentrations compared to those parameters computed from influent pH. The ability to accurately predict pH in serial reuse systems is critical to their design and effective operation. A trial and error solution to the alkalinity–pH system was used to estimate important water quality parameters in serial reuse systems. Transfer of oxygen and carbon dioxide across the air–water interface, at overflow weirs, and impacts of substrate-attached algae and suspended bacteria were modeled. Gas transfer at the weirs was much greater than transfer across the air–water boundary. This simulation model can rapidly estimate influent and effluent concentrations of dissolved oxygen, carbon dioxide, and un-ionized ammonia as a function of water temperature, elevation, water flow, and weir type. The accuracy of the estimates strongly depends on assumed pollutional loading rates and gas transfer at the weirs. The current simulation model is based on mean daily loading rates; the impacts of daily variation loading rates are discussed. Copies of the source code and executable program are available free of charge.

Aquacultural Engineering