Measuring and modelling temporal trends of 226Ra in waters of a Spanish estuary affected by the phosphate industry

The presence and temporal evolution (1990-2001) of (226)Ra in a tidal estuary affected by the phosphate industry has been investigated. Water samples collected in the course of four separate sampling campaigns were analysed for (226)Ra content using a gas flow proportional counter following Ba coprecipitation. Two (226)Ra sources have been identified: direct discharges from the industrial complex and run-off from a phosphogypsum pile. Although activity levels are similar, or even higher, than those found in other environments affected by the phosphate industry, there has been a general decrease in contamination since direct discharges ceased in 1998 due to new regulations from the EU. However, sediments are now acting as a source of Ra to the water column due to redissolution processes. A numerical model of the estuary has been developed to describe quantitatively the experimental results. The model solves the hydrodynamics and the dispersion equation of (226)Ra including interactions with sediments. Model results are, in general, in good agreement with observations.     

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.     

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.     

Interaction of Surface Waves in a Basin with Floating Broken Ice

The multiscale method is used to obtain asymptotic expansions up to the quantities of the third order for the elevations of the surface of the basin and the velocity potential of motion of liquid particles in the wave disturbances formed in the process of nonlinear interaction of periodic running waves of the first and second harmonics in a homogeneous ideal incompressible liquid of constant finite depth covered with broken ice. The dependences of the amplitude-phase structure of disturbances on the ice thickness, depth of the basin, and the parameters of interacting harmonics are investigated. We estimate the error of evaluation of the characteristics of the formed vertical displacement of the surface of the basin and nonlinear mass transfer introduced by neglecting the curvature of the wave profile in the expression for the velocity potential in deducing the kinematic and dynamic surface boundary conditions for nonlinear approximations.     

Mineral assemblages occurring around hydrocarbon vents in the northern Gulf of Mexico

Magnesian calcites are the most abundant authigenic minerals associated with hydrocarbon vents at 25 sites, in water depths ranging from 100 to 600 m in the Green Canyon area and about 2200 m in the Alaminos Canyon area on the Continental Slope of the northern Gulf of Mexico. The most frequently encountered magnesian calcites have 10–15 mol% MgCO₃ and the apparent structural disorder revealed by XRD peak widths increases with Mg substitution. There are no systematic variations in Mg content with respect to water depth or geographic location. The calcite saturation state of the precipitating fluid is primarily determined by the nature of the fluids escaping from the vents, not the ambient seawater.     

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.