From mush to eruption in centuries: assembly of the super-sized Oruanui magma body

The magmatic systems that give rise to voluminous crystal-poor rhyolite magma bodies can be considered to operate on two contrasting timescales: Those governed by longer-term processes by which a magma acquires its chemical and isotopic characteristics (e.g., fractional crystallisation and assimilation), and those operating at shorter timescales during the physical accumulation of the melt-dominant magma body that finally erupts. We explore the compositional and textural relationships between amphibole and orthopyroxene crystals from the 25.4 ka, 530 km³ (magma) Oruanui eruption products (Taupo volcano, New Zealand) to investigate how processes related to the physical assembly of the pre-eruptive magma body are represented in the crystal record. Over 90 % of orthopyroxenes from the volumetrically dominant high-SiO₂ (>74 wt%) rhyolite pumices record textural evidence for a significant disequilibrium event (partial dissolution ± resorption of cores and interiors) prior to the growth of 40–500 μm thick rim zones. This dissolution/regrowth history of orthopyroxene is recorded in the chemistry of co-crystallising amphiboles as a prominent inflection in the concentrations of Mn and Zn, two elements notably enriched in orthopyroxene relative to amphibole. Textural and chemical features, linked with in situ thermobarometric estimates, indicate that a major decompression event preceded the formation of the melt-dominant body. The decompression event is inferred to represent the extraction of large volumes of melt plus crystals from the Oruanui crystal mush/source zone at pressures of 140–300 MPa (~6–12 km depth). Orthopyroxene underwent partial dissolution during ascent before reestablishing in the melt-dominant magma body at pressures of 90–140 MPa (~3.5–6 km). We model Fe–Mg diffusion across the core-rim boundaries along the crystallographic a or b-axes to constrain the timing of this decompression event, which marked establishment of the melt-dominant magma body. Maximum modelled ages indicate that this event did not begin until ~1,600 years before eruption, consistent with constraints from zircon model-age spectra. Once extraction began, it underwent runaway acceleration with a peak extraction age of ~230 years, followed by an apparent period of stasis of ~60 years prior to eruption. The rapidity of the extraction and accumulation processes implies the involvement of a dynamic driving force which, in the rifted continental arc setting of the Taupo Volcanic Zone, seems likely to be represented by magma-assisted extensional tectonic processes.     

Numerical relationships between magnetic parameters measured in Quaternary sediments and global paleoclimatic proxies

The complexity of most geological and geophysical problems prompts sometimes the use of non linear mathematical methods to handle them. An adaptive neuro fuzzy inference system (ANFIS) that combines fuzzy logic with neural networks, is applied here to study a paleoclimate section from the Quaternary sedimentary fill of the Lake Mucubají (western Venezuela). The purpose of this work is to find a set of numerical relationships that could predict the possible connections between oxygen isotope (δ¹⁸O) values from two different locations in the northern hemisphere (Ammersee in southern Germany and an ice core from the Greenland Ice Core Project — GRIP) and rock-magnetic parameters measured in Mucubají samples (i.e. mass-specific magnetic susceptibility — χ, magnetic remanence S-ratio, mass-specific saturation isothermal remanent magnetization — SIRM and anhysteretic remanent magnetization — ARM). The best inferences in terms of coefficient of determionation R² and the Root Mean-Square Error (RMSE) are obtained using those magnetic data as input that include information about magnetite grain size distributions, e.g., SIRM and ARM in FIS structures [1χ, 4ARM] and [4ARM, 1SIRM]. A comparison between Ammersee and GRIP actual data, as well as their corresponding inferences for the FIS structure [4ARM, 1SIRM], reveals a reasonable good inference of global trends for both records, overlooking the regional and/or local paleoclimate forcings that might have affected Ammersee. A better correlation between global isotope paleoclimate records and magnetic proxies, is perhaps prevented by the role played by local and regional paleoclimate and tectonism in Mucubají. We also argue that the ratio of ARM over SIRM appears to be related in a complex way to the onset and to the end of the Younger Dryas. Our novel approach to the assessment of a specific paleoclimate case study shows the potential of the ANFIS technique in solving problems where traditional univariate and multivariate linear regression methods could prove inadequate.     

Younger Dryas to mid-Holocene environmental history of the lowlands of NW Transylvania, Romania

Pollen, micro-charcoal and total carbon analyses on sediments from the Turbuta palaeolake, in the Transylvanian Basin of NW Romania, reveal Younger Dryas to mid-Holocene environmental changes. The chronostratigraphy relies on AMS ¹⁴C measurements on organic matter and U/Th TIMS datings of snail shells. Results indicate the presence of Pinus and Betula open woodlands with small populations of Picea, Ulmus, Alnus and Salix before 12,000 cal yr BP. A fairly abrupt replacement of Pinus and Betula by Ulmus-dominated woodlands at ca. 11,900 cal. yr BP likely represents competition effects of vegetation driven by climate warming at the onset of the Holocene. By 11,000 cal yr BP, the woodlands were increasingly diverse and dense with the expansion of Quercus, Fraxinus and Tilia, the establishment of Corylus and the decline of upland herbaceous and shrubs taxa. The marked expansion of Quercus accompanied by Tilia between 10,500 and 8000 cal yr BP could be the result of low effective moisture associated with both low elevation of the site and with regional change towards a drier climate. At 10,000 cal yr BP, Corylus spread across the region, and by 8000 cal yr BP it replaced Quercus as a dominant forest constituent, with only little representation of Picea abies. Carpinus became established around 5500 cal yr BP, but it was only a minor constituent in local woodlands until ca. 5000 cal yr BP. Results from this study also indicate that the woodlands in the lowlands of Turbuta were never closed.     

Characteristics of Danish estuaries

We review various aspects of the structure and functioning of Danish estuaries from data collected by the National Monitoring Program and from information in published sources. We present data on the physical, chemical, and biological characteristics of estuaries in Denmark, we evaluate the functioning of these systems as filters and transformers of nutrients and we evaluate the outlook for Danish estuaries in the future. Danish estuarine systems are for the most part shallow (<3 m deep), have short residence times, and tend to be heavily loaded with nutrients primarily from agricultural sources. Total average loads from land per unit watershed area are 112 kg P km^?2 yr^?1 and 2,400 kg N km^?2 yr^?1 during the period 1989–1995. The total phosphorus (TP) load in estuaries has been significantly reduced over the last decade, following implementation of the 1987 Action Plan for the Aquatic Environment (Vandmiljøplan in Danish) that prescribed that nitrogen loads to the total aquatic environment should be reduced by 50% and phosphorus loads by 80%. Reductions in the total nitrogen (TN) load have been more modest. Nutrient loading is one of the primary determinants of estuarine nutrient concentration with 70% of the annual variation in TN concentration and 55% of the annual variation in TP concentration explained by variation in the load. Many Danish estuaries have rich communities of macrophytes and benthic filter feeders, such asMytlis edulis andCiona intestinalis, that can control water column chlorophyll concentrations by their filter feeding activities. Many of the estuaries experience hypoxia and anoxia, especially during warm and calm summer months. Further reductions in nutrient loading are expected following implementation of the Action Plan for the Aquatic Environment II, with predicted improvements in oxygen concentrations and in the functioning of these shallow, dynamic estuarine systems.     

Does stress transmission in forelands depend on structural style? Distinctive stress magnitudes during Sevier thin‐skinned and Laramide thick‐skinned layer‐parallel shortening in the Bighorn Basin (USA) revealed by stylolite and calcite twinning paleopiezometry

The Sheep Mountain‐Little Sheep Mountain Anticlines, Bighorn Basin (USA) formed as basement‐cored Laramide structures in the formerly undeformed foreland of the thin‐skinned Sevier orogen. We take advantage of the well‐constrained microstructural network there to reconstruct differential stress magnitudes that prevailed during both Sevier and Laramide layer‐parallel shortening (LPS), before the onset of large‐scale folding. Differential stress magnitudes determined from tectonic stylolites are compared and combined to previous stress estimates from calcite twinning paleopiezometry in the same formations. During stress loading related to LPS, differential stress magnitudes transmitted from the distant Sevier thin‐skinned orogen into the sedimentary cover of the Bighorn basin (11–43 MPa) are higher than the differential stress magnitudes accompanying the early stage of LPS related to the thick‐skinned Laramide deformation (2–19 MPa). This study illustrates that the tectonic style of an orogen affects the transmission of early orogenic stress into the stable continental interior.     

Buried and Submerged Greek Archaeological Coastal Structures and Artifacts as Gauges to Measure Late Holocene Seafloor Subsidence off Calabria,Italy

This synthesis integrates recently acquired archaeological and geological data with earlier documented observations to shed light on the subsidence of ancient Greek coastal facilities in southern Italy. These are now positioned between former shorelines and inner shelf sectors at five Calabrian margin localities. Submergence of coastal to inner shelf facilities has resulted in part from sea‐level rise by about 2 m associated with glacio‐hydro‐isostatic factors since archaic to classic Greek time. This phenomenon alone, however, does not explain the wide variation of measured subsidence rates from site‐to‐site. The marked lowering of coastal site substrates by seismo‐tectonic activity (including extensional fault motion), stratal readjustments at depth, and compaction of underlying sediment sequences is significant. Four of the subsided facilities are positioned near emerged Calabrian areas where prevailing Holocene average annual land uplift rates range to ˜1.0 mm/yr; at the fifth, near Hipponion, terrains have risen by nearly 2 mm/yr. In marked contrast, submerged and/or buried structures record the following late Holocene long‐term average rates of coastal margin subsidence: Sybaris‐Thuri on the Taranto Gulf margin (˜0.5–1.0 mm/yr); Hipponion‐Vibo Valentia along the Tyrrhenian coast (˜0.8 to ˜3.2 mm/yr); and Locri‐Epizefiri, Kaulonia, and Capo Colonna on Calabria's Ionian margin (˜1.6, ˜1.6–2.4, and ˜4.0 mm/yr, respectively). © 2012 Wiley Periodicals, Inc.     

Deconvolving the δO seawater component from subseasonal coral δO and Sr/Ca at Rarotonga in the southwestern subtropical Pacific for the period 1726 to 1997

To reconstruct oceanographic variations in the subtropical South Pacific, 271-year long subseasonal time series of Sr/Ca and δ¹⁸O were generated from a coral growing at Rarotonga (21.5°S, 159.5°W). In this case, coral Sr/Ca appears to be an excellent proxy for sea surface temperature (SST) and coral δ¹⁸O is a function of both SST and seawater δ¹⁸O composition (δ¹⁸O_sw). Here, we focus on extracting the δ¹⁸O_sw signal from these proxy records. A method is presented assuming that coral Sr/Ca is solely a function of SST and that coral δ¹⁸O is a function of both SST and δ¹⁸O_sw. This method separates the effects of δ¹⁸O_sw from SST by breaking the instantaneous changes of coral δ¹⁸O into separate contributions by instantaneous SST and δ¹⁸O_sw changes, respectively. The results show that on average δ¹⁸O_sw at Rarotonga explains ∼39% of the variance in δ¹⁸O and that variations in SST explains the remaining ∼61% of δ¹⁸O variance. Reconstructed δ¹⁸O_sw shows systematic increases in summer months (December-February) consistent with the regional pattern of variations in precipitation and evaporation. The δ¹⁸O_sw also shows a positive linear correlation with satellite-derived estimated salinity for the period 1980 to 1997 (r = 0.72). This linear correlation between reconstructed δ¹⁸O_sw and salinity makes it possible to use the reconstructed δ¹⁸O_sw to estimate the past interannual and decadal salinity changes in this region. Comparisons of coral δ¹⁸O and δ¹⁸O_sw at Rarotonga with the Pacific decadal oscillation index suggest that the decadal and interdecadal salinity and SST variability at Rarotonga appears to be related to basin-scale decadal variability in the Pacific.     

A comparison of estimated and calculated effective porosity

 Effective porosity in solute-transport analyses is usually estimated rather than calculated from tracer tests in the field or laboratory. Calculated values of effective porosity in the laboratory on three different textured samples were compared to estimates derived from particle-size distributions and soil–water characteristic curves. The agreement was poor and it seems that no clear relationships exist between effective porosity calculated from laboratory tracer tests and effective porosity estimated from particle-size distributions and soil–water characteristic curves. A field tracer test in a sand-and-gravel aquifer produced a calculated effective porosity of approximately 0.17. By comparison, estimates of effective porosity from textural data, moisture retention, and published values were approximately 50–90% greater than the field calibrated value. Thus, estimation of effective porosity for chemical transport is highly dependent on the chosen transport model and is best obtained by laboratory or field tracer tests. Received, March 1997 · Revised, August 1997 · Accepted, August 1997     

Synchrotron IR study of hydrous ringwoodite (γ-Mg2SiO4) up to 30 GPa

High-pressure synchrotron infrared (IR) absorption spectra were collected between 650 and 4,000 cm⁻¹ at ambient temperature for hydrous Mg-ringwoodite (γ-Mg₂SiO₄) up to 30 GPa. The main feature in the OH⁻ stretching region is an extremely broad band centred at 3,150 cm⁻¹. The hydrogen bond is strong for most protons and the most probable site for protonation is the tetrahedral edge. With increasing pressure, this band shifts downward while decreasing its integrated intensity until disappearance at a pressure of 25 GPa. Only one band at 2,450 cm⁻¹ and an absorption plateau persist with a maximum wavenumber of 3,800 cm⁻¹. This behaviour is reversible upon pressure release. We interpret this as a second-order phase transition occurring in hydrated Mg-ringwoodite at high pressure (beyond ∼ 25 GPa). This result is compatible with the observation by Kleppe et al. (Phys Chem Miner 29:473–476, 2002a) who suggested the presence of Si–O–Si linkages and/or partial increase in the coordination of Si. Beyond the phase transition, the protons are delocalized and their environment on the ringwoodite structure is probably quite different from that at low pressure. Data obtained in situ at high pressures and temperatures are needed to better understand the effect of protonation on the structure and to better constrain this phase transition.     

Orthophosphate and biotite chemistry from orthopyroxene-bearing migmatites from California and South India: The role of a fluid-phase in the evolution of granulite-facies migmatites

Migmatites from Cone Peak, California, USA and the Satnur-Sangam road, Southern Karnataka, India contain coarser grained orthopyroxene-bearing leucosomes with subordinate biotite in finer grained hornblende-biotite-pyroxene-bearing hosts. At both localities the leucosomes are enriched in quartz and feldspar and have a higher ratio of pyroxene to hornblende + biotite compared to the host rocks. Biotite grains in leucosomes along the Satnur-Sangam road are concentrated at the margins of orthopyroxene grains and have lower abundances of Ti, Fe, and Cl and a higher abundance of F than biotite grains from the host rock. Fluorapatite grains in all rocks from both localities contain monazite inclusions similar to those produced experimentally by metasomatically induced dissolution and reprecipitation. Some fluorapatite grains at both localities are partially rimmed by allanite. The only compositional differences found between fluorapatite grains in the leucosomes and host rocks were higher concentrations of Cl in grains in leucosomes from Cone Peak. The mineralogies of the rocks suggest that the leucosomes formed by dehydration melting reactions that consumed feldspar, quartz, hornblende, and biotite and produced orthopyroxene. Allanite rims at the margins of fluorapatite grains may have formed by the later retrogression of monazite rims formed by incongruent dissolution of fluorapatite in the melt. Biotite grains at the margins of orthopyroxene crystals in the leucosomes from the Satnur-Sangam road apparently formed by retrogression of orthopyroxene upon the solidification of the anatectic melt. A similar high-grade retrogression did not affect orthopyroxene crystals at Cone Peak, indicating that H₂O was removed from the crystallizing leucosomes probably in a low H₂O activity fluid. Compositional differences between the paleosome and neosomes at Cone Peak are best explained by metasomatic interaction with concentrated brines while elevated Cl concentrations in fluorapatites in the leucosome suggest interaction with a Cl-bearing fluid. Brines may have been responsible for an exchange of elements between the host rock along the Satnur-Sangam road and zones of melt generation now marked by leucosomes, but fluid flow appears to have been less vigorous than at Cone Peak.