Episodic weak stratification is a persistent and important feature of the epilimnion of Lake Opeongo. Field studies were conducted in Lake Opeongo in 2009 and 2010 to assess the effect of the epilimnetic weak thermal stratification on turbulent mixing and ecological processes. Near surface thermoclines (as defined by dT/dz > 0.2 ̊C m-1 between 1 m and 2 m) occur approximately 20 % of the time and correspond to periods of high values of gradient Richardson number (Rig). Extended periods of near surface stratification (> 1 hr), account for more than 80% of the stratified period. We compare these findings with previous observations from the Experimental Lakes Area and discuss the biological implications of near surface thermoclines. The persistent weak temperature stratification that characterizes the entire epilimnion of Lake Opeongo, acts to reduce the magnitude of turbulent mixing. During these stratified periods, the values of Rig, increase, with a corresponding decrease in the rates of dissipation of turbulent kinetic energy, (ε ), the turbulence activity parameter, I = ε/νN2, and vertical eddy diffusivity (Kz). Mixing during cold fronts occur over time scales of minutes to hours, which work to erode diurnal thermoclines. The direct implications of weak stratification on aquatic organisms are also assessed. The presence of weak stratification also allows for the existence of internal waves within the epilimnion. A key observation in this thesis is the relationship between enhanced small-scale spatial variability in zooplankton distribution and the presence of internal waves in the weakly stratified epilimnion. To quantify this physical-biological coupling, we compare variance of isotherm displacement and gradient Richardson number (Rig) with small-scale spatial distributions of zooplankton. For smaller size ranges of zooplankton (284-450 μm), we find that spatial variability is statistically greatest for the highest values of variability of isotherm displacement. As vertical velocities estimated from wave characteristics are faster than swimming speeds of small zooplankton, these organisms become passively advected by the internal waves leading to increased spatial variability.

Many deglaciated watersheds possess complex longitudinal profiles and spatially variable sediment sources derived from the last glaciation. Accordingly, downstream patterns in fluvial forms and characteristics of the channel and floodplain may diverge from conventional ‘graded’ watershed models where channel slope decreases downstream. Within a small, deglaciated watershed, five distinct fluvial forms were identified and it was found that the watershed is organized in a spatially variable, but generally inverted arrangement with stable, low-energy forms dominating the upper watershed and more dynamic, high-energy forms in the lower watershed. As a consequence of this inverted and variable downstream succession, and the spatially erratic organization of glacial sediment sources, downstream trends in channel and floodplain characteristics are poorly defined, and in many cases, deviate from expected trends. This is most true for sedimentological variables where multiple co-dominant controls exist, while morphological and hydrological variables better conform to expected downstream trends.

Many deglaciated watersheds possess complex longitudinal profiles and spatially variable sediment sources derived from the last glaciation. Accordingly, downstream patterns in fluvial forms and characteristics of the channel and floodplain may diverge from conventional ‘graded’ watershed models where channel slope decreases downstream. Within a small, deglaciated watershed, five distinct fluvial forms were identified and it was found that the watershed is organized in a spatially variable, but generally inverted arrangement with stable, low-energy forms dominating the upper watershed and more dynamic, high-energy forms in the lower watershed. As a consequence of this inverted and variable downstream succession, and the spatially erratic organization of glacial sediment sources, downstream trends in channel and floodplain characteristics are poorly defined, and in many cases, deviate from expected trends. This is most true for sedimentological variables where multiple co-dominant controls exist, while morphological and hydrological variables better conform to expected downstream trends.

Master of Science / Department of Geology / Saugata Datta / More than 75 million people in the Bengal Delta of eastern India and Bangladesh are exposed to drinking water with dangerously high arsenic (As) concentrations; the worst case of environmental poisoning in human history. Despite recognition of dangers posed to chronic exposure to drinking water with elevated As, its biogeochemical cycle is inadequately constrained in groundwater flow systems due to its complex redox chemistry and microbially-mediated transformations. Arsenic concentrations in Bengal Delta sediments are comparable to global averages, but its highly heterogeneous spatial distribution (on scales of meters to kilometers) in sediments and groundwaters is poorly understood. Though many research efforts have targeted understanding this heterogeneity in Bangladesh, less work has been done in eastern India. Murshidabad (23°56.355‘N, 88°16.156‘E), an eastern district in West Bengal, India, where groundwaters are highly As-affected (~4000 μg/l), was chosen as our study area. Research objectives were: (1) characterize sediment cores (mineralogically, geochemically) and groundwaters (hydrochemically, isotopically) in areas with contrasting As concentrations—west (low-As) and east (high-As) of river Bhagirathi, a major distributary of Ganges flowing through the heart of Murshidabad; (2) describe and understand the extent of spatial variability, laterally and vertically, of dissolved As concentrations in shallow (< 60 m) aquifers, comparing sediment core chemistry to water chemistry; (3) identify source(s) of aquifer recharge and (4) role(s) of inorganic carbon within the aquifer to understand the bioavailability and mobilization of As from sediments to groundwaters. Mineralogical differences between high-As (grey) and low-As (orange-brown) sediments, were the presence of greater amounts of micas, Fe- and Mg-rich clays, amphiboles, carbonates, and apatite in high-As sediments; these were virtually absent from low-As sediments. In high-As areas, As was associated with amorphous and poorly-crystalline Fe-oxyhydroxide phases and labile (specifically-sorbed) phases, especially where Fe(II):Fe[subscript]T was high in the sediments. High-As groundwaters had high As(III):As[subscript]T, iron, bicarbonate, phosphate, and ammonium, and low concentrations of chloride and sulfate. Dry season precipitation was probably the main source of aquifer recharge; lighter values of [superscript]13C in dissolved inorganic carbon resulted from oxidation of natural organic matter. This study points to an idea that both microbially-mediated oxidation-reduction and competitive ion-exchange processes occurring in shallow aquifers of Murshidabad drive As mobilization and sequestration by aquifer sediments.

Doctor of Philosophy / Department of Agronomy / Gerard J. Kluitenberg / The Parsons soil has a sharp increase in clay content from the upper teens in the A horizon to the mid fifties in the Bt horizon. The high clay content continues to the parent material resulting in 1.5 m of dense, slowly permeable subsoil over shale residuum. This project was designed to better understand soil-water management needs of this soil. The main objective was to determine a comprehensive hydrologic balance for the claypan soil. Specific objectives were a) to determine effect of tillage management on select water balance components including water storage and evaporation, b) to quantify relationship between soil water status and crop variables such as emergence and yield, and c) to verify balance findings with predictions from a mechanistic model, specifically HYDRUS 1-D. The study utilized three replicates of an ongoing project in Labette County, Kansas in which till and no-till plots had been maintained in a sorghum [Sorghum bicolor (L.) Moench] – soybean [Glycine max (L.) Merr.] rotation since 1995. Both crops are grown each year in a randomized complete block design. The sorghum plots were equipped with Time Domain Reflectometry (TDR) probes to measure A horizon water content and neutron access tubes for measurement of water throughout the profile. Precipitation, evaporation, and perched water depth were determined at the field scale. Drainage was estimated as negligible after performing hydraulic conductivity measurements on the clayey subsoil. Runoff was determined as the residual in this water balance. Cumulative differences in the hydrologic balances as a result of tillage management were found to be minimal over an entire growing season. However, tillage treatment differences were seen in early season evaporation, surface water content, and the resulting residual runoff values. The chisel-disk treatments had greater evaporation leading to reduced runoff when compared with no-till. There was interaction between tillage treatment and time for surface water content measurements. No effect of tillage treatment was found for whole-profile water content. Crop variables were unaffected by tillage other than the first days emergence, and first days tillering being greater for chisel-disk treatments. No correlation between stored water and crop variables could be found. All aspects of field measurement were well supported by the predictions of the HYDRUS 1-D model.

Current estimates of terrestrial carbon (C) fluxes overlook explicit hydrological controls. In this research project, a spatially explicit hydro-ecological model, BEPS-TerrainLab V2.0 was further developed to improve our understanding of the non-linearities associated with various hydro-ecological processes. A modeling study was conducted in a humid boreal ecosystem in north central Quebec, Canada. The sizes and nature of various ecosystem-C-pools were comprehensively reconstructed under a climate change and disturbance scenario prior to simulation in order to ensure realistic biogeochemical modeling. Further, several ecosystem processes were simulated and validated using field measurements for two years. A sensitivity analysis was also performed. After gaining confidence in the model’s ability to simulate various hydrologically controlled ecophysiological and biogeochemical processes and having understood that topographically driven sub-surface baseflow is the main process determining the soil moisture regime in humid boreal ecosystem, its influence on ecophysiological and biogeochemical processes were investigated. Three modeling scenarios were designed that represent strategies that are currently used in ecological models to represent hydrological controls. These scenarios were: 1) Explicit, where realistic lateral water routing was considered 2) Implicit, where calculations were based on a bucket-modeling approach 3) NoFlow, where the lateral sub-surface flow was turned off in the model. In general, the Implicit scenario overestimated GPP, ET and NEP, as opposed to the Explicit scenario. The NoFlow scenario underestimated GPP and ET but overestimated NEP. The key processes controlling the differences were due to the combined effects of variations in plant physiology, photosynthesis, heterotrophic respiration, autotrophic respiration and nitrogen mineralization; all of which occurred simultaneously in different directions, at different rates, affecting the spatio-temporal distribution of terrestrial C-sources or sinks (NEP). From these results it was clear that lateral water flow does play a significant role in the net terrestrial C distribution and it was discovered that non-explicit forms of hydrological representations underestimate the sizes of terrestrial C-sources rather than C-sinks. The scientific implication of this work demonstrates that regional or global scale terrestrial C estimates could have significant errors if proper hydrological constraints are not considered for modeling ecological processes due to large topographic variations of the Earth’s surface.

I studied how geology, land use and rainfall, correlated with peak flow responses in 110 headwater stream sites during a drought year. Highest discharges were observed in the most developed catchments and in the most poorly drained soils, but specific responses were variable depending on both geology and land disturbance. Redundancy analysis indicated that both surficial geology and land disturbance were important predictors of discharge and that rainfall was in general a poor predictor of discharge. I conclude that responses of headwater streams to individual storms are unpredictable from data generated using GIS, but increased peak flows occur associated with human development, mitigated by surficial geology. The headwater streams that are most vulnerable to flow alterations occur on poorly drained soils, and where urbanization tends to concentrate. Much greater attention to managing water is required if further degradation of stream ecosystems is to be prevented from our future land use.

Upland environments are important sources of mercury (Hg) to downstream wetlands and water bodies. Hydrology is instrumental in facilitating Hg transport within, and export from watersheds. Two complementary studies were conducted to assess the role hydrological processes play in controlling Hg mobility and transport in forested uplands. A field study compared runoff and Hg fluxes from three, replicate hillslope plots during two contrasting spring snowmelt periods, in terms of snowpack depth and timing. Hillslope Hg fluxes were predominately flow-driven. The melting of soil frost significantly delayed a large portion of the Hg flux later into the spring following a winter with minimal snow accumulation. A microcosm laboratory study using a stable Hg isotope tracer applied to intact soil cores investigated the relative controls of soil moisture and precipitation on Hg mobility. Both hydrologic factors control the mobility of contemporary Hg; with greatest Hg flushing from dry soils under high-flow conditions.

Master of Science / Department of Computing and Information Sciences / Mitchell L. Neilsen / In recent years, headline news has been overwhelmed with stories about dam and levee failures including the 2005 levee breaches in New Orleans and the 2006 Kaloko Damfailure in Hawaii that resulted in seven deaths. Since 2000, state and federal agencies have reported 92 dam failures in the United States to the National Performance of Dams Program. Incidents such as these have brought both national and worldwide attention to the need for improved flood warning systems and breach prediction tools for dam embankment and levee failures. (G. J. Hanson, 2010) IESIMBA 2012 is an upgraded version of SIMBA, which has been upgraded from VB6 to C#.NET. The Microsoft Windows-based SIMplified Breach Analysis software (SIMBA) was developed by the USDA Agricultural Research Service in cooperation with Kansas State University. The software was developed for the purpose of analyzing internal erosion, earth embankment breach test data and extending the understanding of the underlying physical processes of breach of an overtopped earth embankment. It is a research tool that is modified routinely to test the sensitivity of the output to various sub-models and assumptions. This software is a test version for use in validation testing of the simplified breach model based on stress and mass failure driven headcut movement. It runs under Microsoft Windows 98SE, Windows 2000, NT, XP, or Vista. The following Input Screens are used to guide the user through development of input data sets.  Model Properties , Dam Profile , Structure Table, Spillway Rating and Hydrograph Data After an input data set has been entered, the data is saved and simulation can be performed on the data stored in memory at any time by selecting Build option. Input and output files are stored in a fixed ASCII text format. The results of the simulation can be viewed in graphical format which are of interest to the researchers at Oklahoma State University, Stillwater by selecting View option.

Little is known about the enzyme-substrate interactions occurring during the dechlorination of chloroform (CF) and 1,1,1-trichloroethane (1,1,1-TCA) by the enrichment culture containing Dehalobacters, hereafter called DHB-CF/MEL. Compound specific isotope analysis (CSIA) is used to investigate the factors which may affect the isotope fractionation observed for CF and 1,1,1-TCA dechlorination. This thesis reports the first isotope enrichment factors observed for CF biodegradation at -27.5‰ ± 0.9‰, thus providing fundamental information for comparing isotope enrichment factors observed during trichlorinated alkane degradation by DHB-CF/MEL. The thesis also reports how the presence of CF and 1,1,1-TCA influences isotope fractionation and explores the possible influence of substrate inhibition on isotope fractionation during 1,1,1-TCA dechlorination. The data suggests that substrate inhibition during 1,1,1-TCA dechlorination by DHB-CF/MEL may not affect carbon isotope fractionation. The results suggest that CSIA is a promising monitoring tool even for the simultaneous biodegradation of CF and 1,1,1-TCA at different 1,1,1-TCA starting concentration.