Modeling fracture porosity development using simple growth laws

A model of porosity development has been developed to investigate general relationships between simple fracture aperture growth laws and fracture porosity in evolved fracture arrays in aquifers. The growth of fracture apertures in two-dimensional orthogonal arrays with initially spatially uncorrelated lognormal aperture distributions has been studied, where aperture growth rate is proportional to an exponent of the flow rate through each fracture. The evolved arrays show geometrical phase changes as a function of the aperture growth rate exponent, e, and the standard deviation of the initial aperture distribution, sigma(z). Low values of e and sigma(z) lead to bimodal aperture distributions, where apertures parallel to flow are preferentially enlarged. At moderate values of e and sigma(z), there is a transition to a regime of more complex geometries consisting of networks of channel-like structures of preferentially enlarged apertures. At larger values of e, array-spanning channel-like paths of preferentially enlarged apertures develop, where the tortuosity of the channel-like paths is a linear function of sigma(z). Following an initial growth phase, during which dynamically stable aperture configurations develop, arrays undergo simple amplification. The geometry of the evolved aperture fields is diverse and they can be highly complex; consequently, parameterization and prediction of their evolution in terms of the initial aperture distributions and growth rate laws is not trivial.     

Biogeochemistry of sulfur and iron in Thioploca-colonized surface sediments in the upwelling area off central chile

The biogeochemistry of sedimentary sulfur was investigated on the continental shelf off central Chile at water depths between 24 and 88 m under partial influence of an oxygen minimum zone. Dissolved and solid iron and sulfur species, including the sulfur intermediates sulfite, thiosulfate, and elemental sulfur, were analyzed at high resolution in the top 20 cm. All stations were characterized by high rates of sulfate reduction, but only the sediments within the Bay of Concepción contained dissolved sulfide. Due to advection and/or in-situ reoxidation of sulfide, dissolved sulfate was close to bottom water values. Whereas the concentrations of sulfite and thiosulfate were mostly in the submicromolar range, elemental sulfur was by far the dominant sulfur intermediate. Although the large nitrate- and sulfur-storing bacteria Thioploca were abundant, the major part of S⁰ was located extracellularly. The distribution of sulfur species and dissolved iron suggests the reaction of sulfide with FeOOH as an important pathway for sulfide oxidation and sulfur intermediate formation. This is in agreement with the sulfur isotope composition of co-existing elemental sulfur and iron monosulfides. In the Bay of Concepción, sulfur isotope data suggest that pyrite formation proceeds via the reaction of FeS with polysulfides or H₂S. At the shelf stations, on the other hand, pyrite was significantly depleted in ³⁴S relative to its potential precursors FeS and S⁰. Isotope mass balance considerations suggest further that pyritization at depth includes light sulfide, potentially originating from bacterial sulfur disproportionation. The δ³⁴S-values of pyrite down to −38‰ vs. V-CDT are among the lightest found in organic-rich marine sediments. Seasonal variations in the sulfur isotope composition of dissolved sulfate indicated a dynamic non-steady-state sulfur cycle in the surface sediments. The ¹⁸O content of porewater sulfate increased with depth at all sites compared to the bottom water composition due to intracellular isotope exchange reactions during microbial sulfur transformations.     

The bailer test: a simple effective pumping test for assessing borehole success

A new pumping test method has been designed around the practical requirements of those working in rural water-supply projects in developing countries. The bailer test needs only simple equipment and can be completed in under 1 h. The test involves removing 20–50 bails from a borehole over a 10-min period and then measuring the recovery. The test is analysed using large-diameter-well analysis (which accounts for well storage) and is appropriate for testing low productivity aquifers (transmissivity 0.1–10 m² d^–1), where water levels are shallow (<20 m depth). The bailer test was developed and trialled in a rural water supply project in Nigeria where it was found to predict similar transmissivity to 5-h constant rate tests (r²=0.9). Analysis of the test was further simplified to provide guidelines for field staff such as community health workers. The likelihood of a borehole sustaining a handpump for 250 people can be indicated by measuring the maximum drawdown and time for 50 and 75% recovery from a bailer test and comparing to a simple table. This simplified test is now being used in this and other rural water-supply projects, and has been modified to indicate whether a borehole can sustain higher yields for small-scale irrigation.     

Nannofossils: the smoking gun for the Canarian hotspot

The origin of volcanism in the Canary Islands has been a matter of controversy for several decades. Discussions have hinged on whether the Canaries owe their origin to seafloor fractures associated with the Atlas Mountain range or to an underlying plume or hotspot of superheated mantle material. However, the debate has recently come to a conclusion following the discovery of nannofossils preserved in the products of the 2011–2012 submarine eruption at El Hierro, which tell us about the age and growth history of the western‐most island of the archipelago. Light coloured, pumice‐like ‘floating rocks’ were found on the sea surface during the first days of the eruption and have been shown to contain fragments of pre‐island sedimentary strata. These sedimentary rock fragments were picked up by ascending magma and transported to the surface during the eruption, and remarkably retained specimens of pre‐island Upper Cretaceous to Pliocene calcareous nannofossils (e.g. coccolithophores). These marine microorganisms are well known biostratigraphical markers and now provide crucial evidence that the westernmost and youngest island in the Canaries is underlain by the youngest sediment relative to the other islands in the archipelago. This finding supports an age progression for the onset of volcanism at the individual islands of the archipeligo. Importantly, as fracture‐related volcanism is known to produce non‐systematic age‐distributions within volcanic alignments, the now‐confirmed age progression corroberates to the relative motion of the African plate over an underlying mantle plume or hotspot as the cause for the present‐day Canary volcanism.     

Defining trace-element alteration halos to skarn deposits hosted in heterogeneous carbonate rocks: Case study from the Cu–Zn Antamina skarn deposit,Peru

The presence of geochemical anomalies, defining haloes around hydrothermal ore deposits, can be used to vector towards mineralization, or identify ore bodies buried at depth. Several important types of ore deposits, including skarn deposits, are often hosted within carbonate-rich sedimentary rocks. Identifying anomalous trace-element concentrations in carbonate rocks is complicated by variable lithology (i.e. siliciclastic component) and volume loss during hydrothermal alteration. In this study of the world-class Antamina skarn deposit in Peru, we use the ratio of metals:immobile elements (e.g. La, Al₂O₃) to differentiate genuine and false geochemical anomalies in limestones and marbles surrounding the skarn deposit. Unaltered limestones are used to define threshold values for metal:immobile element ratios (through use of the median value ± 2 median absolute deviations). Genuine anomalies are identified when metal concentrations exceed those predicted using median + 2 median absolute deviations. In addition, comparison of “four acid” and lithium-borate fusion analytical techniques reveals that the lower cost four-acid techniques give reliable results. Our approach can be used to identify geochemical anomalies and halos related to hydrothermal alteration of carbonate-rich rocks, which have direct application to skarn deposits.     

Natural iodine in a clay formation: Implications for iodine fate in geological disposals

Iodine is one of the most problematic radioisotopes in the context of nuclear waste geological disposal due to its high mobility. Considerable effort has been dedicated to the measurement of its potential retardation during diffusive transport leading to conflicting results, from no retardation to significant retardation, leading in turn to considerable debate. The present study aims at providing new insights into this aspect of the iodine problem by careful quantification of iodine reservoirs in the Callovian-Oxfordian (COx) clay rock taken here as model material for these studies. The present study confirmed the ubiquitous presence of iodine at 1-5 mg kg⁻¹ level in the COx clayey formation. The iodide concentration level in the porewater is also confirmed at a value in the range ∼20-40 μmol L⁻¹, i.e. higher than the expected range of radio-iodine concentration in the far-field of the storage. Surprisingly, most of the iodine was found not to be associated with organic matter but rather in an inorganic form associated with carbonate minerals. This result has potentially significant implications for the fate of radio-iodine. In undisturbed far-field conditions, most natural iodine would not be accessible for isotopic exchange with radioactive iodine, reducing the effective Kd to negligible values. During laboratory experiments, good monitoring of the geochemical parameters (at least the Eh, pH, P_CO2, [Ca] and [Mg]) is mandatory to avoid iodine-bearing carbonate precipitation and to enable rigorous interpretation of the iodide diffusion/retention experiments.     

Effective anisotropic velocity model from surface monitoring of microseismic events

We develop a methodology to obtain a consistent velocity model from calibration shots or microseismicity observed on a buried array. Using a layered 1D isotropic model derived from checkshots as an initial velocity model, we invert P‐wave arrival times to obtain effective anisotropic parameters with a vertical axis of symmetry (VTI). The nonlinear inversion uses iteration between linearized inversion for anisotropic parameters and origin times or depths, which is specific to microseismic monitoring. We apply this technique to multiple microseismic events from several treatments within a buried array. The joint inversion of selected events shows a largely reduced RMS error indicating that we can obtain robust estimates of anisotropic parameters, however we do not show improved source locations. For joint inversion of multiple microseismic events we obtained Thomsen anisotropic parameters ε of 0.15 and δ of 0.05, which are consistent with values observed in active seismic surveys. These values allow us to locate microseismic events from multiple hydraulic fracture treatments separated across thousands of metres with a single velocity model. As a result, we invert the effective anisotropy for the buried array region and are able to provide a more consistent microseismicity mapping for past and future hydraulic fracture stimulations.     

Contrasting pyroclast density spectra from subaerial and submarine silicic eruptions in the Kermadec arc: implications for eruption processes and dredge sampling

Pyroclastic deposits from four caldera volcanoes in the Kermadec arc have been sampled from subaerial sections (Raoul and Macauley) and by dredging from the submerged volcano flanks (Macauley, Healy, and the newly discovered Raoul SW). Suites of 16–32?mm sized clasts have been analyzed for density and shape, and larger clasts have been analyzed for major element compositions. Density spectra for subaerial dry-type eruptions on Raoul Island have narrow unimodal distributions peaking at vesicularities of 80–85%, whereas ingress of external water (wet-type eruption) or extended timescales for degassing generate broader distributions, including denser clasts. Submarine-erupted pyroclasts show two different patterns. Healy and Raoul SW dredge samples and Macauley Island subaerial-emplaced samples are dominated by modes at ~80–85%, implying that submarine explosive volcanism at high eruption rates can generate clasts with similar vesicularities to their subaerial counterparts. A minor proportion of Healy and Raoul SW clasts also show a pink oxidation color, suggesting that hot clasts met air despite 0.5 to >1?km of intervening water. In contrast, Macauley dredged samples have a bimodal density spectrum dominated by clasts formed in a submarine-eruptive style that is not highly explosive. Macauley dredged pyroclasts are also the mixed products of multiple eruptions, as shown by pumice major-element chemistry, and the sea-floor deposits reflect complex volcanic and sedimentation histories. The Kermadec calderas are composite features, and wide dispersal of pumice does not require large single eruptions. When coupled with chemical constraints and textural observations, density spectra are useful for interpreting both eruptive style and the diversity of samples collected from the submarine environment.     

Physico‐chemical features of lake Pupuke,Auckland

Summary

Lake Pupuke, situated in the North Island of New Zealand near Auckland (Latitude 36° 53'S, Longitude 174° 52'E), is a small, circular body of water (maximum width 1.2 km, maximum depth 55 m). From November 1966 to November 1967 observations on its physico‐chemical factors were made twice monthly.

The range of temperatures in the year was 12.0–24.2°c at the surface and 11.3–14.5°c at 40 m. A thermocline farmed between early October and early January, when it became strongly developed at about 17 m and persisted until late April. Between early May and late June the stratification was breaking down and between mid July and mid September the lake was isothermal.

The values for light transmitted per metre through the water were high from mid August to mid December (53–76% per metre), very low from mid December until mid February (23–49% per metre), high in March (75% per metre) and relatively constant over the winter months (53–64% per metre). The amount of light transmitted depended closely on the density of phytoplankton. Changes in the depth at which a Secchi disc disappeared paralleled the light transmission values, and the mean Secchi depth was 2.2 m, range 1.0–5.2 m.

Oxygen concentrations down to 30 m were highest at all depths from mid November to early January (112–134% saturation). From early January to late April oxygen was strongly stratified (epilimnion 55–100% saturation, hypolimnion 24–77% saturation). During winter oxygen values gradually increased at all depths.

The range of pH was 7.8–9.5 units at the surface and 7.8–8.6 units at 30 m. From early January to late April pH values were strongly stratified but from mid July to mid August they were nearly uniform. The range of alkalinity (mainly HCO₃ ^? ions) was 76.0–80.8 ppm CaCO₃ and values showed stratification in summer. Carbonate ion was present from early October to mid June, but free carbon dioxide was only found in trace amounts.

The range of nitrate ion concentrations was 0.05–0.31 ppm; nitrate became obviously depleted in summer. Only one nitrite value was significant, 0.25 ppm. other nitrite values never exceeded 0.10 ppm. The concentration of ammoniacal nitrogen built up to a summer maximum of 0.95 ppm in the hypolimnion (annual range 0–0.95 ppm); values for albuminoid ammonia (0.016–1.03 ppm of nitrogen) varied with the quantity of living matter present.

The annual mean for total ion concentration (3.30 meq./1) was high and the order of ions was Na⁺>Mg²⁺>Ca²⁺>K⁺ and HCO₃ ^?>CI^?>SO₄ ^2?. Such dissolved salts may have arrived as airborne particles, dissolved in rain, or as sea spray.     

Variation in sunscreen compounds (mycosporine‐like amino acids) for marine species along a gradient of ultraviolet radiation transmission within doubtful sound,New Zealand

We examined the response of four species of New Zealand marine algae (Ecklonia radiata, Apophlaea lyallii, Rhodymenia spp., Ulva lactuca) and a sea urchin (Evechinus chloroticus) to spatial variation in ultraviolet radiation (UV‐R) by examining the concentration of UV‐R absorbing compounds known as mycosporine‐like amino acids (MAAs). The purpose was to understand how, and the degree to which, local marine species could potentially respond to any future increases in incident UV‐R in the New Zealand marine environment. The research was undertaken in Doubtful Sound, where we observed a gradient of water column UV‐R transmission along the 40 km length of the fiord. We examined spatial differences in MAAs along the UV‐B gradient in the macrophytes and temporal changes in MAAs in sea urchin gonads. Among the algae, thallus MAA concentrations (nmol mg^–1 protein) ranged from 12.5 to 87.8 in E. radiata, from 433.1 to 1446.4 in A. lyallii, 12.7 to 103.4 in Rhodymenia spp., but were not detected in U. lactuca. For E. chloroticus, gonadal MAA concentrations ranged from 83.9 to 224.3 nmol mg^–1 protein spatially, and over the year from 1.85 to 14.12 nmol mg^–1 dry weight (DW) depending on site and gametogenic cycle. Laboratory manipulations indicated that concentrations of MAAs in E. chloroticus gonads and eggs are influenced by diet. MAA concentration could be correlated with UV‐B intensities in two of the algal species. E. chloroticus MAA concentrations could also be correlated with UV‐B transmission, which we concluded was a reflection of the greater ingestion and accumulation of MAA‐rich macrophytes at those sites where higher ambient UV‐R induced greater MAA concentrations to occur in the algae. Given this, we suggest that one response of marine species to increases in UV‐B would be an increase in the synthesis and/or accumulation of MAAs for photoautotrophs and a dietary accumulation of those MAAs in E. chloroticus, an important herbivore in this system.