Geochemistry of coral from Papua New Guinea as a proxy for ENSO ocean–atmosphere interactions in the Pacific Warm Pool

A Porites sp. coral growing offshore from the Sepik and Ramu Rivers in equatorial northern Papua New Guinea has yielded an accurate 20-year history (1977–1996) of sea surface temperature (SST), river discharge, and wind-induced mixing of the upper water column. Depressions in average SSTs of about 0.5–1.0 °C (indicated by coral Sr/Ca) and markedly diminished freshwater runoff to the coastal ocean (indicated by coral δ¹⁸O, δ¹³C and UV fluorescence) are evident during the El Niño – Southern Oscillation (ENSO) events of 1982–1983, 1987 and 1991-1993. The perturbations recorded by the coral are in good agreement with changes in instrumental SST and river discharge/precipitation records, which are known to be diagnostic of the response of the Pacific Warm Pool ocean–atmosphere system to El Niño. Consideration of coastal ocean dynamics indicates that the establishment of northwest monsoon winds promotes mixing of near-surface waters to greater depths in the first quarter of most years, making the coral record sensitive to changes in the Asian–Australian monsoon cycle. Sudden cooling of SSTs by 1°C following westerly wind episodes, as indicated by the coral Sr/Ca, is consistent with greater mixing in the upper water column at these times. Furthermore, the coral UV fluorescence and oxygen isotope data indicate minimal contribution of river runoff to surface ocean waters at the beginning of most years, during the time of maximum discharge. This abrupt shift in flood-plume behaviour appears to reflect the duration and magnitude of northwest monsoon winds, which tend to disperse flood plume waters to a greater extent in the water column when wind-mixing is enhanced. Our results suggest that a multi-proxy geochemical approach to the production of long coral records should provide comprehensive reconstructions of tropical paleoclimate processes operating on interannual timescales.     

Pyroclastic Flow Hazard at Volcán Citlaltépetl

Volcán Citlaltépetl (Pico de Orizaba) with an elevation of 5,675 m is the highest volcano in North America. Its most recent catastrophic events involved the production of pyroclastic flows that erupted approximately 4,000, 8,500, and 13,000 years ago. The distribution of mapped deposits from these eruptions gives an approximate guide to the extent of products from potential future eruptions. Because the topography of this volcano is constantly changing computer simulations were made on the present topography using three computer algorithms: energy cone, FLOW2D, and FLOW3D. The Heim Coefficient (), used as a code parameter for frictional sliding in all our algorithms, is the ratio of the assumed drop in elevation (H) divided by the lateral extent of the mapped deposits (L). The viscosity parameter for the FLOW2D and FLOW3D codes was adjusted so that the paths of the flows mimicked those inferred from the mapped deposits. We modeled two categories of pyroclastic flows modeled for the level I and level II events. Level I pyroclastic flows correspond to small but more frequent block-and-ash flows that remain on the main cone. Level II flows correspond to more widespread flows from catastrophic eruptions with an approximate 4,000-year repose period. We developed hazard maps from simulations based on a National Imagery and Mapping Agency (NIMA) DTED-1 DEM with a 90 m grid and a vertical accuracy of ±30 m. Because realistic visualization is an important aid to understanding the risks related to volcanic hazards we present the DEM as modeled by FLOW3D. The model shows that the pyroclastic flows extend for much greater distances to the east of the volcano summit where the topographic relief is nearly 4,300 m. This study was used to plot hazard zones for pyroclastic flows in the official hazard map that was published recently.     

Building and Testing an <Emphasis Type="Italic">a priori</Emphasis> Geophysical Model for Western Eurasia and North Africa

We construct and evaluate a new three-dimensional model of crust and upper mantle structure in Western Eurasia and North Africa (WENA) extending to 700 km depth and having 1° parameterization. The model is compiled in an a priori fashion entirely from existing geophysical literature, specifically, combining two regionalized crustal models with a high-resolution global sediment model and a global upper mantle model. The resulting WENA1.0 model consists of 24 layers: water, three sediment layers, upper, middle, and lower crust, uppermost mantle, and 16 additional upper mantle layers. Each of the layers is specified by its depth, compressional and shear velocity, density, and attenuation (quality factors, Q _P and Q _S ). The model is tested by comparing the model predictions with geophysical observations including: crustal thickness, surface wave group and phase velocities, upper mantle _n velocities, receiver functions, P-wave travel times, waveform characteristics, regional 1-D velocities, and Bouguer gravity. We find generally good agreement between WENA1.0 model predictions and empirical observations for a wide variety of independent data sets. We believe this model is representative of our current knowledge of crust and upper mantle structure in the WENA region and can successfully be used to model the propagation characteristics of regional seismic waveform data. The WENA1.0 model will continue to evolve as new data are incorporated into future validations and any new deficiencies in the model are identified. Eventually this a priori model will serve as the initial starting model for a multiple data set tomographic inversion for structure of the Eurasian continent.     

Amoeboid olivine aggregates and related objects in carbonaceous chondrites: records of nebular and asteroid processes

Amoeboid olivine aggregates (AOAs) are the most common type of refractory inclusions in CM, CR, CH, CV, CO, and ungrouped carbonaceous chondrites Acfer 094 and Adelaide; only one AOA was found in the CB_b chondrite Hammadah al Hamra 237 and none were observed in the CB_a chondrites Bencubbin, Gujba, and Weatherford. In primitive (unaltered and unmetamorphosed) carbonaceous chondrites, AOAs consist of forsterite (Fa_<2), Fe, Ni-metal (5-12 wt% Ni), and Ca, Al-rich inclusions (CAIs) composed of Al-diopside, spinel, anorthite, and very rare melilite. Melilite is typically replaced by a fine-grained mixture of spinel, Al-diopside, and ±anorthite; spinel is replaced by anorthite. About 10% of AOAs contain low-Ca pyroxene replacing forsterite. Forsterite and spinel are always ¹⁶O-rich (δ^17,18O∼−40‰ to −50‰), whereas melilite, anorthite, and diopside could be either similarly ¹⁶O-rich or ¹⁶O-depleted to varying degrees; the latter is common in AOAs from altered and metamorphosed carbonaceous chondrites such as some CVs and COs. Low-Ca pyroxene is either ¹⁶O-rich (δ^17,18O∼−40‰) or ¹⁶O-poor (δ^17,18O∼0‰). Most AOAs in CV chondrites have unfractionated (∼2-10×CI) rare-earth element patterns. AOAs have similar textures, mineralogy and oxygen isotopic compositions to those of forsterite-rich accretionary rims surrounding different types of CAIs (compact and fluffy Type A, Type B, and fine-grained, spinel-rich) in CV and CR chondrites. AOAs in primitive carbonaceous chondrites show no evidence for alteration and thermal metamorphism. Secondary minerals in AOAs from CR, CM, and CO, and CV chondrites are similar to those in chondrules, CAIs, and matrices of their host meteorites and include phyllosilicates, magnetite, carbonates, nepheline, sodalite, grossular, wollastonite, hedenbergite, andradite, and ferrous olivine.Our observations and a thermodynamic analysis suggest that AOAs and forsterite-rich accretionary rims formed in ¹⁶O-rich gaseous reservoirs, probably in the CAI-forming region(s), as aggregates of solar nebular condensates originally composed of forsterite, Fe, Ni-metal, and CAIs. Some of the CAIs were melted prior to aggregation into AOAs and experienced formation of Wark-Lovering rims. Before and possibly after the aggregation, melilite and spinel in CAIs reacted with SiO and Mg of the solar nebula gas enriched in ¹⁶O to form Al-diopside and anorthite. Forsterite in some AOAs reacted with ¹⁶O-enriched SiO gas to form low-Ca pyroxene. Some other AOAs were either reheated in ¹⁶O-poor gaseous reservoirs or coated by ¹⁶O-depleted pyroxene-rich dust and melted to varying degrees, possibly during chondrule formation. The most extensively melted AOAs experienced oxygen isotope exchange with ¹⁶O-poor nebular gas and may have been transformed into magnesian (Type I) chondrules. Secondary mineralization and at least some of the oxygen isotope exchange in AOAs from altered and metamorphosed chondrites must have resulted from alteration in the presence of aqueous solutions after aggregation and lithification of the chondrite parent asteroids.     

Neoproterozoic and Cambro-Ordovician magmatism in the Variscan Kłodzko Metamorphic Complex (West Sudetes,Poland): new insights from U/Pb zircon dating

The Kodzko Metamorphic Complex (KMC) in the Central Sudetes consists of meta-sedimentary and meta-igneous rocks metamorphosed under greenschist to amphibolite facies conditions. They are comprised in a number of separate tectonic units interpreted as thrust sheets. In contrast to other Lower Palaeozoic volcano-sedimentary successions in the Sudetes, the two uppermost units (the Orla-Googowy unit and the Kodzko Fortress unit) of the KMC contain meta-igneous rocks with supra-subduction zone affinities. The age of the KMC was previously assumed to be Early Palaeozoic–Devonian, based on biostratigraphic findings in the lowermost tectonic unit. Our geochronological study focused on the magmatic rocks from the two uppermost tectonic units, exposed in the SW part of the KMC. Two orthogneiss samples from the Orla-Googowy unit yielded ages of 500.4±3.1 and 500.2±4.9 Ma, interpreted to indicate the crystallization age of the granitic precursors. A plagioclase gneiss from the same tectonic unit, intimately interlayered with metagabbro, provided an upper intercept age of 590.1±7.2 Ma, which is interpreted as the time of igneous crystallization. From the topmost Kodzko Fortress unit, a metatuffite was studied, which contains a mixture of genetically different zircon grains. The youngest ²⁰⁷Pb/²⁰⁶Pb ages, which cluster at ca. 590-600 Ma, are interpreted to indicate the maximum depositional age for this metasediment. The results of this study are in accord with a model that suggests a nappe structure for the KMC, with a Middle Devonian succession at the base and Upper Proterozoic units at structurally higher levels. It is suggested here that the KMC represents a composite tectonic suture that juxtaposes elements of pre-Variscan basement, intruded by the Lower Ordovician granite, against a Middle Palaeozoic passive margin succession. The new ages, combined with the overall geochemical variation in the KMC, indicate the existence of rock assemblages representing a Gondwana active margin. The recognition of Neoproterozoic subduction-related magmatism provides additional arguments for the hypothesis that equivalents of the Teplá-Barrandian domain are exposed in the Central Sudetes.     

Picophytoplankton and bacterioplankton in the Mississippi River plume and its adjacent waters

Picoplankton abundance and distribution in the Mississippi River plume and its adjacent waters were studied during two cruises in April (high discharge) and October (low discharge) 2000 using flow cytometry. Concentrations of photosynthetic picoplankton,Synechococcus and picoeukaryotes were low in the turbid plume water but high in the coastal waters—i.e., the green waters resulting from mixing of river and oceanic waters. In this region, three types ofSynechococcus, characterized by their phycoerythrin chromophore composition, were found:Synechococcus cells with a low phycourobilin to phycoerythrobilin ratio (PUB:PEB) occurred throughout the region and dominated the totalSymechococcus abundance during both seasons; high PUB:PEB cells, which are the dominant strains in the open or blue ocean, occurred only at the outer shelf stations; and PEB-onlySynechococcus were abundant in most of the surveyed area during april, but were not observed during October.Prochlorococcus cyanobacteria only occurred at the oceanic stations, but extended farther inshore in October compared to April. This was a consequence of the reduced discharge and plume size during October. Picophytoplankton were a less important component of total phytoplankton biomass in the turbid river water and more important in the oligotrophic Gulf water. Seasonally, the contribution of picophytoplankton to total phytoplankton biomass in the surveyed area was higher during low discharge in October than during high discharge in April, even though the spring 2000 river discharge was unusually low and might not present a typical high discharge scenario. The abundance of heterotrophic bacteria was weakly correlated to chlorophylla (chla) concentration, but better correlated to picophytoplankton biomass. A higher proportion of High DNA bacteria occurred in the river-impacted regions during both seasons, with the ratio of High DNA bacteria to Low DNA bacteria significantly higher in April.     

Using analytic element models to delineate drinking water source protection areas

Since 1999, Ohio EPA hydrogeologists have used two analytic element models (AEMs), the proprietary software GFLOW and U.S. EPA's WhAEM, to delineate protection areas for 535 public water systems. Both models now use the GFLOW2001 solution engine, integrate well with Geographic Information System (GIS) technology, have a user-friendly graphical interface, are capable of simulating a variety of complex hydrogeologic settings, and do not rely upon a model grid. These features simplify the modeling process and enable AEMs to bridge the gap between existing simplistic delineation methods and more complex numerical models. Ohio EPA hydrogeologists demonstrated that WhAEM2000 and GFLOW2000 were capable of producing capture zones similar to more widely accepted models by applying the AEMs to eight sites that had been previously delineated using other methods. After the Ohio EPA delineated protection areas using AEMs, more simplistic delineation methods used by other states (volumetric equation and arbitrary fixed radii) were applied to the same water systems to compare the differences between various methods. GIS software and two-tailed paired t-tests were used to quantify the differences in protection areas and analyze the data. The results of this analysis demonstrate that AEMs typically produce significantly different protection areas than the most simplistic delineation methods, in terms of total area and shape. If the volumetric equation had been used instead of AEMs, Ohio would not have protected 265 km2 of critical upgradient area and would have overprotected 269 km2 of primarily downgradient land. Since an increasing number of land-use restrictions are being tied to drinking water protection areas, this analysis has broad policy implications.     

Conventional and combined pump-and-treat systems under nonuniform background flow

We investigate the performance of vertical hydraulic barriers in combination with extraction wells for the partial hydraulic isolation of contaminated aquifer areas. The potential advantage of such combinations compared to a conventional pump-and-treat system has already been demonstrated in a previous study. Here we extend the scope of the performance analysis to the impact of uncertainty in the regional flow direction as well as to highly heterogeneous aquifer transmissivity distributions. In addition, two new well-barrier scenarios are proposed and analyzed. The hydraulic efficiency of the scenarios is rated based on the expected (mean) reduction of the pumping rate that is required to achieve downgradient contaminant capture. The uncertain spatial distribution of aquifer transmissivity is considered by means of unconditioned Monte Carlo simulations. The significance of uncertain background flow conditions is incorporated by computing minimized pumping rates for deviations of the regional flow direction up to 30 degrees from a normative base case. The results give an answer on how pumping rates have to be changed for each barrier-well combination in order to achieve robust systems. It is exposed that in comparison to installing exclusively wells, the barrier-supported approach generally yields savings in the (average) pumping rate. The particular efficiency is shown to be highly dependent on the interaction of variance and integral scale of transmissivity distribution, well and barrier position, as well as direction of background flow.