Source mechanism of volcanic tremor: fluid-driven crack models and their application to the 1963 kilauea eruption

We propose a model for the mechanism of magma transport based on a fluid-filled tensile crack driven by the excess pressure of fluid. Such a transport mechanism can generate seismic waves by a succession of jerky crack extensions, if the fracture strength of rock varies in space, or if there is a difference between the dynamic and static values of the critical stress intensity factor. We also find that the opening and closing of a narrow channel connecting two fluid-filled cracks may be a source of seismic waves. Using the finite-difference method, we calculated the vibration of dry and fluid-filled cracks generated by: (1) a jerky extension at one end or at both ends and (2) a jerky opening of a narrow channel connecting two cracks. We then calculated the far-field and near-field radiation from these vibrating cracks. The spectra show peaked structures, but interestingly, most high-frequency peaks are only present in the near-field and cannot be transmitted to the far-field. The spectral features described above are often observed for volcanic tremors and in some cases for seismic signals associated with hydraulic fracturing experiments.We first consider as a model of volcanic tremor randomly occurring jerky crack extensions, and derive a formula relating the tremor amplitude to the excess pressure in the magma, the incremental area in each extension, and the frequency of extensions. These parameters are also constrained by other observations, such as the rate of magma flow.Our model was tested quantitatively against observations made in one of the best-described case histories of volcanic tremor: the October 5–6, 1963 Kilauea flank eruption. We found that a single, long crack extending from the summit to the eruptive site cannot explain the observations. The model of a steadily expanding crack ran into difficulties when quantitative comparisons were made with observations. The extension of crack area needed to explain the amplitude of volcanic tremor should accompany a large increase in tremor period which was not observed.Our second model is a chain of cracks connected by narrow channels which open and close. The length of each crack is around 1 km, the channel area connecting neighboring cracks is about 10³m², and the channel opens jerkily with the magmatic excess pressure of about 20 bars. The frequency of jerky opening of each channel is about once in 15 seconds. The channel is closed after each jerky opening, as soon as magma is moved through the channel.     

Time to talk knowledge

As part of the secondary school qualification repackaging, the Ministry of Education has requested the geography education community define what constitutes their “critical body of knowledge.” In this article, I discuss three knowledge‐oriented themes focused on: (a) recent assessment trends in geography; (b) the predictive value of NCEA Level 3 geography with 100‐level success and (c) the nature of knowledge as expressed by the New Zealand Curriculum and assessment. Taken together, these themes suggest a need for a serious conversation about knowledge in school geography, as a pre‐cursor to the re‐design of the geography standards matrix.     

U–Pb and Th–Pb dating of apatite by LA-ICPMS

Apatite is a common U- and Th-bearing accessory mineral in igneous and metamorphic rocks, and a minor but widespread detrital component in clastic sedimentary rocks. U–Pb and Th–Pb dating of apatite has potential application in sedimentary provenance studies, as it likely represents first cycle detritus compared to the polycyclic behavior of zircon. However, low U, Th and radiogenic Pb concentrations, elevated common Pb and the lack of a U–Th–Pb apatite standard remain significant challenges in dating apatite by LA-ICPMS, and consequently in developing the chronometer as a provenance tool.This study has determined U–Pb and Th–Pb ages for seven well known apatite occurrences (Durango, Emerald Lake, Kovdor, Mineville, Mud Tank, Otter Lake and Slyudyanka) by LA-ICPMS. Analytical procedures involved rastering a 10 μm spot over a 40 × 40 μm square to a depth of 10 μm using a Geolas 193 nm ArF excimer laser coupled to a Thermo ElementXR single-collector ICPMS. These raster conditions minimized laser-induced inter-element fractionation, which was corrected for using the back-calculated intercept of the time-resolved signal. A Tl–U–Bi–Np tracer solution was aspirated with the sample into the plasma to correct for instrument mass bias. External standards (Ple?ovice and 91500 zircon, NIST SRM 610 and 612 silicate glasses and STDP5 phosphate glass) along with Kovdor apatite were analyzed to monitor U–Pb, Th–Pb, U–Th and Pb–Pb ratiosCommon Pb correction employed the ²⁰⁷Pb method, and also a ²⁰⁸Pb correction method for samples with low Th/U. The ²⁰⁷Pb and ²⁰⁸Pb corrections employed either the initial Pb isotopic composition or the Stacey and Kramers model and propagated conservative uncertainties in the initial Pb isotopic composition. Common Pb correction using the Stacey and Kramers (1975) model employed an initial Pb isotopic composition calculated from either the estimated U–Pb age of the sample or an iterative approach. The age difference between these two methods is typically less than 2%, suggesting that the iterative approach works well for samples where there are no constraints on the initial Pb composition, such as a detrital sample. No ²⁰⁴Pb correction was undertaken because of low ²⁰⁴Pb counts on single collector instruments and ²⁰⁴Pb interference by ²⁰⁴Hg in the argon gas supply.Age calculations employed between 11 and 33 analyses per sample and used a weighted average of the common Pb-corrected ages, a Tera–Wasserburg Concordia intercept age and a Tera–Wasserburg Concordia intercept age anchored through common Pb. The samples in general yield ages consistent (at the 2σ level) with independent estimates of the U–Pb apatite age, which demonstrates the suitability of the analytical protocol employed. Weighted mean age uncertainties are as low as 1–2% for U- and/or Th-rich Palaeozoic–Neoproterozoic samples; the uncertainty on the youngest sample, the Cenozoic (31.44 Ma) Durango apatite, ranges from 3.7–7.6% according to the common Pb correction method employed. The accurate and relatively precise common Pb-corrected ages demonstrate the U–Pb and Th–Pb apatite chronometers are suitable as sedimentary provenance tools. The Kovdor carbonatite apatite is recommended as a potential U–Pb and Th–Pb apatite standard as it yields precise and reproducible ²⁰⁷Pb-corrected, ²³²Th–²⁰⁸Pb, and common Pb-anchored Tera–Wasserburg Concordia intercept ages.     

Is rainfall intensity significant in the rainfall–runoff process within tropical rainforests of northeast Queensland? The Hewlett regression analyses revisited

Following the statistical analyses of long‐term rainfall‐runoff records from research basins in humid temperate latitudes, Hewlett and co‐workers extended the global challenge to disprove their findings that rainfall intensity was non‐significant. This paper responds to Hewlett's challenge as no preceding analyses have involved forested basins in a tropical cyclone‐prone area. Based on a 7 year rainfall‐runoff record, quickflow (QF), peak flow (QP) and quickflow response ratios (QRR) were regressed as dependent variables against rainfall parameters (intensity, P_i, amount, P), storm duration, D and antecedent flow, I. These data sets were categorised into total streamflow (Q) classes and stratified into three seasons, (monsoon, post‐monsoon and dry) for forested and cleared catchments. Where rainfall variable collinearity met acceptable levels, the addition of P_i to regression models including P, D, I contributed up to 9% and 66% of the respective variations in quickflow and peak flow. For the highest Q storm classes (monsoon), P_i alone accounted for up to 67% and 91% of the variation in Q_F and Q_P respectively and was the dominant influence on Q_P for all seasons. The very high rainfall intensities experienced in the monsoon season is a causal factor why these results differ from those of other research drainage basins. Surprisingly, P_i continued to have a significant influence on Q_F for dry season classes when less‐intense rainfall occurs. Further the results were similar for both catchments across all seasons. P was the dominant independent variable affecting Q_F above a threshold Q of 50 mm (monsoon), as rainfall contributes directly to saturation overland flow and return flow under saturated conditions. Further although QRR increased with increasing Q for each season, the regression results for that parameter were poor possibly due to the non‐linearity of the rainfall‐runoff relationship. Copyright © 2010 John Wiley & Sons, Ltd.     

Ar/Ar thermochronology of the Kampa Dome, southern Tibet: Implications for tectonic evolution of the North Himalayan gneiss domes

Structural and thermochronological studies of the Kampa Dome provide constraints on timing and mechanisms of gneiss dome formation in southern Tibet. The core of Kampa Dome contains the Kampa Granite, a Cambrian orthogneiss that was deformed under high temperature (sub-solidus) conditions during Himalayan orogenesis. The Kampa Granite is intruded by syn-tectonic leucogranite dikes and sills of probable Oligocene to Miocene age. Overlying Paleozoic to Mesozoic metasedimentary rocks decrease in peak metamorphic grade from kyanite + staurolite grade at the base of the sequence to unmetamorphosed at the top. The Kampa Shear Zone traverses the Kampa Granite — metasediment contact and contains evidence for high-temperature to low-temperature ductile deformation and brittle faulting. The shear zone is interpreted to represent an exhumed portion of the South Tibetan Detachment System. Biotite and muscovite ⁴⁰Ar/³⁹Ar thermochronology from the metasedimentary sequence yields disturbed spectra with 14.22 ± 0.18 to 15.54 ± 0.39 Ma cooling ages and concordant spectra with 14.64 ± 0.15 to 14.68 ± 0.07 Ma cooling ages. Petrographic investigations suggest disturbed samples are associated with excess argon, intracrystalline deformation, mineral and fluid inclusions and/or chloritization that led to variations in argon systematics. We conclude that the entire metasedimentary sequence cooled rapidly through mica closure temperatures at  14.6 Ma. The Kampa Granite yields the youngest biotite ⁴⁰Ar/³⁹Ar ages of  13.7 Ma immediately below the granite–metasediment contact. We suggest that this age variation reflects either varying mica closure temperatures, re-heating of the Kampa Granite biotites above closure temperatures between 14.6 Ma and 13.7 Ma, or juxtaposition of rocks with different thermal histories. Our data do not corroborate the “inverse” mica cooling gradient observed in adjacent North Himalayan gneiss domes. Instead, we infer that mica cooling occurred in response to exhumation and conduction related to top-to-north normal faulting in the overlying sequence, top-to-south thrusting at depth, and coeval surface denudation.     

Natural aggregates in the performance and durability of concrete: chemical characteristics

This article continues from the earlier feature in Geology Today (2011, v.27, n.4) and looks at the influence of chemical reactivity within concrete due to the presence of unsound aggregates or deleterious material incorporated in the aggregate used in the concrete mix. Some of the more disruptive problems from chemical unsoundness come from the presence of chlorides, which cause rusting in reinforced concrete or, sulphate attack from ground waters. Both of these reactions are relatively common in Britain. Alkali–silica reactions with aggregates and cement is not too prevalent in Britain but can be very prevalent and disruptive elsewhere in the world.     

Guidelines on the use of structure-from-motion photogrammetry in geomorphic research

As a topographic modelling technique, structure-from-motion (SfM) photogrammetry combines the utility of digital photogrammetry with a flexibility and ease of use derived from multi-view computer vision methods. In conjunction with the rapidly increasing availability of imagery, particularly from unmanned aerial vehicles, SfM photogrammetry represents a powerful tool for geomorphological research. However, to fully realize this potential, its application must be carefully underpinned by photogrammetric considerations, surveys should be reported in sufficient detail to be repeatable (if practical) and results appropriately assessed to understand fully the potential errors involved. To deliver these goals, robust survey and reporting must be supported through (i) using appropriate survey design, (ii) applying suitable statistics to identify systematic error (bias) and to estimate precision within results, and (iii) propagating uncertainty estimates into the final data products. © 2019 John Wiley & Sons, Ltd.     

River stabilisation due to changing climate and vegetation during the late Quaternary in western Tasmania, Australia

The Stanley River in western Tasmania, Australia, contains sub-fossil rainforest logs within the channel and floodplain. Of the more than 85 radiocarbon dates obtained, all but 3 date from 17 ka to the present and permit an interpretation of fluvial and related environmental changes over this period. Particular attention is focused on the interactive relationship between the river and its riparian rainforest. Following the Last Glacial Maximum, the Stanley River was a laterally active gravel-load system reworking most of its valley floor in the upstream reaches. With ameliorating conditions at the end of the Pleistocene, climate became less seasonal and flow regimes less energetic. Huon pines already present in the catchment, re-asserted themselves in the form of dense tree cover along the river banks and floodplains with basal floodplain deposition shifting from gravels to coarse sands and granules. By about 3.5 ka, a further change in climate reduced stream discharges substantially. As a result the channel reduced in size, transported finer sediment, became laterally stable, and the floodplain accreted with overbank deposits of sand and silt. Huon pines falling into the channel formed obstructions of woody debris, some surviving for 2 ka. These have reduced stream power and boundary shear stress, further contributing to channel stability. Generational sequences of Huon pines on the river banks, some extending back 1–2 ka, are additional evidence of this stability. Since the Pleistocene, changing climate and the re-establishment of dense riparian rainforest appear to have stabilised the river channels and floodplains of western Tasmania.     

Water Table Uncertainties due to Uncertainties in Structure and Properties of an Unconfined Aquifer

We consider two sources of geology‐related uncertainty in making predictions of the steady‐state water table elevation for an unconfined aquifer. That is the uncertainty in the depth to base of the aquifer and in the hydraulic conductivity distribution within the aquifer. Stochastic approaches to hydrological modeling commonly use geostatistical techniques to account for hydraulic conductivity uncertainty within the aquifer. In the absence of well data allowing derivation of a relationship between geophysical and hydrological parameters, the use of geophysical data is often limited to constraining the structural boundaries. If we recover the base of an unconfined aquifer from an analysis of geophysical data, then the associated uncertainties are a consequence of the geophysical inversion process. In this study, we illustrate this by quantifying water table uncertainties for the unconfined aquifer formed by the paleochannel network around the Kintyre Uranium deposit in Western Australia. The focus of the Bayesian parametric bootstrap approach employed for the inversion of the available airborne electromagnetic data is the recovery of the base of the paleochannel network and the associated uncertainties. This allows us to then quantify the associated influences on the water table in a conceptualized groundwater usage scenario and compare the resulting uncertainties with uncertainties due to an uncertain hydraulic conductivity distribution within the aquifer. Our modeling shows that neither uncertainties in the depth to the base of the aquifer nor hydraulic conductivity uncertainties alone can capture the patterns of uncertainty in the water table that emerge when the two are combined.