Komsomolskaya diamondiferous eclogites: evidence for oceanic crustal protoliths

The Komsomolskaya kimberlite is one of numerous (>1,000) kimberlite pipes that host eclogite xenoliths on the Siberian craton. Eclogite xenoliths from the adjacent Udachnaya kimberlite pipe have previously been geochemically well characterized; however, data from surrounding diamond-bearing kimberlite pipes from the center of the craton are relatively sparse. Here, we report major- and trace-element data, as well as oxygen isotope systematics, for mineral separates of diamondiferous eclogite xenoliths from the Komsomolskaya kimberlite, suggesting two distinct subgroups of a metamorphosed, subducted oceanic crustal protolith. Using almandine contents, this suite can be divided into two subgroups: group B1, with a high almandine component (>20 mol%) and group B2, with a low almandine component (<20 mol%). Reconstructed REE profiles for B1 eclogites overlap with typical oceanic basalts and lack distinct Eu anomalies. In addition, elevated oxygen isotope values, which are interpreted to reflect isotopic exchange with seawater at low temperatures (<350 °C), are consistent with an upper-oceanic crustal protolith. Reconstructed REE profiles for B2 eclogites are consistent with oceanic gabbros and display distinct Eu anomalies, suggesting a plagioclase-rich cumulate protolith. In contrast to B1, B2 eclogites do not display elevated oxygen isotope values, suggesting an origin deep within the crustal pile, where little-to-no interaction with hydrothermal fluids has occurred. Major-element systematics were reconstructed based on mineral modes; group B1 eclogites have higher MgO wt% and lower SiO₂ wt%, with respect to typical oceanic basalts, reflecting a partial melting event during slab subduction. Calculated residues from batch partial melt modeling of a range of Precambrian basalts overlap with group B1 trace-element chemistry. When taken together with the respective partial melt trajectories, these melting events are clearly linked to the formation of Tonalite–Trondhjemite–Granodiorite (TTG) complexes. As a result, we propose that many, if not all, diamondiferous eclogite xenoliths from Komsomolskaya represent mantle ‘restites’ that preserve chemical signatures of Precambrian oceanic crust.     

Distribution and Abundance of Arsenic in the Soils and Sediments of O‘ahu,Hawai‘i

Lack of high-spatial-resolution soil and sediment arsenic data for Hawai‘i has generated substantial disagreement between researchers and regulators regarding the magnitude of natural levels of arsenic in Hawai‘i and rendered difficult the defining of areas of anthropogenically elevated arsenic. Our earlier research into the occurrence of arsenic in terrestrial and marine environments revealed widely disparate concentrations of arsenic with no apparent spatial pattern. To better understand the distribution and abundance of arsenic in soils and sediments of O‘ahu, we collected an additional 64 samples at locations chosen to represent different environments with varying degrees of human impact. We found surface arsenic values that ranged from 0.28 to 740 ppm with a median concentration of 8.1 ppm, which is above the global median of 5 ppm and US soil median of 5.2 ppm. Higher concentrations of arsenic (up to 913 ppm) were encountered at depth in soil cores. The median arsenic in streambed sediments from one of our earlier studies of 6.1 ppm was comparable to the conterminous US median of 6.3 ppm; however, we encountered arsenic concentrations as high as 43.9 ppm (median = 8.60 ppm, n = 75) in marine sediments in recent work off the leeward coast of O‘ahu. Overall, arsenic in the soils and sediments of O‘ahu is elevated relative to world and national values, but there still is no readily discernible pattern in the distribution of arsenic to explain these elevated values.     

Combining long- and short-term probabilistic volcanic hazard assessment with cost-benefit analysis to support decision making in a volcanic crisis from the Auckland Volcanic Field,New Zealand

By using BET_VH, we propose a quantitative probabilistic hazard assessment for base surge impact in Auckland, New Zealand. Base surges resulting from phreatomagmatic eruptions are among the most dangerous phenomena likely to be associated with the initial phase of a future eruption in the Auckland Volcanic Field. The assessment is done both in the long-term and in a specific short-term case study, i.e. the simulated pre-eruptive unrest episode during Exercise Ruaumoko, a national civil defence exercise. The most important factors to account for are the uncertainties in the vent location (expected for a volcanic field) and in the run-out distance of base surges. Here, we propose a statistical model of base surge run-out distance based on deposits from past eruptions in Auckland and in analogous volcanoes. We then combine our hazard assessment with an analysis of the costs and benefits of evacuating people (on a 1 × 1-km cell grid). In addition to stressing the practical importance of a cost-benefit analysis in creating a bridge between volcanologists and decision makers, our study highlights some important points. First, in the Exercise Ruaumoko application, the evacuation call seems to be required as soon as the unrest phase is clear; additionally, the evacuation area is much larger than what is recommended in the current contingency plan. Secondly, the evacuation area changes in size with time, due to a reduction in the uncertainty in the vent location and increase in the probability of eruption. It is the tradeoff between these two factors that dictates which cells must be evacuated, and when, thus determining the ultimate size and shape of the area to be evacuated.     

Application of surface-wave tomography to mineral exploration: a case study from Siilinjärvi,Finland

In order to assess the feasibility and validity of surface-wave tomography as a tool for mineral exploration, we present an active seismic three-dimensional case study from the Siilinjärvi mine in Eastern Finland. The aim of the survey is to identify the formation carrying the mineralization in an area south of the main pit, which will be mined in the future. Before acquiring the data, we performed an accurate survey design to maximize data coverage and minimize the time for deployment and recollection of the equipment. We extract path-averaged Rayleigh-wave phase-velocity dispersion curves by means of a two-station method. We invert them using a computationally efficient tomographic code which does not require the computation of phase-velocity maps and inverts directly for one-dimensional S-wave velocity models. The retrieved velocities are in good agreement with the data from a borehole in the vicinity, and the pseudo three–dimensional S-wave velocity volume allows us to identify the geological contact between the formation hosting most of the mineralization and the surrounding rock. We conclude that the proposed method is a valid tool, given the small amount of equipment used and the acceptable amount of time required to process the data.     

Estimating groundwater evapotranspiration rates using diurnal water-table fluctuations in a semi-arid riparian zone

In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ET_G) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ET_G, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ET_G to potential evapotranspiration (PET), K _c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K _c, in semi-arid rangelands. Before vegetation senescence, K _c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K _c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ET_G is primarily from estimates of specific yield (S _y), which is difficult to quantify in heterogeneous sediments. ET_G values may be more reliable because the range of acceptable S _y is smaller than K _c and S _y does not change with vegetation type or soil moisture.     

Detecting the stepwise propagation of the Eastern Sicily thrust belt (Italy): insight from thermal and thermochronological constraints

An extensive dataset of vitrinite reflectance, FTIR parameters on organic matter, illite content in mixed layers illite‐smectite, apatite fission tracks and U‐Th/He dating has been used to reconstruct the stepwise propagation of the Eastern Sicily fold‐and‐thrust belt during Late Palaeogene and Neogene times. The results indicate that the fold‐and‐thrust belt is divisible into two levels of thermal maturity. These levels consist of a less evolved level of thermal maturity that records limited sedimentary burial and minor heating, and a more evolved level of thermal maturity that indicates tectonic burial and exhumation at different times. Deformation and exhumation of shallowly buried units are linked to wedge forward propagation by low‐angle thrusts, whereas the evolution of deeply buried units is associated with tectonic imbrications by duplex formation and steep thrusts. The two tectonic styles alternate during evolution of the fold‐and‐thrust belt under low erosion rates.     

Multi-Modal Monitoring of Slip Along Frictional Discontinuities

Seismic wave transmission and digital image correlation (DIC) are employed to study slip processes along frictional discontinuities. A series of biaxial compression experiments are performed on gypsum specimens with non-homogeneous contact surfaces. The specimens are composed of two blocks with perfectly mated contact surfaces with a smooth surface with low frictional strength on the upper half and a rough surface with high frictional strength on the lower half. Compressional, P, and shear, S, wave pulses were transmitted through the discontinuity while digital images of the specimen surface were acquired during the test. A distinct peak in the amplitude of transmitted wave occurs prior to the peak shear strength and is considered a “precursor” to the failure. Precursors indicate that slip initiates from the smooth surface and extends to the rough surface as the shear load is increased. From the DIC data, slip is identified as a jump in the displacement field along the fracture that initiates from the smooth surface and propagates to the rough surface. Precursors are associated with an increase in the rate of slip across the discontinuity and are a measure of the reduction in the fracture shear stiffness.