Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Searching for the Stonyhurst discs

Projects to drill into and recover samples from active faults are progressing in two widely separated locations: off the coast of Japan and in California.     

Helium isotopes in some historical lavas from Mount Vesuvius

³He/⁴He ratios in lavas erupted during the last 360 years at Mt. Vesuvius are between 2.2 and 2.7 R_A (R_A = atmospheric ratio of 1.39 × 10⁻⁶), and are among the lowest values measured in young volcanic rocks. They are also identical to values measured in summit crater fumaroles sampled during 1987–1991. This agreement indicates that the ³He/⁴He ratio in the crater fumaroles faithfully tracks the magmatic value. The relatively low and uniform ³He/⁴He ratio in the lavas reflects either a mantle source enriched in (U + Th)/³He, or a mixture of magmatic and crustal components.     

Trends in the aggregated rate of pre-eruptive volcano-tectonic seismicity at Kilauea volcano, Hawaii

Accelerating rates of volcano-tectonic (VT) earthquakes are commonly observed during volcanic unrest. Understanding the repeatability of their behaviour is essential to evaluating their potential to forecast eruptions. Quantitative eruption forecasts have focused on changes in precursors over intervals of weeks or less. Previous studies at basaltic volcanoes in frequent eruption, such as Kilauea in Hawaii and Piton de La Fournaise on Réunion, suggest that VT earthquake rates tend to follow a power-law acceleration with time about 2 weeks before eruption, but that this trend is often obscured by random fluctuations (or noise) in VT earthquake rate. These previous studies used a stacking procedure, in which precursory sequences for several eruptions are combined to enhance the signal from an underlying acceleration in VT earthquake rate. Such analyses assume a common precursory trend for all eruptions. This assumption is tested here for the 57 eruptions and intrusions recorded at Kilauea between 1959 and 1984. Applying rigorous criteria for selecting data (e.g. maximum depth; restricting magnitudes to be greater than the completeness magnitude, 2.1), we find a much less pronounced increase in the aggregate rate of earthquakes than previously reported. The stacked trend is also strongly controlled by the behaviour of one particular pre-eruptive sequence. In contrast, a robust signal emerges among stacked VT earthquake rates for a subset of the eruptions and intrusions. The results are consistent with two different precursory styles at Kilauea: (1) a small proportion of eruptions and intrusions that are preceded by accelerating rates of VT earthquakes over intervals of weeks to months and (2) a much larger number of eruptions that show no consistent increase until a few hours beforehand. The results also confirm the importance of testing precursory trends against data that have been filtered according to simple constraints on the spatial distribution and completeness magnitude of the VT earthquakes.     

Continental-scale patterns in soil geochemistry and mineralogy: Results from two transects across the United States and Canada

In 2004, the US Geological Survey (USGS) and the Geological Survey of Canada (GSC) initiated a pilot study that involved collection of more than 1500 soil samples from 221 sites along two continental transects across Canada and the United States. The pilot study was designed to test and refine protocols for a soil geochemical survey of North America. The two transects crossed a wide array of soil parent materials, soil ages, climatic conditions, landforms, land covers and land uses. Sample sites were selected randomly at approximately 40-km intervals from a population defined as all soils of the continent. At each site, soils representing 0 to 5 cm depth, and the O, A, and C horizons, if present, were collected and analyzed for their near-total content of over 40 major and trace elements. Soils from 0–5 cm depth were also collected for analysis of organic compounds. Results from the transects confirm that soil samples collected at a 40-km spacing reveal coherent, continental- to subcontinental-scale geochemical and mineralogical patterns that can be correlated to aspects of underlying soil parent material, soil age and climate influence. The geochemical data also demonstrate that at the continental-scale the dominance of any of these major factors that control soil geochemistry can change across the landscape. Along both transects, soil mineralogy and geochemistry change abruptly with changes in soil parent materials. However, the chemical influence of a soil’s parent material can be obscured by changing climatic conditions. For the transects, increasing precipitation from west to east and increasing temperature from north to south affect both soil mineralogy and geochemistry because of climate effects on soil weathering and leaching, and plant productivity. Regional anomalous metal concentrations can be linked to natural variations in soil parent materials, such as high Ni and Cr in soils developed on ultramafic rocks in California or high P in soils formed on weathered Ordovician limestones in central Kentucky. On local scales, anomalous metal concentrations recognized in soil profiles, such as high P in soils from animal confinement sites, are consistent with local anthropogenic disturbances. At a larger scale, the distribution of Hg across the west to east transect demonstrates that it can be difficult to distinguish between natural or anthropogenic contributions and that many factors can contribute to an element’s spatial distribution.Only three samples in a subset of seventy-three 0–5 cm depth soil samples from the north to south transect had organochlorine pesticides values above the method detection limit, apparently related to historic usage of the pesticides DDT and dieldrin.     

The volcanic and geothermally active Campi Flegrei caldera: an integrated multidisciplinary image of its buried structure

The Campi Flegrei caldera in southern Italy is one of the greatest geohazard areas on Earth. Evidence of an active magmatic and geothermal system is provided by ongoing ground uplift, with volcano-tectonic and long-period (LP) seismicity, the persistent degassing of ~1500 tonnes of CO₂ per day, the presence of hot fumaroles at temperatures of 90–150 °C, brine-rich aquifers (with total dissolved solids up to 33 g l^?1) and high thermal gradients in the crust (with temperatures reaching 420 °C at 3,050 m b.s.l.). Since the 1940s, more than 100 exploratory boreholes have been drilled in the area to depths of 80–3,100 m by the Azienda Geologica Italiana Petroli (AGIP) and the Società Anonima Forze Endogene Napoletane (SAFEN). To date, however, no systematic reanalysis of the drilling data has been carried out, and the buried volcanic structure has not been updated using the most recent scientific results and previous findings. By integrating unpublished data from the AGIP and SAFEN reports with published information from geological, volcanological, petrological, petrophysical and geophysical studies, this paper presents an improved picture of the Campi Flegrei caldera that will be useful for volcanic hazard assessment and mitigation in the Naples area and for future research planning. The results suggest that intra-caldera activity has been influenced by how the magmatic system at depths greater than about 4 km has determined the transfer of magma, volatiles, and heat to the overlying geothermal system and, ultimately, to the surface. In particular, intriguing is that the most volcanically active central-eastern sector of the caldera, which is subject to intense bradyseismic ground movement and gas emission, coincides with a structurally delimited subsurface rock volume characterized by an uprising of the 100 °C isotherm, a deep water supply to the shallower aquifer, the early disappearance of secondary calcite, LP seismicity and high seismic S-wave attenuation. In this area, we also document evidence of repeated injection at depths of c. 1.5–3.0 km of isolated and small-volume batches of magma, where occurred their crystallization and degassing. Shallow intrusions and degassing of magma are thus identified as two of the key processes that drive unrest in Campi Flegrei.     

Precursors to dyke-fed eruptions at basaltic volcanoes: insights from patterns of volcano-tectonic seismicity at Kilauea volcano,Hawaii

To investigate the physical controls on volcano-tectonic (VT) precursors to eruptions and intrusions at basaltic volcanoes, we have analyzed the spatial and temporal patterns of VT earthquakes associated with 34 eruptions and 23 dyke intrusions that occurred between 1960 and 1983 at Kilauea, in Hawaii. Eighteen of the 57 magmatic events were preceded by an acceleration of the mean rate of VT earthquakes located close to the main shallow magma reservoir. Using a maximum-likelihood technique and the Bayesian Information Criterion for model preference, we demonstrate that an exponential acceleration is preferred over a power-law acceleration for all sequences. These sequences evolve over time-scales of weeks to months and are consistent with theoretical models for the approach to volcanic eruptions based on the growth of a population of fractures in response to an excess magma pressure. Among the remaining 40 magmatic events, we found a significant correlation between swarms of VT earthquakes located in the mobile south-flank of Kilauea and eruptions and intrusions. The behaviour of these swarms suggests that at least some of the magmatic events are triggered by transient episodes of elevated rates of aseismic flank movement, which could explain why many eruptions and intrusions are not preceded by longer-term precursory signals. In none of the 57 cases could a precursory sequence be used to distinguish between the approach to an eruption or an intrusion, so that, even when a precursory sequence is recognized, there remains an empirical chance of about 40% (24 intrusions from 57 magmatic events) of issuing a false alarm for an imminent eruption.     

Pilot studies for the North American Soil Geochemical Landscapes Project – Site selection, sampling protocols, analytical methods, and quality control protocols

In 2004, the US Geological Survey (USGS) and the Geological Survey of Canada sampled and chemically analyzed soils along two transects across Canada and the USA in preparation for a planned soil geochemical survey of North America. This effort was a pilot study to test and refine sampling protocols, analytical methods, quality control protocols, and field logistics for the continental survey. A total of 220 sample sites were selected at approximately 40-km intervals along the two transects. The ideal sampling protocol at each site called for a sample from a depth of 0–5 cm and a composite of each of the O, A, and C horizons. The <2-mm fraction of each sample was analyzed for Al, Ca, Fe, K, Mg, Na, S, Ti, Ag, As, Ba, Be, Bi, Cd, Ce, Co, Cr, Cs, Cu, Ga, In, La, Li, Mn, Mo, Nb, Ni, P, Pb, Rb, Sb, Sc, Sn, Sr, Te, Th, Tl, U, V, W, Y, and Zn by inductively coupled plasma-mass spectrometry and inductively coupled plasma-atomic emission spectrometry following a near-total digestion in a mixture of HCl, HNO₃, HClO₄, and HF. Separate methods were used for Hg, Se, total C, and carbonate-C on this same size fraction. Only Ag, In, and Te had a large percentage of concentrations below the detection limit. Quality control (QC) of the analyses was monitored at three levels: the laboratory performing the analysis, the USGS QC officer, and the principal investigator for the study. This level of review resulted in an average of one QC sample for every 20 field samples, which proved to be minimally adequate for such a large-scale survey. Additional QC samples should be added to monitor within-batch quality to the extent that no more than 10 samples are analyzed between a QC sample. Only Cr (77%), Y (82%), and Sb (80%) fell outside the acceptable limits of accuracy (% recovery between 85 and 115%) because of likely residence in mineral phases resistant to the acid digestion.A separate sample of 0–5-cm material was collected at each site for determination of organic compounds. A subset of 73 of these samples was analyzed for a suite of 19 organochlorine pesticides by gas chromatography. Only three of these samples had detectable pesticide concentrations. A separate sample of A-horizon soil was collected for microbial characterization by phospholipid fatty acid analysis (PLFA), soil enzyme assays, and determination of selected human and agricultural pathogens. Collection, preservation and analysis of samples for both organic compounds and microbial characterization add a great degree of complication to the sampling and preservation protocols and a significant increase to the cost for a continental-scale survey. Both these issues must be considered carefully prior to adopting these parameters as part of the soil geochemical survey of North America.     

The evolution of lava flow-fields: observations of the 1981 and 1983 eruptions of Mount Etna,Sicily

The eruptions of Mount Etna in 1981 on the north flank and 1983 on the south flank of the volcano were of strikingly different character. The former was a short duration, high effusion rate eruption producing for the most part a simple flow-field; the latter was of relatively long duration and low effusion rate, producing a compound flow-field of overlapping flows.Despite the differences between the eruptive behaviour of these two events and the way in which the flow-field developed, both the flow-fields achieved about the same maximum length. This is considered fortuitous. The evidence suggests that the main 1981 flow stopped because the lava supply ceased and was thus volume controlled. The 1983 flow-field had a more complex history of branching, but in this case it appears that, for the longest individual flow, cooling played an important role in controlling the maximum extent of the flows.