Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.     

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.     

The practice of DEM stream burning revisited

Stream burning is a common flow enforcement technique used to correct surface drainage patterns derived from digital elevation models (DEM). The technique involves adjusting the elevations of grid cells that are coincident with the features of a vector hydrography layer. This paper focuses on the problematic issues with common stream burning practices, particularly the topological errors resulting from the mismatched scales of the hydrography and DEM data sets. A novel alternative stream burning method is described and tested using five DEMs of varying resolutions (1 to 30 arc‐seconds) for an extensive area of southwestern Ontario, Canada. This TopologicalBreachBurn method uses total upstream channel length (TUCL) to prune the vector hydrography layer to a level of detail that matches the raster DEM grid resolution. Network pruning reduces the occurrence of erroneous stream piracy caused by the rasterization of multiple stream links to the same DEM grid cell. The algorithm also restricts flow within individual stream reaches, further reducing erroneous stream piracy. In situations where two vector stream features occupy the same grid cell, the new tool ensures that the larger stream, designated by higher TUCL, is given priority. TUCL‐based priority minimizes the impact of the topological errors that occur during the stream rasterization process on modeled regional drainage patterns. The test data demonstrated that TopologicalBreachBurn produces highly accurate and scale‐insensitive drainage patterns and watershed boundaries. The drainage divides of four large watersheds within the study region that were delineated from the TopologicalBreachBurn‐processed DEMs were found to be highly accurate when compared with the official watershed boundaries, even at the coarsest grid resolutions, with Kappa index of agreement values ranging from 0.952 to 0.921. The corresponding Kappa coefficient values for a traditional stream burning method (FillBurn) ranged from 0.953 to 0.490, demonstrating a significant decrease in mapping accuracy at coarser DEM grid resolutions. Copyright © 2015 John Wiley & Sons, Ltd.     

Measuring the significance of a divide to local drainage patterns

This article presents a framework for estimating a new topographic attribute derived from digital elevation models (DEMs) called maximum branch length (B _max). Branch length is defined as the distance travelled along a flow path initiated at one grid cell to the confluence with the flow path passing through a second cell. B _max is the longest branch length measured for a grid cell and its eight neighbours. The index provides a physically meaningful method for assessing the relative significance of drainage divides to the dispersion of materials and energy across a landscape, that is, it is a measure of ‘divide size’. B _max is particularly useful for studying divide network structure, for mapping drainage divides, and in landform classification applications. Sensitivity analyses were performed to evaluate the robustness of estimates of B _max to the algorithm used to estimate flow lengths and the prevalence of edge effects resulting from inadequate DEM extent. The findings suggest that the index is insensitive to the specific flow algorithm used but that edge effects can result in significant underestimation along major divides. Edge contamination can, however, be avoided by using an appropriately extensive DEM.     

On the origins of GEMS grains

From their birth as condensates in the outflows of oxygen-rich evolved stars, processing in interstellar space, and incorporation into disks around new stars, amorphous silicates predominate in most astrophysical environments. Amorphous silicates were a major building block of our Solar System and are prominent in infrared spectra of comets. Anhydrous interplanetary dust particles (IDPs) thought to derive from comets contain abundant amorphous silicates known as GEMS (glass with embedded metal and sulfides) grains. GEMS grains have been proposed to be isotopically and chemically homogenized interstellar amorphous silicate dust. We evaluated this hypothesis through coordinated chemical and isotopic analyses of GEMS grains in a suite of IDPs to constrain their origins. GEMS grains show order of magnitude variations in Mg, Fe, Ca, and S abundances. GEMS grains do not match the average element abundances inferred for ISM dust containing on average, too little Mg, Fe, and Ca, and too much S. GEMS grains have complementary compositions to the crystalline components in IDPs suggesting that they formed from the same reservoir. We did not observe any unequivocal microstructural or chemical evidence that GEMS grains experienced prolonged exposure to radiation.We identified four GEMS grains having O isotopic compositions that point to origins in red giant branch or asymptotic giant branch stars and supernovae. Based on their O isotopic compositions, we estimate that 1-6% of GEMS grains are surviving circumstellar grains. The remaining 94-99% of GEMS grains have O isotopic compositions that are indistinguishable from terrestrial materials and carbonaceous chondrites. These isotopically solar GEMS grains either formed in the Solar System or were completely homogenized in the interstellar medium (ISM). However, the chemical compositions of GEMS grains are extremely heterogeneous and seem to rule out this possibility. Based on their solar isotopic compositions and their non-solar elemental compositions we propose that most GEMS grains formed in the nebula as late-stage non-equilibrium condensates.     

Heavitree Quartzite,a Neoproterozoic (ca 800–760 Ma), high‐energy,tidally influenced,ramp association,Amadeus Basin,central Australia

The Neoproterozoic Heavitree Quartzite is widespread in the Amadeus Basin and has correlatives in all of the major central Australian intracratonic basins. The origin of the formation is enigmatic, not only because of its widespread sheet‐like distribution and uniformity of composition, but also because intense silicification makes facies studies difficult. Recently discovered exposures at the eastern end of the basin are relatively free of diagenetic quartz allowing a detailed study of sedimentary structures and an understanding of the depositional architecture of the formation. The formation, which consists largely of pale‐tan or white quartzose sandstone interbedded with rare laminated mudstone and conglomerate intervals, was deposited in at least four depositional sequences. The sheet‐like nature of the sandstone results from an abundant supply of sediments deposited in a high‐energy, open, shelf‐like environment on a regionally subsiding, low‐gradient ramp. Environmental settings switched both laterally and temporally between sand waves deposited by reversing tidal flow and higher velocity unidirectional currents involving dunes and plane beds. In the early stages of deposition, mud‐dominated, tidal‐flat environments alternated with higher energy, sand‐dominated, tidally influenced settings. However, in the later stages of deposition a major eustatic sea‐level fall moved base‐level basinwards, earlier sediments were reworked by streams to form a ravinement surface, gravel was carried well into the basin and fines largely disappeared from the environment. Gravel deposition was followed by a return to high‐energy, tidally influenced deposits involving large sand waves or dunes. Towards the top of the formation sand waves deposited by reversing tidal currents gradually decline and are eventually replaced by dunes deposited by unidirectional current flow. The transition to the shallow‐marine, anoxic rocks of the Bitter Springs Formation is gradational in response to increased accommodation in a ramp setting which lacked a clearly defined shelf break. The Heavitree Quartzite was probably deposited as a direct response to the events surrounding the assembly and breakup of Rodinia, in particular peneplanation during regional uplift in response to a rising mantle plume followed by broad regional subsidence as the plume decayed prior to the breakup of the supercontinent. The large supply of quartz sand resulted from peneplanation associated with the rising plume and the lack of soil‐stabilising vascular plants, an environmental setting with no modern analogue. The ultimate disposition of fines is not known but, given the environment of deposition, it is likely that they were removed during peneplanation and bypassed the sag basin completely.     

Climate Change and the Disappearance of Malaria from England

As the world's climate continues to change there is concern that this may cause malaria to spread to new areas. Here we examine whether past changes in temperature, in addition to social changes, may have affected malaria in England. Our analysis indicates that cold summers experienced during the 1800s may have contributed to the disappearance of this disease from the country. As our summers become progressively warner it is unlikely that malaria will become firmly re-established on these shores, although the possibility of future outbreaks should not be ignored.     

Mineralogy,petrology, chronology,and exposure history of the Chelyabinsk meteorite and parent body

Three masses of the Chelyabinsk meteorite have been studied with a wide range of analytical techniques to understand the mineralogical variation and thermal history of the Chelyabinsk parent body. The samples exhibit little to no postentry oxidation via Mössbauer and Raman spectroscopy indicating their fresh character, but despite the rapid collection and care of handling some low levels of terrestrial contamination did nonetheless result. Detailed studies show three distinct lithologies, indicative of a genomict breccia. A light‐colored lithology is LL5 material that has experienced thermal metamorphism and subsequent shock at levels near S4. The second lithology is a shock‐darkened LL5 material in which the darkening is caused by melt and metal‐troilite veins along grain boundaries. The third lithology is an impact melt breccia that formed at high temperatures (~1600 °C), and it experienced rapid cooling and degassing of S₂ gas. Portions of light and dark lithologies from Chel‐101, and the impact melt breccias (Chel‐102 and Chel‐103) were prepared and analyzed for Rb‐Sr, Sm‐Nd, and Ar‐Ar dating. When combined with results from other studies and chronometers, at least eight impact events (e.g., ~4.53 Ga, ~4.45 Ga, ~3.73 Ga, ~2.81 Ga, ~1.46 Ga, ~852 Ma, ~312 Ma, and ~27 Ma) are clearly identified for Chelyabinsk, indicating a complex history of impacts and heating events. Finally, noble gases yield young cosmic ray exposure ages, near 1 Ma. These young ages, together with the absence of measurable cosmogenic derived Sm and Cr, indicate that Chelyabinsk may have been derived from a recent breakup event on an NEO of LL chondrite composition.     

Arsenic attenuation in tailings at a former Cu–W–As mine,SW Finland

Nearly half a century after mine closure, release of As from the Ylöjärvi Cu–W–As mine tailings in groundwater and surface water run-off was observed. Investigations by scanning electron microscopy (SEM), electron microprobe analysis (EMPA), synchrotron-based micro-X-ray diffraction (μ-XRD), micro-X-ray absorption near edge structure (μ-XANES) and micro-extended X-ray absorption fine structure (μ-EXAFS) spectroscopy, and a sequential extraction procedure were performed to assess As attenuation mechanisms in the vadose zone of this tailings deposit. Results of SEM, EMPA, and sequential extractions indicated that the precipitation of As bearing Fe(III) (oxy)hydroxides (up to 18.4 wt.% As₂O₅) and Fe(III) arsenates were important secondary controls on As mobility. The μ-XRD, μ-XANES and μ-EXAFS analyses suggested that these phases correspond to poorly crystalline and disordered As-bearing precipitates, including arsenical ferrihydrite, scorodite, kaňkite, and hydrous ferric arsenate (HFA). The pH within 200 cm of the tailings surface averaged 5.7, conditions which favor the precipitation of ferrihydrite. Poorly crystalline Fe(III) arsenates are potentially unstable over time, and their transformation to ferrihydrite, which contributes to As uptake, has potential to increase the As adsorption capacity of the tailings. Arsenic mobility in tailings pore water at the Ylöjärvi mine will depend on continued arsenopyrite oxidation, dissolution or transformation of secondary Fe(III) arsenates, and the As adsorption capacity of Fe(III) (oxy)hydroxides within this tailings deposit.