Link between SSZ ophiolite formation, emplacement and arc inception, Northland, New Zealand: U–Pb SHRIMP constraints; Cenozoic SW Pacific tectonic implications

New U–Pb age-data from zircons separated from a Northland ophiolite gabbro yield a mean ²⁰⁶Pb/²³⁸U age of 31.6 ± 0.2 Ma, providing support for a recently determined 28.3 ± 0.2 Ma SHRIMP age of an associated plagiogranite and  29–26 Ma ⁴⁰Ar/³⁹Ar ages (n = 9) of basalts of the ophiolite. Elsewhere, Miocene arc-related calc-alkaline andesite dikes which intrude the ophiolitic rocks contain zircons which yield mean ²⁰⁶Pb/²³⁸U ages of 20.1 ± 0.2 and 19.8 ± 0.2 Ma. The ophiolite gabbro and the andesites both contain rare inherited zircons ranging from 122–104 Ma. The Early Cretaceous zircons in the arc andesites are interpreted as xenocrysts from the Mt. Camel basement terrane through which magmas of the Northland Miocene arc lavas erupted. The inherited zircons in the ophiolite gabbros suggest that a small fraction of this basement was introduced into the suboceanic mantle by subduction and mixed with mantle melts during ophiolite formation.

We postulate that the tholeiitic suite of the ophiolite represents the crustal segment of SSZ lithosphere (SSZL) generated in the southern South Fiji Basin (SFB) at a northeast-dipping subduction zone that was initiated at about 35 Ma. The subduction zone nucleated along a pre-existing transform boundary separating circa 45–20 Ma oceanic lithosphere to the north and west of the Northland Peninsula from nascent back arc basin lithosphere of the SFB. Construction of the SSZL propagated southward along the transform boundary as the SFB continued to unzip to the southeast. After subduction of a large portion of oceanic lithosphere by about 26 Ma and collision of the SSZL with New Zealand, compression between the Australian Plate and the Pacific Plate was taken up along a new southwest-dipping subduction zone behind the SSZL. Renewed volcanism began in the oceanic forearc at 25 Ma producing boninitic-like, SSZ and within-plate alkalic and calc-alkaline rocks. Rocks of these types temporally overlap ophiolite emplacement and subsequent Miocene continental arc construction.     

Geochemistry and magmatic properties of eruption episodes from Haroharo linear vent zone, Okataina Volcanic Centre, New Zealand during the last 10 kyr

Post-10 ka rhyolitic eruptions from the Haroharo linear vent zone, Okataina Volcanic Centre, have occurred from several simultaneously active vents spread over 12 km. Two of the three eruption episodes have tapped multiple compositionally distinct homogeneous magma batches. Three magmas totalling ~8 km³ were erupted during the 9.5 ka Rotoma episode. The most evolved Rotoma magma (SiO₂=76.5–77.9 wt%, Sr=96–112 ppm) erupted from a southeastern vent, and is characterised by a cummingtonite-dominant mineralogy, a temperature of 739±14°C, and fO₂ of NNO+0.52±0.11. The least evolved (SiO₂=75.0–76.4 wt%, Sr=128–138 ppm, orthopyroxene+ hornblende-dominant) Rotoma magma erupted from several vents, and was hotter (764±18°C) and more reduced (NNO+0.40±0.13). The ~11 km³ Whakatane episode occurred at 5.6 ka and also erupted three magmas, each from a separate vent. The most evolved (SiO₂=73.3–76.2 wt%, Sr=88–100 ppm) Whakatane magma erupted from the southwestern (Makatiti) vent and is cummingtonite-dominant, cool (745±11°C), and reduced (NNO+0.34±0.08). The least evolved (SiO₂=72.8–74.1 wt%, Sr=132–134 ppm) magma was erupted from the northeastern (Pararoa) vent and is characterised by an orthopyroxene+ hornblende-dominant mineralogy, temperature of 764±18°C, and fO₂ of NNO+0.40±0.13. Compositionally intermediate magmas were erupted during the Rotoma and Whakatane episodes are likely to be hybrids. A single ~13 km³ magma erupted during the intervening 8.1 ka Mamaku episode was relatively homogeneous in composition (SiO₂=76.1–76.8 wt%, Sr=104–112 ppm), temperature (736±18°C), and oxygen fugacity (NNO+0.19±0.12). Some of the vents tapped a single magma while others tapped several. Deposit stratigraphy suggests that the eruptions alternated between magmas, which were often simultaneously erupted from separate vents. Both effusive and explosive activity alternated, but was predominantly effusive (>75% erupted as lava domes and flows). The plumbing systems which fed the vents are inferred to be complex, with magma experiencing different conditions in the conduits. As the eruption of several magmas was essentially concurrent, the episodes were likely triggered by a common event such as magmatic intrusion or seismic disturbance.     

Uncertainties in projected impacts of climate change on European agriculture and terrestrial ecosystems based on scenarios from regional climate models

The uncertainties and sources of variation in projected impacts of climate change on agriculture and terrestrial ecosystems depend not only on the emission scenarios and climate models used for projecting future climates, but also on the impact models used, and the local soil and climatic conditions of the managed or unmanaged ecosystems under study. We addressed these uncertainties by applying different impact models at site, regional and continental scales, and by separating the variation in simulated relative changes in ecosystem performance into the different sources of uncertainty and variation using analyses of variance. The crop and ecosystem models used output from a range of global and regional climate models (GCMs and RCMs) projecting climate change over Europe between 1961–1990 and 2071–2100 under the IPCC SRES scenarios. The projected impacts on productivity of crops and ecosystems included the direct effects of increased CO₂ concentration on photosynthesis. The variation in simulated results attributed to differences between the climate models were, in all cases, smaller than the variation attributed to either emission scenarios or local conditions. The methods used for applying the climate model outputs played a larger role than the choice of the GCM or RCM. The thermal suitability for grain maize cultivation in Europe was estimated to expand by 30–50% across all SRES emissions scenarios. Strong increases in net primary productivity (NPP) (35–54%) were projected in northern European ecosystems as a result of a longer growing season and higher CO₂ concentrations. Changing water balance dominated the projected responses of southern European ecosystems, with NPP declining or increasing only slightly relative to present-day conditions. Both site and continental scale models showed large increases in yield of rain-fed winter wheat for northern Europe, with smaller increases or even decreases in southern Europe. Site-based, regional and continental scale models showed large spatial variations in the response of nitrate leaching from winter wheat cultivation to projected climate change due to strong interactions with soils and climate. The variation in simulated impacts was smaller between scenarios based on RCMs nested within the same GCM than between scenarios based on different GCMs or between emission scenarios.     

A numerical simulation of residual circulation in Tampa Bay. Part I: Low-frequency temporal variations

The residual (time-average) salinity and circulation in a numerical ocean model of the Tampa Bay estuary are shown to experience significant temporal variation under realistic forcing conditions. A version of the Estuarine Coastal Ocean Model developed for Tampa Bay with 70 by 100 horizontal grid points and 11 sigma levels is examined for the years 2001–2003. Model output variables are averaged over the entire time of the simulation to generate long-term residual fields. The residual axial current is found to be dominated by the buoyancy-driven baroclinic circulation with an outflow (southwestward) at the surface and to the sides of the shipping channel, and an inflow (northeastward) usually occurring subsurface within or above the shipping channel. Averages over 30 d are used to examine variations in the residual fields. During the simulation the average surface salinity near the head of Tampa Bay varies with the freshwater inflow, from 12‰ to 33%. At the bay mouth salinity varies from 30%. to 36%.. A localized measure of the baroclinic circulation in the shipping channel indicates the residual circulation can vary strongly, attaining a magnitude triple the long-term mean value. The baroclinic circulation can be disrupted, going to near zero or even reversing, when the buoyancy-driven flow is weak and the surface winds are to the northeast. Three time periods, representing different environmental conditions, are chosen to examine these results in detail. A scaling argument indicates the relative strength of buoyancy versus wind as ΔρgH²(LC _Dω²)⁻¹, where δρ is head-to-mouth density difference across the bay,g is gravitational acceleration,H is depth,L is bay length,C _D is the surface wind drag coefficient, andw is wind speed. Tampa Bay is usually in the buoyancy dominated regime. The importance of winds in the weak-buoyancy case is demonstrated in an additional simulation without wind stress.     

Infrared imagery of a turbulent intrusion in a stratified environment

A time sequence of airborne infrared imagery provides a unique view of phenomena associated with a turbulent tidal intrusion into a stratified bay. During flood tide, cooler water from the Strait of Juan de Fuca is observed to penetrate Sequim Bay (Washington, U.S.V.) as a turbulent jet. After separating from the shoreline, the jet collapses into the stratified middle part of the bay, forming a mushroom-shaped head consisting of a semicircular plunge front and areas of recirculating flow. As the plunge front advances into the estuary, a set of nonlinear internal waves emerges and propagates toward the relatively stagnant southern part of the bay, where they are a potential source of vertical mixing. This range of phenomena is expected based on laboratory studies, but has not been seen previously in a natural setting.     

Impulsively generated MHD waves and their detectability in solar coronal loops

Impulsively generated magnetohydrodynamic waves in solar coronal loops, with arbitrary plasma , are studied numerically by a flux-corrected transport algorithm. Numerical results show that the total reflection which occurs in the region of low Alfvén speed leads to trapped fast kink magnetosonic waves. These waves propagate along the slab and exhibit periodic, quasi-periodic, and decay phases. As a consequence of the difference in wave propagation speeds, the time signatures of the slow magnetosonic waves are delayed in time in comparison to the time signatures of the fast magnetosonic and Alfvén waves. An interaction between the waves can generate a longer lasting and complex quasi-periodic phase of the fast wave. We discuss also the observational detectability of such MHD waves in optical, radio, and soft X-ray wavelenghts.     

Paleolimnological significance of spinose and non-spinose morphs of Pyrgophorus hibbardi (Leonard and Franzen, 1944)

The basal portion of the Ogallala Formation (=Laverne Formation) (Lower Pliocene) Beaver County, Oklahoma, contains an interesting assemblage of non-marine fossil molluscs that include both spinose and non-spinose forms of the aquatic gastropod species Pyrgophorus hibbardi. The origin and paleolimnological significance of the spinose morph has been a source of much conjecture that has influenced environmental reconstructions of this assemblage. In one hypothesis the spinose forms of P. hibbardi are assumed to be associated with brackish water conditions by analogy with some populations of a related hydrobiid Potamopyrgus jenkinsi. To test the hypothesis that the spinose forms lived under different water conditions than the non-spinose morphs, we analyzed 10 specimens each of the two varieties for stable oxygen and carbon isotope ratios in the shell aragonite.The mean isotope ratios for the smooth and spinose morphs show no significant difference (oxygen: t = 0.28, df = 18, P (T t) 0.78 n.s.; carbon: t = 0.96, df = 18, P (T t) 0.35 n.s). We conclude that the lack of a statistically significant difference between the means of the oxygen and carbon isotope values for the smooth and spinose morphs suggests that the two forms lived in waters having similar isotope signatures. The considerable range in oxygen isotope values recorded by both morphs of P. hibbardi, including values as high as 5–6, suggest that both morphs were associated with waters which were periodically evaporatively enriched in ¹⁸O.