Production and nutrient dynamics of aSyringodium filiforme Kütz. Seagrass bed in Indian River Lagoon,Florida

A year-long analysis of the characteristics of the seagrassSyringodium filiforme and the associated dynamics of the nutrient pool in the sediment pore water was done to assess co-variation. Changes in seagrass growth rate and standing stock throughout the year were accompanied by seasonal changes in the nutrient pools. The link between plant production and morphometrics and the sediment nutrient pool was found to be predominantly physiological, with the plant balancing the ability to photosynthesize with the nutrients needed for maintaining production. Measurements of whole plant growth for this seagrass, rather than the more typical leaf growth measurements, show that the production of new shoots and rhizome elongation for these plants represents as substantial amount of growth that usually goes unmeasured. Further, these whole plant growth measures demonstrate the rapid lateral rhizome spread of this species, exceeding one meter per plant per year. The primary cause of seasonal variation in the yearly seagrass cycle was investigated. Correlation analysis supported the hypothesis that the major factor controlling seasonal variation in this seagrass was light. During the peak growing season, however, growth was not regulated by light but by nitrogen. Depletion of the sediment ammonium pool and reduction in pore water ammonium relative to adsorbed ammonium, as well as changes in N content of seagrass leaves, support our hypothesis of peak growing season nitrogen limitation. Our results forSyringodium filiforme in terrigenous sediments are in contrast to our recent findings of phosphorus limitation in this same species occurring in carbonate sediments.     

Effects of two hurricanes onSyringodium filiforme, manatee grass, within the Loxahatchee River estuary, Southeast Florida

In September 2004, the Loxahatchee River Estuary was affected by Hurricanes Frances and Jeanne, which resulted in a monthly rainfall record of 610 mm and abnormally high freshwater discharges to the system. The occurrence, density, and biomass ofSyringodium filiforme in the Loxahatchee River Estuary declined significantly following the September 2004 storms based on 15 mo of pre-hurricane monitoring and 12 mo of post-hurricane monitoring. Throughout posthurricane monitoring,S. filiforme showed no sign of recovery, thoughHalophila johnsonii increased considerably during the post-hurricane period. Freshwater discharges resulting from the September 2004 hurricanes lowered minimum daily salinity values to near zero and increased standard deviation of daily salinity values to 11‰. Extremely low minimum daily salinity values and high daily salinity fluctuations likely resulted in the observed decline ofS. filiforme. We advise the use of minimum daily salinity values when assessing seagrass habitat suitability or when modeling the effects of alternative water management scenarios.     

Biodeposition by a fouling community in the Indian River,Florida

Biodeposition rates were studied for a fouling community with a biomass of 6–10 kg per m² dry wt including shells in which the barnacle Balanus eburneus was a dominant species. The fouling community filtered Indian River lagoon water containing 2–15 mg per 1 mud-size particles and deposited them as sand-size fecal pellets. Measurements of the fecal pellet flux by sediment traps indicated seasonal variations between 16.7 and 74.8 g per m² per day. A significant correlation was found between fecal pellet flux and temperature (r=0.90; p<0.001). The average flux of fecal pellet deposition was four times greater than the average flux of suspended particle settling without biological influence. Suspended sediment concentration did not significantly affect the rate of biodeposition. Annual biodeposition was 18 kg per m².     

A comparison of mercury in estuarine fish between Florida Bay and the Indian River Lagoon, Florida, USA

Concentrations of mercury (Hg) in fish were compared between two Florida estuaries, the Indian River Lagoon and Florida Bay. The objective was to determine if differences in Hg concentration exist and to attempt to relate those differences to sources of Hg. Five hundred and thirteen estuarine fish were collected and analyzed for Hg concentration. Fish species collected were black drum, bluefish, bonnethead shark, common snook, crevalle jack, gafftopsail catfish, gray snapper, Mayan cichlid, pompano, red drum, sheepshead, southern flounder, spadefish, and spotted seatrout. Analysis of variance of species-specific Hg data among the three defined regions of eastern and western Florida Bay and the Indian River Lagoon substantiated regional differences. Proximity to known anthropogenic sources of Hg appeared to be a significant factor in the distribution of Hg concentration among the fish collected. Sufficient numbers of crevalle jack, gray snapper, and spotted seatrout were collected to permit statistical analysis among regions. Hg concentrations in all three of these species from eastern Florida Bay were higher than those collected in the other two areas. A major fraction of the estuarine fish collected in eastern Florida Bay exceeded one or more State of Florida or U.S. Food and Drug Administration fish consumption health advisory criteria. In general, fish from western Florida Bay contained less Hg than those from the Indian River Lagoon, and fish from the Indian River contained less Hg than those from eastern Florida Bay. Crevalle jack from all areas and spotted seatrout from Florida Bay were placed on a consumption advisory in Florida. Detailed study of Florida Bay food web dynamics and Hg biogeochemical cycling is recommended to better understand the processes underlying the elevated Hg levels in fish from eastern Florida Bay. This information may be vital in the formulation of appropriate strategies in the ongoing South Florida restoration process.     

The demise and recovery of seagrass in the northern Indian River Lagoon,Florida

Seagrass both disappeared and recovered within 4 yr in one region of northern Indian River Lagoon (IRL). For the specific area referred to as Turnbull Bay, a relatively pristine area of the IRL, over 100 ha of seagrass completely disappeared from 1996 to 1997 and then recovered by 2000. Based on lagoon-wide mapping from aerial photographs taken every 2–3 years since 1986, coverage of seagrass in Turnbull Bay declined from 124 ha in 1989 to 34 ha by 1999 and increased to 58 ha in 2003. Bi-annual monitoring of fixed seagrass transects tells a more detailed story. Species composition along the Turnbull transect shifted fromHalodule wrightii toRuppia maritima beginning in 1995, and macroalgal abundance increased. By the summer of 1997, seagrass completely disappeared along the transect, as well as in most of the surrounding areas in Turnbull Bay; macroalgae covered much of the sediment surface. No significant water quality changes were detected. Light attenuation and suspended solid values did increase after the seagrass disappeared. Porewater sulfide concentrations, taken after all the grass was gone in 1997, were high (2,000 μM), but did improve by 1998 (1,200 μM). Seagrass recovery was rapid and occurred in the reverse sequence of species loss. Seedlings ofR. maritima were the first colonizers, then patches ofH. wrightii appeared. In 2000,Halophila engelmannii returned in the deeper water (>0.6m). By the summer of 2000, the beds had completely recovered. We conclude that this demise was a natural event caused by a long-term buildup of seagrass biomass and a thick (10–15 cm) layer of organic detritus and ooze. We surmise that such a crash and subsequent recovery may be a natural cycle of decline and recovery within this semirestricted, poorly-flushed area. The frequency of this cycle remains uncertain.     

Monitoring and Modeling of Syringodium filiforme (Manatee Grass) in Southern Indian River Lagoon

Decreased salinity and submarine light associated with hurricanes of 2004?C2005 impacted seagrass habitats in the Florida coastal zone. A combination of salinities ??20 and light attenuation ??1.5?m^?1 resulting from the freshwater discharge in 2005 were among the drivers for a widespread decrease in the coverage and biomass of Syringodium filiforme (manatee grass) in 2006. These observations provided an opportunity to develop and apply a modeling framework to simulate responses of S. filiforme to variable water quality. The framework connects water column variables to field monitoring of seagrass abundance and salinity growth response experiments. The base model was calibrated with macrophyte abundance observed in southern Indian River Lagoon (IRL) from 2002 to 2007 and tested against shoot data from a different time (1997?C2002) and nearby location in the IRL. Model shoot biomass (gC?m^?2) was similar to field observations (r ²?=?0.70) while responding to monthly seasonal fluctuations in salinity and light throughout the 6-year simulations. Field and model results indicated that S. filiforme growth and survival were sensitive to, and increased with, rising salinity throughout 2007. This modeling study emphasizes that discharge, salinity, and submarine light are inter-dependent variables affecting South Florida seagrass habitats on seasonal to inter-annual time scales.     

Responses of common fouling organisms in the Indian River,Florida, to various predation and disturbance intensities

The response of some common fouling organisms to increased predation by scraping and decreased predation by caging is described. Substrate coverage by many colonial forms such as Perophora viridis and Diplosoma macdonaldi was not affected by changes in predation intensity whereas coverage by many solitary forms such as Spirorbis sp. and Styela plicata increased when predation was reduced. These differences in responses of colonial and solitary species may be because solitary species reproduce sexually and resettle newly opened space whereas colonial species can rapidly expand adjacent colonies into newly cleared space.     

Setting seagrass depth, coverage, and light targets for the Indian River Lagoon system, Florida

Seagrass protection and restoration in Florida’s Indian River Lagoon system (IRLS) is a mutual goal of state and federal programs. These programs require, the establishment of management targets indicative of seagrass recovery and health. We used three metrics related to seagrass distribution: areal coverage, depth limit, and light requirement. In order to account for the IRLS’s spatial heterogeneity and temporal variability, we developed coverage and depth limit targets for each of its 19 segments. Our method consisted of two steps: mapping the union of seagrass coverages from all availabe mapping years (1943, 1986, 1989, 1992, 1994, 1996, and 1999) to delineate wherever seagrass had been mapped and determining the distribution of depth limits based on 5,615 depth measurements collected on or very near the deep-edge boundary of the union coverage. The frequency distribution of depth limits derived from the union coverage, along with the median (50th percentile) and maximum (95th percentile) depth limits, serve as the seagrass depth targets for each segment. The median and maximum depth targets for the IRLS vary among segments from 0.8 to 1.8 and 1.2 to 2.8 m, respectively.Halodule wrightii is typically the dominant seagrass species at the deep-edge of IRLS grass beds. We set light requirement targets by using a 10-yr record of light data (1990–1999) and the union coverage depth limit distributions from the most temporally stable seagrass segments. The average annual light requirement, based on the medians of the depth limit distributions, is 33 ± 17% of the subsurface light. The minimum annual light requirement, based on of the 95th percentile of the depth distributions, is 20 ± 14%; the minimum growing season light requirement (March to mid September) is essentially the same (20 ± 13%). Variation in depth limits and light requirements, is probably due to factors other than light that influence the depth limit of seagrasses (e.g., competition, physical disturbance). The methods used in this study are robust when applied to large or long-term data sets and can be applied to other estuaries where grass beds are routinely monitored and mapped.     

Notes on the occurrence,distribution, and biology of elasmobranch fishes in the Indian River lagoon system,Florida

Fifteen species of elasmobranchs, eight sharks and seven rays, have been recorded with reasonable certainty from the Indian River lagoon system on the central east coast of Florida. We collected four shark and six ray species during a three and one-half year study of the northern portion of the lagoon system. Five of these appear to be year-round residents, and the remainder utilize the area only at restricted times of the year or as a nursery ground. The most abundant resident species areDasyatis sayi, D. sabina, andCarcharhinus leucas. Pristis pectinata, once a common resident species, has been extirpated from the lagoons. The distribution of ocean inlets and salinity appear to be major factors affecting diversity and numbers of elasmobranchs in the Indian River system.     

Tidal and nontidal flushing of Florida’s Indian River Lagoon

The flushing of Florida’s Indian River Lagoon is investigated as a response to tidal and low-frequency lagoon-shelf exchanges in the presence of freshwater gains and losses. A one-dimensional computer model uses the continuity equation to convert water-level variations into both advective transport within the lagoon and lagoon-shelf exchanges. The model also incorporates transport by longitudinal diffusion. Flushing is quantified by calculating the 50% renewal time, R₅₀, for each of 16 segments. R₅₀ is calculated for tidal exchanges enhanced by 0–30 cm nontidal fluctuations in coastal sea level, then for a range of rainfall rates. In both series of simulations, results suggest that in the northern sub-basin, R₅₀ increases dramatically with distance from the inlet due to relatively weak tidal and nontidal exchanges. A 50% renewal occurs in about one tidal cycle just inside Sebastian Inlet; at the northern end of the northern sub-basin, R₅₀ is over 230 d, and only coastal sea-level variations on the order of 30 cm and/or dry season rainfall rates decrease R₅₀ to less than 1 yr. R₅₀ is 1 wk or less throughout the central and southern sub-basins, where lagoon-shelf exchanges occur through two inlets. Simulations involving seasonal variations in precipitation and evaporation indicate that maximum and minimum rates of freshwater input lead minimum and maximum salinities by time periods on the order of 2–3 wk for the lagoon as a whole and in the northern sub-basin. The central and southern sub-basins respond in 1–2 wk.     

The disruption and recovery of fish communities in the Indian River Lagoon, Florida, following two hurricanes in 2004

During the summer of 2004, four hurricanes (Charley, Frances, Ivan, and Jeanne) affected Florida between August 13 and September 27. Two storms (Frances: category 2 and Jeanne: category 3) made landfall in the southern portion of the Indian River Lagoon (IRL) on the east-central coast of Florida. The presence of Florida Fish and Wildlife Conservation Commission's long-term fisheries monitoring program in the IRL provided a unique opportunity to examine the effects of large tropical events on estuarine fish communities. Increased sampling efforts to monitor the effects of tropical disturbances on the fish community within the IRL and one of its major tributaries (St. Sebastian River) were initiated within days after the passing of the last hurricane (Jeanne). The objectives of the study were to characterize changes to the composition of the fish community within the lagoon and river immediately after the passage of two hurricanes, and to examine the recovery of the fish communities. Analyses indicated that immediately after the last hurricane passed, community diversity within the estuary decreased following these storms due to the absence of many marine species, whereas the fish community within the St. Sebastian River shifted to one containing a greater percentage of freshwater species. Recovery of the community structure to pre-hurricane conditions was evident within several weeks following the last hurricane, and by mid December 2004 (ca. 3 mo after the last storm), there was little difference between the pre-hurricane and post-hurricane fish communities.     

Climatic Trends and Temporal Patterns of Phytoplankton Composition, Abundance, and Succession in the Indian River Lagoon, Florida, USA

This paper describes the results of 10 years of water quality monitoring in the Indian River Lagoon Florida, with special emphasis on the relationships between trends in climatic conditions and the distribution, composition, and abundance of the phytoplankton community. The Indian River Lagoon, which spans 220 km of Florida’s east coast, is a region of particular concern because of the rapid rate of human development throughout the region and the hydrologically restricted character of the lagoon, which heightens the potential for algal bloom. Water sampling was carried out on a monthly to twice-monthly basis at six sites located in the northern and central lagoon. The 10-year study included both extended periods of below and above average rainfall. A number of ecologically distinct regions exist within the lagoon, which differ considerably in water exchange properties and watershed inputs. The northern lagoon is characterized by longer water residence times, lower phosphorus concentrations, higher nitrogen concentrations, and more stable salinity conditions than the central lagoon. Mean phytoplankton biovolumes were substantially higher at the sites in the northern lagoon than at the sites in the central lagoon, and algal blooms were more common and intense in the former region. Inter-annual patterns of phytoplankton biovolume were also different in the northern and central lagoon. In the northern lagoon, phytoplankton biovolumes were lowest during the drought period, from the autumn of 1998 to the spring of 2001. By contrast, algal bloom events in the central lagoon were not only less frequent but were not tied to periods of high rainfall. The most widespread and common bloom formers were the potentially toxic dinoflagellate Pyrodinium bahamense var. bahamense and two centric diatoms, Dactyliosolen fragilissimus and Cerataulina pelagica. Many of the biovolume peaks observed over the study period were attributable to these three species. The results of time series modeling of phytoplankton dynamics further highlighted the disparities between the two regions of the lagoon in terms of the suite of parameters that best predict the observed trends in the biomass of phytoplankton. Overall, the outcome of this initial modeling effort in the Indian River Lagoon suggests that time series approaches can help define the factors that influence phytoplankton dynamics.     

The impacts of the 2004 hurricanes on hydrology, water quality, and seagrass in the Central Indian River Lagoon, Florida

Between August 14 and September 26, 2004, four tropical weather systems (Charley, Frances, Ivan, and Jeanne) affected the central Indian River Lagoon (IRL). The central IRL received a prodigious amount of rainfall for the 2 mo, between 72 and 83 cm, which is a once-in-50-yr rainfall event. High stream discharges were generated that, combined with wind-suspended sediments, significantly reduced salinities and water transparency. In September, salinities among central IRL segments dropped from 30 psu or more to ≤15 psu, color increased from a low of 10 pcu to ≥100 pcu, and turbidity increased from ≤3 NTU up to 14 NTU. Evidence of the hurricanes' physical effects on seagrasses (burial, no scour) was limited to just one of the more than 25 sites inspected. Within 2 to 3 mo following the hurricane period, most parameters related to water transparency returned to or showed improvement over their prehurricane (February–July 2004) levels. Unseasonably low salinities (<20 psu) and moderately high color (>20 pcu) were observed through spring 2005, largely attributable to a relatively long residence time and a wetter-than-average spring season in 2005. By the end of the study period (July 2006), the central IRL generally showed a continuation of two opposite seagrass trends—an increase in depthlimit coverage but a decline in coverage density—that began before 2004. Also, within a limited reach of the central IRL, there was a temporary shift in species composition in summer 2005 (Ruppia maritima increased asHalodule wrightü decreased). It is likely that the persistently low salinities (not color) in 2004–2005 affected the species composition and coverage density. This study reveals that seagrasses are resilient to the acute effects of hurricanes and underscores the need to reduce chronic, an thropogenic effects on seagrasses.     

Relationships between phytoplankton standing crop and physical, chemical, and biological gradients in the Suwannee River and plume region, U.S.A.

The Suwannee River (USA) is an amber stained, nutrient rich, blackwater river which flows into relatively clear oceanic waters resulting in the formation of a coastal region with unique physical, chemical, and biological gradients. The intent of this study was to describe the spatial and temporal variability of phytoplankton as it relates to these gradients. Ten stations along a transect ranging from 5 km up river to 31 km offshore, were sampled during four different flow regimes. All four sampling periods included in our study of the Suwannee River and plume region exhibited a similar pattern of phytoplankton abundance; low phytoplankton biomass in the Suwannee River and offshore stations with an area of elevated biomass seaward of the Suwannee River outflow. The results of our analysis of light and nutrient limitation in the region support the hypothesis that this spatial pattern of phytoplankton abundance is strongly influenced by color dependent light limitation in the river and outflow area, combined with nutrient limitation offshore. Our results suggest that both light and nutrient availability control abundance and composition of phyto plankton in this coastal area.     

Seasonal variation of gold and arsenic in biogeochemical samples from a disseminated gold deposit in the Northern Cordillera

An orientation survey was conducted over the Quartz Mountain, Oregon, hot-spring type disseminated gold deposit to address three questions critical to successful exploration of the area: What is the relationship between bedrock geology and the trace-element content of trees and soils; do these relationships change seasonally; are these relationships similar in the two tree species which discontinuously blanket the area?Twig, needle and wood samples were collected four times during the year from both ponderosa pine (Pinus ponderosa) and white fir (Abies concolor). Soils were collected once. All samples were analyzed for Au and As and the wood samples were analyzed also for Sb.Arsenic was the best pathfinder element, having little analytical, spatial, or seasonal variation in the twig and needle samples. Anomalous levels were 150 μg/g (ppm) As in pine and 30 μg/g in fir. Gold showed analytical variation of 20–30%, sample duplicate variation up to 90%, spatial changes, and seasonal variation ranging from winter-summer background [10–20 ng/g (ppb) Au] to anomalies of 300–400 ng/g (ppb) Au in pine and 100 μg/g (ppb) Au in fir in the spring and fall. Antimony could not be completely evaluated as a pathfinder element because it was only determined in wood samples.Needles and twigs from both species proved to be viable sample media for exploration of the area. Wood was not a suitable medium because of low, erratic values, perhaps due in part to analytical technique. Needles had 2 to 20 times higher As concentrations than did twigs. Twigs had a slightly higher Au content [25 ng/g (ppb)] than did needles. The pine samples were higher in both elements than were the fir samples. The data, normalized by species, could be treated as one homogeneous data base.The soil Au and As data outlined the mineralization clearly with a central Au anomaly [100 ng/g (ppb) and greater] surrounded by As anomalies [100 μg/g (ppm) and greater] over a distance of 594 m. Neither spring nor fall tree Au anomalies were coincident with the soil Au anomaly. Consistent throughout the year, tree As anomalies coincided with the soil As anomalies, but covered a smaller area. Both the Au and As anomalies in trees appeared to be related to bedrock contact zones rather than to the soil the trees were growing in.     

Seasonal variation of the salinity in the Zuari estuary,Goa, India

The annual salt budget of the Zuari is examined. The characteristics of the estuary differ markedly from the low run off season during November–May to the heavy run off period of the southwest monsoon from June to October. During November–May the estuary is vertically mixed and the two processes controlling the transport of salt are run off induced advective transport out of the estuary, and tidally induced diffusive transport into the estuary. The magnitude of the latter is about 20% larger, leading to a salinity rise in the estuary. The diffusion coefficient has been estimated to be 233 ± 101 m²/sec. With the onset of the southwest monsoon, the run off increases dramatically, and the estuary loses about 75% of its salt during the first two months of the season. About 2/3 of this loss is recovered in the next two months when the run off decreases. Because the estuary is partially stratified during June–October, gravitational circulation is expected to play a role in addition to tidal diffusion and run off. The magnitude of its contribution has, however, not yet been determined.     

黄河源区多年冻土活动层土壤水文过程季节变异分析

为从整体上认识多年冻土活动层土壤水文过程季节变异特性,以黄河源区巴颜喀拉山北坡冻土剖面为例,结合大气降水、冻土土壤水分、冻土层上水的野外观测,采用HYDRUS-1D软件冻融模块进行模拟分析,分析冻融作用对活动层土壤水文过程的影响,研究结果表明:(1)冻土层上水位与土壤水热之间存在着相互影响、相互作用的关系,依据活动层土壤温度变化,基于冻融过程,多年冻土活动层土壤水分与冻土层上水位可划分为冻结稳定、快速融化、融化稳定和快速冻结4个阶段。(2)降雨入渗是坡面尺度下活动层土壤水文过程的主要驱动力,活动层冻融锋面是主要限制性因素,受冻融过程影响,冻结期降雨减少,土壤冻结,土壤储水能力下降,土壤水分下渗停止,坡面侧向流动减弱,土壤水分和冻土层上水位处于下降趋势;融化期降雨增多,土壤融化,土壤储水能力上升,土壤水分下渗强烈,坡面侧向流动增强,土壤水分和冻土层上水位处于上升趋势。(3)受坡面地形影响,上坡活动层厚度大于下坡,上坡冻融锋面变化较下坡平缓,上坡土壤水分和冻土层上水位的变化幅度相对下坡较为平缓,而上坡土壤水分相对下坡含量较低,下坡冻土层上水位相对稳定。     

长江干流颗粒有机碳及其同位素组成的季节性输送特征

2007年夏季和冬季对长江干流悬浮物进行了季节性采样,系统分析了悬浮物颗粒有机碳含量及稳定碳同位素的组成,研究了其空间分布以及季节变化特征。结果显示,长江干流夏季颗粒有机碳含量在0．4％～1．3％,冬季含量在0．7%-2．2％。冬季和夏季颗粒有机碳平均δ^13C值分别为-24．74‰和～24．83‰,季节性差异不大,在...     

Alongshore variation in the rip current hazard at Pensacola Beach, Florida

Many drowning and near drownings at Pensacola Beach, Florida are attributed to rip currents, the strong seaward-flowing currents that extend from the shoreline to the line of breakers and sometimes beyond. While surf forecasts assume that the rip hazard is uniform alongshore and that the (erosion) rips are ephemeral features, evidence is presented to suggest that the rip hazard at Pensacola Beach is not uniform alongshore. Rather the rip current “hotspots” develop as a consequence of an alongshore variation in the surf similarity parameter and nearshore state on the order of ~1,450 m. The variation is forced by transverse ridges on the inner shelf that force wave refraction and focusing at the ridge crests. This creates a more dissipative, rhythmic bar and beach morphology at the ridges and rougher surf. Between ridges, where wave heights and periods are smaller and the outermost bar is forced closer to the shoreline, the nearshore is in a (more reflective) bar and rip state during red flag conditions. Drownings between 2000 and 2009 are shown to be clustered between transverse ridges and in the years following a hurricane or tropical storm (2000–2003 and 2005–2008) when the bar and rip morphology first develops as the shore face recovers. This continues until the innermost bar attaches to the beach face unless the bar system is reset by another tropical storm or hurricane. It is argued that the rip hazard is dependent on the alongshore covariation of the environmental forcing with the individual and group behavior in both time and space, even on what appears to be a relatively uniform beach environment.     

Submarine groundwater discharge is an important net source of light and middle REEs to coastal waters of the Indian River Lagoon, Florida, USA

Porewater (i.e., groundwater) samples were collected from multi-level piezometers across the freshwater-saltwater seepage face within the Indian River Lagoon subterranean estuary along Florida’s (USA) Atlantic coast for analysis of the rare earth elements (REE). Surface water samples for REE analysis were also collected from the water column of the Indian River Lagoon as well as two local rivers (Eau Gallie River, Crane Creek) that flow into the lagoon within the study area. Concentrations of REEs in porewaters from the subterranean estuary are 10-100 times higher than typical seawater values (e.g., Nd ranges from 217 to 2409 pmol kg⁻¹), with submarine groundwater discharge (SGD) at the freshwater-saltwater seepage face exhibiting the highest REE concentrations. The elevated REE concentrations for SGD at the seepage face are too high to be the result of simple, binary mixing between a seawater end-member and local terrestrial SGD. Instead, the high REE concentrations indicate that geochemical reactions occurring within the subterranean estuary contribute substantially to the REE cycle. A simple mass balance model is used to investigate the cycling of REEs in the Indian River Lagoon and its underlying subterranean estuary. Mass balance modeling reveals that the Indian River Lagoon is approximately at steady-state with respect to the REE fluxes into and out of the lagoon. However, the subterranean estuary is not at steady-state with respect to the REE fluxes. Specifically, the model suggests that the SGD Nd flux, for example, exported from the subterranean estuary to the overlying lagoon waters exceeds the combined input to the subterranean estuary from terrestrial SGD and recirculating marine SGD by, on average, ∼100 mmol day⁻¹. The mass balance model also reveals that the subterranean estuary is a net source of light REEs (LREE) and middle REEs (MREE) to the overlying lagoon waters, but acts as a sink for the heavy REEs (HREE). Geochemical modeling and statistical analysis further suggests that this fractionation occurs, in part, due to the coupling between REE cycling and iron redox cycling within the Indian River Lagoon subterranean estuary. The net SGD flux of Nd to the Indian River Lagoon is ∼7-fold larger than the local effective river flux to these coastal waters. This previously unrecognized source of Nd to the coastal ocean could conceivably be important to the global oceanic Nd budget, and help to resolve the oceanic “Nd paradox” by accounting for a substantial fraction of the hypothesized missing Nd flux to the ocean.