Phytoplankton production in two east coast estuaries: Photosynthesis-light functions and patterns of carbon assimilation in Chesapeake and Delaware Bays

Chesapeake Bay is a large and productive estuary that has received close scrutiny in recent years because of indications that its water quality and biota have been damaged by man's activities. Data on primary production for the estuary as a whole, however, are surprisingly sparse. We describe here the distribution of photosynthetic carbon assimilation by phytoplankton in Chesapeake Bay, and relate productivity patterns to hydrographic characteristics of the estuary. Between March 1982 and April 1983, a series of four cruises was conducted on Chesapeake Bay, and two cruises on the urbanized Delaware Bay for comparison. The upper Chesapeake and Delaware were highly turbid with high concentrations of suspended particulate matter and dissolved inorganic nutrients. Low chlorophyll concentrations were usually found in these areas of high turbidity, despite the abundance of nutrients, suggesting light limitation. Application of Wofsy's (1983) model of phytoplanton growth confirmed this suggestion. Chlorophyll and productivity maxima usually occurred seaward of the turbidity maxima where light penetration increased and suffient nutrients were present to support active phytoplankton growth. Further seaward of the chlorophyll maxima in the Chesapeake, the photic zone depth increased, concentrations of nutrients decreased, and phytoplankton biomass decreased, suggesting that nutrient availability, rather than light, controlled phytoplankton growth in the lower portion of the estuary. In contrast to the Chesapeake, Delaware Bay was more turbid, had generally higher nutrient concentrations, and was lower in phytoplankton productivity. The chlorophyll maxima and region of rapid phytoplankton growth occurred further toward the lower estuary and shelf regions in Delaware Bay because the high turbidity extended further seaward. Nutrients were never depleted at the shelf end of the estuary sufficiently to retard phytoplankton growth. Photosynthesis-irradiance (P-I) curves from simulated in situ and constant intensity incubations showed a strong correlation of the light-limited slope (a^B) with the light-saturated rate ( ) on each cruise. Spatial variations in corresponded to patterns of phytoplankton abundance, as did integral production (PP) and carbon-based growth rates (μ_C, μ_m), and photosynthetic parameters varied significantly with temperature.     

A calderas field at the Marifil Formation, new volcanogenic interpretation, Norpatagonian Massif, Argentina

Major volcanogenic structures within the Marifil Formation allowed to determine the location of a claderas field that spreads along 400 kilometers. This field has at least three large calderas with diameters near to 100 kilometers, recognizable by regional distribution of associated dikes and rhyolite lava flows present at the caldera border that separates a monotonous ignimbritic plateau from a depressed (150 meters) inner volcanic collapsed caldera. The Marifil ignimbritic plateau cover more than 50,000 square kilometers with thicknesses that reach 800 meters. Associated with these plateau ignimbrites there is an important fluorite mineralization.     

Quantification of karst aquifer discharge components during storm events through end-member mixing analysis using natural chemistry and stable isotopes as tracers

Karst aquifer components that contribute to the discharge of a water supply well in the Classical Karst (Kras) region (Italy/Slovenia) were quantitatively estimated during storm events. Results show that water released from storage within the epikarst may comprise as much as two-thirds of conduit flow in a karst aquifer following rainfall. Principal components analysis (PCA) and end-member mixing analysis (EMMA) were performed using major ion chemistry and the stable isotopes of water (δ¹⁸O, δ²H) and of dissolved inorganic carbon (δ¹³C_DIC) to estimate mixing proportions among three sources: (1) allogenic river recharge, (2) autogenic recharge, and (3) an anthropogenic component stored within the epikarst. The sinking river most influences the chemical composition of the water-supply well under low-flow conditions; however, this proportion changes rapidly during recharge events. Autogenic recharge water, released from shallow storage in the epikarst, displaces the river water and is observed at the well within hours after the onset of precipitation. The autogenic recharge end member is the second largest component of the well chemistry, and its contribution increases with higher flow. An anthropogenic component derived from epikarstic storage also impacts the well under conditions of elevated hydraulic head, accounting for the majority of the chemical response at the well during the wettest conditions.     

Late Quaternary Bison diminution on the Great Plains of North America: evaluating the role of human hunting versus climate change

Body size changes of Bison and mortality age structure data document the effects of climate-driven environmental change and human hunting pressure on large mammals in North America. Morphometric and mortality data are drawn from 58 archaeological and 9 paleontological localities dating between 37,000 and 250 calBP. Proxy information on body size is based on measurements recorded on 901 adult calcanei (os calcis) and 1026 humeri. In addition, published mortality profiles from 24 archaeological faunal assemblages spanning the last 14,000 years were used to approximate the age structure of bison populations. These results suggest that dramatic diminution in bison body size occurred in several short bursts, rather than a continuous gradual decline. These periods of rapid size reduction correlate with times of ecological reorganization, when aridity-driven changes in grasslands decreased forage quality and availability. Mortality age data indicate that the decrease in body size occurred in a context where there was no evidence for a progressively severe juvenile bias in bison populations. Overall, it appears that the changes in body size were a reaction to environmental conditions rather than the result of human predation pressure.     

Renewed uplift at the Yellowstone Caldera measured by leveling surveys and satellite radar interferometry

 A first-order leveling survey across the northeast part of the Yellowstone caldera in September 1998 showed that the central caldera floor near Le Hardy Rapids rose 24±5 mm relative to the caldera rim at Lake Butte since the previous survey in September 1995. Annual surveys along the same traverse from 1985 to 1995 tracked progressive subsidence near Le Hardy Rapids at an average rate of –19±1 mm/year. Earlier, less frequent surveys measured net uplift in the same area during 1923–1976 (14±1 mm/year) and 1976–1984 (22±1 mm/year). The resumption of uplift following a decade of subsidence was first detected by satellite synthetic aperture radar interferometry, which revealed approximately 15 mm of uplift in the vicinity of Le Hardy Rapids from July 1995 to June 1997. Radar interferograms show that the center of subsidence shifted from the Sour Creek resurgent dome in the northeast part of the caldera during August 1992 to June 1993 to the Mallard Lake resurgent dome in the southwest part during June 1993 to August 1995. Uplift began at the Sour Creek dome during August 1995 to September 1996 and spread to the Mallard Lake dome by June 1997. The rapidity of these changes and the spatial pattern of surface deformation suggest that ground movements are caused at least in part by accumulation and migration of fluids in two sill-like bodies at 5–10 km depth, near the interface between Yellowstone's magmatic and deep hydrothermal systems. Received: 30 November 1998 / Accepted: 16 April 1999     

A digital event-recording ocean bottom seismometer capsule

The ocean bottom seismometer capsule contains a 1 Hz. vertical seismometer and triggerable or programmable digital recording system. The output of the seismometer is continuously digitized at a preselected rate of 64, 128, or 256 samples/sec. The digital data words are mixed with a time code and synchronization characters, serialized and passed through a 1536 sample shift register which acts as a delay line. The serial output bits are then encoded and recorded on a SONY TC800B tape recorder which is turned on when a seismic event occurs. The event trigger occurs when the seismic signal jumps to 8 times the time averaged input signal. A memory may be programmed to run the recorder on a schedule so that small amplitude signals from refraction shots are sure to be recorded. Data are recovered using the same recorder for playback and a decoder which provides an analog output for field data interpretation or a digital output for computer analysis. An acoustic transponder allows precise ranges between the capsule and ship to be determined. In addition, commands for the capsule to release or to transmit diagnostic data may be given from the surface ship. The capsule falls freely to the ocean bottom. After a predetermined time or when a release command is received, it is released from a 68 kg steel tripod and floats to the surface. A dual timer and explosive bolt system is used to increase recovery reliability.The first capsules were designed and constructed between October 1972 and October 1973. Good results were obtained from 38 out of 43 launchings made on six expeditions in 1974, 1975, and 1976. Four capsules have been lost.     

Bond length and radii variations in fluoride and oxide molecules and crystals

Molecular orbital calculations completed on fluoride molecules containing first and second row cations have generated bond lengths, R, that match those observed for coordinated polyhedra in crystals to within ～0.04 Å, on average. The calculated bond lengths and those observed for fluoride crystals can be ranked with the expression R=Kp ^?0.22, where p=s/r, s is the Pauling strength of the bond, r is the row number of the cation and K=1.34. The exponent -0.22 (≈ -2/9) is the same as that observed for oxide, nitride and sulfide molecules and crystals. Bonded radii for the fluoride anion, obtained from theoretical electron density maps, increase linearly with bond length. Those calculated for the cations as well as for the fluoride anion match calculated promolecule radii to within ～0.03 Å, on average, suggesting that the electron density distributions in the vicinity of the minima along the bond paths possess a significant atomic component despite bond type. Bonded radii for Si and O ions provided by experimental electron density maps measured for the oxides coesite, danburite and stishovite match those calculated for a series of monosilicic acid molecules. The resulting radii increase with bond length and coordination number with the radius of the oxide ion increasing at a faster rate than that of the Si cation. The oxide ion within danburite exhibits several distinct radii, ranging between 0.9 and 1.2 Å, rather than a single radius with each exhibiting a different radius along each of the nonequivalent bonds with B, Si and Ca. Promolecule radii calculated for the coordinated polyhedra in danburite match procrystal radii obtained in a structure analysis to within 0.002 Å. The close agreement between these two sets of radii and experimentally determined bonded radii lends credence to Slater's statement that the difference between the electron density distribution observed for a crystal and that calculated for a procrystal (IAM) model of the crystal “would be small and subtle, and very hard to determine by examination of the total charge density.”     

Cross-shelf exchange in a model of the Ross Sea circulation and biogeochemistry

Transport of warm, nutrient-rich Circumpolar Deep Water (CDW) onto Antarctic continental shelves and coastal seas has important effects on physical and biological processes. The present study investigates the locations of this transport and its dynamics in the Ross Sea with a high-resolution three-dimensional numerical model. The model circulation is forced by daily wind stress along with heat and salt fluxes calculated from atmospheric climatologies by bulk formulae. All surface fluxes are modified by an imposed climatological ice cover. Waters under the Ross Ice Shelf are not included explicitly, but their effect on temperature and salinity is imposed in a buffer zone at the southern end of the model domain. A simple nutrient uptake is calculated based on the climatological chlorophyll distribution and Monod uptake kinetics.Model circulation is strongly affected by bottom topography, due to weak stratification, and agrees with schematics of the general flow and long-term current measurements except near the southern boundary. The sea-surface temperature is similar to satellite estimates except that the warmest simulated temperatures are slightly higher than observations. There is a significant correlation between the curvature of the shelf break and the transport across the shelf break. A momentum term balance shows that momentum advection helps to force flow across the shelf break in specific locations due to the curvature of the bathymetry (that is, the isobaths curve in front of the flow). For the model to create a strong intrusion of CDW onto the shelf, it appears two mechanisms are necessary. First, CDW is driven onto the shelf at least partially due to momentum advection and the curvature of the shelf break; then, the general circulation on the shelf takes the CDW into the interior.     

The moon’s permanent magnetic field: a cratered-shell model

As part of our study of the larger-scale remanent magnetic field of the Moon, we have examined the effects of cratering in an otherwise spherically symmetrical shell magnetized by a concentric dipolar magnetic fieldH _o to an intensity of magnetizationc H _o, wherec is a constant. In our initial model, we assume that the material excavated from the craters is distributed with random orientation and thus does not contribute to the remanent dipole momentM _g . We further assume that the mare fill does not contribute significantly toM _g . We choose the magnetizing dipole momentM _o and the constantc such that the magnitude of the productcH _o ≃ 3 × 10⁻⁴Г at the outer surface of the shell in the equatorial plane of the dipole. This value of the intensity of remanent magnetization was chosen to be within the range 10⁻⁷−10⁻³Г’; the intensities of thermo-remanent magnetization exhibited by Apollo samples. Finally, we use the locations and diameters of the 10 largest craters on the Moon and the depth-to-diameter ratios of Pike’s formulation to model approximately the excavation of the magnetized shell. The remanent dipole momentM _g was calculated for each of three orthogonal orientations of the magnetizing dipoleM _o. The three magnitudes ofM _g fall in the range 4 × 10¹⁸−1 × 10¹⁹Г cm³ which is close to the upper limit of 10¹⁹Г cm³ estimated forM _g from the field measurements obtained with the Apollo subsatellites. Further, the distribution of the craters is such as to produce a significant transverse component ofM _g with acute angles between the spin axis andM _g in the range 51°–77°.