Permeability and degassing of dome lavas undergoing rapid decompression: An experimental determination

The gas permeability of volcanic rocks may influence various eruptive processes. The transition from a quiescent degassing dome to rock failure (fragmentation) may, for example, be controlled by the rocks permeability, in as much as it affects the speed by which a gas overpressure in vesicles is reduced in response to decompression. Using a modified shock-tube-based fragmentation bomb (Alidibirov and Dingwell 1996a,b; Spieler et al. 2003a), we have measured unsteady-state permeability at a high initial pressure differential. Following sudden decompression above the rock cylinder, pressurized gas flows through the sample. Two pressure transducers record the pressure signals above and below the sample. A transient 1D filtration code has been developed to calculate permeability using the experimental decay curve of the lower pressure transducer. Additionally an analytical steady-state method to achieve permeability is presented as an alternative to swiftly predict the sample permeability in a sufficiently precise manner. Over 100 permeability measurements have been performed on samples covering a wide range of porosity. The results show a general positive relationship between porosity and permeability with a high data scatter. Our preferred interpretation of the results is a combination of two different, but overlapping effects. We propose that at low porosities, gas escape occurs predominantly through microcracks or elongated micropores and therefore could be described by simplified forms of Kozeny–Carman relations (Carman 1956) and fracture flow models. At higher porosities, the influence of vesicles becomes progressively stronger as they form an increasingly connected network. Therefore, a model based on the percolation theory of fully penetrable spheres is used, as a first approximation, to describe the permeability-porosity trend. In the data acquired to date it is evident, that in addition to the porosity control, the samples bubble size, shape and distribution strongly influence the permeability. This leads to a range of permeability values up to 2.5 orders of magnitude at a given porosity.     

A magnetotelluric study of the Damara Belt in Namibia: 1. Regional scale conductivity anomalies

The Namibian margin is dominated by the late Proterozoic to early Cambrian fold belts of the Damara Orogen, which wrap around and separate the Congo and Kalahari Cratons. This mosaic of relatively ‘soft’ fold belts and ‘hard’ cratons apparently controlled the path for the opening of the South Atlantic in the early Cretaceous. The continents split along the coast-parallel fold belts of the Damara Orogen while the inland fold belt (Damara Belt) was effected by extension and widespread igneous intrusion but never developed to the rift stage. This paper is concerned with the interpretation of magnetotelluric (MT) data along a 200 km NW-SE profile across the Damara Belt in NW Namibia. The regional, two-dimensional electrical resistivity model and the induction vector data exhibit three distinctive zones: (i) a generally very resistive upper crust which is typical for the granites and metasediments of the Damara Belt, (ii) two subvertical conductors in upper to mid-crustal levels which correlate with major tectonic zone boundaries and (iii) a highly conductive middle to lower crust in the southern part of the profile. The geometry of the conductive structures could reflect a regional shear system in which upper crustal listric faults pass into a detachment zone in the middle crust. We interpret the high electrical conductivity in terms of graphite (or other forms of mineralization) enrichment along the shear planes. This zone of crustal weakness may have originated in Pre-Damara times and had probably experienced several episodes of crustal reactivation before the intrusion of basaltic dike swarms during the Cretaceous rifting and magmatism associated with the opening of the South Atlantic.     

Bofedal response to climate variability,local management,and water extraction:A case study of Chucuyo,Northern Chile

Andean high-altitude wetlands are important ecosystems that serve a range of social-ecological functions.In the Andes,bofedales,a specific type of peat-producin...     

Quantifying groundwater–surface water interactions in a proglacial valley,Cordillera Blanca,Peru

A myriad of downstream communities and industries rely on streams fed by both groundwater discharge and glacier meltwater draining the Cordillera Blanca, Northern Peruvian Andes, which contains the highest density of glaciers in the tropics. During the dry season, approximately half the discharge in the region's proglacial streams comes from groundwater. However, because of the remote and difficult access to the region, there are few field methods that are effective at the reach scale to identify the spatial distribution of groundwater discharge. An energy balance model, Rhodamine WT dye tracing, and high‐definition kite‐borne imagery were used to determine gross and net groundwater inputs to a 4‐km reach of the Quilcay River in Huascaran National Park, Peru. The HFLUX computer programme ( http://hydrology.syr.edu/hflux.html ) was used to simulate the Quilcay River's energy balance using stream temperature observations, meteorological measurements, and kite‐borne areal photography. Inference from the model indicates 29% of stream discharge at the reach outlet was contributed by groundwater discharge over the study section. Rhodamine WT dye tracing results, coupled with the energy balance, show that approximately 49% of stream water is exchanged (no net gain) with the subsurface as gross gains and losses. The results suggest that gross gains from groundwater are largest in a moraine subreach but because of large gross losses, net gains are larger in the meadow subreaches. These insights into pathways of groundwater–surface water interaction can be applied to improve hydrological modelling in proglacial catchments throughout South America. Copyright © 2016 John Wiley & Sons, Ltd.     

Calculations of the seasonal loadings of North Sea sediments with persistent aromatic organochlorines

The content and distribution of ΣPCB (as the sum of 24 individual components), HCB and p,p′-DDE was determined in sediments taken from the North Sea during two seasons. The mass-balance of aromatic organochlorine compounds in the upper 2 cm of sediments is calculated and dicussed in relation to pollutant dispersal.     

Cosmogenic isotope burial dating of fluvial sediments from the Lower Rhine Embayment,Germany

Cosmogenic isotope burial dating, using ¹⁰Be and ²⁶Al, was applied to Plio–Pleistocene fluvial successions from the Lower Rhine Embayment, Germany. The approach consists of three principal steps: (1) measurement of cosmogenic nuclides in depth profiles, (2) modelling of hypothetical nuclide concentrations based on a first-order conceptualisation of the geological context and the principal succession of depositions and subsequent erosional and burial phases, and (3) using parameter estimation to identify values for the a priori unknown model parameters (burial age, initial nuclide concentrations, terrace erosion rates) that result in minimal disagreement between hypothetical and measured nuclide concentrations.The Late Pliocene Kieseloolite Formation was dated to 3650 ± 1490 ka and the Early Pleistocene Waalre Formation to 900 ± 280 ka. The unconformably overlying Upper Terrace Formation revealed ages of 740 ± 210 ka and 750 ± 250 ka for the two different sites. These findings are in good agreement with independent age control derived by bio-, magneto-, and litho-stratigraphy. Furthermore, identifiability and uncertainty analysis reveal a relationship between burial depth and sensitivity of isotope concentrations to burial age and erosion rate. Our results indicate that using shallower buried samples would enable a considerably more robust estimation of the burial age and the terrace erosion rate. Uncertainties arose mainly from nuclide measurements, and not from the uncertainties derived from modelling or insufficient knowledge of nuclide production and decay properties.     

Rare earth abundances in Deccan trap basalts

Twenty one basaltic flows of known mineralogical and chemical composition from the Deccan of western India have been analyzed for nine rare-earth elements (REE), together with Hf, Ta and Th by instrumental neutron activation methods. The flows, which are tholeiites and transitional alkali basalts, are relatively enriched in the Light REE and have total REE abundances systematically related to their bulk petrochemistry. The results are compared with REE abundances of basalts from other provinces.     

Application of the EnKF and Localization to Automatic History Matching of Facies Distribution and Production Data

The performance of the ensemble Kalman filter (EnKF) for continuous updating of facies location and boundaries in a reservoir model based on production and facies data for a 3D synthetic problem is presented. The occurrence of the different facies types is treated as a random process and the initial distribution was obtained by truncating a bi-Gaussian random field. Because facies data are highly non-Gaussian, re-parameterization was necessary in order to use the EnKF algorithm for data assimilation; two Gaussian random fields are updated in lieu of the static facies parameters. The problem of history matching applied to facies is difficult due to (1) constraints to facies observations at wells are occasionally violated when productions data are assimilated; (2) excessive reduction of variance seems to be a bigger problem with facies than with Gaussian random permeability and porosity fields; and (3) the relationship between facies variables and data is so highly non-linear that the final facies field does not always honor early production data well. Consequently three issues are investigated in this work. Is it possible to iteratively enforce facies constraints when updates due to production data have caused them to be violated? Can localization of adjustments be used for facies to prevent collapse of the variance during the data-assimilation period? Is a forecast from the final state better than a forecast from time zero using the final parameter fields?To investigate these issues, a 3D reservoir simulation model is coupled with the EnKF technique for data assimilation. One approach to enforcing the facies constraint is continuous iteration on all available data, which may lead to inconsistent model states, incorrect weighting of the production data and incorrect adjustment of the state vector. A sequential EnKF where the dynamic and static data are assimilated sequentially is presented and this approach seems to have solved the highlighted problems above. When the ensemble size is small compared to the number of independent data, the localized adjustment of the state vector is a very important technique that may be used to mitigate loss of rank in the ensemble. Implementing a distance-based localization of the facies adjustment appears to mitigate the problem of variance deficiency in the ensembles by ensuring that sufficient variability in the ensemble is maintained throughout the data assimilation period. Finally, when data are assimilated without localization, the prediction results appear to be independent of the starting point. When localization is applied, it is better to predict from the start using the final parameter field rather than continue from the final state.