A reappraisal of the 1835 eruption of Cosigüina and its atmospheric impact

In the great January 1835 eruption of Cosigüina volcano, Nicaragua, andesitic magma and lithic material were erupted over a period of at least three days. Proximal facies consist of clastogenic lava, scoria-fall, and lithic ash-fall produced by phreatomagmatic to vulcanian or plinian activity, together with surge deposits and lithic block-falls. Pyroclastic flow deposits covered some flanks of the volcano and entered the sea in the Gulf of Fonseca. Little record exists of the distal ash-fall, thus the total bulk volume erupted can only be roughly constrained to 2.9–5.6 km³. Furthermore, the amount of juvenile material is thought to be small. A recent study of volatiles in 1835 scoria suggests sulfur release from the magma was negligible. This reappraisal indicates that the Cosigüina eruption probably had little global climatic impact. Despite its violent nature, the magnitude of the eruption was modest. The eruption occurred too late to initiate the Northern Hemisphere cooling trend form 1828–1836. Dry fogs and other atmospheric optical phenomena usually observed after eruptions that contribute significantly to the stratospheric aerosol burden were not recorded after 1835.     

Symmetrical and segmented variation of physical and geochemical characteristics of the central american volcanic front

The regional variation of physical and geochemical characteristics of Central American volcanoes occurs in two fundamentally different patterns. The first pattern is symmetrical about Nicaragua. Crustal thickness, silica contents of mafic lavas and volcanic edifice heights are lowest in Nicaragua and increase smoothly toward Costa Rica to the south and Guatemala to the north. Magma density is maximum in Nicaragua and decreases smoothly outward. The regional variation in crustal thickness is just enough so that magma densities, calculated at appropriate Moho pressures, are the same at the base of the crust throughout the region. This is consistent with magma ponding at the base of the crust. The bulk compositions of Central American basalts show the same symmetrical variation. Suites of Nicaraguan basalts plotted in pseudo-ternary CMAS projections indicate large olivine and plagioclase primary-phase volumes. Toward Costa Rica and Guatemala the olivine and plagioclase fields inferred from suites of basaltic lavas are smaller, which is consistent with fractionation at increasing depth.The second pattern is the segmentation of the volcanic front and the plate margin in general. The segmentation strongly affects the spacing and size of volcanic centers. At segment boundaries volcanic centers are generally small and unusually widely spaced. Toward segment interiors volcano spacing and size increase systematically. The LIL element contents of lavas strongly reflect this pattern. For lavas with similar silica contents the larger the volcano, the higher the LIL element contents. The relationships between segmentation, volcano spacing and volcano size are compatible with diapiric rise of magma accumulated in narrow ribbons near the upper surface of the underthrust slab. The relationship between volcano volume and LIL element content is qualitatively in agreement with an open-system fractionation model.     

Decametric radio measurement of Jupiter's rotation period

A total of 26 measurements of Jupiter's 12-year average rotation period were made at frequencies of 18, 20, and 22.2 MHz at observatories in Florida and Chile. An improved method was employed in which histograms of occurrence probability vs central meridian longitude obtained at the same frequency and observatory during apparitions about 12 years (one Jovian year) apart were cross correlated. The longitude shift giving maximum cross correlation was used to correct the initially assumed rotation period value. The mean of the measurements is 9 hr 55 min 29.689 sec, with a standard deviation of the mean of 0.005 sec. This is about 0.02 sec, or 4 standard deviations, less than the System III (1965) value. The measurements indicate that the rotation period was not changing (linearly) at a rate in excess of 0.03 sec/yr. If the synoptic monitoring program is continued through the next maximum of the jovicentric declination of the Earth (D_E), we will probably be able to detect a rate of change in rotation period as small as 0.002 sec/yr. This accuracy might be sufficient to reveal a secular drift in Jupiter's magnetic field.     

Book Review: Geostatistics for Engineers and Earth Scientists,by Ricardo A. Olea. Kluwer Academic Publishers,Boston, 1999, 303 p., $135.00 (U.S.), hardcover,ISBN 0-7923-8523-3

Mathematical Geosciences -     

Weathering of Oil in a Surficial Aquifer

The composition of crude oil in a surficial aquifer was determined in two locations at the Bemidji, MN, spill site. The abundances of 71 individual hydrocarbons varied within 16 locations sampled. Little depletion of these hydrocarbons (relative to the pipeline oil) occurred in the first 10 years after the spill, whereas losses of 25% to 85% of the total measured hydrocarbons occurred after 30 years. The C_6‐30 n‐alkanes, toluene, and o‐xylene were the most depleted hydrocarbons. Some hydrocarbons, such as the n‐C_10–24 cyclohexanes, tri‐ and tetra‐ methylbenzenes, acyclic isoprenoids, and naphthalenes were the least depleted. Benzene was detected at every sampling location 30 years after the spill. Degradation of the oil led to increases in the percent organic carbon and in the δ ¹³C of the oil. Another method of determining hydrocarbon loss was by normalizing the total measured hydrocarbon concentrations to that of the most conservative analytes. This method indicated that the total measured hydrocarbons were depleted by 47% to 77% and loss of the oil mass over 30 years was 18% to 31%. Differences in hydrocarbon depletion were related to the depth of the oil in the aquifer, local topography, amount of recharge reaching the oil, availability of electron acceptors, and the presence of less permeable soils above the oil. The results from this study indicate that once crude oil has been in the subsurface for a number of years there is no longer a “starting oil concentration” that can be used to understand processes that affect its fate and the transport of hydrocarbons in groundwater.     

Hydraulic Conductivity Calibration of Logging NMR in a Granite Aquifer,Laramie Range,Wyoming

In granite aquifers, fractures can provide both storage volume and conduits for groundwater. Characterization of fracture hydraulic conductivity (K) in such aquifers is important for predicting flow rate and calibrating models. Nuclear magnetic resonance (NMR) well logging is a method to quickly obtain near-borehole hydraulic conductivity (i.e., K_NMR) at high-vertical resolution. On the other hand, FLUTe flexible liner technology can produce a K profile at comparable resolution but requires a fluid driving force between borehole and formation. For three boreholes completed in a fractured granite, we jointly interpreted logging NMR data and FLUTe K estimates to calibrate an empirical equation for translating borehole NMR data to K estimates. For over 90% of the depth intervals investigated from these boreholes, the estimated K_NMR are within one order of magnitude of K_FLUTe. The empirical parameters obtained from calibrating the NMR data suggest that “intermediate diffusion” and/or “slow diffusion” during the NMR relaxation time may occur in the flowing fractures when hydraulic aperture are sufficiently large. For each borehole, “intermediate diffusion” dominates the relaxation time, therefore assuming “fast diffusion” in the interpretation of NMR data from fractured rock may lead to inaccurate K_NMR estimates. We also compare calibrations using inexpensive slug tests that suggest reliable K_NMR estimates for fractured rock may be achieved using limited calibration against borehole hydraulic measurements.     

Petrology,geochronology,and geophysical characterization of Mesoproterozoic rocks in central Illinois,USA

The Precambrian basement rocks of the Eastern Granite-Rhyolite Province(EGRP)in central Illinois(midcontinent region of North America)exhibit a complex history of early volcanism,granite emplacement,and intrusion of mafic rocks.A comprehensive suite of dedicated petrographic analyses,geophysical logs,and drill core from four basement-penetrating wells,two-dimensional and three-dimensional seismic reflection data,and U-Pb age data from the Illinois Basin–Decatur Project(IBDP)and Illinois Carbon Capture Storage(ICCS)Project site provide new constraints for interpreting the Precambrian basement of the Illinois Basin.These new data reveal the basement to be compositionally and structurally complex,having typical EGRP felsic volcanic rocks intruded by the first reported gabbro in the Precambrian basement in Illinois.Zircons(n?29)from rhyolite give a U-Pb weighted mean average age of 14679 Ma.Zircons(n?3)from a gabbro dike that intrudes the rhyolite yield a concordia age of 107312 Ma,which corresponds to Grenville-age extension and represents the first Grenville-age rock in Illinois and in the EGRP.A high-resolution three-dimensional seismic reflection volume,coincident with the four wells,provides a context for interpreting the petrological data and implies a high degree of heterogeneity for basement rocks at the IBDP–ICCS site,as also shown by the drill cores.The occurrence of Grenville-age gabbro is related to a prominent bowl-like structure observed on local two-dimensional seismic reflection profiles and the three-dimensional volume that is interpreted as a deep-seated mafic sill complex.Furthermore,heterogeneities such as the brecciated EGRP rhyolite and later gabbro intrusion observed in the basement lithology at the IBDP–ICCS may reflect previously unknown distal elements of the 1.1 Ga Midcontinent Rift in the EGRP and more likely Grenville-age extension.     

Production regimes in four Eastern Boundary Current systems

High productivity (maxima 3 gCm⁻² day⁻¹) of the Eastern Boundary Currents (EBCs), i.e. the California, Peru-Humboldt, Canary and Benguela Currents, is driven by a combination of local forcing and large-scale circulation. The characteristics of the deep water brought to the surface by upwelling favorable winds depend on the large-scale circulation patterns. Here we use a new hydrographic and nutrient climatology together with satellite measurements of the wind vector, sea-surface temperature (SST), chlorophyll concentration, and primary production modeled from ocean color to quantify the meridional and seasonal patterns of upwelling dynamics and biological response. The unprecedented combination of data sets allows us to describe objectively the variability for small regions within each current and to characterize the governing factors for biological production. The temporal and spatial environmental variability was due in most regions to large-scale circulation, alone or in combination with offshore transport (local forcing). The observed meridional and seasonal patterns of biomass and primary production were most highly correlated to components representing large-scale circulation. The biomass sustained by a given nutrient concentration in the Atlantic EBCs was twice as large as that of the Pacific EBCs. This apparent greater efficiency may be due to availability of iron, physical retention, or differences in planktonic community structure.     

Developing a 150 ka luminescence chronology for the barrier dunes of the southern Cape, South Africa

The southern coastline of South Africa exhibits extensive aeolianites and spectacular barrier dunes; the ages and palaeoenvironmental significance of which have, until recently, received limited attention. We present preliminary results for samples taken from a single 6.5 m section in the seaward barrier dune near to the town of Wilderness. The reliability of these ages has been investigated using a series of dose recovery and pulsed irradiation experiments, as well as through analyses of single aliquot and single grain D_e distributions. Results indicate that the SAR protocol is providing accurate and precise D_e determinations and that adequate sensitivity correction is being made. The application of pulsed irradiations does not result in significant changes in the dose–response curves or averaged D_e's. Single grain D_e distributions are broad, with significant levels of overdispersion, which may reflect beta dose heterogeneity effects.

Peaks in aeolian activity date to 157–154 and 140–125 ka, with a slower rate of deposition recorded from 115 to 85 ka. Brief phases at 40 and 21 ka are also identified. Some ages, notably those from OIS 3 and OIS 2, are surprising from a geomorphological perspective, given the likely distance of the dune from the shoreline at such times, although they are not inconsistent with evidence from the east coast of South Africa. The results reveal multiple phases of Late Quaternary barrier dune accumulation, and suggest that aeolian activity in this region has been episodic and potentially more complex than previously thought.     

Re‐shaping hysteretic behaviour—spectral analysis and design equations for semi‐active structures

Semi‐active dampers offer significant capability to reduce dynamic wind and seismic structural response. A novel resetable device with independent valve control laws that enables semi‐active re‐shaping of the overall structural hysteretic behaviour has been recently developed, and a one‐fifth scale prototype experimentally validated. This research statistically analyses three methods of re‐shaping structural hysteretic dynamics in a performance‐based seismic design context. Displacement, structural force, and total base‐shear response reduction factor spectra are obtained for suites of ground motions from the SAC project. Results indicate that the reduction factors are suite invariant. Resisting all motion adds damping in all four quadrants and showed 40–60% reductions in the structural force and displacement at the cost of a 20–60% increase in total base‐shear. Resisting only motion away from equilibrium adds damping in quadrants 1 and 3, and provides reductions of 20–40%, with a 20–50% increase in total base‐shear. However, only resisting motion towards equilibrium adds damping in quadrants 2 and 4 only, for which the structural responses and total base‐shear are reduced 20–40%. The spectral analysis results are used to create empirical reduction factor equations suitable for use in performance based design methods, creating an avenue for designing these devices into structural applications. Overall, the reductions in both response and base‐shear indicate the potential appeal of this semi‐active hysteresis sculpting approach for seismic retrofit applications—largely due to the reduction of the structural force and overturning demands on the foundation system. Copyright © 2006 John Wiley & Sons, Ltd.