Foraminiferal assemblages in Biscayne Bay,Florida, USA: Responses to urban and agricultural influence in a subtropical estuary

This study assessed foraminiferal assemblages in Biscayne Bay, Florida, a heavily utilized estuary, interpreting changes over the past 65 years and providing a baseline for future comparisons. Analyses of foraminiferal data at the genus level revealed three distinct biotopes. The assemblage from the northern bay was characterized by stress-tolerant taxa, especially Ammonia, present in low abundances (∼2.0 × 10³ foraminifers/gram) though relatively high diversity (∼19 genera/sample). The southwestern margin of the bay was dominated by Ammonia and Quinqueloculina, an assemblage characterized by the lowest diversities (∼12 genera/sample) and highest abundances (∼1.1 × 10⁴ foraminifers/gram), influenced by both reduced salinity and elevated organic-carbon concentrations. A diverse assemblage of smaller miliolids and rotaliids (∼26 genera/sample) characterized the open-bay assemblage, which also had a significant component (∼10%) of taxa that host algal endosymbionts. In the past 65 years, populations of symbiont-bearing taxa, which are indicators of normal-marine conditions, have decreased while stress-tolerant taxa, especially Ammonia spp., have increased in predominance.     

A nutrient loading budget for Biscayne Bay, Florida

The water quality in Biscayne Bay has been significantly affected by past and continuing coastal and watershed development. The nutrient concentrations in the Bay have been dramatically changed by the conversion of natural creeks and sheet flow freshwater inputs to rapid and episodic canal inputs from the large and rapidly expanding Miami metropolitan area. This study is an evaluation of nutrient loadings to Biscayne Bay for 1994-2002 from canal, atmospheric, and groundwater sources. Dissolved inorganic nitrogen (DIN, as nitrate, nitrite, and ammonium) and total phosphorus (TP) loadings by the canals were influenced by their geographic locations relative to discharge amount, watershed land use, stormwater runoff, and proximity to landfills. Annual budgets showed that canals contributed the bulk of N loading to the bay as 1687.2 metric ton N yr(-1) (88% total load). Direct atmospheric DIN load for Biscayne Bay was only 231.7 ton N yr(-1), based on surface area. Of the canal DIN load, nitrate+nitrite (NO(x)(-)) loading (1294.5 ton N yr(-1)) made up a much greater proportion than that of ammonium (NH(4)(+), 392.6 ton N yr(-1)). In the urbanized north and central Bay, canal DIN load was evenly split between NO(x)(-) and NH(4)(+). However, in the south, 95% of the DIN load was in the form of NO(x)(-), reflecting the more agricultural land use. Contrary to N, canals contributed the only 66% of P load to the bay (27.5 ton P yr(-1)). Atmospheric TP load was 14 ton Pyr(-1). In the north, canal P load dominated the budget while in the south, atmospheric load was almost double canal load. Groundwater inputs, estimated only for the south Bay, represented an important source of N and P in this zone. Groundwater input of N (141 ton N yr(-1)) was about equal to atmospheric load, while P load (5.9 ton P yr(-1)) was about equal to canal load.     

Spatial patterning of water quality in Biscayne Bay, Florida as a function of land use and water management

An objective classification analysis was performed on a water quality data set from 25 sites collected monthly during 1994-2003. The water quality parameters measured included: TN, TON, DIN, NH4+, NO3-, NO2-, TP, SRP, TN:TP ratio, TOC, DO, CHL A, turbidity, salinity and temperature. Based on this spatial analysis, Biscayne Bay was divided into five zones having similar water quality characteristics. A robust nutrient gradient, driven mostly by dissolved inorganic nitrogen, from alongshore to offshore in the main Bay, was a large determinant in the spatial clustering. Two of these zones (Alongshore and Inshore) were heavily influenced by freshwater input from four canals which drain the South Dade agricultural area, Black Point Landfill, and sewage treatment plant. The North Bay zone, with high turbidity, phytoplankton biomass, total phosphorus, and low DO, was affected by runoff from five canals, the Munisport Landfill, and the urban landscape. The South Bay zone, an embayment surrounded by mangrove wetlands with little urban development, was high in dissolved organic constituents but low in inorganic nutrients. The Main Bay was the area most influenced by water exchange with the Atlantic Ocean and showed the lowest nutrient concentrations. The water quality in Biscayne Bay is therefore highly dependent of the land use and influence from the watershed.     

Evaluation of time-space distributions of submarine ground water discharge

Submarine ground water discharge (SGD) rates were measured continuously by automated seepage meters to evaluate the process of ground water discharge to the ocean in the coastal zone of Suruga Bay, Japan. The ratio of terrestrial fresh SGD to total SGD was estimated to be at most 9% by continuous measurements of electrical conductivity of SGD. Semidiurnal changes of SGD due to tidal effects and an inverse relation between SGD and barometric pressure were observed. Power spectrum density analyses of SGD, sea level, and ground water level show that SGD near shore correlated to ground water level changes and SGD offshore correlated to sea level changes. SGD rates near the mouth of the Abe River are smaller than those elsewhere, possibly showing the effect of the river on SGD. The ratio of terrestrial ground water discharge to the total discharge to the ocean was estimated to be 14.7% using a water balance method.     

Occurrence of IRGAROL 1051 in coastal waters from Biscayne Bay,Florida, USA

Surface water samples from marinas, commercial ports and open bay areas collected from Biscayne Bay and the Miami River, Florida, USA, were analyzed for the occurrence of IRGAROL 1051 by GC/MS. The anifouling boosting herbicide was found in 80% (46/57) of the samples collected between March 1999 and September 2000. Concentrations within the bay range between non-detected (<1 ppt) and 61 ppt (ng/L) and were generally low compared with levels reported in European or Japanese waters. Aside from the elevated concentrations observed along the Miami River South Fork (61 ppt), the highest concentrations observed in the bay corresponded to marinas with high density of pleasure craft and restricted water circulation. In contrast, occurrence of IRGAROL 1051 along the commercial port or the cruise line terminal was generally lower (<1-2.2 ppt). Concentrations around Coconut Grove Marina were consistently higher (5-12 ppt) than the rest of the bay waters during the whole period of time surveyed.     

Stream bottom resistivity tomography to map ground water discharge

This study investigates the effectiveness of direct current electrical resistivity as a tool for assessing ground water/surface water interactions within streams. This research has shown that patterns of ground water discharge can be mapped at the meter scale, which is important for understanding stream water quality and ecosystem function. Underwater electrical resistivity surveys along a 107-m stream section within the Burd Run Watershed in South Central Pennsylvania identified three resistivity layers: a resistive (100 to 400 Ωm) surface layer corresponding to the streambed sediments, a conductive (20 to 100 Ωm) middle layer corresponding to residual clay sediments, and a resistive (100 to 450 Ωm) bottom layer corresponding to the carbonate bedrock. Tile probing to determine the depth to the bedrock and resistivity test box analysis of augered sediment samples confirmed these interpretations of the resistivity data. Ground water seeps occurred where the resistivity data showed that the residual clays were thinnest and bedrock was closest to the streambed. Plotting the difference in resistivity between two surveys, one conducted during low-stage and the other during high-stage stream conditions, showed changes in the conductivity of the pore fluids saturating the sediments. Under high-stream stage conditions, the top layer showed increased resistivity values for sections with surface water infiltration but showed nearly constant resistivity in sections with ground water seeps. This was expressed as difference values less than 50 Ωm in the area of the seeps and greater than 50 Ωm change for the streambed sediments saturated by surface water. Thus, electrical resistivity aided in characterizing ground water discharge zones by detecting variations in subsurface resistivity under high- and low-stream stage conditions as well as mapping subsurface heterogeneities that promote these exchanges.     

Geological controls on submarine groundwater discharge in Long Bay,South Carolina (USA)

A combination of geophysical methods including continuous electrical resistivity and high-resolution Chirp sub-bottom profiling were utilized to characterize geologic controls on pore fluid salinity in the nearshore of Long Bay, SC. Resistivity values ranged from less than 1 Ω m to greater than 40 Ω m throughout the bay. Areas of elevated electrical resistivity suggest the influence of relatively fresher water on pore water composition. Geophysical evidence alone does not eliminate all ambiguity associated with lithological and porosity variations that may also contribute to electrical structure of shallow marine sediments. The anomalous field is of sufficient magnitude that lithological variation alone does not control the spatial distribution of elevated electrical resistivity zones. Geographical distribution of electrical anomalies and structures interpreted from nearby sub-bottom profiles indicates abrupt changes in shallow geologic units control preferential pathways for discharge of fresh water into the marine environment. Shore parallel resistivity profiles show dramatic decreases in magnitude with increasing distance from shore, suggesting a significant portion of the terrestrially driven fresh SGD in Long Bay is occurring via the surficial aquifer within a few hundred meters of shore.     

Michigan basin regional ground water flow discharge to three Great Lakes

Ground water discharge to the Great Lakes around the Lower Peninsula of Michigan is primarily from recharge in riparian basins and proximal upland areas that are especially important to the northern half of the Lake Michigan shoreline. A steady-state finite-difference model was developed to simulate ground water flow in four regional aquifers in Michigan's Lower Peninsula: the Glaciofluvial, Saginaw, Parma-Bayport, and Marshall aquifers interlayered with the Till/"red beds," Saginaw, and Michigan confining units, respectively. The model domain was laterally bound by a continuous specified-head boundary, formed from lakes Michigan, Huron, St. Clair, and Erie, with the St. Clair and Detroit River connecting channels. The model was developed to quantify regional ground water flow in the aquifer systems using independently determined recharge estimates. According to the flow model, local stream stages and discharges account for 95% of the overall model water budget; only 50% enters the lakes directly from the ground water system. Direct ground water discharge to the Great Lakes' shorelines was calculated at 36 m3/sec, accounting for 5% of the overall model water budget. Lowland areas contribute far less ground water discharge to the Great Lakes than upland areas. The model indicates that Saginaw Bay receives only approximately 1.13 m3/sec ground water; the southern half of the Lake Michigan shoreline receives only approximately 2.83 m3/sec. In contrast, the northern half of the Lake Michigan shoreline receives more than 17 m3/sec from upland areas.     

PRO-GRADE: GIS toolkits for ground water recharge and discharge estimation

PRO-GRADE is an ESRI ArcGIS 9.2 plug-in package that consists of two separate toolkits: (1) the p attern r ecognition o rganizer for g eographic i nformation s ystem (PRO-GIS) and (2) the g round water r echarge a nd d ischarge e stimator for GIS (GRADE-GIS). PRO-GIS is a collection of several existing image-processing algorithms into one user interface to offer the flexibility to extract spatial patterns according to the user's needs. GRADE-GIS is a ground water recharge and discharge estimation interface using a mass balance method that requires only hydraulic conductivity, water table, and bedrock elevation data for simulating two-dimensional steady-state unconfined aquifers. PRO-GRADE was developed to assist ongoing assessments of the water resources in Illinois and Wisconsin, and is being used to assist several ground water resource studies in several locations in the United States. The advantage of using PRO-GRADE is to enable fast production of initial recharge and discharge maps that can be further enhanced by using a follow-up ground water flow model with parameter estimation codes. PRO-GRADE leverages ArcGIS to provide a computer-assisted framework to support expert judgment in order to efficiently select alternative recharge and discharge maps that can be used as (1) guidelines for field study planning and decision making; (2) initial conditions for numerical simulation; and (3) screening for alternative model selection and prediction/parameter uncertainty evaluation. In addition, PRO-GRADE allows for more easy and rapid correlation of those maps with other hydrologically relevant geospatial data.     

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d⁻¹. ³H/³He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d⁻¹) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.     

Using Rn to estimate submarine groundwater discharge (SGD) and the associated nutrient fluxes into Xiangshan Bay,East China Sea

Continuous radon (²²²Rn) monitoring was conducted at two stations (site A and site B) with different perpendicular distance from the shoreline in Xiangshan Bay, East China Sea. Based on a ²²²Rn balance model (various sources and sinks of ²²²Rn in coastal water), the average rate of SGD was estimated to be 0.69 cm/day and 0.23 cm/day for site A and site B, respectively. The results from a nutrient analysis of the groundwater indicate that the associated nutrients fluxes loading through the SGD pathway were 4.27 × 10⁶ mol/day for DIN, 2.24 × 10⁴ mol/day for DIP and 1.82 × 10⁶ mol/day for DSi, respectively, which were comparable to or even higher than the levels observed in the local streams. Therefore, adequate attention should be paid to the importance of SGD as one source of nutrients during the eutrophication control process in this area.     

Ecosystem under pressure: ballast water discharge into Galveston Bay, Texas (USA) from 2005 to 2010

Ballast water exchange processes facilitate the dispersal and unnatural geographic expansion of phytoplankton, including harmful algal bloom species. From 2005 to 2010, over 45,000 vessels (≈ 8000 annually) travelled across Galveston Bay (Texas, USA) to the deep-water ports of Houston (10th largest in the world), Texas City and Galveston. These vessels (primarily tankers and bulkers) discharged ≈ 1.2 × 10(8) metrictons of ballast water; equivalent to ≈ 3.4% of the total volume of the Bay. Over half of the ballast water discharged had a coastwise origin, 96% being from US waters. Galveston Bay has fewer non-indigenous species but receives a higher volume of ballast water discharge, relative to the highly invaded Chesapeake and San Francisco Bays. Given the magnitude of shipping traffic, the role of Galveston Bay, both as a recipient and donor region of non-indigenous phytoplankton species is discussed here in terms of the invasibility risk to this system by way of ballast water.     

Evaluating the source and seasonality of submarine groundwater discharge using a radon-222 pore water transport model

Pore water radon (²²²Rn) distributions from Indian River Lagoon, Florida, are characterized by three zones: a lower zone where pore water ²²²Rn and sediment-bound radium (²²⁶Ra) are in equilibrium and concentration gradients are vertical; a middle zone where ²²²Rn is in excess of sediment-bound ²²⁶Ra and concentration gradients are concave-downward; and an upper zone where ²²²Rn concentration gradients are nearly vertical. These ²²²Rn data are simulated in a one-dimensional numerical model including advection, diffusion, and non-local exchange to estimate magnitudes of submarine groundwater discharge components (fresh or marine). The numerical model estimates three parameters, fresh groundwater seepage velocity, irrigation intensity, and irrigation attenuation, using two Monte Carlo (MC) simulations that (1) ensure the minimization algorithm converges on a global minimum of the merit function and the parameter estimates are consistent within this global minimum, and (2) provide 90% confidence intervals on the parameter estimates using the measured ²²²Rn activity variance. Model estimates of seepage velocities and discharge agree with previous estimates obtained from numerical groundwater flow models and seepage meter measurements and show the fresh water component decreases offshore and varies seasonally by a factor of nine or less. Comparison between the discharge estimates and precipitation patterns suggests a mean residence time in unsaturated and saturated zones on the order of 5 to 7 months. Irrigation rates generally decrease offshore for all sampling periods. The mean irrigation rate is approximately three times greater than the mean seepage velocity although the ranges of irrigation rates and seepage velocities are the same. Possible mechanisms for irrigation include density-driven convection, wave pumping, and bio-irrigation. Simulation of both advection and irrigation allows the separation of submarine groundwater discharge into fresh groundwater and (re)circulated lagoon water.     

Interactions of river discharge and tidal modulation in a tropical estuary, NE Brazil

Observations of thermohaline properties and currents were undertaken in the Curimataú River estuary (6°18′S), Rio Grande do Norte state (RN), Brazil, during consecutive neap–spring tidal cycles in the austral autumn rainy season. Highly asymmetric neap tide along channel velocities (−0.4 to 0.9 m s⁻¹) and highly stratified conditions were generated by an increase of the buoyancy energy from the freshwater input (R _iE≈5.6). During the spring-tidal cycle the river discharge decreased and the longitudinal velocity components were higher, less asymmetrical (−0.8 to 1.1 m s⁻¹) and semidiurnal, associated with moderately stratified conditions (R _iE≈0.1) due to the increase of the kinetic tidal energy forcing mechanism. The overall salinity variation from surface to bottom during two tidal cycles was from 20.5 to 36.3 and 29 to 36.7 in the neap and spring tide experiments, respectively; in the last experiment, the tropical water (TW) mass intrusion was enhanced. The net salt transport reversed from down to up estuary during the neap and spring tide experiments, respectively, varied from 6.0 to –2.0 kg m⁻¹ s⁻¹, an indication of changes in the main forcing of the estuary dynamics. Evaluation of a classical steady analytical model, in comparison with nearly steady experimental vertical profiles of velocity, shows an agreement classifiable as reasonably fair.     

Transport through Tavernier Creek: a primary flushing pathway for Northeast Florida Bay

Current meter data from a 411-day study are used to characterize the movement of water through Tavernier Creek, a tidal channel in the Upper Florida Keys that connects northeastern Florida Bay with the narrow continental shelf on the Atlantic Ocean side of the Keys. The record reveals active tidal and nontidal exchanges. Strongest flood and ebb current speeds commonly reach 50 cm s⁻¹. Low-frequency exchanges are highly coherent with the across-shelf component of local wind stress over time scales in excess of 2.5 days. Bay-shelf exchanges are investigated in four ways. Current measurements made while a drogue was tracked from one end of the creek to the other provide a relationship that can be applied to the time series of current meter data. Results suggest that ocean water reaches the bay end of the creek on 92% of the floods, and bay water reaches the ocean end on 94% of the ebbs. The Eulerian tidal excursion calculated from the amplitude of the M ₂ tidal constituent is 1.42 times the length of the channel, and half-tidal cycle Eulerian displacements are commonly 1.5 times the length of the channel. Salinity measurements over a 165-day period document the arrival of bay and ocean water at a study site at about the midpoint of Tavernier Creek. Results suggest that the creek becomes completely flushed after about 450,000 m³ of water have entered from either end. Histograms of ebb and flood volume transports indicate that half-tidal cycle transports are commonly between 800,000 and 1,100,000 m³. The long-term movement of water through Tavernier Creek is a net outflow from Florida Bay. Results support the idea that Tavernier Creek serves as an effective conduit for exchanging bay and ocean water, and especially for draining the northeast corner of Florida Bay.     

Arsenic speciation and distribution in the Carnon estuary following the acute discharge of contaminated water from a disused mine

In February 1991 Wheal Jane, a tin mine in S.W. England, was closed, the mine drainage pumps were removed and underground water levels were left to rise. By January 1992 the build up of water in the mine was sufficient to cause the Nangiles adit to burst, and contaminated mine water containing high levels of arsenic overflowed into the Carnon River. The river in turn drained into an estuary system which reaches the sea at Falmouth. A contingency plan was introduced, with water being pumped from the mine and treated before discharge through a tailings dam. This brought the discharge from Nangiles under control by late February 1992, although flow recommenced from mid-April until early July 1992.The immediate impact of the discharge and treatment operation on the distribution and speciation of arsenic in the river was the presence of high concentrations of dissolved arsenite, with a level of 240 μg As l⁻¹ being measured below the mine tailings dam. Elevated concentrations of arsenite, typically 9 μg As l⁻¹, persisted in the saline region of the Carnon estuary. By the following July the levels of arsenite had dropped significantly in the river but the levels in Restronguet Creek remained essentially unchanged. In comparison with the arsenite concentrations, the arsenate levels were comparatively low in the whole river/estuary system in the February following the initial discharge. A significant increase was however found in the saline region in the following July.Methylated arsenic species were absent from the river but were present in the biologically productive saline region during the warmer summer sampling period. The July concentrations of dimethylarsenic ranged from 1.0 to 2.1 μg As l⁻¹ and monomethylarsenic from 0.2 to 0.5 μg As l⁻¹.     

Modelling residence‐time response to freshwater input in Apalachicola Bay,Florida, USA

Residence time of an estuary can be used to estimate the rate of removal of freshwater and pollutants from river inflow. In this study, a calibrated three‐dimensional hydrodynamic model was used to determine residence time in response to the change of freshwater input in Apalachicola Bay. The bay is about 40 km long and 7 km wide, with an average 3 m water depth. Through hydrodynamic model simulations, the spatial and temporal salinity and the total freshwater volume in the bay were calculated. Then the freshwater fraction method was used to estimate the residence time. Results indicate that the residence time in Apalachicola Bay typically ranges between 3 and 10 days for the daily freshwater input ranging from 177 m³/s to 4561 m³/s. Regression analysis of model results shows that an exponential regression equation can be used to correlate the estuarine residence time to changes of freshwater input. Copyright © 2002 John Wiley & Sons, Ltd.     

Distribution and characteristics of suspended particulate matter in a heavily eutrophic estuary, Tokyo Bay, Japan

Bulk suspended particulate matter (SPM), chlorophyll a (Chl a), ignition loss, particulate organic carbon (POC), organic nitrogen, inorganic phosphorus (PIP), and organic phosphorus were investigated in the heavily eutrophic Arakawa River estuary, Japan. Chl a was high (approximately 35 microg l(-1)) in summer and low (approximately 6.7 microg l(-1)) in winter and autumn. POC from living phytoplankton accounted for approximately 34% and approximately 70% of total POC during low- and high-biomass seasons, respectively. During the low-biomass season, detrital POC distribution was conservative, and less reactive, land-derived materials mainly composed particulate organic materials (POM), but complex mixing of land-derived POM and autochthonous planktonic detritus caused nonconservative detrital POC behavior during the high-biomass season. PIP concentration in SPM decreased with increasing salinity, likely by desorption of soluble orthophosphate (ortho-P). The ortho-P released from SPM, 56% of the ortho-P input from the Arakawa River to the bay, was a significant potential source of biologically available phosphorus causing eutrophication of coastal environments.     

Investigation of submarine groundwater discharge to tidal rivers: Evidence for regional and local scale seepage

Fluxes of submarine groundwater discharge (SGD) were investigated into two tidal rivers on the north and south shore of Long Island, NY, during July 2015. Ground‐based handheld thermal infrared (TIR) imagery, combined with direct push‐point piezometer sampling, documented spatially heterogeneous small‐scale intertidal seepage zones. Pore waters were relatively fresh and enriched in nitrogen (N) within these small‐scale seeps. Pore waters sampled just 20 cm away, outside the boundary of the ground‐based TIR‐located seepage zone, were more saline and lower in N. These ground‐based TIR‐identified seeps geochemically represented the terrestrial fresh groundwater endmember, whereas N in pore waters sampled outside of the TIR‐identified seeps was derived from the remineralization of organic matter introduced into the sediment by tidal seawater infiltration. A ²²²Rn (radon‐222) time‐series was used to quantify fresh SGD‐associated N fluxes using the N endmembers sampled from the ground‐based TIR pore water profiles. N fluxes were up‐scaled to groundwater seepage zones identified from high‐resolution airborne TIR imagery using the two‐dimensional size of the airborne TIR surface water anomalies, relative to the N flux from the time‐series sampling location. Results suggest that the N load from the north‐shore tidal river to Long Island Sound is underrepresented by at least 1.6–3.6%, whereas the N load from SGD to a south‐shore tidal river may be up to 9% higher than previous estimates. These results demonstrate the importance of SGD in supplying nutrients to the lower reaches of tidal rivers and suggest that N loads in other tidal river environments may be underestimated if SGD is not accounted for.