Regional-scale, fully coupled modelling of stream-aquifer interaction in a tropical catchment

The planning and management of water resources in the Pioneer Valley, north-eastern Australia requires a tool for assessing the impact of groundwater and stream abstractions on water supply reliabilities and environmental flows in Sandy Creek (the main surface water system studied). Consequently, a fully coupled stream-aquifer model has been constructed using the code MODHMS, calibrated to near-stream observations of watertable behaviour and multiple components of gauged stream flow. This model has been tested using other methods of estimation, including stream depletion analysis and radon isotope tracer sampling. The coarseness of spatial discretisation, which is required for practical reasons of computational efficiency, limits the model's capacity to simulate small-scale processes (e.g., near-stream groundwater pumping, bank storage effects), and alternative approaches are required to complement the model's range of applicability. Model predictions of groundwater influx to Sandy Creek are compared with baseflow estimates from three different hydrograph separation techniques, which were found to be unable to reflect the dynamics of Sandy Creek stream-aquifer interactions. The model was also used to infer changes in the water balance of the system caused by historical land use change. This led to constraints on the recharge distribution which can be implemented to improve model calibration performance. (c) 2006 Elsevier B.V. All rights reserved.

Swash-aquifer interaction in the vicinity of the water table exit point on a sandy beach

The coupling of sandy beach aquifers with the swash zone in the vicinity of the water table exit point is investigated through simultaneous measurements of the instantaneous shoreline (swash front) location, pore pressures and the water table exit point. The field observations reveal new insights into swash-aquifer coupling not previously gleaned from measurements of pore pressure only. In particular, for the case where the exit point is seaward of the observation point, the pore pressure response is correlated with the distance between the exit point and the shoreline in that when the distance is large the rate of pressure drop is fast and when the distance is small the rate decreases. The observations expose limitations in a simple model describing exit point dynamics which is based only on the force balance on a particle of water at the sand surface and neglects subsurface pressures. A new modified form of the model is shown to significantly improve the model-data comparison through a parameterization of the effects of capillarity into the aquifer storage coefficient. The model enables sufficiently accurate predictions of the exit point to determine when the swash uprush propagates over a saturated or a partially saturated sand surface, potentially an important factor in the morphological evolution of the beach face. Observations of the shoreward propagation of the swash-induced pore pressure waves ahead of the runup limit shows that the magnitude of the pressure fluctuation decays exponentially and that there is a linear increase in time lags, behavior similar to that of tidally induced water table waves. The location of the exit point and the intermittency of wave runup events is also shown to be significant in terms of the shore-normal energy distribution. Seaward of the mean exit point location, peak energies are small because of the saturated sand surface within the seepage face acting as a "rigid lid'' and limiting pressure fluctuations. Landward of the mean exit point the peak energies grow before decreasing landward of the maximum shoreline position.

Diagenesis and aquifer characteristics of the Sherwood Sandstone Group in North Nottinghamshire

The Sherwood Sandstone Group forms an important aquifer in Eastern England, which in North Nottinghamshire comprises the Nottingham Castle and Lenton Sandstone Formations. The aquifer is formed by an alluvial red-bed sequence dominated by medium-coarse grained sandstones which are texturally immature to submature and have only been subjected to shallow burial diagenesis. These sandstones reached the mature stage of the meso diagenetic regime, and four stages are recognized in their diagenetic history depending upon the physical/chemical processes prevailing and the subsequent effect on porosity and permeability. Stage "One" represents changes including dissolution of unstable silicates, clay replacement, red colouration and precipitation of authigenic minerals (quartz, feldspar, illite, l/S, kaolinite, dolomite, ferroan calcite, calcite). The net result of these changes was porosity reduction. Stage "Two" included changes due to mechanical compaction which resulted in minor porosity reduction. Stage "Three" was the main phase of secondary porosity enhancement. Stage "Four" represents changes taking place in the present groundwater where porosity and permeability may have been increased by dissolution and partly reduced by kaolinite precipitation. Porosity measured by water-resaturation and Hg-injection gave average values of 25.63% and 24.85% respectively. The results are comparable and showed marked correlation especially in highly porous/permeable rocks. Porosity measurements from photomicrographs were markedly offset from laboratory results. Horizontal Kw ranged between 1.43 x 10-5 and 1.13 x 10-1 mm/sec, with an average of  1.68 x 10-2 mm/sec. The estimated KHg ranged between 7.29 x 10-6 and 6.99 x 10-2 mm/sec with an average of 1.47 x 10-2 mm/sec. Both results are significantly correlated for highly porous/permeable rocks. The hydraulic properties are highly dependent upon the diagenetic properties (as most of the pores present are of secondary origin) as well as the pore size distribution. The chemistry of these groundwaters indicates that they are under-saturated with respect to dolomite, calcite, K-feldspar, l/S clay, and montmorillonite. The precipitation of kaolinite,and to a lesser extent illite, is favoured in the present groundwater regime.

Quantitative diatom-based reconstruction of paleoenvironmental conditions in Florida Bay and Biscayne Bay, U.S.A.

The spatial and temporal distribution of modern diatom assemblages in surface sediments, on the most dominant macrophytes, and in the water column at 96 locations in Florida Bay, Biscayne Bay and adjacent regions were examined in order to develop paleoenvironmental prediction models for this region. Analyses of these distributions revealed distinct temporal and spatial differences in assemblages among the locations. The differences among diatom assemblages living on subaquatic vegetation and sediments, and in the water column were significant. Because concentrations of salts, total phosphorus (WTP), total nitrogen (WTN) and total organic carbon (WTOC) are partly controlled by water management in this region, diatom-based models were produced to assess these variables. Discriminant function analyses showed that diatoms can also be successfully used to reconstruct changes in the abundance of diatom assemblages typical for different habitats and life habits. ^ To interpret paleoenvironmental changes, changes in salinity, WTN, WTP and WTOC were inferred from diatoms preserved in sediment cores collected along environmental gradients in Florida Bay (4 cores) and from nearshore and offshore locations in Biscayne Bay (3 cores). The reconstructions showed that water quality conditions in these estuaries have been fluctuating for thousands of years due to natural processes and sea-level changes, but almost synchronized shifts in diatom assemblages occurred in the mid-1960’s at all coring locations (except Ninemile Bank and Bob Allen Bank in Florida Bay). These alterations correspond to the major construction of numerous water management structures on the mainland. Additionally, all the coring sites (except Card Sound Bank, Biscayne Bay and Trout Cove, Florida Bay) showed decreasing salinity and fluctuations in nutrient levels in the last two decades that correspond to increased rainfall in the 1990’s and increased freshwater discharge to the bays, a result of increased freshwater deliveries to the Everglades by South Florida Water Management District in the 1980’s and 1990’s. Reconstructions of the abundance of diatom assemblages typical for different habitats and life habits revealed multiple sources of diatoms to the coring locations and that epiphytic assemblages in both bays increased in abundance since the early 1990’s. ^

Spatio-temporal patterns and nutrient status of macroalgae in a heavily managed region of Biscayne Bay, Florida, USA

The coastal bays of South Florida are located downstream of the Florida Everglades, where a comprehensive restoration plan will strongly impact the hydrology of the region. Submerged aquatic vegetation communities are common components of benthic habitats of Biscayne Bay, and will be directly affected by changes in water quality. This study explores community structure, spatio-temporal dynamics, and tissue nutrient content of macroalgae to detect and describe relationships with water quality. The macroalgal community responded to strong variability in salinity; three distinctive macroalgal assemblages were correlated with salinity as follows: (1) low-salinity, dominated by Chara hornemannii and a mix of filamentous algae; (2) brackish, dominated by Penicillus capitatus, Batophora oerstedii, and Acetabularia schenckii; and (3) marine, dominated by Halimeda incrassata and Anadyomene stellata. Tissue-nutrient content was variable in space and time but tissues at all sites had high nitrogen and N:P values, demonstrating high nitrogen availability and phosphorus limitation in this region. This study clearly shows that distinct macroalgal assemblages are related to specific water quality conditions, and that macroalgal assemblages can be used as community-level indicators within an adaptive management framework to evaluate performance and restoration impacts in Biscayne Bay and other regions where both freshwater and nutrient inputs are modified by water management decisions.

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.

Determining Spatial and Temporal Inputs of Freshwater, Including Submarine Groundwater Discharge, to a Subtropical Estuary Using Geochemical Tracers, Biscayne Bay, South Florida

Geochemical mixing models were used to decipher the dominant source of freshwater (rainfall, canal discharge, or groundwater discharge) to Biscayne Bay, an estuary in south Florida. Discrete samples of precipitation, canal water, groundwater, and bay surface water were collected monthly for 2 years and analyzed for salinity, stable isotopes of oxygen and hydrogen, and Sr2+/Ca2+ concentrations. These geochemical tracers were used in three separate mixing models and then combined to trace the magnitude and timing of the freshwater inputs to the estuary. Fresh groundwater had an isotopic signature (δ 18O = −2.66‰, δD −7.60‰) similar to rainfall (δ 18O = −2.86‰, δD = −4.78‰). Canal water had a heavy isotopic signature (δ 18O = −0.46‰, δD  = −2.48‰) due to evaporation. This made it possible to use stable isotopes of oxygen and hydrogen to separate canal water from precipitation and groundwater as a source of freshwater into the bay. A second model using Sr2+/Ca2+ ratios was developed to discern fresh groundwater inputs from precipitation inputs. Groundwater had a Sr2+/Ca2+ ratio of 0.07, while precipitation had a dissimilar ratio of 0.89. When combined, these models showed a freshwater input ratio of canal/precipitation/groundwater of 37%:53%:10% in the wet season and 40%:55%:5% in the dry season with an error of ±25%. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for 1–2% of the total fresh and saline water input.

Determination of vertical and horizontal pathways of injected fresh wastewater into a deep saline aquifer (Florida, USA) using natural chemical tracers

Two deep-well injection sites in south Florida, USA, inject an average of 430 million liters per day (MLD) of treated domestic fresh wastewater into a deep saline aquifer 900 m below land surface. Elevated levels of NH3 (highest concentration 939 µmol) in the overlying aquifer above ambient concentrations (concentration less than 30 µmol) were evidence of the upward migration of injected fluids. Three pathways were distinguished based on ammonium, chloride and bromide ratios, and temperature. At the South District Wastewater Treatment Plant, the tracer ratios showed that the injectate remained chemically distinct as it migrated upwards through rapid vertical pathways via density-driven buoyancy. The warmer injectate (mean 28°C) retained the temperature signal as it vertically migrated upwards; however, the temperature signal did not persist as the injectate moved horizontally into the overlying aquifers. Once introduced, the injectate moved slowly horizontally through the aquifer and mixed with ambient water. At the North District Wastewater Treatment Plant, data provide strong evidence of a one-time pulse of injectate into the overlying aquifers due to improper well construction. No evidence of rapid vertical pathways was observed at the North District Wastewater Treatment Plant.

Impact of Late Holocene climate variability and anthropogenic activities on Biscayne Bay (Florida, U.S.A.): Evidence from diatoms

Shallow marine ecosystems are experiencing significant environmental alterations as a result of changing climate and increasing human activities along coasts. Intensive urbanization of the southeast Florida coast and intensification of climate change over the last few centuries changed the character of coastal ecosystems in the semi-enclosed Biscayne Bay, Florida. In order to develop management policies for the Bay, it is vital to obtain reliable scientific evidence of past ecological conditions. The long-term records of subfossil diatoms obtained from No Name Bank and Featherbed Bank in the Central Biscayne Bay, and from the Card Sound Bank in the neighboring Card Sound, were used to study the magnitude of the environmental change caused by climate variability and water management over the last ~ 600 yr. Analyses of these records revealed that the major shifts in the diatom assemblage structures at No Name Bank occurred in 1956, at Featherbed Bank in 1966, and at Card Sound Bank in 1957. Smaller magnitude shifts were also recorded at Featherbed Bank in 1893, 1942, 1974 and 1983. Most of these changes coincided with severe drought periods that developed during the cold phases of El Niño Southern Oscillation (ENSO), Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO), or when AMO was in warm phase and PDO was in the cold phase. Only the 1983 change coincided with an unusually wet period that developed during the warm phases of ENSO and PDO. Quantitative reconstructions of salinity using the weighted averaging partial least squares (WA-PLS) diatom-based salinity model revealed a gradual increase in salinity at the three coring locations over the last ~ 600 yr, which was primarily caused by continuously rising sea level and in the last several decades also by the reduction of the amount of freshwater inflow from the mainland. Concentration of sediment total nitrogen (TN), total phosphorus (TP) and total organic carbon (TOC) increased in the second half of the 20th century, which coincided with the construction of canals, landfills, marinas and water treatment plants along the western margin of Biscayne Bay. Increased magnitude and rate of the diatom assemblage restructuring in the mid- and late-1900s, suggest that large environmental changes are occurring more rapidly now than in the past.

Ecology and distribution of diatoms in Biscayne Bay, Florida (USA): Implications for bioassessment and paleoenvironmental studies

The spatial and temporal distribution of planktonic, sediment-associated and epiphytic diatoms among 58 sites in Biscayne Bay, Florida was examined in order to identify diatom taxa indicative of different salinity and water quality conditions, geographic locations and habitat types. Assessments were made in contrasting wet and dry seasons in order to develop robust assessment models for salinity and water quality for this region. We found that diatom assemblages differed between nearshore and offshore locations, especially during the wet season when salinity and nutrient gradients were steepest. In the dry season, habitat structure was primary determinant of diatom assemblage composition. Among a suite of physicochemical variables, water depth and sediment total phosphorus (STP) were most strongly associated with diatom assemblage composition in the dry season, while salinity and water total phosphorus (TP) were more important in the wet season. We used indicator species analysis (ISA) to identify taxa that were most abundant and frequent at nearshore and offshore locations, in planktonic, epiphytic and benthic habitats and in contrasting salinity and water quality regimes. Because surface water concentrations of salts, total phosphorus, nitrogen (TN) and organic carbon (TOC) are partly controlled by water management in this region, diatom-based models were produced to infer these variables in modern and retrospective assessments of management-driven changes. Weighted averaging (WA) and weighted averaging partial least squares (WA-PLS) regressions produced reliable estimates of salinity, TP, TN and TOC from diatoms (r2 = 0.92, 0.77, 0.77 and 0.71, respectively). Because of their sensitivity to salinity, nutrient and TOC concentrations diatom assemblages should be useful in developing protective nutrient criteria for estuaries and coastal waters of Florida.