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.     

Corrections to the improved lunar ephemeris

Based primarily upon the formation of new conditional equations using analytical partial derivatives of the moon's mean elements, meridian circle observations of the moon from 1952–67 have been examined to determine corrections to the constants of lunar theory and to the fundamental coordinate system (FK4). With certain exceptions, the new corrections are in agreement with those published earlier by the author. Systematic corrections to FK4 are surprisingly large, although in agreement with some other recent determinations. New corrections to the lunar ephemeris, resulting from the discussion, are also presented.     

Libration effects for retrograde satellites in the restricted three-body problem

Numerical explorations of the restricted problem have shown that for stable large nonperiodic retrograde satellite orbits, the motion can be decomposed into a fast reference motion and a slow libration aroundB ₂ We study here this libration in the circular plane Hill's case, for which the reference motion is elliptic. We establish the equations of motion for the coordinates of the centre of this ellipse. We find two integrals of motion: the first is the semi-major axis of the ellipse; the second is essentially Jacobi's integral, translated into the new coordinates. We give a formula for the period of the libration and we find its limiting value for small libration amplitudes. A numerical verification gives very good agreement for all these results.     

The fusion crust of the Winchcombe meteorite: A preserved record of atmospheric entry processes

Fusion crusts form during the atmospheric entry heating of meteorites and preserve a record of the conditions that occurred during deceleration in the atmosphere. The fusion crust of the Winchcombe meteorite closely resembles that of other stony meteorites, and in particular CM2 chondrites, since it is dominated by olivine phenocrysts set in a glassy mesostasis with magnetite, and is highly vesicular. Dehydration cracks are unusually abundant in Winchcombe. Failure of this weak layer is an additional ablation mechanism to produce large numbers of particles during deceleration, consistent with the observation of pulses of plasma in videos of the Winchcombe fireball. Calving events might provide an observable phenomenon related to meteorites that are particularly susceptible to dehydration. Oscillatory zoning is observed within olivine phenocrysts in the fusion crust, in contrast to other meteorites, perhaps owing to temperature fluctuations resulting from calving events. Magnetite monolayers are found in the crust, and have also not been previously reported, and form discontinuous strata. These features grade into magnetite rims formed on the external surface of the crust and suggest the trapping of surface magnetite by collapse of melt. Magnetite monolayers may be a feature of meteorites that undergo significant degassing. Silicate warts with dendritic textures were observed and are suggested to be droplets ablated from another stone in the shower. They, therefore, represent the first evidence for intershower transfer of ablation materials and are consistent with the other evidence in the Winchcombe meteorite for unusually intense gas loss and ablation, despite its low entry velocity.     

A seasonal climate model of the atmospheres of the giant planets at the voyager encounter time: I. Saturn's stratosphere

A radiative seasonal model which incorporates a multilayer radiative transfer treatment at wave-lengths longward of 7 μm is presented and applied to Saturn's stratosphere. Opacities due to H₂-He, CH₄, C₂H₂, and C₂H₆ are included. Season-dependent insolation is shown to produce a strong hemispheric asymmetry decreasing with depth at the Voyager encounter times, and seasonal amplitudes of 30°K at the poles are predicted in the high stratosphere. The ring-modulated dependence of the insolation and the orbital eccentricity are shown to have a significant effect. Calculations agree closely with the Voyager 1 and 2 radio occultation ingress profiles recorded at 76°S and 36.5°S for CH₄/H₂ = 3.5 + 1.4/? 1.0 × 10^?3;the estimated errors include modeling systematic errors and uncertainties in the occultations profiles. The possible role of aerosols in the stratospheric heating is analyzed. The Voyager 2 egress profile recorded at 31°S cannot be reproduced by calculations. Some constraints on the C₂H₂ and C₂H₆ abundances are derived. The upper portion of the occultation profiles (p < 3mbar) can be matched for C₂H₂/H₂ = 1.0 + 1.3/?0.6 × 10^?7, C₂H₆/H₂ = 1.5 + 1.8/?0.9 × 10^?6 at 76°S and C₂H₂/H₂ = 4 + 6/?4 × 10^?8, C₂H₆/H₂ = 6 + 9/?6 × 10^?7 at 36.5°N. At the northern occultation latitude, the discrepancy with the concentrations derived from analysis of IRIS spectra by R. Courtin, D. Gautier, A. Marten, B. Bézard, and R. Hanel (1984, Astrophys. J.287) can be explained by a sharp variation of the mixing ratios of these gases with altitude in the upper stratosphere. Other interpretations are discussed.     

Physiological and behavioral responses of Bathynerita naticoidea (Gastropoda: Neritidae) and Methanoaricia dendrobranchiata (Polychaeta: Orbiniidae) to hypersaline conditions at a brine pool cold seep

Bathynerita naticoidea (Gastropoda: Neritidae) and Methanoaricia dendrobranchiata (Polychaeta: Orbiniidae) are two of the most abundant invertebrates associated with cold‐seep mussel beds in the Gulf of Mexico. At the methane seep known as Brine Pool NR‐1 (27 °43.415 N, 91 °16.756 W; 650 m depth), which is surrounded by a broad band of mussels (Bathymodiolus childressi), these species have distinctly different patterns of abundance, with the gastropod being found mostly at the outer edge of the mussel bed (average density in November 2003: 817 individuals·m⁻² in outer zone, 20·m⁻² in inner zone) and the polychaete being found almost exclusively near the inner edge (average density in November 2003: 3155 individuals·m⁻² in inner zone, 0·m⁻² in outer zone), adjacent to the brine pool itself. The salinity of the brine pool exceeds 120, so we hypothesized that M. dendrobranchiata should be more tolerant of high salinities than B. naticoidea. The opposite proved to be true. The gastropods were capable of withstanding salinities at least as high as 85, whereas the polychaetes died at salinities higher than 75. Both species were osmoconformers over the range of salinities (35–75) tested. Behavioral responses of B. naticoidea to salinities of 50, 60, and 70 were investigated in inverted vertical haloclines. Gastropods generally did not enter water of salinity greater than 60, but tolerated short periods at 60. Behavioral avoidance of brine should limit the vertical distribution of B. naticoidea in the inner zone to the top 2.5–5 cm of the mussel bed. Behavior is also a likely (though unproven) mechanism for controlling horizontal distribution of this species across the mussel bed. Methanoaricia dendrobranchiata can tolerate short excursions into the brine, but probably avoids hypersaline conditions by aggregating on the tops of the mussels.     

Impact of grazing by benthic eukaryotic organisms on the nitrogen sediment–water exchange in the North Sea

The metabolic inhibitor cycloheximide was used to estimate the influence of primarily unicellular eukaryotes (heterotrophic protozoa) on nutrient recycling in different types of sediments in the North Sea. Fluxes of dissolved inorganic nitrogen across the sediment–water interface were measured in undisturbed sediment cores (controls) and compared to fluxes in sediment cores with cycloheximide added. If eukaryotes play an important role in nutrient recycling, one would expect to find lower nitrogen sediment–water effluxes in cores with cycloheximide due to the inactivation of eukaryotes. This important role hypothesised for eukaryotes was not generally observed: Only at four of the nineteen stations were ammonium effluxes significantly higher in controls than in cores with cycloheximide, and at five stations nitrate effluxes were significantly higher in the controls than in the cores with cycloheximide. Eukaryotic activity apparently contributed to the sediment–water exchange of ammonium through mineralisation of organic matter, nitrification and the subsequent release of ammonium and nitrate at these stations. At most other stations no differences were obtained between controls and cores with cycloheximide. This suggests that bacteria were the most important nutrient mineralisers at these stations at the time of the cruises.     

A multiscaling‐based intensity–duration–frequency model for extreme precipitation

Rainfall intensity–duration–frequency (IDF) curves are a standard tool in urban water resources engineering and management. They express how return levels of extreme rainfall intensity vary with duration. The simple scaling property of extreme rainfall intensity, with respect to duration, determines the form of IDF relationships. It is supposed that the annual maximum intensity follows the generalized extreme value (GEV) distribution. As well known, for simple scaling processes, the location parameter and scale parameter of the GEV distribution obey a power law with the same exponent. Although, the simple scaling hypothesis is commonly used as a suitable working assumption, the multiscaling approach provides a more general framework. We present a new IDF relationship that has been formulated on the basis of the multiscaling property. It turns out that the GEV parameters (location and scale) have a different scaling exponent. Next, we apply a Bayesian framework to estimate the multiscaling GEV model and to choose the most appropriate model. It is shown that the model performance increases when using the multiscaling approach. The new model for IDF curves reproduces the data very well and has a reasonable degree of complexity without overfitting on the data.