Development of a mathematical and computer model to assess water quality impacts of hazardous compounds from minor gasoline spills in the inter-tidal zone of the Miami River, Florida

The objective of this study was to develop a model to predict transport and fate of gasoline components of environmental concern in the Miami River by mathematically simulating the movement of dissolved benzene, toluene, xylene (BTX), and methyl-tertiary-butyl ether (MTBE) occurring from minor gasoline spills in the inter-tidal zone of the river. Computer codes were based on mathematical algorithms that acknowledge the role of advective and dispersive physical phenomena along the river and prevailing phase transformations of BTX and MTBE. Phase transformations included volatilization and settling. ^ The model used a finite-difference scheme of steady-state conditions, with a set of numerical equations that was solved by two numerical methods: Gauss-Seidel and Jacobi iterations. A numerical validation process was conducted by comparing the results from both methods with analytical and numerical reference solutions. Since similar trends were achieved after the numerical validation process, it was concluded that the computer codes algorithmically were correct. The Gauss-Seidel iteration yielded at a faster convergence rate than the Jacobi iteration. Hence, the mathematical code was selected to further develop the computer program and software. The model was then analyzed for its sensitivity. It was found that the model was very sensitive to wind speed but not to sediment settling velocity. ^ A computer software was developed with the model code embedded. The software was provided with two major user-friendly visualized forms, one to interface with the database files and the other to execute and present the graphical and tabulated results. For all predicted concentrations of BTX and MTBE, the maximum concentrations were over an order of magnitude lower than current drinking water standards. It should be pointed out, however, that smaller concentrations than the latter reported standards and values, although not harmful to humans, may be very harmful to organisms of the trophic levels of the Miami River ecosystem and associated waters. This computer model can be used for the rapid assessment and management of the effects of minor gasoline spills on inter-tidal riverine water quality. ^

Sediment Dynamics of a Shallow Hypereutrophic lake: Lake Jesup, Florida, USA

Improved knowledge of sediment dynamics within a lake system is important for understanding lake water quality. This research was focused on an assessment of the vertical sediment flux in Lake Jesup, a shallow (1.3 m average depth) hypereutrophic lake of central Florida. Sediment dynamics were assessed at varying time scales (daily to weekly) to understand the transport of sediments from external forces; wind, waves, precipitation and/or runoff. Four stations were selected within the lake on the basis of water depth and the thicknesses of unconsolidated (floc) and consolidated sediments. At each of these stations, a 10:1 (length to diameter) high aspect ratio trap (STHA) was deployed to collect particulate matter for a one to two week period. The water and sediment samples were collected and analyzed for total carbon (TC), total phosphorus (TP) and total nitrogen (TN). Mass accumulation rates (MAR) collected by the traps varied from 77 to 418 g m-2 d-1 over seven deployments. TN, TP and TC sediment concentrations collected by the traps were consistently higher than the sediments collected by coring the lake bottom and is most likely associated with water column biomass. A yearly nutrient budget was determined from August 2009 to August 2010 with flux calculated as 2,033,882 mt yr-1.

Forecasting Responses of Seagrass Distributions to Changing Water Quality Using Monitoring Data

Extensive data sets on water quality and seagrass distributions in Florida Bay have been assembled under complementary, but independent, monitoring programs. This paper presents the landscape-scale results from these monitoring programs and outlines a method for exploring the relationships between two such data sets. Seagrass species occurrence and abundance data were used to define eight benthic habitat classes from 677 sampling locations in Florida Bay. Water quality data from 28 monitoring stations spread across the Bay were used to construct a discriminant function model that assigned a probability of a given benthic habitat class occurring for a given combination of water quality variables. Mean salinity, salinity variability, the amount of light reaching the benthos, sediment depth, and mean nutrient concentrations were important predictor variables in the discriminant function model. Using a cross-validated classification scheme, this discriminant function identified the most likely benthic habitat type as the actual habitat type in most cases. The model predicted that the distribution of benthic habitat types in Florida Bay would likely change if water quality and water delivery were changed by human engineering of freshwater discharge from the Everglades. Specifically, an increase in the seasonal delivery of freshwater to Florida Bay should cause an expansion of seagrass beds dominated by Ruppia maritima and Halodule wrightii at the expense of the Thalassia testudinum-dominated community that now occurs in northeast Florida Bay. These statistical techniques should prove useful for predicting landscape-scale changes in community composition in diverse systems where communities are in quasi-equilibrium with environmental drivers.

Historic primary producer communities linked to water quality and hydrologic changes in the northern Everglades

The northern Everglades Water Conservation Areas have experienced recent ecological shifts in primary producer community structure involving marl periphyton mats and dense Typha-dominated macrophyte stands. Multiple investigations have identified phosphorus (P) as a driver of primary producer community structure, but effects of water impoundment beginning in the 1950s and changes in water hardness [e.g., (CaCO3)] have also been identified as a concern. In an effort to understand pre-1950, primary producer community structure and identify community shifts since 1950, we measured pigment proxies on three sediment cores collected in Water Conservation Area-2A (WCA-2A) along a phosphorus enrichment gradient. Photosynthetic pigments, sediment total phosphorus content (TP), organic matter, total organic carbon and nitrogen were used to infer historic primary producer communities and changes in water quality and hydrology regulating those communities. Excess 210Pb was used to establish historic dates for the sediment cores. Results indicate the northern area of WCA-2A increased marl deposition and increased algal abundance ca. 1920. This increase in (presumably) calcareous periphyton before intensive agriculture and impoundment suggest canal-derived calcium inputs and to some extent early drainage effects played a role in initiating this community shift. The northern area community then shifted to Typha dominance around 1965. The areas to the south in WCA-2A experienced increased marl deposition and algal abundance around or just prior to 1950s impoundment, the precise timing limited by core age resolution. Continued increases in algal abundance were evident after 1950, coinciding with impoundment and deepening of canals draining into WCA-2A, both likely increasing water mineral and nutrient concentrations. The intermediate site developed a Typha-dominated community ca. 1995 while the southern-most core site WCA-2A has yet to develop Typha dominance. Numerous studies link sediment TP >650 mg P/kg to marsh habitat degradation into Typha-dominance. The northern and intermediate cores where Typha is currently support this previous research by showing a distinct shift in the sediment record to Typha dominance corresponding to sediment TP between 600 and 700 mg P/kg. These temporal and spatial differences are consistent with modern evidence showing water-column gradients in mineral inputs (including Ca, carbonates, and phosphorus) altering primary producer community structure in WCA-2A, but also suggest hydroperiod has an effect on the mechanisms regulating periphyton development and Typha dominance.

Compilation of downward flux observations from sediment trap deployments in the Atlantic Ocean - Contribution to EURO-BASIN's Data integration

We provide a compilation of downward fluxes (total mass, POC, PON, BSiO2, CaCO3, PIC and lithogenic/terrigenous fluxes) from over 6000 sediment trap measurements distributed in the Atlantic Ocean, from 30 degree North to 49 degree South, and covering the period 1982-2011. Data from the Mediterranean Sea are also included. Data were compiled from different sources: data repositories (BCO-DMO, PANGAEA), time series sites (BATS, CARIACO), published scientific papers and/or personal communications from PI's. All sources are specifed in the data set. Data from the World Ocean Atlas 2009 were extracted to provide each flux observation with contextual environmental data, such as temperature, salinity, oxygen (concentration, AOU and percentage saturation), nitrate, phosphate and silicate.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

Fragilariopsis kerguelensis images from sediment core PS1768-8

Background: Light microscopic analysis of diatom frustules is widely used both in basic and applied research, notably taxonomy, morphometrics, water quality monitoring and paleo-environmental studies. In these applications, usually large numbers of frustules need to be identified and / or measured. Although there is a need for automation in these applications, and image processing and analysis methods supporting these tasks have previously been developed, they did not become widespread in diatom analysis. While methodological reports for a wide variety of methods for image segmentation, diatom identification and feature extraction are available, no single implementation combining a subset of these into a readily applicable workflow accessible to diatomists exists. Results: The newly developed tool SHERPA offers a versatile image processing workflow focused on the identification and measurement of object outlines, handling all steps from image segmentation over object identification to feature extraction, and providing interactive functions for reviewing and revising results. Special attention was given to ease of use, applicability to a broad range of data and problems, and supporting high throughput analyses with minimal manual intervention. Conclusions: Tested with several diatom datasets from different sources and of various compositions, SHERPA proved its ability to successfully analyze large amounts of diatom micrographs depicting a broad range of species. SHERPA is unique in combining the following features: application of multiple segmentation methods and selection of the one giving the best result for each individual object; identification of shapes of interest based on outline matching against a template library; quality scoring and ranking of resulting outlines supporting quick quality checking; extraction of a wide range of outline shape descriptors widely used in diatom studies and elsewhere; minimizing the need for, but enabling manual quality control and corrections. Although primarily developed for analyzing images of diatom valves originating from automated microscopy, SHERPA can also be useful for other object detection, segmentation and outline-based identification problems.

d13C and d15N stable isotopes values, total nitrogen and C/N ratios of sediment core GeoB4801-1

Rising stable nitrogen isotope ratios (d15N) in dated sediment records of the German Bight/SE North Sea track river-induced coastal eutrophication over the last 2 centuries. Fully exploiting their potential for reconstructions of pristine conditions and quantitative analysis of historical changes in the nitrogen cycle from these sediment records requires knowledge on processes that alter the isotopic signal in non-living organic matter (OM) of sinking particles and sediments. In this study, we analyze the isotopic composition of particulate nitrogen (PN) in the water column during different seasons, in surface sediments, and in sediment cores to assess diagenetic influences on the isotopic composition of OM. Amino acid (AA) compositions of suspended matter, surface sediments, and dated cores at selected sites of the German Bight serve as indicators for quality and degradation state of PN. The d15N of PN in suspended matter had seasonal variances caused by two main nitrate sources (oceanic and river) and different stages of nitrate availability during phytoplankton assimilation. Elevated d15N values (> 20 per mil) in suspended matter near river mouths and the coast coincide with a coastal water mass receiving nitrate with elevated isotope signal (d15N > 10 per mil) derived from anthropogenic input. Particulate nitrogen at offshore sites fed by oceanic nitrate having a d15N between 5 and 6 per mil had low d15N values (< 2 per mil), indicative of an incipient phytoplankton bloom. Surface sediments along an offshore-onshore transect also reflect the gradient of low d15N of nitrate in offshore sites to high values near river mouths, but the range of values is smaller than between the end members listed above and integrates the annual d15N of detritus. Sediment cores from the coastal sector of the gradient show an increasing d15N trend (increase of 2.5 per mil) over the last 150 years. This is not related to any change in AA composition and thus reflects eutrophication. The d15N signals from before AD 1860 represent a good estimation of pre-industrial isotopic compositions with minimal diagenetic overprinting. Rising d13C in step with rising d15N in these cores is best explained by increasing productivity caused by eutrophication.

Pollen counts of sediment core MFM6

In recent years enormous success has been achieved in varve counting of the Eifel maar lakes, but a detailed correlation with the biostratigraphy has been missing. In this paper, we present new palynological results of the Lateglacial sequences from Holzmaar Lake and Meerfelder Maar Lake based on annually laminated sediments. In particular, the Meerfelder Maar has great potential, because, in contrast to the Holzmaar, the sequence between the Ulmener Maar Tephra (11 000 varve years BP) and the Laacher See Tephra (12 880 varve years BP) including the Younger Dryas is undisturbed and complete. Therefore, we currently use the Meerfelder Maar chronology (Brauer et al., 1999b) as an independent varve calendar for the biostratigraphy of the Lateglacial. The palynological signals of both maar lakes are in good agreement and can easily be correlated with one another and with type sections/type regions in northwestern Germany and Jutland. The sequences of the Eifel maar lakes have the quality of hypostratotypes with regional biozones based on an absolute time scale.

Pollen analysis of sediment profile WO1 from Wollingster See

A palynological study of a 15 m sediment core from the centre of Lake Wollingst (water depth 14,5 m) is presented. The pollen record shows 3 lateglacial thermomers, called Meiendorf, Bölling, Alleröd and the early holocene Friesland-Thermomer. The succession of forest vegetation taking place on the lake surroundings during the Holocene was typical for older moraine soils which are poor in nutrients: forest vegetation started with birch and pine, followed by hazel, oak and elm in the Boreal and by alder, lime and ash-tree in the Atlantic. Beech and hornbeam reached the area during Subboreal. However, due to the poor soils they spread out only after the Iron Age. With the deforestation during the medieval time the lake lost its character of a primeval forest lake. Lake Wollingst was oligotrophic since its origin at the end of the Pleniglacial. After medieval forest-clearing the lake has changed its quality of water particularly in connection with hemp- and flax-rotting. The modem sediments in this profile are completely disturbed. They contain reworked material, a lot of blue-green algae and remains of Bosmina longirostris indicating eutrophic conditions.

Particle flux measured on deep sea sediment trap VP-2_trap (Appendix A2.7)

A 17 month record of vertical particle flux of dry weight, carbonate and organic carbon were 25.8, 9.4 and 2.4g/m**2/y, respectively. Parallel to trap deployments, pelagic system structure was recorded with high vertical and temporal resolution. Within a distinct seasonal cycle of vertical particle flux, zooplankton faecal pellets of various sizes, shapes and contents were collected by the traps in different proportions and quantities throughout the year (range: 0-4,500 10**3/m**2/d). The remains of different groups of organisms showed distinct seasonal variations in abundance. In early summer there was a small maximum in the diatom flux and this was followed by pulses of tinntinids, radiolarians, foraminiferans and pteropods between July and November. Food web interactions in the water column were important in controlling the quality and quantity of sinking materials. For example, changes in the population structure of dominant herbivores, the break-down of regenerating summer populations of microflagellates and protozooplankton and the collapse of a pteropod dominated community, each resulted in marked sedimentation pulses. These data from the Norwegian Sea indicate those mechanisms which either accelerate or counteract loss of material via sedimentation. These involve variations in the structure of the pelagic system and they operatè on long (e.g. annual plankton succession) and short (e.g. the end of new production, sporadic grazing of swarm feeders) time scales. Connecting investigation of the water column with a high resolution in time in parallel with drifting sediment trap deployments and shipboard experiments with the dominant zooplankters is a promising approach for giving a better understanding of both the origin and the fate of material sinking to the sea floor.

In situ characterization of sediment oxygen demand across the spectrum of non-resuspension to full resuspension with varying bed shear stress

Sediment oxygen demand (SOD) can be a significant oxygen sink in various types of water bodies, particularly slow-moving waters with substantial organic sediment accumulation. In most settings where SOD is a concern, the prevailing hydraulic conditions are such that the impact of sediment resuspension on SOD is not considered. However, in the case of Bubbly Creek in Chicago, Illinois, the prevailing slack water conditions are interrupted by infrequent intervals of very high flow rates associated with pumped combined sewer overflow (CSO) during intense hydrologic events. These events can cause resuspension of the highly organic, nutrient-rich bottom sediments, resulting in precipitous drawdown of dissolved oxygen (DO) in the water column. While many past studies have addressed the dependence of SOD on near-bed velocity and bed shear stress prior to the point of sediment resuspension, there has been limited research that has attempted to characterize the complex and dynamic phenomenon of resuspended-sediment oxygen demand. To address this issue, a new in situ experimental apparatus referred to as the U of I Hydrodynamic SOD Sampler was designed to achieve a broad range of velocities and associated bed shear stresses. This allowed SOD to be analyzed across the spectrum of no sediment resuspension associated with low velocity/ bed shear stress through full sediment resuspension associated with high velocity / bed shear stress. The current study split SOD into two separate components: (1) SODNR is the sediment oxygen demand associated with non-resuspension conditions and is a surface sink calculated using traditional methods to yield a value with units (g/m2/day); and (2) SODR is the oxygen demand associated with resuspension conditions, which is a volumetric sink most accurately characterized using non-traditional methods and units that reflect suspension in the water column (mg/L/day). In the case of resuspension, the suspended sediment concentration was analyzed as a function of bed shear stress, and a formulation was developed to characterize SODR as a function of suspended sediment concentration in a form similar to first-order biochemical oxygen demand (BOD) kinetics with Monod DO term. The results obtained are intended to be implemented into a numerical model containing hydrodynamic, sediment transport, and water quality components to yield oxygen demand varying in both space and time for specific flow events. Such implementation will allow evaluation of proposed Bubbly Creek water quality improvement alternatives which take into account the impact of SOD under various flow conditions. Although the findings were based on experiments specific to the conditions in Bubbly Creek, the techniques and formulations developed in this study should be applicable to similar sites.