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.     

A nitrate and silicate budget in the equatorial Pacific Ocean: a coupled physical–biological model study

A coupled physical–biological model was developed to simulate the low-silicate, high-nitrate, and low-chlorophyll (LSHNLC) conditions in the equatorial Pacific Ocean and used to compute a detailed budget in the Wyrtki box (5°N–5°S, 180–90°W) for the major sources and cycling of nitrogen and silicon in the equatorial Pacific. With the incorporation of biogenic silicon dissolution, NH₄ regeneration from organic nitrogen and nitrification of ammonia in the model, we show that silicon recycling in the upper ocean is less efficient than nitrogen. As the major source of nutrients to the equatorial Pacific, the Equatorial Undercurrent provides slightly less Si(OH)₄ than NO₃ to the upwelling zone, which is defined as 2.5°N–2.5°S. As a result, the equatorial upwelling supplies less Si(OH)₄ than NO₃ into the euphotic zone in the Wyrtki box, having a Si/N supply ratio of about 0.85 (2.5 vs. 2.96 mmolm⁻² day⁻¹). More Si(OH)₄ than NO₃ is taken up with a Si/N ratio of 1.17 (2.72 vs. 2.33 mmolm⁻² day⁻¹) within the euphotic zone. The difference between upwelling supply and biological uptake is balanced by nutrient regeneration and horizontal advection. Excluding regeneration, the net silicate and nitrate uptakes are nearly equal (1.76 vs. 1.84 mmolm⁻² day⁻¹). However, biogenic silica export production is slightly higher than organic nitrogen (1.74 vs. 1.59 mmolm⁻² day⁻¹) following a 1.1 Si/N ratio. In the central equatorial Pacific, low silicate concentrations limit diatom growth; therefore non-diatom new production accounts for most of the new production. Higher silicate supply in the east maintains elevated diatom growth rates and new production associated with diatoms dominate upwelling zone. In contrast, the new production associated with small phytoplankton is nearly constant or decreases eastward along the equator. The total new production has a higher rate in the east than in the west, following the pattern of surface silicate. This suggests that silicate regulates the diatom production, total new production, and thereby carbon cycle in this area. The modeled mean primary production is 48.4 mmolCm⁻² day⁻¹, representing the lower end of direct field measurements, while new production is 15.0 mmolCm⁻² day⁻¹, which compares well with previous estimates.     

Inter-annual variability in larval supply to populations of three invertebrate taxa in the northern California Current

We investigated sources of inter-annual variability in larval supply to crab and sea urchin populations at Bodega Head and Point Reyes in northern California. During the spring and summer upwelling seasons of the years 1992 through 1997 we monitored the weekly settlement rates of nine species of crabs and two species of sea urchins. As observed in previous studies, daily values of alongshore windstress, temperature and salinity provided evidence for the poleward flow of relatively warm, low salinity water from south of Point Reyes, an apparent retention zone, during upwelling relaxation events. In years dominated by these events (1992, 1993, 1995 and 1996) we observed that alongshore windstress, temperature and salinity were coherent and temperature was significantly correlated with cancrid crab settlement. During these years the magnitude of cancrid crab settlement and the fraction of cancrid crabs relative to other crab species settling were high. Over four years of concurrent sampling there was consistently greater cancrid crab settlement at the Point Reyes site, within the retention zone, than at Bodega Head. Settlement of non-cancrid crabs (porcellanids, grapsids, pagurids and majids) was not as closely linked to intra-annual patterns of upwelling and relaxation, possibly due to the shorter seasonal availability of larvae allowing for the influence of fewer relaxation events. Settlement of this group among years was positively correlated with environmental indicators of strong seasonal upwelling; high salinity, Bakun upwelling index and low temperature. Sea urchin settlement events were observed in June and July of 1992, 1994 and 1997 during warming periods when salinity and temperature were increasing and alongshore windstress was low. Across the six years of the study, we found that cancrid crab larvae had a more even seasonal availability than larvae of non-cancrid species, which settled in greatest numbers during the early portion of the upwelling season. Sea urchins settled in greatest numbers during the later part of the upwelling season. Together these patterns demonstrate the taxon-specific way that inter-annual variability in larval supply is forced by the coincidence of larval availability with favorable physical transport mechanisms.     

Multi-bioindicator study of a small estuary in Vendée (France)

Surficial sediment (about 5 mm) was scraped off in 15 stations selected in the Vie Estuary, a small estuary on the French Atlantic coast. The samples, prepared using the suitable methods, were used for comparing the behaviour of periphytic algae, sporopollinic and other palynologic material, ostracoda and foraminifera for proposing a synthetic evaluation of the information provided by each of them. Dinoflagellate cysts, which are found in an almost monospecific assemblage in low-salinity water, are dispersed by tidal currents, but the location of tintinnids in only the lower reaches shows a limited inward transport. Sporopollinic material is clearly related to the adjacent vegetation, indicating a local origin, and then little influence of transport on the microbenthos, even if fragments of freshwater crustaceans and of Pediastrum alga are found down to the mouth in the sporopollinic deposits. This limited transport of microbenthos allows the use of periphytic algae, benthic ostracoda and benthic foraminifera as indicators of local environmental conditions. It was then possible to determine the extent of marine influence in the estuary, to discriminate between the main channel and adjacent basins, and between different stages of evolution of the basins, from a recently dug basin to a eutrophicated one.     

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.     

The relationship between the amplitudinal characteristics of the high-frequency spectral components of wind-generated waves and the wind velocity over the sea

The relationship between the RMS amplitudes of the wind wave spectral components and the wind speed has been studied at ten frequencies in the band of 0.65–23 Hz. To measure the parameters of the high-frequenci waves, a resistance elevation wave gauge was operated, which was deployed in the Black See on an oceanographic platform near Katsively. The correlation between the wave amplitudes and the wind velocity at high frequencies of 5–23 Hz, corresponding to gravitation-capillary ripples, was found to reach a value of 0.8. At lower frequencies of 0.65–4.3 Hz, corresponding to short gravity waves, it dropped to 0.5–0.7. The response of spectral components to the wind speed variations in the gravity-capillary range is higher than in the range of short gravity waves. The results obtained differ from Phillips' idea about a saturated range for the frequency form of the spectrum of high-frequency gravity waves, since a linear dependence of the spectral amplitudes on the wind speed is established at a wind of force 1–8.Translated by Mikhail M. Trufanov.     

Heated diver decompression shelter: The “womb concept”

A hybrid thermal protection method using waste heat from a surface-mounted outboard motor is shown to create a warm “micro-climate” environment for divers. The effects of surface heater capacities, water flow rates, shelter volume and shelter insulation on micro-climate temperatures are characterized. During long, cold-water decompression stops this method offers a reliable, low-cost alternative to surface-supplied hot water suits or diver-carried heating systems. An added bonus for divers using closed-circuit breathing apparatus is prolonged durations of their carbon dioxide scrubbers when surrounded by the warm water “micro-climate”. Closed-circuit and open-circuit options of this diver decompression shelter concept are evaluated.     

Numerical experiments with data assimilation in the simplest dynamic-stochastic model for the ocean's synoptic variability closed at the level of the second moments

The results of numerical experiments with data assimilation in a non-linear model of the synoptic dynamics of the world's oceans are reported. The model's response to the data assimilation is studied at different values of the relative error in determining the initial fields and assimilation discreteness. The results of application of the algorithms of optimum filtration and optimum interpolation are compared. The advantages of the filtration algorithm employed are discussed only for the case of small errors in determining the initial fields and low assimilation discreteness. The possibility of deterioration of the prediction estimates in data assimilation with discreteness less than the typical time of variability of the average fields is noted.Translated by Mikhail M. Trufanov.     

The bottom currents in the area of abyssal hills in the north-east tropical Pacific Ocean

The results of simultaneous measurements of the bottom (6.25 and 35 m above the bottom) currents, deep currents, and surface currents made at three points in the north-east tropical Pacific Ocean are given. The bottom intensification of the current velocity is revealed in a layer of 35–25 m above the bottom. The estimation of the thickness of the bottom boundary layer (BBL) indicates that the velocity intensification is observed over the boundary layer upper border. A 10-day long benthic storm with a maximum measured velocity of 13 cm/s was revealed 6 m above the bottom. As was found, the origin of the benthic storm is associated with the penetration of an anticyclonic eddy down to the bottom.Translated by Mikhail M. Trufanov.