Geology and eruptive history of the Rabaul Caldera area, Papua New Guinea

Rabaul Caldera is the most recently active (1937–1943) of four adjoining volcanic centres aligned north-south through the northern extremity of eastern New Britain. Geological mapping after the 1983–1985 Rabaul seismic and deformation crisis has partially revealed a long and complex eruption history dominated by numerous explosive eruptions, the largest accompanied by caldera collapse. The oldest exposed eruptives are the basaltic pre-caldera cone Tovanumbatir Lavas K/Ar dated at 0.5 Ma. The dacitic Rabaul Quarry Lavas exposed in the caldera wall and K/Ar dated at 0.19 Ma, are overlain by a sequence of dacitic and andesitic pyroclastic flow and fall deposits. Uplifted coral reef limestones, interbedded within the pyroclastic sequence on the northeast coast, suggest that explosive eruptions in the Rabaul area had commenced prior to the 0.125 Ma last interglacial high sea level stand. The pyroclastic sequence includes the large Boroi Ignimbrites and Malaguna Pyroclastics both ⁴⁰Ar/³⁹Ar dated at about 0.1 Ma, and the Barge Tunnel Ignimbrite ⁴⁰Ar/³⁹Ar dated at around 0.04 Ma. Few reliable ages exist for the many younger eruptives. These include Holocene ignimbrites of the latest caldera-forming eruptions—the Raluan Pyroclastics variously dated (¹⁴C) at either about 3500 or 7000 yr B.P., and the ca. 1400 yr B.P. Rabaul Pyroclastics. At least eight intracaldera eruptions have occurred since the 1400 yr B.P. collapse, building small pyroclastic and lava cones within the caldera.A major erosional episode is evident as a widespread unconformity in the upper pyroclastic stratigraphy at Rabaul. Lacking relevant radiometric ages, this episode is assumed to have occurred during last glaciation low sea levels and is here arbitarily dated at ca. ?20 ka. At least five, possibly nine, significant ignimbrite eruptions have occurred at Rabaul during the last ?20 ka. The new eruptive history differs considerably from that previously published, which considered ignimbrite eruption and caldera collapse to have first occurred at 3500 yr B.P.Rabaul volcanism has been dominated by two main types: (a) basaltic and basaltic andesite cone building eruptions; and (b) dacitic, and rarely andesitic or rhyolitic, plinian/ignimbrite eruptions of both high- and low-aspect ratio types. The 1400 yr B.P. Rabaul Ignimbrite is a type example of a low-aspect ratio, high-energy, and potentially very damaging eruption. Fine vitric ash deposits, common in the Rabaul pyroclastic sequence, demonstrate the frequent modification of eruptions by external water probably related to early caldera lakes or bays. Interbedding of these fine ashes with plinian pumice lapilli beds suggests that many early eruptions occurred from multiple vents, located in both wet and dry areas.     

Fluid convection and mass transfer in porous sandstones—a theoretical model

When porous, fluid-saturated, bodies of geologic dimensions are subjected to temperature gradients on the order of the normal geothermal gradient (25°C/km) it is shown that large (km) scale eddy currents will spontaneously arise and persist in these bodies. The velocity of these fluid currents is on the order of 1 meter per year, but taken over a period of several million years, the mass flux accompanying the fluid flow is large enough to produce significant porosity changes. It is shown that diffusive mass transfer is generally a negligible component of the total mass flux when characteristic distances are greater than grain diameters and it is concluded that most post-depositional reservoir cementation is due to slowly circulating fluids. In terms of a broad perspective, reservoir diagenesis can be regarded as a low heat flow analogue of hydrothermal alteration, occupying the low fluid flow regime of the spectrum, and passing into the hydrothermal regime as the heat flow and fluid velocities increase.     

Reaction relationships of chrome-spinels in igneous rocks — further evidence from the layered intrusions of Rhum and Mull,Inner Hebrides,Scotland

Compositions of chrome-spinels in peridotites from the layered igneous rocks of Rhum and of the Ben Buie intrusion, Mull, show evidence of post-depositional reaction relationships with the cumulus olivine and/or the intercumulus liquid. Some seam-forming spinels in the Ben Buie intrusion are highly aluminous, more so than those in the Rhum intrusion. Occasional zoned spinels, in both intrusions, show enrichment in Al towards their rims. The evidence suggests that the spinel reaction trend involving Cr-Al exchange proceeds in the direction of Al-enrichment. Equilibrium between chrome-spinel and cumulus olivine appears to have been attained rarely, thus limiting the application of the olivine-spinel geothermometer.     

Electronic structures of sulfide minerals — Theory and experiment

The sulfide minerals exhibit a rich diversity in sturctural chemistry and in electrical, magnetic and other physical properties. Models based on molecular orbital theory and incorporating some elements of band theory can be developed to describe the diverse valence electron behavior in these minerals. Qualitative models can be proposed on the basis of observed properties, and the models can be tested and refined using experimental data from X-ray emission and X-ray photoelectron spectroscopy and quantum mechanical calculations performed on cluster units which form the basic building blocks of the crystals. This approach to chemical bonding in sulfide minerals is illustrated for binary non-transition metal sulfides (ZnS, CdS, HgS, PbS), binary transition metal sulfides (FeS₂, CoS₂, NiS₂, CuS₂ ZnS₂) and more complex sulfides (CuFeS₂, Cu₂S, Ag₂S, CuS, Co₃S₄, CuCo₂S₄, Fe₃S₄). The relationship between qualitative and quantitative theories is reviewed with reference to the pyrite-marcasite-arsenopyrite-loellingite series of minerals. Application of the models to understanding structure-determining principles, relative stabilities, solid solution limits and properties such as color, reflectance and hardness are discussed.     

Detection Time for Plausible Changes in Annual Precipitation, Evapotranspiration, and Streamflow in Three Mississippi River Sub-Basins

We use diagnostic studies of off-line variable infiltration capacity (VIC) model simulations of terrestrial water budgets and 21st-century climate change simulations using the parallel climate model (PCM) to estimate the time required to detect predicted changes in annual precipitation (P), evapotranspiration (E), and discharge (Q) in three sub-basins of the Mississippi River Basin. Time series lengths on the order of 50–350 years are required to detect plausible P, E, and Q trends in the Missouri, Ohio, and Upper Mississippi River basins. Approximately 80–160, 50, and 140–350 years, respectively, are needed to detect the predicted P, E, and Q trends with a high degree of statistical confidence. These detection time estimates are based on conservative statistical criteria (α = 0.05 and β = 0.10) associated with low probability of both detecting a trend when it is not occurring (Type I error) and not detecting a trend when it is occurring (Type II error). The long detection times suggest that global-warming-induced changes in annual basin-wide hydro-climatic variables that may already be occurring in the three basins probably cannot yet be detected at this level of confidence. Furthermore, changes for some variables that may occur within the 21st century might not be detectable for many decades or until the following century – this may or may not be the case for individual recording station data. The long detection times for streamflow result from comparatively low signal-to-noise ratios in the annual time series. Finally, initial estimates suggest that faster detection of acceleration in the hydrological cycle may be possible using seasonal time series of appropriate hydro-climatic variables, rather than annual time series.     

Solute and isotope constraint of groundwater recharge simulation in an arid environment,Abu Dhabi Emirate,United Arab Emirates

A simple, well-constrained simulation of solute increases in a downgradient direction was used in the shallow unconfined aquifer of eastern Abu Dhabi Emirate (United Arab Emirates). The simulation indicates that the observed exponential increase in solute concentrations results from a combination of upward transport of solutes from underlying mudstones and evaporites, and groundwater losses by evaporation. Groundwater recharge and discharge flux in unconfined regional aquifers in arid regions commonly are difficult to estimate because there are few constraints on the flux of water lost or gained from the system. Total dissolved solids (TDS) and deuterium isotopes (δ²H) in groundwater are used to constrain estimated fluxes to the shallow aquifer of eastern Abu Dhabi Emirate. Vertical upward transport of solutes from underlying mudstones and evaporites accounts for solute increases along approximately the first 80 km of the simulated flow path, but a combination of upward solute transport and evaporation is necessary to explain observed solute concentrations beyond 80 km. Mobilization and transport of solutes in the unsaturated zone by recharging precipitation is not a significant factor.

---

Résumé Une simulation simple et bien contrainte de l’augmentation des solutés dans la direction de l’écoulement a été réalisée sur l’aquifère phréatique et libre de l’Est de l’émirat d’Abu Dhabi (émirats Arabes Unis). La simulation indique que l’augmentation observée de solutés, de forme exponentielle, est le produit de la combinaison du transport ascendant de solutés provenant des mudstones et des évaporites sous-jacents, et de l’évaporation des eaux souterraines. La recharge des eaux souterraines et les flux de vidange dans les aquifères régionaux libres des régions arides sont habituellement difficiles à estimer du fait du peu de contraintes agissant sur la perte ou le gain d’eau dans le système. La charge dissoute totale (TDS, en anglais) et les isotopes deutérium (δ²H) des eaux souterraines sont utilisés pour contraindre l’estimation des flux en direction de l’aquifère phréatique de l’Est de l’émirat d’Abu Dhabi. Le transport vertical ascendant des solutés à partir des mudstones et des évaporites fait parti des processus qui concentrent les solutés approximativement le long des 80 premiers kilomètres du trajet des écoulements simulés, mais une combinaison du transport ascendant des solutés et de l’évaporation est nécessaire pour expliquer les concentrations observées en soluté au delà des 80 km. La mobilisation et le transport des solutés par les précipitations efficaces dans la zone non saturée ne représentent pas des facteurs significatifs.

---

Resumen Se utilizó una simulación simple y muy restringida, de incrementos de soluto en la dirección del flujo, en el acuífero libre poco profundo del Emirato oriental de Abu Dhabi, (Emiratos árabes Unidos). La simulación indica que el aumento exponencial observado en las concentraciones del soluto, resulta de una combinación de transporte ascendente de solutos desde las evaporitas y lodolitas subyacentes, y por pérdidas del agua subterránea por evaporación. Los flujos de recarga y descarga de agua subterránea, en los acuíferos libres regionales en regiones áridas, normalmente son difíciles estimar, porque hay algunas restricciones en el flujo de agua perdido o ganado por el sistema. Los sólidos disueltos totales (TDS) y los isótopos de deuterio (δ²H) en el agua subterránea, se usan para forzar los flujos estimados hacia el acuífero poco profundo del Emirato oriental de Abu Dhabi. El transporte ascendente vertical de solutos desde las lodolitas y evaporitas subyacentes, involucra los aumentos del soluto a lo largo de aproximadamente los primeros 80 Km. de la dirección de flujo simulada, pero es necesaria una combinación de transporte ascendente del soluto y la evaporación para explicar las concentraciones del soluto observadas más allá de 80 Km. La movilización y transporte de solutos en la zona no saturada, por recarga debida a precipitación, no es un factor significante.

     

Large-Eddy Simulation Of Neutral Turbulent Flow Over Rough Sinusoidal Ridges

The requirements for a credible large-eddy simulation of neutral, turbulent flow over hills with an aerodynamically rough surface are discussed, in order to select a suitable case for simulation. As well as providing adequate resolution within the dynamically important inner region, obtaining a realistic upstream or undisturbed mean velocity profile is also of critical importance. A distributed drag canopy formulation has been introduced to the model to allow it to obtain such a profile while resolving very close to the rough surface. Simulations have then been performed of flow over ridges of varying heights. The results from the steepest case, which is just on the verge of separation, are compared with wind-tunnel observations. It is shown that the large-eddy simulation results are in much better agreement with the experimental data than are the results from a simple first-order mixing-length closure model. An encouraging lack of sensitivity of the simulation results to numerical resolution is also demonstrated.     

Accuracy of ozonesonde measurements in the troposphere

Measurements of atmospheric ozone at concentrations typical of the free troposphere have been compared for ECC sondes and a UV-absorption photometer, using a Bendix chemiluminescent analyser as a transfer standard. Comparisons were conducted in the laboratory and in the tropospheric part of the atmosphere. It was found that the measurements agreed to within 4% provided that the background current for the ECC sonde was measured before exposure to ozone in the preparation procedure, and was assumed to be constant throughout the sonde flight. These results confirm those of earlier experiments and mean that the methods currently used to correct for the background current in the troposphere need to be revised.     

A comparative allometric study of the morphometry of the gills of an alkalinity adapted cichlid fish,Oreochromis alcalicus grahami,of Lake Magadi,Kenya

A morphometric analysis of the gills ofOreochromis alcalicus grahami has been carried out on specimens from ecologically distinct lagoons and a water-holding tank of Lake Magadi, a highly alkaline salt lake situated in a volcanically active region of the southern part of the Great Rift Valley in Kenya. The data were compared with those fromOreochromis niloticus, a close relative that lives in fresh water and with data from other fresh water and marine fish. Our primary goal was to identify the possible adaptive features which enable the fish to survive in an environment characterized by severely fluctuating levels of oxygen, a condition exacerbated by factors such as high temperature, alkalinity and osmolarity. The specimens ofO. a. grahami from the south-western lagoons of the lake had gills better adapted for gas exchange with a body mass specific diffusing capacity for oxygen which was about 2 times greater than that of the gills of the specimens from the fish spring lagoons and 2.5 times that of those from the water-holding tanks. Some parameters of the gills ofO. a. grahami, e.g. the gill filament length and number of gill filaments are significantly greater than those ofO. niloticus but the number of secondary lamellae, area of secondary lamellae and the diffusing capacity of the gills are similar in the two species. Compared with most other fish, the gills ofO. a. grahami appear to be particularly well adapted for gas exchange especially by having a thin water-blood barrier. Perhaps in no other extant fish have the gills had to be so exquisitely designed to meet environmental extremes and regulate complex and at times conflicting functions such as gas exchange, iono-regulation, acid-base balance and nitrogenous waste excretion as inO. a. grahami