NPZ Models of Plankton Dynamics: Their Construction, Coupling to Physics, and Application

Nutrient-phytoplankton-zooplankton (NPZ) models have been in use in oceanography for at least three decades, and are still a common research tool. Given the discoveries of the last two decades, particularly concerning the role of bacteria in the plankton, there are questions as to whether NPZ models can still be supported as a useful tool in planktonic research. Here I review the construction of NPZ models, and some of the physical platforms they have been coupled to. I then discuss the applications of NPZ-physical models, and conclude that they still constitute an important and viable research tool, provided that the questions being explored are clearly stated. This revised version was published online in August 2006 with corrections to the Cover Date.     

Casualties in dry bulk shipping (1963–1996)

In recent years, bulk carriers have been identified with high risks of catastrophic structural failure and foundering, and with heavy loss of human life. This study, based upon Lloyd's of London casualty records, identified four risk factors that had significant, independent effect upon the likelihood of a bulk carrier foundering. The risk of foundering increased with the age of the ship, and was related to the ship's flag of registration. Most importantly, however, increased risks were found for heavy cargoes of iron ore and scrap steel or iron, and for trading routes to the Far East and from Europe to North America. Additional safety measures, in particular regarding ship design and high-risk trades, may well be necessary to reduce the high casualty rates.     

Light Requirements for Growth and Survival of Eelgrass (<Emphasis Type="Italic">Zostera marina</Emphasis> L.) in Pacific Northwest (USA) Estuaries

We developed light requirements for eelgrass in the Pacific Northwest, USA, to evaluate the effects of short- and long-term reductions in irradiance reaching eelgrass, especially related to turbidity and overwater structures. Photosynthesis-irradiance experiments and depth distribution field studies indicated that eelgrass productivity was maximum at a photosynthetic photon flux density (PPFD) of about 350–550 μmol quanta m⁻² s⁻¹. Winter plants had approximately threefold greater net apparent primary productivity rate at the same irradiance as summer plants. Growth studies using artificial shading as well as field monitoring of light and eelgrass growth indicated that long-term survival required at least 3 mol quanta m⁻² day⁻¹ on average during spring and summer (i.e., May-September), and that growth was saturated above about 7 mol quanta m⁻² day⁻¹. We conclude that non-light-limited growth of eelgrass in the Pacific Northwest requires an average of at least 7 mol quanta m⁻² day⁻¹ during spring and summer and that long-term survival requires a minimum average of 3 mol quanta m⁻² day⁻¹.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Production of organic molecules in the outer solar system by proton irradiation: Laboratory simulations

In the past few years considerable attention has been given to the determination of likely compounds that could account for the various colors observed in the outer solar system: and to possible formation mechanisms for these compounds. Many experiments have been done using electrical discharges (Chadha, M. S., et al., 1971, Icarus15, 39) and ultraviolet light (Khare, B. N., and Sagan, C., 1973, Icarus20, 311) on mixtures of CH₄, NH₃, and H₂S, which are most likely the dominant minor constituents of the atmospheres of Jupiter, Saturn, Titan, and possibly the other satellites early in their histories. Colored polymers, usually brownish-red, have been produced in these experiments. With the passage of Pioneer 10 around Jupiter, there is another source of energy worthy of consideration, energetic protons (and electrons). Preliminary experiments to investigate the formation of colored polymers and other interesting molecules by the irradiation of gas mixtures by protons are discussed. Two to four Mev protons were used, with corresponding beam fluxes (as measured at 6R_J from the planet) equivalent to approximately 80 Earth years at Jupiter per hour of exposure. As in the other types of experiments, colored polymers have been produced. An important feature of this work is the presence or absence of absorption at 5 μm in the different materials produced; Titan is quite dark at this wavelength and Io is fairly bright. Such features may provide criteria for accepting or rejecting various materials produced in these experiments as reasonable coloring agents for the outer solar system.     

Point-connected rigid bodies in a topological tree

Identical equations of motion are shown to emerge for a system ofn+1 rigid bodies all interconnected byn points, each of which is common to two bodies, by means of each of the following derivation procedures, all of which employ a kinematical identity developed by Hooker and Margulies: The Hooker-Margulies/Hooker equations; Kane's quasicoordinate formulation of D'Alembert's principle; the combination of Lagrange's generalized coordinate equations and Lagrange's quasicoordinate equations; and the combination of Lagrange's generalized coordinate equations and the vector rotational equationM=H applied to the total system and resolved into a vector basis fixed in a reference body of the system. Thus the previously published Hooker-Margulies/Hooker equations are shown to be the natural result of several derivation procedures other than the Newton-Euler method originally used, provided that the central kinematical identity of the original derivation of Hooker and Margulies is employed.     

Incorporating solid solutions in reactive transport equations using a kinetic discrete-composition approach

A new approach is proposed for incorporating solid solution reactions into mass conservation equations describing reaction paths in both closed and open systems. The method is applicable to problems involving advective, dispersive, and diffusive transport in a porous medium. By representing the continuously variable solid solution composition with a discrete set of stoichiometric solids that span composition space, combined with a kinetic formulation of their rates of reaction, a self-determining spatial and temporal evolution of the solid solution concentration and composition is obtained. It is demonstrated that equilibrium of an aqueous solution with a stoichiometric solid derived from a solid solution corresponds to equilibrium of the solid solution itself if and only if equilibrium of the stoichiometric solid is stable. One advantage of this approach is that it is unnecessary to introduce any additional compositional variables to represent the solid solution. Discretization may be over the entire range of composition space, or over some subset depending on the system. A major consequence of the kinetic discrete-composition solid solution representation is that modeling solid solutions is similar to modeling pure mineral phases with the exception of a weighting factor applied to reaction rates of stoichiometric solids corresponding to a common solid solution. With this approach, precipitation leads to a discrete zonation of the solid solution that approximates the continuous variation in composition expected for the actual solid solution. The approach is demonstrated for a hypothetical ideal and non-ideal binary solid solution A_xB_1−xC for a reaction path formulation and reactive transport involving advection and diffusion.     

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).