Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study

A Lagrangian particle method embedded within a 2-D finite element code, is used to study the transport and ocean–estuary exchange processes in the well-mixed Great Bay Estuarine System in New Hampshire, USA. The 2-D finite element model, driven by residual, semi-diurnal and diurnal tidal constituents, includes the effects of wetting and drying of estuarine mud flats through the use of a porous medium transport module. The particle method includes tidal advection, plus a random walk model in the horizontal that simulates sub-grid scale turbulent transport processes. Our approach involves instantaneous, massive [O(500,000)] particle releases that enable the quantification of ocean–estuary and inter-bay exchanges in a Markovian framework. The effects of the release time, spring–neap cycle, riverine discharge and diffusion strength on the intra-estuary and estuary–ocean exchange are also investigated.The results show a rather dynamic interaction between the ocean and the estuary with a fraction of the exiting particles being caught up in the Gulf of Maine Coastal Current and swept away. Three somewhat different estimates of estuarine residence time are calculated to provide complementary views of estuary flushing. Maps of residence time versus release location uncover a strong spatial dependency of residence time within the estuary that has very important ramifications for local water quality. Simulations with and without the turbulent random walk show that the combined effect of advective shear and turbulent diffusion is very effective at spreading particles throughout the estuary relatively quickly, even at low (1 m²/s) diffusivity. The results presented here show that a first-order Markov Chain approach has applicability and a high potential for improving our understanding of the mixing processes in estuaries.     

Modeling bed-load transport of coarse sediments in the Great Bay Estuary, New Hampshire

Current, sea level and bed-load transport are investigated in the Lower Piscataqua River section of the Great Bay Estuary, New Hampshire, USA—a well-mixed and geometrically complex system with low freshwater input, having main channel tidal currents ranging between 0.5 and 2 m s⁻¹. Current and sea level forced by the M₂M₄M₆ tides at the estuarine mouth are simulated by a vertically averaged, non-linear, time-stepping finite element model. The hydrodynamic model uses a fixed boundary computational domain and accounts for flooding–drying of tidal flats by making use of a groundwater component. Inertia terms are neglected in comparison with pressure gradient and bottom friction terms, which is consistent with the observed principal dynamic balance for this section of the system. The accuracy of hydrodynamic predictions in the study area is demonstrated by comparison with four tidal elevation stations and two cross-section averaged current measurements. Simulated current is then used to model bed-load transport in the vicinity of a rapidly growing shoal located in the main channel of the lower system. Consisting of coarse sand and gravel, the shoal must be dredged every five to eight years. Two approaches are taken—an Eulerian parametric method in which nodal bed-load flux vectors are averaged over the tidal cycle and a Lagrangian particle tracking approach in which a finite number of sediment particles are released and tracked. Both methods yield pathways and accumulations in agreement with the observed shoal formation and the long-term rate of sediment accumulation in the shoal area.     

Groundwater control of mangrove surface elevation: Shrink and swell varies with soil depth

We measured monthly soil surface elevation change and determined its relationship to groundwater changes at a mangrove forest site along Shark River, Everglades National Park, Florida. We combined the use of an original design, surface elevation table with new rod-surface elevation tables to separately track changes in the mid zone (0–4 m), the shallow root zone (0–0.35 m), and the full sediment profile (0–6 m) in response to site hydrology (daily river stage and daily groundwater piezometric pressure). We calculated expansion and contraction for each of the four constituent soil zones (surface [accretion and erosion; above 0 m], shallow zone [0–0.35 m], middle zone [0.35–4 m], and bottom zone [4–6]) that comprise the entire soil column. Changes in groundwater pressure correlated strongly, with changes in soil elevation for the entire profile (Adjusted R² = 0.90); this relationship was not proportional to the depth of the soil profile sampled. The change in thickness of the bottom soil zone accounted for the majority (R² = 0.63) of the entire soil profile expansion and contraction. The influence of hydrology on specific soil zones and absolute elevation change must be considered when evaluating the effect of disturbances, sea level rise, and water management decisions on coastal wetland systems.     

Long‐term streamflow trends in Hawai'i and implications for native stream fauna

Climate change has fundamentally altered the water cycle in tropical islands, which is a critical driver of freshwater ecosystems. To examine how changes in streamflow regime have impacted habitat quality for native migratory aquatic species, we present a 50‐year (1967–2016) analysis of hydrologic records in 23 unregulated streams across the five largest Hawaiian Islands. For each stream, flow was separated into direct run‐off and baseflow and high‐ and low‐flow statistics (i.e., Q10 and Q90) with ecologically important hydrologic indices (e.g., frequency of flooding and low flow duration) derived. Using Mann–Kendall tests with a running trend analysis, we determined the persistence of streamflow trends through time. We analysed native stream fauna from ~400 sites, sampled from 1992 to 2007, to assess species richness among islands and streams. Declines in streamflow metrics indicated a general drying across the islands. In particular, significant declines in low flow conditions (baseflows), were experienced in 57% of streams, compared with a significant decline in storm flow conditions for 22% of streams. The running trend analysis indicated that many of the significant downward trends were not persistent through time but were only significant if recent decades (1987–2016) were included, with an average decline in baseflow and run‐off of 10.90% and 8.28% per decade, respectively. Streams that supported higher native species diversity were associated with moderate discharge and baseflow index, short duration of low flows, and negligible downward trends in flow. A significant decline in dry season flows (May–October) has led to an increase in the number of no‐flow days in drier areas, indicating that more streams may become intermittent, which has important implications for mauka to makai (mountain to ocean) hydrological connectivity and management of Hawai'i's native migratory freshwater fauna.     

Dynamics and role of the Durban cyclonic eddy in the KwaZulu-Natal Bight ecosystem

The semi-permanent Durban Eddy is a mesoscale, lee-trapped, cold-core cyclonic circulation that occurs off the east coast of South Africa between Durban in the north and Sezela, some 70 km to the south. When present, strong north-eastward countercurrents reaching 100 cm s^–1 are found inshore. It is hypothesised that the cyclone is driven by the strong south-westward flowing Agulhas Current offshore of the regressing shelf edge near Durban. Analysis of ADCP data and satellite imagery shows the eddy to be present off Durban approximately 55% of the time, with an average lifespan of 8.6 days, and inter-eddy periods of 4 to 8 days. After spin-up the eddy breaks loose from its lee position and propagates downstream on the inshore boundary of the Agulhas Current. The eddy is highly variable in occurrence, strength and downstream propagation speeds. There is no detectable seasonal cycle in eddy occurrence, with the Natal Pulse causing more variability than any seasonal signal. A thermistor array deployed in the eddy centre, together with ship CTD data, indicates upward doming of the thermal structure in the eddy core associated with cooler water and nutrients being moved higher in the water column, stimulating primary production. Together with the use of satellite imagery, our findings indicate a second mechanism of upwelling, viz. divergent upwelling in the northern limb of the eddy. Satellite-tracked surface drifters released in the eddy demonstrated the potential for nutrient-rich eddy water to be transported northwards along the inshore regions of the KwaZulu-Natal (KZN) Bight, thus contributing to the functioning of the bight ecosystem, as well as southwards along the KZN and Transkei coasts – both by the eddy migrating downstream and by eddy water being recirculated into the inshore boundary of the Agulhas Current itself.     

Sharks caught in the KwaZulu-Natal bather protection programme,South Africa. 14. The smooth hammerhead shark Sphyrna zygaena (Linnaeus)

The current study provides long-term catch-rate, biological and feeding data for smooth hammerhead sharks, Sphyrna zygaena, caught in South Africa’s KwaZulu-Natal bather protection programme. In total, 2 512 S. zygaena were caught in net installations between 1978 and 2014, and 72 S. zygaena were caught on drumlines between 2007 and 2014. There was no significant log-linear year trend in the net catch rate over time (slope = 0.0054, t = 1.808, p = 0.07). However, there was a significant temporal increase in mean size of the captured sharks (slope = 0.0012, t = 3.502, p < 0.001). A quasi-Poisson generalised additive mixed model showed that increasing latitude, winter months, colder sea temperatures and the deployment of drumlines all had a significant positive effect on the catch rate of sharks in nets. The size frequency of the catch was unimodal, with significantly more females caught in the nets and more males on the drumlines. The majority (93.1%) of all sharks caught were immature and measured between 80 and 120 cm precaudal length. Teleosts and cephalopods dominated the sharks’ diet in terms of all dietary indices. The prey species consumed indicate that immature S. zygaena are feeding primarily within the pelagic zone of shallow coastal habitats.     

Vulnerabilities,competition and rights in a context of climate change toward equitable water governance in Peru's Rio Santa Valley

As glacial retreat changes the hydrology of Peru's Rio Santa, water demand is growing, pollution is worsening, and competition for water among economic sectors, political jurisdictions and upstream and downstream water users is intensifying. The vulnerability of highland communities, food producers, and poor urban neighborhoods in the Santa watershed in the face of these changes is magnified by inequities in water governance, giving rise to water conflict. Peru's new water regime defines water as an economic good and seeks to centralize control over water. This article analyzes implications of this regime for ensuring equity and managing conflict. It concludes that Peru's water regime is more likely to address equity issues when faced with concerted citizen action.     

Recent accretion in mangrove ecosystems based on137Cs and210Pb

Accretion rates were measured in fringe and basin mangrove forests in river and tidally dominated sites in Terminos Lagoon, Mexico, and a basin mangrove forest in Rookery Bay, Florida, USA. Accretion rates were determined using the radionuclides²¹⁰Pb and¹³⁷Cs. Consolidation-corrected accretion rates for the Rookery Bay cores, ranged from 1.4 to 1.7 mm yr^?1, with an average rate of 1.6 mm yr^?1. Rates at the Mexico sites ranged from 1.0 to 4.4 mm yr^?1, with an average of 2.4 mm yr^?1. Determination of rates in these mangrove forests was greatly affected by the consolidation corrections which decreased the apparent accretion rate by over 50% in one case. Accretion rates at basin sites compare favorably with a reported 1.4 to 1.6 mm yr^?1 rate of sea-level rise, indicating little or no subsidence at inland locations. Accretion rates in fringe sites are generally greater than basin sites, indicating greater subsidence rates in these sediments over longer time intervals.     

A regional model for studies of atmosphere-ice-ocean interaction in the western Arctic

Summary Asa step in the development of a fully coupled regional model of the atmosphere-ice-ocean system, atmospheric and sea ice models have been adapted to a western Arctic domain centered on the Bering Strait. Lateral boundary conditions derived from operational analyses drive the models through simulations on grids having horizontal resolutions of 21 km and 7 km. Sensitivities to the presence of sea ice are large after only 48 hours, by which time the surface temperatures in the Bering and Chukchi Seas are 10–15°C higher without sea ice than with sea ice. The temperatures, in turn, modify the fields of sea level pressure, surface wind and precipitation. By influencing the surface wind stress through the static static stability, the surface state feeds back to the surface momentum exchange, ice/ocean transport, and the rate of formation of new ice. The results also show a resolution-dependence of the surface winds, precipitation rates and new ice formation rates, particularly in areas in which the coastal configuration and topography are spatially complex. The experiments will be augmented by the implementation of an ocean model on the same grids.With 12 Figures