Geochemical determination of the fate and transport of injected fresh wastewater to a deep saline aquifer

Deep well injection into non-potable saline aquifers of treated domestic wastewater has been used in Florida for decades as a safe and effective alternative to ocean outfall disposal. The objectives of this study were to determine the fate and transport of injected wastewater at two deep well injection sites in Miami Dade County, Florida, USA. Detection of ammonium in the Middle Confining units of the Floridan aquifer above the injection zone at both sites has been interpreted as evidence of upward migration of injected wastewater, posing a risk to underground sources of drinking water. Historical water quality data, including ammonia, chloride, temperature, and pH from existing monitoring wells at both sites from 1983 to 2008, major ions collected monthly from 2006 and 2008, and a synoptic sampling event for stable isotopes, tritium, and dissolved gases in 2008, were used to determine the source of ammonium in groundwater and possible migration pathways. Geochemical modeling was used to determine possible effects of injected wastewater on native water and aquifer matrix geochemistry. Injected wastewater was determined to be the source of elevated ammonium concentrations above ambient water levels, based on the results of major ion concentrations, tritium, dissolved noble gases and 15N isotopes analyses. Various possible fluid migration pathways were identified at the sites. Data for the south site suggest buoyancy-driven vertical pathways to overlying aquifers bypassing the confining units, with little mixing of injected wastewater with native water as it migrated upward. Once it is introduced into an aquifer, the injectate appeared to migrate advectively with the regional groundwater flow. Geochemical modeling indicated that CO 2-enriched injected wastewater allowed for carbonate dissolution along the vertical pathways, enhancing permeability along these flowpaths. At the north site, diffusive upward flow through the confining units or offsite vertical pathways were determined to be possible, however no evidence was detected for any on-site confining unit bypass pathway. No evidence was observed at either site of injected wastewater migration to the Upper Floridan aquifer, which is used as a municipal water supply and for aquifer storage and recovery.

Arsenic transport by zebrafish aquaglyceroporins

Background Arsenic is one of the most ubiquitous toxins and endangers the health of tens of millions of humans worldwide. It is a mainly a water-borne contaminant. Inorganic trivalent arsenic (AsIII) is one of the major species that exists environmentally. The transport of AsIII has been studied in microbes, plants and mammals. Members of the aquaglyceroporin family have been shown to actively conduct AsIII and its organic metabolite, monomethylarsenite (MAsIII). However, the transport of AsIII and MAsIII in in any fish species has not been characterized. Results In this study, five members of the aquaglyceroporin family from zebrafish (Danio rerio) were cloned, and their ability to transport water, glycerol, and trivalent arsenicals (AsIII and MAsIII) and antimonite (SbIII) was investigated. Genes for at least seven aquaglyceroporins have been annotated in the zebrafish genome project. Here, five genes which are close homologues to human AQP3, AQP9 and AQP10 were cloned from a zebrafish cDNA preparation. These genes were namedaqp3, aqp3l, aqp9a, aqp9b and aqp10 according to their similarities to the corresponding human AQPs. Expression of aqp9a, aqp9b, aqp3, aqp3l and aqp10 in multiple zebrafish organs were examined by RT-PCR. Our results demonstrated that these aquaglyceroporins exhibited different tissue expression. They are all detected in more than one tissue. The ability of these five aquaglyceroporins to transport water, glycerol and the metalloids arsenic and antimony was examined following expression in oocytes from Xenopus leavis. Each of these channels showed substantial glycerol transport at equivalent rates. These aquaglyceroporins also facilitate uptake of inorganic AsIII, MAsIII and SbIII. Arsenic accumulation in fish larvae and in different tissues from adult zebrafish was studied following short-term arsenic exposure. The results showed that liver is the major organ of arsenic accumulation; other tissues such as gill, eye, heart, intestine muscle and skin also exhibited significant ability to accumulate arsenic. The zebrafish larvae also accumulate considerable amounts of arsenic. Conclusion This is the first molecular identification of fish arsenite transport systems and we propose that the extensive expression of the fish aquaglyceroporins and their ability to transport metalloids suggests that aquaglyceroporins are the major pathways for arsenic accumulation in a variety of zebrafish tissues. Uptake is one important step of arsenic metabolism. Our results will contribute to a new understanding of aquatic arsenic metabolism and will support the use of zebrafish as a new model system to study arsenic associated human diseases.

Influence of hydrobiogeochemistry on transport of chromium through manganese -containing sediments: Experimental and modeling approaches

Chromium (Cr) is a metal of particular environmental concern, owing to its toxicity and widespread occurrence in groundwater, soil, and soil solution. A combination of hydrological, geochemical, and microbiological processes governs the subsurface migration of Cr. Little effort has been devoted to examining how these biogeochemical reactions combine with hydrologic processes influence Cr migration. This study has focused on the complex problem of predicting the Cr transport in laboratory column experiments. A 1-D reactive transport model was developed and evaluated against data obtained from laboratory column experiments. ^ A series of dynamic laboratory column experiments were conducted under abiotic and biotic conditions. Cr(III) was injected into columns packed with β-MnO 2-coated sand at different initial concentrations, variable flow rates, and at two different pore water pH (3.0 and 4.0). In biotic anaerobic column experiments Cr(VI) along with lactate was injected into columns packed with quartz sand or β-MnO2-coated sand and bacteria, Shewanella alga Simidu (BrY-MT). A mathematical model was developed which included advection-dispersion equations for the movement of Cr(III), Cr(VI), dissolved oxygen, lactate, and biomass. The model included first-order rate laws governing the adsorption of each Cr species and lactate. The equations for transport and adsorption were coupled with nonlinear equations for rate-limited oxidation-reduction reactions along with dual-monod kinetic equations. Kinetic batch experiments were conducted to determine the reduction of Cr(VI) by BrY-MT in three different substrates. Results of the column experiments with Cr(III)-containing influent solutions demonstrate that β-MnO2 effectively catalyzes the oxidation of Cr(III) to Cr(VI). For a given influent concentration and pore water velocity, oxidation rates are higher, and hence effluent concentrations of Cr(VI) are greater, at pH 4 relative to pH 3. Reduction of Cr(VI) by BrY-MT was rapid (within one hour) in columns packed with quartz sand, whereas Cr(VI) reduction by BrY-MT was delayed (57 hours) in presence of β-MnO 2-coated sand. BrY-MT grown in BHIB (brain heart infusion broth) reduced maximum amount of Cr(VI) to Cr(III) followed by TSB (tryptic soy broth) and M9 (minimum media). The comparisons of data and model results from the column experiments show that the depths associated with Cr(III) oxidation and transport within sediments of shallow aquatic systems can strongly influence trends in surface water quality. The results of this study suggests that carefully performed, laboratory column experiments is a useful tool in determining the biotransformation of redox-sensitive metals even in the presence of strong oxidant, like β-MnO2. ^

Examining Seasonally Pulsed Detrital Transport in the Coastal Everglades Using a Sediment Tracing Technique

The Comprehensive Everglades Restoration includes plans to restore freshwater delivery to Taylor Slough, a shallow drainage basin in the Southern Everglades, ultimately resulting in increased freshwater flow to the downstream Taylor River estuary. The effect of altered hydrologic regime on the transport dynamics of flocculent, estuarine detritus is not well understood. We utilized a paramagnetic sediment tracer to examine detrital transport in three Taylor River pond/creek pairs during early wet versus late wet transition season estuarine flow conditions. Flux of floc tracer was greatest in the downstream direction during all observations, and was most pronounced during the early wet season, coincident with shallower water depth and faster discharge from northern Taylor River. Floc tracer was more likely to move upriver during the late wet/dry season. We observed a floc tracer transport velocity of approximately 1.74 to 1.78 m/day across both seasonal hydrologic conditions. Tracer dynamics were also surprisingly site-dependent, which may highlight the importance of channel geomorphology in regulating hydrologic and sediment transport conditions. Our data suggest that restoration of surface water delivery to Taylor River will influence downstream loading of detritus material into riverine ponds. These detrital inputs have the potential to enhance ecosystem primary productivity and/or secondary productivity.

Simulating everglades national park hydrology and phosphorus transport under existing and future scenarios using numerical modeling

The Florida Everglades has a long history of anthropogenic changes which have impacted the quantity and quality of water entering the system. Since the construction of Tamiami Trail in the 1920's, overland flow to the Florida Everglades has decreased significantly, impacting ecosystems from the wetlands to the estuary. The MIKE Marsh Model of Everglades National Park (M3ENP) is a numerical model, which simulates Everglades National Park (ENP) hydrology using MIKE SHE/MIKE 11software. This model has been developed to determine the parameters that effect Everglades hydrology and understand the impact of specific flow changes on the hydrology of the system. ^ As part of the effort to return flows to the historical levels, several changes to the existing water management infrastructure have been implemented or are in the design phase. Bridge construction scenarios were programed into the M3ENP model to review the effect of these structural changes and evaluate the potential impacts on water levels and hydroperiods in the receiving Northeast Shark Slough ecosystem. These scenarios have shown critical water level increases in an area which has been in decline due to low water levels. Results from this work may help guide future decisions for restoration designs. ^ Excess phosphorus entering Everglades National Park in South Florida may promote the growth of more phosphorus-opportunistic species and alter the food chain from the bottom up. Two phosphorus transport methods were developed into the M3ENP hydrodynamic model to determine the factors affecting phosphorus transport and the impact of bridge construction on water quality. Results showed that while phosphorus concentrations in surface waters decreased overall, some areas within ENP interior may experience an increase in phosphorus loading which the addition of bridges to Tamiami Trail. Finally, phosphorus data and modeled water level data was used to evaluate the spectral response of Everglades vegetation to increasing phosphorus availability using Landsat imagery.^

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

A Modeling Approach to Determine the Impacts of Land Use and Climate Change Scenarios on the Water Flux of the Upper Mara River

With the flow of the Mara River becoming increasingly erratic especially in the upper reaches, attention has been directed to land use change as the major cause of this problem. The semi-distributed hydrological model Soil and Water Assessment Tool 5 (SWAT) and Landsat imagery were utilized in the upper Mara River Basin in order to 1) map existing field scale land use practices in order to determine their impact 2) determine the impacts of land use change on water flux; and 3) determine the impacts of rainfall (0%, ±10% and ±20%) and air temperature variations (0% and +5%) based on the Intergovernmental Panel on Climate Change projections on the water flux of the 10 upper Mara River. This study found that the different scenarios impacted on the water balance components differently. Land use changes resulted in a slightly more erratic discharge while rainfall and air temperature changes had a more predictable impact on the discharge and water balance components. These findings demonstrate that the model results 15 show the flow was more sensitive to the rainfall changes than land use changes. It was also shown that land use changes can reduce dry season flow which is the most important problem in the basin. The model shows also deforestation in the Mau Forest increased the peak flows which can also lead to high sediment loading in the Mara River. The effect of the land use and climate change scenarios on the sediment and 20 water quality of the river needs a thorough understanding of the sediment transport processes in addition to observed sediment and water quality data for validation of modeling results.

Physical oceanography, nutrients, and d18O measured on water bottle samples during POLARSTERN cruise ARK-XXII/2

Extremely low summer sea-ice coverage in the Arctic Ocean in 2007 allowed extensive sampling and a wide quasi-synoptic hydrographic and d18O dataset could be collected in the Eurasian Basin and the Makarov Basin up to the Alpha Ridge and the East Siberian continental margin. With the aim of determining the origin of freshwater in the halocline, fractions of river water and sea-ice meltwater in the upper 150 m were quantified by a combination of salinity and d18O in the Eurasian Basin. Two methods, applying the preformed phosphate concentration (PO*) and the nitrate-to-phosphate ratio (N/P), were compared to further differentiate the marine fraction into Atlantic and Pacific-derived contributions. While PO*-based assessments systematically underestimate the contribution of Pacific-derived waters, N/P-based calculations overestimate Pacific-derived waters within the Transpolar Drift due to denitrification in bottom sediments at the Laptev Sea continental margin. Within the Eurasian Basin a west to east oriented front between net melting and production of sea-ice is observed. Outside the Atlantic regime dominated by net sea-ice melting, a pronounced layer influenced by brines released during sea-ice formation is present at about 30 to 50 m water depth with a maximum over the Lomonosov Ridge. The geographically distinct definition of this maximum demonstrates the rapid release and transport of signals from the shelf regions in discrete pulses within the Transpolar Drift. The ratio of sea-ice derived brine influence and river water is roughly constant within each layer of the Arctic Ocean halocline. The correlation between brine influence and river water reveals two clusters that can be assigned to the two main mechanisms of sea-ice formation within the Arctic Ocean. Over the open ocean or in polynyas at the continental slope where relatively small amounts of river water are found, sea-ice formation results in a linear correlation between brine influence and river water at salinities of about 32 to 34. In coastal polynyas in the shallow regions of the Laptev Sea and southern Kara Sea, sea-ice formation transports river water into the shelf's bottom layer due to the close proximity to the river mouths. This process therefore results in waters that form a second linear correlation between brine influence and river water at salinities of about 30 to 32. Our study indicates which layers of the Arctic Ocean halocline are primarily influenced by sea-ice formation in coastal polynyas and which layers are primarily influenced by sea-ice formation over the open ocean. Accordingly we use the ratio of sea-ice derived brine influence and river water to link the maximum in brine influence within the Transpolar Drift with a pulse of shelf waters from the Laptev Sea that was likely released in summer 2005.

Polonium-210 and Lead-210 activities measured on 17 water bottle profiles and 50 surface water samples during POLARSTERN cruise ARK-XXII/2

Polonium-210 and Lead-210 have been measured in the water column and on suspended particulate matter during the POLARSTERN cruise ARK-XXII/2. The data have been submitted to Pangaea following a Polonium-Lead intercalibration exercise organized by GEOTRACES, where the AWI lab results range within the data standard deviation from 10 participating labs. Polonium-210 and Lead-210 in the ocean can be used to identify the sources and sinks of suspended matter. In seawater, Polonium-210 (210Po) and Lead-210 (210Pb) are produced by stepwise radioactive decay of Uranium-238. 210Po (138 days half life) and 210Pb (22.3 years half life) have high affinities for suspended particles. Those radionuclides are present in dissolved form and adsorbed onto particles. Following adsorption onto particle surfaces, 210Po especially is transported into the interior of cells where it bonds to proteins. In this way, 210Po also accumulates in the food chain. 210Po is therefore considered to be a good tracer for POC, and traces particle export over a timescale of month. 210Pb (22.3 years half life) adsorbs preferably onto structural components of cells, biogenic silica and lithogenic particles, and is therefore a better tracer more rapidly sinking matter. Our goal during ARK XXII/2 was to trace pathways of particulate and dissolved matter leaving the Siberian Shelf. The pathways of particulate and dissolved matter will be followed by the combined use of 210Po and 234Th as a tracer pair (and perhaps 210Pb) for particle flux (Cai, P.; Rutgers van der Loeff, MM (2008) doi:10.1594/PANGAEA.708354). This information gathered from the water column will be complemented with the results of the 210Po-210Pb study in sea ice (Camara-Mor, P, Instituto de Ciencias del Mar-SCIC, Barcelona, Spain) to provide a more thorough picture of particle transport from the shelf to the open sea and from surface to depth.

Sediment and pore water data of five cores from the continental slope off Uruguay

Assessing frequency and extent of mass movement at continental margins is crucial to evaluate risks for offshore constructions and coastal areas. A multidisciplinary approach including geophysical, sedimentological, geotechnical, and geochemical methods was applied to investigate multistage mass transport deposits (MTDs) off Uruguay, on top of which no surficial hemipelagic drape was detected based on echosounder data. Nonsteady state pore water conditions are evidenced by a distinct gradient change in the sulfate (SO4**2-) profile at 2.8 m depth. A sharp sedimentological contact at 2.43 m coincides with an abrupt downward increase in shear strength from approx. 10 to >20 kPa. This boundary is interpreted as a paleosurface (and top of an older MTD) that has recently been covered by a sediment package during a younger landslide event. This youngest MTD supposedly originated from an upslope position and carried its initial pore water signature downward. The kink in the SO4**2- profile approx. 35 cm below the sedimentological and geotechnical contact indicates that bioirrigation affected the paleosurface before deposition of the youngest MTD. Based on modeling of the diffusive re-equilibration of SO4**2- the age of the most recent MTD is estimated to be <30 years. The mass movement was possibly related to an earthquake in 1988 (approx. 70 km southwest of the core location). Probabilistic slope stability back analysis of general landslide structures in the study area reveals that slope failure initiation requires additional ground accelerations. Therefore, we consider the earthquake as a reasonable trigger if additional weakening processes (e.g., erosion by previous retrogressive failure events or excess pore pressures) preconditioned the slope for failure. Our study reveals the necessity of multidisciplinary approaches to accurately recognize and date recent slope failures in complex settings such as the investigated area.

Sulfate and methane concentrations measured in pore water of sediment core GeoB13809-1 and GeoB13849-1

Several previous studies have shown that submarine mass-movements can profoundly impact the shape of pore water profiles. Therefore, pore water geochemistry and diffusion models were proposed as tools for identifying and dating recent (max. several thousands of years old) mass-transport deposits (MTDs). In particular, sulfate profiles evidentially indicate transient pore water conditions generated by submarine landslides. After mass-movements that result in the deposition of sediment packages with distinct pore water signatures, the sulfate profiles can be kink-shaped and evolve into the concave and linear shape with time due to molecular diffusion. Here we present data from the RV METEOR cruise M78/3 along the continental margin off Uruguay and Argentina. Sulfate profiles of 15 gravity cores are compared with the respective acoustic facies recorded by a sediment echosounder system. Our results show that in this very dynamic depositional setting, non-steady state profiles occur often, but are not exclusively associated with mass-movements. Three sites that show acoustic indications for recent MTDs are presented in detail. Where recent MTDs are identified, a geochemical transport/reaction model is used to estimate the time that has elapsed since the perturbation of the pore water system and, thus, the timing of the MTD emplacement. We conclude that geochemical analyses are a powerful complementary tool in the identification of recent MTDs and provide a simple and accurate way of dating such deposits.