Original Research ARTICLE
Radiolytic Hydrogen Production in the Subseafloor Basaltic Aquifer
- 1Graduate School of Oceanography, University of Rhode Island, Narragansett, RI, USA
- 2Department of Earth and Environment, Boston University, Boston, MA, USA
Hydrogen (H2) is produced in geological settings by dissociation of water due to radiation from radioactive decay of naturally occurring uranium (238U, 235U), thorium (232Th) and potassium (40K). To quantify the potential significance of radiolytic H2 as an electron donor for microbes within the South Pacific subseafloor basaltic aquifer, we use radionuclide concentrations of 43 basalt samples from IODP Expedition 329 to calculate radiolytic H2 production rates in basement fractures. The samples are from three sites with very different basement ages and a wide range of alteration types. U, Th, and K concentrations vary by up to an order of magnitude from sample to sample at each site. Comparison of our samples to each other and to the results of previous studies of unaltered East Pacific Rise basalt suggests that significant variations in radionuclide concentrations are due to differences in initial (unaltered basalt) concentrations (which can vary between eruptive events) and post-emplacement alteration. However, there is no clear relationship between alteration type and calculated radiolytic yields. Local maxima in U, Th, and K produce hotspots of H2 production, causing calculated radiolytic rates to differ by up to a factor of 80 from sample to sample. Fracture width also greatly influences H2 production, where microfractures are hotspots for radiolytic H2 production. For example, H2 production rates normalized to water volume are 190 times higher in 1 μm wide fractures than in fractures that are 10 cm wide. To assess the importance of water radiolysis for microbial communities in subseafloor basaltic aquifers, we compare electron transfer rates from radiolysis to rates from iron oxidation in subseafloor basalt. Radiolysis appears likely to be a more important electron donor source than iron oxidation in old (>10 Ma) basement basalt. Radiolytic H2 production in the volume of water adjacent to a square cm of the most radioactive SPG basalt may support as many as 1500 cells.
The oceanic basement contains the largest aquifer on Earth. Its fractured rock contains nearly 2% of Earth’s total volume of seawater (Johnson and Pruis, 2003). Although the extent of life and microbial activity in oceanic basement is not well known, a variety of evidence suggests that microbes reside within the aquifer (Cowen et al., 2003; Edwards et al., 2012; Jungbluth et al., 2013; Lever et al., 2013; Orcutt et al., 2013). Fisk et al. (1998) and Staudigel et al. (2008) report weathering textures suggestive of microbial alteration in subseafloor basaltic glass. Microorganisms have been found in fluid flowing through 3.5 million year old basalt on the Juan de Fuca ridge flank (Cowen et al., 2003). DNA and isotopic signatures of mineral alteration provide evidence of microbes and microbial activity in ridge-flank basalt (Lever et al., 2013).
Physical and chemical properties limit microbial habitability of the oceanic basement. For example, habitability within subseafloor basalt is constrained by availability of electron donors [e.g., organic carbon, ferrous iron (Fe2+), and hydrogen (H2)] and electron acceptors [e.g., oxygen (O2), nitrate (NO3-), and sulfate (SO42-)] (Madigan et al., 2000; Bach and Edwards, 2003; D’Hondt et al., 2004). A number of studies have suggested that microbial life in igneous-rock aquifers may be supported by oxidation of electron donors native to the rock or produced by water–rock interactions (Pedersen, 1993; Stevens and McKinley, 1995; Kelley et al., 2001; Chapelle et al., 2002; Bach and Edwards, 2003; Edwards et al., 2005). For example, basement basalt has high concentrations of reduced elements, specifically iron (Fe) and sulfur (S) (Bach and Edwards, 2003). Oxidation of these elements with O2 or NO3- in the seawater provides energy that microorganisms might utilize. Water–rock interactions that can produce electron donors in the form of molecular H2 include serpentinization (Kelley et al., 2001) and radiolysis of water due to radioactive decay of radionuclides within the rock (Pedersen, 1993; Lin et al., 2005a,b; Blair et al., 2007; D’Hondt et al., 2009, 2015; Edwards et al., 2011; Lollar et al., 2014). Water radiolysis within subseafloor basalt is the focus of this study.
Water radiolysis is the decomposition of water molecules by ionizing radiation produced during the decay of radioactive elements (Debierne, 1914; Le Caër, 2011). The principal radioactive elements that produce ionizing radiation in basalt are uranium (238U and 235U), thorium (232Th), and potassium (40K), which collectively emit alpha (α), beta (β), and gamma (γ) radiation as they and their daughter nuclides decay. Transfer of energy from this radiation excites and ionizes water molecules, producing several chemical species: eaq-, HO•, H•, HO2•, H3O+, OH-, H2O2, and H2 (Spinks and Woods, 1990; Le Caër, 2011). The distribution and rate of formation of these products depends on the linear energy transfer (LET) of the radiation (the amount of energy deposited by the radiation along its path). Low-LET radiation (γ-rays and β particles) ionizes water discretely along the radiation path. High-LET radiation (e.g., α particles) deposits energy densely along the particles’ track. The radiolytic yields of the radicals decline with LET. The radicals are short-lived and highly reactive; they recombine in the radiation track to produce stable decomposition products (H2 and H2O2). Production of H2 and H2O2 increases with increasing LET (Pastina and LaVerne, 2001; Le Caër, 2011).
For subseafloor environments, we are particularly interested in the production of the reductant H2. Many organisms catabolically utilize H2, including methanogens, sulfate-reducers, iron reducers, and nitrate reducers (Fang and Zhang, 2011). There is evidence that some of these organisms, specifically sulfate reducers and methanogens, are active in subseafloor basalt (Lever et al., 2013). Radiolysis undoubtedly occurs in subseafloor basalt, as both water and radiation are present. Edwards et al. (2012) suggested that in the old and relatively weathered basaltic basement of the South Pacific Gyre (SPG), radiolytic H2 may be the dominant electron donor.
The potential of radiolytic H2 as an electron donor has been studied for other environments, such as continental crust (Pedersen, 1993; Lin et al., 2005a,b) and deep-sea sediment (Blair et al., 2007). In these environments, water radiolysis supplies H2 that may support microbial communities. In deep-sea sediment, H2 may be especially important to microbial communities where organic carbon availability is very low (Blair et al., 2007; D’Hondt et al., 2009).
Here, we focus on calculating the rate of radiolytic H2 production in oceanic basement. More specifically, we use our recently published water radiolysis model to calculate H2 production rates in fractures of South Pacific subseafloor basalt (Dzaugis et al., 2015). Our model more accurately quantifies radiolytic H2 production in fractured hard rock then previous models, which are more appropriate for physically homogenous environments where solid particles are small compared to the distance that the ionizing radiation travels, such as deep-sea sediment.
To determine variations in radionuclide concentrations of relatively old SPG basalt, we analyzed basalt samples from three SPG sites with basement ages of circa 13 to 100 Ma. Our samples include a wide variety of alteration types. To identify effects of basement age, initial composition, and alteration on radiolysis rates, we compare our H2 production calculations for these samples to calculations for which we use previously published U, Th, and K data from young basalt samples on the East Pacific Rise (EPR) (Gale et al., 2013). In addition, we quantify the effect of fracture width on the H2 production rates. We assess the possible significance of radiolytic H2 production for microbial communities in subseafloor basaltic aquifers by comparing the electron transfer flux from radiolytic H2 to the electron transfer flux from iron oxidation in the basalt. Finally, we estimate the number of cells that might be supported by radiolysis in the basement.
Materials and Methods
Integrated Ocean Drilling Program (IODP) Expedition 329 collected subseafloor basalt samples at three sites (Sites U1365, U1367, U1368; Figure 1). We analyzed 43 samples from these sites to determine the range of U, Th, and K concentrations in SPG basalt.
FIGURE 1. Site locations. Colored dots mark site locations: Site U1365 (yellow), Site U1367 (orange), Site U1368 (red). The site localities are superimposed on a bathymetry map with the East Pacific Rise in light gray. Black lines indicate basement age (20 Myr contours).
Site U1365 is located in the western part of the SPG, where the seafloor is approximately 5695 m below sea level (Expedition 329 Scientists, 2011a). At this location, basement is mainly composed of massive lava flows that are approximately 100 Ma (Dunlea et al., 2015a). Recovery of cored basalt was 74.6% at this site. Throughout the cored sequence, alteration extent of the recovered basalt varies between 2 and 95%, with alteration mainly associated with veins and vesicles near flow margins and breccias (Expedition 329 Scientists, 2011a). We analyzed 17 samples from Site U1365 for this study.
The second site, U1367, is east of Site U1365, with seafloor approximately 4288 m below sea level (Expedition 329 Scientists, 2011b). The basement is younger at this site, with an approximate age of 33.5 Ma. The uppermost basement is mostly composed of fractured pillow fragments with a small thin flow at the base of the recovered basalt. Due to the fractured nature of the basalt cored at this site, core recovery was low (11.2%). The extent of alteration in recovered basalt from this site varies between 2 and 25%. Most basalt alteration at Site U1367 occurs around vein-filled fractures of chilled margins (Expedition 329 Scientists, 2011b). We analyzed five samples from Site U1367.
Site U1368 is the eastern-most site of Expedition 329. Seafloor here is 3740 m below sea level (Expedition 329 Scientists, 2011c). Basement at Site U1368 is approximately 13.5 Ma. Core recovery was 27.6%. The recovered portion of the basement is dominantly composed of pillow basalt. Basalt alteration at this site is mainly in veins, vesicles, along chilled margins, and within volcaniclastic breccias. The recovered basalt varies in alteration extent from 2 to 60% (Expedition 329 Scientists, 2011c). There are 21 samples from Site U1368 used in this study.
Our 43 samples are characterized by wide ranges of alteration extent and alteration type. We assign samples to the following rock types based on macroscopic visual appearance; brown halos, dark gray halos, carbonate veins, iron oxyhydroxide (Fe) stained, breccia, mixed alteration, and minimally altered. The term brown halo refers to all halos that vary from red to brown in color. Similarly, dark gray halo refers to halos ranging from very dark gray to dark green in color. Fe staining includes samples with iron oxyhydroxide staining of secondary minerals or filling of veins and vesicles. Iron oxyhydroxides leave behind a bright red-orange color. We define samples as mixed-alteration when they are visibly altered but without a dominant single alteration type, such as samples with veins and vesicles filled with carbonate, clays and iron oxyhydroxides and have secondary mineral emplacement. Finally, we classify all aphyric samples without visibly altered regions as minimally altered.
Where possible, we analyzed samples of different visibly altered regions within the same basaltic rock to compare the radioactive element concentrations of the different alteration zones. For example, we took two samples from a piece of light gray basalt with a brown halo. For all analyses of alteration halos and carbonate veins, we separated each alteration type from the rest of the rock. Other alteration categories included a mixture of altered rock and background basalt, as they could not be fully separated. We analyzed the 43 samples for their U, Th, and K concentrations.
U, Th, K Measurements
We used well-established protocols to measure concentrations of U, Th, and K in the 43 samples by inductively coupled plasma-emission spectrometry (ICP-ES) and inductively coupled plasma-mass spectrometry (ICP-MS; Dunlea et al., 2015b). We analyzed these samples at Boston University with a VG PlasmaQuad Excell ICP-MS for U and Th concentrations, and a Jobin-Yvon (JY) Ultima-C ICP-ES for K concentrations. Based on replicate analysis, U and Th precision was 2 and 1%, respectively, of their measured values. The K measurements were within 1% of the measured value. To assess analytical accuracy, we analyzed BHVO-2 Standard Reference Material independently from our calibrations. The measured values agree with the reported accepted values within the analytical precision (Jochum et al., 2005). We list all U, Th, and K data in Supplementary Table S1.
Radiolytic H2 Production Model
Our radiolysis model calculates production rates in water near a radionuclide-containing solid (Dzaugis et al., 2015). Previous models, such as those used in sediment (Blair et al., 2007), are not applicable to basalt, because studies of basalt cannot assume homogenous porosity or grain/crystal size smaller than the stopping distance of the radiation. Dzaugis et al. (2015) gives a detailed description of the model that we apply to our samples. In this section, we describe the inputs and assumptions that we used to calculate radiolytic H2 production rates in fractures.
There are four dominant parent radionuclides in basalt: 238U, 235U, 232Th, and 40K. When these nuclides decay, α, β, and/or γ radiation is emitted. Each type of radiation has different properties that affect radiolytic yield. For example, α particles travel short distances (10s of μm) but have high initial energies and produce the most H2 molecules per unit of energy absorbed. Our calculations assume that the entire decay series of 238U, 235U, and 232Th are in secular equilibrium. This assumption is valid for basalt older then 377,000 years (five half-lives of 230Th; LaTourrette et al., 1993). However, younger basalt is characterized by isotopic disequilibrium within the 238U-decay series, because fractionation during partial melting leads to excess 230Th relative to 238U (LaTourrette et al., 1993). Due to this 230Th enrichment, we slightly underestimate radiolytic H2 production rates when we assume secular equilibrium for young basalt.
The distance that α or β radiation travels before losing all of its kinetic energy is called its stopping distance (Spinks and Woods, 1990). Gamma radiation loses energy exponentially with distance and therefore is not assigned a specific stopping distance. Instead, we use the maximum distance traveled by γ-rays to be 10 half-distances, at which point less then 0.1% of their initial energy remains. Stopping distance and γ travel distance depend on the matrix; for example, a 5 MeV α particle will travel about 20 μm in basalt but 40 μm in water. If radiation is emitted from a radionuclide farther from the water interface than the stopping distance, it does not contribute to water radiolysis. Once radiation reaches the basalt-water interface, it continuously ionizes water along its path until it reaches its stopping distance (Le Caër, 2011) or re-enters basalt from the water. Our model incorporates radioactivity, decay energy of each radionuclide, and how each type of radiation attenuates energy along its path. It sums all of the radiation that is absorbed in a fracture. Consequently, the inputs to our model are U, Th, and K concentration, initial energy for all radiation emitted from the isotopes, energy-range relationships for α, β, and γ radiation in basalt and water, H2 yield per unit energy for each type of radiation, and distance the radiation travels through water. For the calculations in this paper, we assume there is one centimeter of basalt on either side of the fracture, unless noted otherwise.
We used the data in Supplementary Table S1 to calculate the radioactivity of our SPG basalt samples for 238U, 235U, 232Th, and 40K. We used published radionuclide data (Gale et al., 2013) to calculate the radioactivity of EPR basalt at its time of formation. We used the program RadDecay (Hacker, 1997) to find the initial energy for all radiation emitted from radionuclides in basalt. We calculated stopping distances of α and β particles and half-distances of γ-rays for basalt and water using the energy-range data from the ASTAR, ESTAR (Berger et al., 2005), and X-ray Attenuation (Hubbell and Seltzer, 2004) programs in the NIST database. To determine stopping distances in basalt, we used energy-range data for borosilicate glass (in the NIST database, borosilicate glass is the material with electron density most similar to oceanic basalt). Using these relationships, we developed equations to calculate travel distance given an initial energy. Supplementary Table S2 summarizes the equations we use in this study.
Using radioactivity, energy and range data, and how many H2 molecules are produced per 100 eV (G-values), we calculated volume-normalized H2 production rates for basalt fractures of widths between 1 um and 1 m. These rates assume (i) basalt of the same composition on both sides of a fracture and (ii) homogeneous distribution of the radionuclides throughout the basalt.
South Pacific H2 Production Rates
We calculated radiolytic H2 production rates as a function of fracture width for all 43 samples. We show the results of these calculations separately for each site in Figure 2. Volume-normalized H2 production decreases greatly as fracture width increases (from 1μm to 1 m in width; Figure 2). For each width, the differences in H2 yields between samples are due to variations in radioactive element concentrations. Site U1365 contains the oldest and generally most altered basalt of the three sites (Expedition 329 Scientists, 2011a). It exhibits the lowest H2 production rates (Figure 2A), while Site U1368, with the youngest basalt age, generally has the highest rates (Figure 2C). The SPG sample with the lowest H2 production rate (approximately 3X lower then other samples) is from Site U1365 with a Th concentration below the detection limit (0.01 ppm). We calculated the radiolytic H2 production rate for this sample using the below detection limit value.
FIGURE 2. Radiolytic H2 production rates calculated from radionuclide concentrations of our SPG samples. Volume-normalized rates are shown as a function of fracture width for (A) Site U1365, (B) Site U1367, and (C) Site U1368. Black lines mark maximum and minimum rates calculated from EPR glass samples (Gale et al., 2013). Note that both axes are log scale.
Radionuclide Compositional Variance
U and K variance in basalt is due to geochemical composition at the time of basalt formation (initial composition) and post-emplacement alteration. In contrast, Th is not significantly modified during post-emplacement alteration (Kelley et al., 2005). We use Th to constrain abundance variations due to initial composition and U/Th ratios to constrain the magnitude of U alteration.
In our samples, Th ranges from below the detection limit (0.01 ppm) to 1.17 ppm. This range is comparable to the entire range of values observed in unaltered basalt glass from the EPR (compositional data from Gale et al., 2013). The U/Th ratios of our samples range from 0.32 to 4.1 (excluding samples with Th below detection limit). This is much greater than ratios observed in unaltered EPR basaltic glass (Gale et al., 2013). However, it is similar to the range in altered material from other studies (e.g., Kelley et al., 2003). About 1/3 of our samples have U/Th ratios that fall within two standard deviations of the mean U/Th ratio in unaltered EPR glass samples. Most of the SPG samples that fall into this category are from the site with youngest basement age, U1368 (11 samples). Only one and two samples at Site U1367 and U1365, respectively, have U/Th ratios indicative of unaltered composition.
H2 Yields and Alteration
We separately describe the calculated radiolytic rates of the three SPG sites because they have different geographic origins, different eruptive histories, and different degrees of deviation from the mean U/Th ratio of unaltered EPR glass. Site U1365 exhibits a 21-fold range in rates, whereas Sites U1367 and U1368 show 7- to 8-fold ranges. Overall, Site U1365 has the lowest rates, while Site U1368 has the highest (Figure 2).
Our samples show no clear relationship between alteration type and H2 production rates. The H2 production rates calculated for all alteration categories, including minimally altered samples, seem to span the entire range; no alteration type exhibits a significantly different range of H2 production rates than the others (Figure 3). There is some separation of alteration types within the individual sites. At Site U1365, iron stained samples tend to have higher rates, while at Site U1368 the mixed alteration samples and the breccia are highest (Figure 3). Two of the three mixed alteration samples at U1368 have U/Th ratios within the unaltered range indicating that the high range may not completely be due to alteration. Many of the H2 production rates for the different samples overlap; the rates calculated for each sample at the basalt-water interface are given in Supplementary Figure S1.
FIGURE 3. Relation of rock alteration to radiolytic H2 production rates. Volume-normalized rates are given for Sites (A) U1365, (B) U1367, (C) U1368. Colors show the alteration type to which we assigned each sample. Note that both axes are log scale. As in Figure 2, rates are given as a function of fracture width.
To further investigate the impact of alteration, we calculate the U due to alteration (Ualt) as
where the subscript meas refers to measured values and is the average ratio in unaltered EPR basalt (0.37 ± 0.08). We then calculate H2 yield due to Ualt and determine its fractional contribution to the total yield. There are 27 samples that have excess U that is significantly (more then two standard deviations) different than zero. For these samples, the H2 yield based on Ualt ranges from 0 to 4.9 × 10-1 ± 1.5 × 10-2nM H2 yr-1, and can contribute up to 85 ± 3% of the radiolytic H2 produced (average of 31 ± 7% for all samples with altered U/Th ratios).
Influence of Fracture Width and Basalt Width on H2 Production Rates
While compositional variation leads to a large range of radiolytic H2 production rates (almost two orders of magnitude within the SPG samples), fracture width has an even greater effect on volume-normalized H2 production (moles per vol. of water per time). The production rate per volume of water decreases as fracture width increases. This decline in volume-normalized rates is most pronounced after α and β particles run out of energy. To illustrate this effect, we calculated volume-normalized production rates for a range of fracture widths that occur in basement basalt (1 μm to 1 m). Volume-normalized H2 production rates are highest in microfractures (<10 μm), regardless of radionuclide concentration, and strongly decrease as fracture width increases (Figure 2). Production rates in 1 μm wide fractures differ by more than three orders of magnitude from rates in 1 m wide fractures. Volume-normalized H2 production rates are highest at the rock–water interface, due to high dose rates. However, if production rate is normalized to the surface area of fractures, it increases with fracture width as more radiation, especially from γ-rays, is absorbed in wider fractures than narrower fractures.
Radiolytic H2 production rates also vary with the thickness of basalt that abuts a fracture (Figure 4). To illustrate this effect, we calculated production rates based on a single SPG sample and three different thicknesses: 1 m, 1 cm and 100 μm. Thickness affects the amount of radiation emitted to water. One meter of basalt is approximately equivalent to an infinite basalt thickness because less than 0.1% of the radiation travels beyond a meter (10 half-distances of γ-rays); the amount of energy that reaches the water approaches its maximum at about a meter of basalt.
FIGURE 4. Effect of basalt thickness on radiolytic H2 production rates. (A) Volume-normalized radiolytic H2 production for 1-m thick basalt using an average SPG composition of our samples, (B) 1-cm thick basalt, and (C) 100-μm thick basalt. Total production rates are shown as black lines. Rates due to α radiation are shown in blue, β radiation in green and γ radiation in red. X-axis is log scale.
This dependence on basalt thickness is clearly illustrated by comparing H2 production by β and γ radiation with 100 μm thick basalt (Figure 4C) to H2 production by β and γ radiation with 1 cm thick or 1 m thick basalt (Figures 4A,B). The decrease in thickness from 1 m to 100 μm of basalt causes a 20% decrease in β-produced H2 and a 98.9% decrease in γ-produced H2. Because the stopping distance for β and γ radiation is much greater than 100 μm in basalt, the change in thickness greatly impacts the number of β particles and γ-rays emitted from the basalt. Despite these effects of basalt thickness on H2 production by β and γ radiation, total radiolytic H2 production rate is not greatly affected by basalt thickness because α radiation is responsible for most radiolytic H2 production in our examples (Figure 4). Total H2 production rate is only 4% less with 1 cm of basalt and 10% less with 100 μm than 1 m of basalt. In all three cases, α-dose rate is the same, because α particles penetrate no more than 50 μm of basalt (25 μm on average). Therefore, in all three scenarios, all the α radiation that can reach fracture water does.
In this section, we discuss (i) factors that impact radiolytic H2 production rates, (ii) their potential importance to electron donor fluxes in subseafloor basalt, and (iii) the number of cells that might be supported by rates of radiolytic H2 production in representative fractures. We first consider the fractured nature of oceanic basement and how both fracture width and basalt thickness affect radiolysis rates. We address how source-melt composition and subsequent low-temperature alteration play a role in determining the distribution of radionuclides in our sample set and rates of radiolytic H2 production. We then compare electron donor fluxes from radiolytic H2 production and from iron oxidation. Finally, we estimate numbers of cells that might be supported by radiolysis in SPG basaltic fractures, based on comparison to per-cell O2 consumption rates in SPG sediment.
Factors that Influence H2 Production
We focused this study on radiolysis in fractures because most of the water in oceanic basement resides in, and travels through, fractures (Fisher, 1998). The three primary categories of porosity in oceanic crust are (i) macroscopic features, such as lithological boundaries (e.g., surfaces of lava flows), voids associated with pillow basalts, and collapse structures; (ii) microcracks which have narrow widths and limited lateral extent, and (iii) vesicles and other primary porosity, which we do not consider since they are typically disconnected and isolated. Of these three categories, the first is most closely related to formation-scale permeability (Fisher, 1998). Fluid flow in the igneous basement is mostly through interconnected fractures in the oceanic basement (Wang, 2004), e.g., through fractures between pillow basalts and along the margins of lava flows. This relationship between fluid flow and interconnected fractures is apparent in the distribution of alteration halos, which are mainly associated with veins in the basalt, indicating that chemical transport is very limited in the low porosity, impermeable basalt matrix (Bach and Edwards, 2003).
Our focus on quantifying radiolytic rates in fractures differs from previous studies of radiolysis in continental crust (e.g., Lin et al., 2005a), which assume porosity to be homogenously distributed throughout the rock. Such an assumption is very problematic for estimating radiolytic H2 production in hard-rock aquifers, where most porosity is in heterogeneously distributed fractures and much of the radiation emitted from radionuclides in the rock will never reach water. As we discussed in Section “Results,” properties that most affect H2 production rates in hard-rock aquifers include fracture width, rock thickness and rock composition. Of these properties, fracture width and rock composition are the most important.
Fracture width greatly influences radiolytic H2 production rates (Figures 2–4) because volume-normalized H2 production rapidly decreases with distance from the rock-water interface. This decrease in H2 production is due to the limited ranges of α and β radiation and increased volume of water. Of the three types of radiation, α radiation has the highest initial energy and highest G-value at 1.2 molecules H2 per 100 eV (Pastina and LaVerne, 2001). However, it also has the shortest range. Consequently, H2 production by α radiation extends only several tens of microns into fracture water. Beta and gamma radiation have lower G-values (0.6 and 0.45 molecules H2 100 eV-1, respectively) than α particles (Kohan et al., 2013; Mustaree et al., 2014), but β particles can travel 1000s of μm in water and γ-rays travel for tens of cm. Consequently, α radiation dominates total radiolytic H2 production near the rock-water interface, while γ-produced H2 is highest in fractures greater than 1 cm (Figure 4). Due to the very high rate of H2 production from α radiation, volume-normalized H2 production rates in 1 μm fractures are on average 190 times higher than fractures 10 cm in width and 1.6 × 103 times higher than in 1 m fractures. In short, microfractures are hotspots of radiolytic H2 production.
To a much lesser extent, total radiolytic H2 production varies with the thickness of the abutting rock (Figure 4). Over a four-order magnitude of change in basalt thickness (from 100 μm to 1 m of basalt), H2 production rates change by only 10% at the rock–water interface. This change is due to differences in the absorbed doses of β and γ radiation. This change is small because α radiation from deeper in the rock does not penetrate to the surface. If the thickness of basalt facing a fracture is less than 50 μm, the rates will drop more dramatically due to a lower absorbed α-dose rate.
Along with microfractures, high concentrations of U, Th, and K create hotspots of radiolytic H2 production. Uranium, Th, and K concentrations differ widely from sample to sample within and between the three sites. To illustrate the range of H2 yields due to initial composition variance, we show a H2 production range for unaltered basaltic EPR glass samples (Gale et al., 2013). The range of calculated H2 production at the basalt-water interface is approximately 10-fold for the group of EPR samples with Th composition within two standard deviations of the mean (∼90% of the samples) (black lines, Figure 2). The range is much greater (93-fold) if we consider all EPR glass measurements reported by Gale et al. (2013), which include measurements of enriched mid-ocean ridge basalt (E-MORB) that have much higher U and Th concentrations. Basement comprised of E-MORB will have higher radiolytic yields than normal ocean basalt.
The 10-fold difference in H2 yields from EPR data suggests that basaltic source-melt composition has a significant effect on radiolytic rates. Given the wide range of compositions exhibited by the EPR samples, we expect variation in radiolytic H2 production rates from one SPG site to another, as they all have different geographic locations and ages of origin (Figure 1). In addition, the basalt at Site U1365 has a very different eruption history than the basalt at Sites U1367 and U1368. The basalt drilled at Site U1365 was accreted during medium to fast spreading and is comprised mostly of large sheet flows. In contrast, the basalts drilled at Sites U1367 and U1368 are predominantly pillow basalt, and were likely produced during slower spreading than at Site U1365 (Zhang et al., 2013; Zhang and Smith-Duque, 2014).
Site U1365 has lower radiolytic H2 production rates than Sites U1367 and U1368, despite having the oldest and most altered basalt (Figure 2). This can be attributed to low initial concentrations, as indicated by Th values. Thorium concentrations of U1365 samples are all below 0.2 ppm. These concentrations are lower than those in the samples from the other localities and fall into the bottom half of EPR Th data from Gale et al. (2013). Even with enrichment of U during alteration (indicated by high U/Th ratios) at Site U1365, the radiolytic rates are low compared to the other localities indicating that they are greatly affected by initial radioisotope concentrations. This variation in source-melt composition is consistent with the results of Zhang and Smith-Duque (2014), who used differences in initial composition to explain geochemical variance between Sites U1365 and U1368. In addition, there is also likely vertical variation at each site due to changes in source composition between eruptive events (e.g., Bergmanis et al., 2007). Zhang et al. (2013) reported some vertical variation at Site U1368. Using U, Th and K from Zhang et al. (2013) results in a 6X difference in radiolytic rates between their lowest normal-MORB sample and the sample they suggest was influenced by an enriched-MORB source. This difference is similar to the sevenfold range in rates calculated from our radioisotope concentration data for this site.
Concentrations of U and K often increase with alteration (e.g., Staudigel et al., 1996; Teagle et al., 1996; Kelley et al., 2003). We use U/Th ratios to constrain how much excess U a sample contains. We don’t calculate excess values for K, because K/Th ratios in fresh basalt are not consistent enough to calculate excess K in this manner. Potassium is typically concentrated in such alteration minerals as smectite and zeolite, as well as K-feldspar (Alt and Teagle, 2003; Bartetzko, 2005). Shipboard logs of natural gamma radiation (NGR) show clear evidence of such K concentration at alteration fronts in Site U1365 basalt (D’Hondt et al., 2011, 2013). However of the three radionuclides, K contributes the least to total H2 production because it only produces β and γ radiation when it decays. Using an average concentration for the SPG samples, K contributes 13% to total production at the basalt–water interface. Uranium and Th decay series account for 65 and 22% of the H2 produced, respectively. Of our SPG basalt samples, K decay dominates H2 production in only one sample (U1365E-3R-4W 25/30), where it accounts for 48% of the total H2 yield.
The U decay series is overwhelmingly the largest contributor to H2 production. Uranium enrichment often occurs in carbonate veins and at oxidation/reduction fronts during basalt alteration (Alt et al., 1992; Farr et al., 2001; Kelley et al., 2003). Twenty-seven of our SPG samples exhibit U/Th ratios that indicate U enrichment, excluding the sample with Th below the detection limit. In these twenty-seven samples, U enrichment ranges up to a 10X increase, with an average U/Th ratio indicating a 2X increase in U. Our samples bracket the typical enrichment value of 5X for altered basalt when compared to unaltered glass (Kelley et al., 2003). The samples with excess U increase H2 yields on average by 29 ± 6%, 60 ± 3%, and 23 ± 10% at Sites U1365, U1367, and U1368, respectively. Therefore, the variation in radiolytic H2 production rates that we calculate for bulk basalt likely accounts for much of the variation that would be seen in the oceanic crust.
As shown by the variety of visual alteration categories and the wide range of U/Th ratios exhibited by our SPG samples, the kind, and extent of alteration vary greatly on a variety of spatial scales. The most intensely altered samples are likely those with most U enrichment; however that does not mean that they have the highest radionuclide concentrations. This can make links between alteration extent and radiolytic yields difficult to clearly identify. Separation of samples by alteration category or by U/Th ratios does not result in obvious patterns, suggesting that there is no simple relationship between radiolytic rates and alteration. Visual inspection of the samples is not enough to fully assess the impact of alteration on radiolytic rates. For example, some of the SPG samples with brown and dark gray alteration halos have U/Th ratios that indicate no U enrichment. Other studies, such as the SPG study of Zhang and Smith-Duque (2014) and the western Pacific study of Kelley et al. (2003) have shown U enrichment especially within these zones. Much closer analysis of mineral and chemical composition would be needed to fully assess effects of alteration on radionuclide concentrations and radiolytic H2 production rates of SPG basalt.
The rates we calculate are based on bulk-rock analyses and may underestimate radiolytic H2 production associated with specific alteration phases. For localized mineral phases with high U, Th, or K concentrations rates will be much higher. Staudigel et al. (2008) suggest that weathering of basaltic glass is the dominant process of chemical exchange between basalt and seawater. Alteration of basaltic glass produces palagonite (Staudigel et al., 2008; Türke et al., 2015), a very porous mineral (∼14–38 wt% H2O; Pauly et al., 2011) with high concentrations of radioactive elements (Türke et al., 2015).
Palagonite is of special interest for studies of radiolytic H2 production in subseafloor basalt because its porous nature and its enriched radionuclide concentrations must result in high radiolytic rates (Türke et al., 2015). Türke et al. (2015) calculated H2 production rates within palagonite rims and show that H2 can accumulate to minimum concentrations needed for hydrogentrophy at North Pond (on the western flank of the Mid-Atlantic Ridge). For their calculations, Türke et al. (2015) used a porosity model based on protocols from Lin et al. (2005a) and Blair et al. (2007). Because palagonite is so porous, these models’ assumptions of fine grain size and homogeneously distributed porosity are broadly relevant for palagonite calculations. However, these models’ assumption that all of the radiation interacts with water is not appropriate for relatively thin palagonite occurrences. Türke et al. (2015) note this issue and calculate yields assuming 1% efficiency. However, there is no explicit way to calculate the actual efficiency using their model. The average H2 production rates calculated by Türke et al. (2015) are about four times higher than average SPG rates in microfractures. Consequently, palagonite and other similarly porous alteration minerals may significantly increase radiolytic H2 production relative to neighboring basalt.
Calculating the global rate of radiolytic H2 production in oceanic basement requires firm constraints on the distribution of fracture widths in the basement, mean concentrations of U, Th and K in the basalt, and porosity and abundance of alteration minerals in the basement. Mapping the horizontal and vertical distribution of radiolytic H2 production rates in subseafloor basalt requires information beyond the mean, including knowledge of horizontal and vertical variation in fracture width and number, radionuclide (U, Th, and K) concentrations, and porosity of alteration phases. Spatial variation in radionuclide concentrations is especially difficult to constrain, because it in turn depends on both (i) initial composition (which can vary from region to region, and even between successive eruptive events at the same location; Bergmanis et al., 2007) and (ii) alteration history (including the extent and kinds of secondary mineralization; Alt, 2004).
Comparison of Electron Donor Fluxes from Radiolytic H2 Production and Iron Oxidation
To assess the potential importance of radiolytic H2 production for microbial communities in subseafloor basalt, we compare our calculated radiolytic H2 production rates to iron oxidation rates calculated from Fe(III)/ΣFe ratios compiled by Bach and Edwards (2003; Figure 5).
FIGURE 5. Comparison of electron transfer rates supportable by Fe oxidation and radiolytic H2 production. Ranges of electron (e-) transfer rates due to Fe oxidation are shown for basalt of two different ages, 1 and 10 Ma shown by the gray and orange shading, respectively. The blue-hatched region indicates the range of electron transfer supportable by radiolytic H2 production.
To compare these rates, we convert them to electron transfer rates per unit volume. For Fe oxidation, one electron is transferred during oxidation of Fe(II) to Fe(III); for H2, two electrons are transferred for each H2 molecule produced. For direct comparison, we assume a cubic meter of basalt with 10% porosity, density of 2950 kg/m3 and 8 wt.% Fe (Bach and Edwards, 2003). We use an initial Fe(III)/ΣFe ratio of 0.15 ± 0.05 for fresh basalt and a final ratio of 0.45 ± 0.15 for basalt greater than 10 Ma (Bach and Edwards, 2003). These values include the variation associated with the average Fe(III)/ΣFe ratios at their respective basalt ages. Given these ratios, we calculate a mean oxidation rate of 0.19 ± 0.1 mol Fe oxidized kg basalt-1 Myr-1 for the first 10 Myr of the basalt’s existence. Considerable variability in the Fe(III)/ΣFe ratios given by Bach and Edwards (2003) indicates that different regions within the basement have higher, or lower, Fe oxidation rates. For H2 calculations, we assume that each unique fracture width makes up the 10% porosity in the basalt. For example, there are one thousand 100 μm fractures or one 10 cm fracture in our calculations to account for the correct water-rock ratio.
Given these assumptions, for oceanic crust younger than 10 Ma, calculated electron transfer rates (mol e- per m3 rock per year) due to Fe oxidation are about three orders of magnitude higher than electron transfer rates associated with radiolytic H2 production (Figure 5). However, after circa 10 Ma, mean Fe(III)/ΣFe ratios appear to plateau, suggesting that all available Fe has been oxidized (Bach and Edwards, 2003). This means that Fe oxidation rates in basalt older than circa 10 Ma are near zero. In this older basalt, radiolytic H2 may be more important than Fe(II) as an electron donor for microbial communities.
This change from Fe(II) to radiolytic H2 as the predominant electron donor is consistent with the palagonite-based conclusion of Türke et al. (2015). In young ridge flanks, Fe oxidation is a main source of energy for microorganisms due to interaction of seawater with fresh basalt in fractures. However, in older crust, after permeability decreases and access to fresh basalt is prevented by alteration phases, the availability of Fe decreases and the dominant electron donor may be radiolytic H2 (Türke et al., 2015).
Numbers of Cells that Might be Supported by Radiolytic H2 Production in Fractures
To estimate how many cells might be supported by oxidation of radiolytic H2, we divide the radiolytic H2 flux for different fracture widths by the net per-cell oxygen reduction rates in SPG sediment (Table 1). As with our comparison to Fe(II) oxidation, we convert both radiolytic H2 production rates and oxygen reduction rates to electron transfer rates. Net per-cell oxygen reduction rates in subseafloor SPG sediment range between 4.2 × 10-17 ± 2.6 × 10-17 mol e- cell-1 yr-1 and 6.8 × 10-15 ± 1.3 × 10-15 mol e- cell-1 yr-1 (D’Hondt et al., 2015). Using these maximum and minimum oxygen reduction rates and the radiolytic electron fluxes (mol e- cm-2 yr-1) for discrete fracture widths in our least and most radioactive basalt samples, we can estimate the number of cells that might be supported by H2 production. If the cells respire at the same rate as the aerobic communities in subseafloor SPG sediment, radiolytic H2 production in the volume of water adjacent to a square cm of the least radioactive SPG basalt might support up to 30 cells. H2 production in the same volume of water adjacent to a square cm of the most radioactive SPG basalt might support up to 1500 cells.
The extent of life in oceanic crust must depend in large part on the availability of electron donors and acceptors. Water radiolysis produces H2, which can be metabolized by microorganisms. Microfractures and local maxima in radionuclide concentrations serve as hotspots for radiolytic H2 production and may also act as hotspots for microbial life. Differences in initial melt composition and low-temperature alteration by seawater affect concentrations and distributions of U, Th, and K within the basement, and consequently also change radiolytic H2 production rates. Our calculations suggest that in young (less than 10 Ma) basalt, oxidizable Fe(II) is a far more abundant electron donor than radiolytic H2. However, in older seafloor, where little Fe is accessible to oxidants in the formation water, radiolytic H2 may be the dominant electron donor. Radiolytic H2 in the water adjacent to a square cm of SPG basaltic fracture may support up to 103 cells if the cells respire at the same rate as net per-cell oxygen reduction in subseafloor SPG sediment.
All authors listed, have made substantial, direct, and intellectual contribution to the work, and approved it for publication.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the reviewers for helpful comments. We thank the National Aeronautics and Space Administration (grant NNX12AD65G) and the U.S. National Science Foundation (grant NSF-OCE-1130735 to SD and NSF OCE1130531 to RM, and through the Center for Dark Energy Biosphere Investigations; grant NSF-OCE-0939564 to SD) for funding this study. Portions of this material are based upon work supported while RM was serving at the National Science Foundation. This is C-DEBI publication number 287.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.00076
Alt, J. C. (2004). “Alteration of the upper oceanic crust: mineralogy, chemistry, and processes,” in Hydrogeology of the Oceanic Lithosphere, eds E. E. Davis and H. Elderfield (New York, NY: Cambridge University Press), 497–595.
Alt, J. C., and Teagle, D. A. (2003). Hydrothermal alteration of upper oceanic crust formed at a fast-spreading ridge: mineral, chemical, and isotopic evidence from ODP Site 801. Chem. Geol. 201, 191–211. doi: 10.1016/S0009-2541(03)00201-8
Alt, J. C., France-Lanord, C., Floyd, P. A., Castillo, P., and Galy, A. (1992). “Low-temperature hydrothermal alteration of Jurassic ocean crust, Site 801,” in Proceeding of the Ocean Drilling Program Scientific Results, Vol. 129, eds R. L. Larson and Y. Lancelot (College Station, TX: Ocean Drilling Program), 415–427.
Bach, W., and Edwards, K. (2003). Iron and sulfide oxidation within the basaltic ocean crust: implications for chemolithoautotrophic microbial biomass production. Geochim. Cosmochim. Acta 67, 3871–3887. doi: 10.1016/S0016-7037(03)00304-1
Berger, M. J., Coursey, J. S., Zucker, M. A., and Chang, J. (2005). ESTAR, PSTAR, and ASTAR: Computer Programs for Calculating Stopping-Power and Range Tables for Electrons, Protons, and Helium Ions (version 1.2.3). Gaithersburg, MD: National Institute of Standards and Technology. Available at: http://physics.nist.gov/Star [accessed on December 30, 2015].
Bergmanis, E. C., Sinton, J., and Rubin, K. H. (2007). Recent eruptive history and magma reservoir dynamics on the southern east pacific rise at 17°30′S. Geochem. Geophys. Geosyst. 8:Q12O06. doi: 10.1029/2007GC001742
Chapelle, F. H., O’Neill, K., Bradley, P. M., Methe, B. A., Ciufo, S. A., Knobel, L. L., et al. (2002). A hydrogen-based subsurface microbial community dominated by methanogens. Nature 415, 312–314. doi: 10.1038/415312a
Cowen, J. P., Giovannoni, S. J., Kenig, F., Johnson, H. P., Butterfield, D., Rappé, M. S., et al. (2003). Fluids from aging ocean crust that support microbial life. Science 299, 120–123. doi: 10.1126/science.1075653
D’Hondt, S., F. Inagaki, C., Zarikian, A., and the IODP Expedition. 329 Scientists (2013). IODP Expedition 329: life and habitability beneath the seafloor of the South Pacific Gyre. Sci. Drill. 15, 4–10.
D’Hondt, S., Inagaki, F., Alvarez Zarikian, C., Abrams, L. J., Dubois, N., Engelhardt, T., et al. (2015). Presence of oxygen and aerobic communities from seafloor to basement in deep-sea sediment. Nat. Geosci. 8, 299–304. doi: 10.1038/NGEO2387
D’Hondt, S. F., Inagaki, C. A. Zarikian, A., and the Expedition Scientists (2011). “South Pacific Gyre Subseafloor Life,” in Proceeding of the Ocean Drilling Program, 329, Tokyo: Integrated Ocean Drilling Program Management International, Inc.
D’Hondt, S., Jørgensen, B. B., Miller, D. J., Batzke, A., Blake, R., Cragg, B. A., et al. (2004). Distributions of microbial activities in deep subseafloor sediments. Science 306, 2216–2221. doi: 10.1126/science.1101155
D’Hondt, S., Spivack, A. J., Pockalny, R., Ferdelman, T. G., Fischer, J. P., Kallmeyer, J., et al. (2009). Subseafloor sedimentary life in the South Pacific Gyre. Proc. Natl. Acad. Sci. U.S.A. 106, 11651–11656. doi: 10.1073/pnas.0811793106
Dunlea, A. G., Murray, R. W., Sauvage, J., Pockalny, R. A., Spivack, A. J., Harris, R. N., et al. (2015a). Cobalt-based age models of pelagic clay in the South Pacific Gyre. Geochem. Geophys. Geosyst. 16, 2694–2710. doi: 10.1002/2015GC005892
Dunlea, A. G., Murray, R. W., Sauvage, J., Spivack, A. J., Harris, R. N., and D’Hondt, S. (2015b). Dust, volcanic ash, and the evolution of the South Pacific Gyre through the Cenozoic. Paleoceanography 30, 1078–1099. doi: 10.1002/2015PA002829
Dzaugis, M. E., Spivack, A. J., and D’Hondt, S. (2015). A quantitative model of water radiolysis and chemical production rates near radionuclide-containing solids. Radiat. Phys. Chem. 115, 127–134. doi: 10.1016/j.radphyschem.2015.06.011
Edwards, K. J., Bach, W., and McCollom, T. M. (2005). Geomicrobiology in oceanography: microbe–mineral interactions at and below the seafloor. Trends Microbiol. 13, 449–456. doi: 10.1016/j.tim.2005.07.005
Edwards, K., Fisher, A., and Wheat, C. (2012). The deep subsurface biosphere in igneous ocean crust: frontier habitats for microbiological exploration. Front. Microbiol. 3:8. doi: 10.3389/fmicb.2012.00008
Expedition 329 Scientists (2011a). “Site U1365,” in Proceeding of the Integrated Ocean Drilling Program Management International, 329, eds S. D’Hondt, F. Inagaki, C. A. Alvarez Zarikian, and the Expedition 329 Scientists (Tokyo: Integrated Ocean Drilling Program Management International, Inc.).
Expedition 329 Scientists (2011b). “Site U1367,” in Proceeding of the Integrated Ocean Drilling Program Management International, 329, eds S. D’Hondt, F. Inagaki, C. A. Alvarez Zarikian, and the Expedition 329 Scientists (Tokyo: Integrated Ocean Drilling Program Management International, Inc.).
Expedition 329 Scientists (2011c). “Site U1368,” in Proceeding of the Integrated Ocean Drilling Program Management International, 329, eds S. D’Hondt, F. Inagaki, C. A. Alvarez Zarikian, and the Expedition 329 Scientists (Tokyo: Integrated Ocean Drilling Program Management International, Inc.).
Hubbell, J. H., and Seltzer, S. M. (2004). Tables of X-Ray Mass Attenuation Coefficients and Mass Energy-Absorption Coefficients (version 1.4). Gaithersburg, MD: National Institute of Standards and Technology. Available at: http://physics.nist.gov/xaamdi [accessed on January 13, 2016].
Jochum, K. P., Nohl, U., Herwig, K., Lammel, E., Stoll, B., and Hofmann, A. W. (2005). GeoReM: a new geochemical database for reference materials and isotopic standards. Geostand. Geoanalyt. Res. 29, 333–338. doi: 10.1111/j.1751-908X.2005.tb00904.x
Jungbluth, S. P., Grote, J., Lin, H.-T., Cowen, J. P., and Rappé, M. S. (2013). Microbial diversity within basement fluids of the sediment-buried juan de fuca ridge flank. ISME 7, 161–172. doi: 10.1038/ismej.2012.73
Kelley, D. S., Karson, J. A., Blackman, D. K., Früh-Green, G. L., Butterfield, D. A., Lilley, M. D., et al. (2001). An off-axis hydrothermal vent field near the mid-atlantic ridge at 30° N. Nature 412, 145–149. doi: 10.1038/35084000
Kohan, L. M., Sanguanmith, S., Meesungnoen, J., Causey, P., Stuart, C. R., and Jay-Gerin, J. (2013). Self-Radiolysis of Tritiated Water. 1. A comparison of the effects of 60Co Γ-Rays and Tritium β-particles on water and aqueous solutions at room temperature. RSC Adv. 3:19282. doi: 10.1039/c3ra42984a
LaTourrette, T. Z., Kennedy, A. K., and Wasserburg, G. J. (1993). Thorium-uranium fractionation by garnet: evidence for a deep source and rapid rise of oceanic basalts. Science 261, 739–742. doi: 10.1126/science.261.5122.739
Lever, M. A., Rouxel, O., Alt, J. C., Shimizu, N., Ono, S., Coggon, R. M., et al. (2013). Evidence for microbial carbon and sulfur cycling in deeply buried ridge flank basalt. Science 339, 1305–1308. doi: 10.1126/science.1229240
Lin, L.-H., Slater, G. F., Lollar, B. S., Lacrampe-Couloume, G., and Onstott, T. C. (2005a). The yield and isotopic composition of radiolytic H2, a potential energy source for the deep subsurface biosphere. Geochim. Cosmochim. Acta 69, 893–903. doi: 10.1016/j.gca.2004.07.032
Lin, L.-H., Hall, J., Lippmann-Pipke, J., Ward, J. A., Sherwood Lollar, B., DeFlaun, M., et al. (2005b). Radiolytic H2 in continental crust: nuclear power for deep subsurface microbial communities. Geochem. Geophys. Geosyst. 6, 7. doi: 10.1029/2004GC000907
Lollar, B. S., Onstott, T. C., Lacrampe-Couloume, G., and Ballentine, C. J. (2014). The contribution of the Precambrian continental lithosphere to global H2 production. Nature 516, 379–382. doi: 10.1038/nature14017
Mustaree, S., Meesungnoen, J., Butarbutar, S. L., Causey, P., Stuart, C. R., and Jay-Gerin, J. P. (2014). Self-radiolysis of tritiated water. 3. The OH scavenging effect of bromide ions on the yield of H2O2 in the radiolysis of water by 60Co γ-rays and tritium β-particles at room temperature. RSC Adv. 4, 43572–43581. doi: 10.1039/C4RA06707J
Orcutt, B. N., Wheat, C. G., Rouxel, O., Hulme, S., Edwards, K. J., and Bach, W. (2013). Oxygen consumption rates in subseafloor basaltic crust derived from a reaction transport model. Nat. Commun. 4:2539. doi: 10.1038/ncomms3539
Staudigel, H., Furnes, H., McLoughlin, N., Banerjee, N. R., Connell, L. B., and Templeton, A. (2008). 3.5 billion years of glass bioalteration: volcanic rocks as a basis for microbial life? Earth Sci. Rev. 89, 156–176. doi: 10.1016/j.earscirev.2008.04.005
Staudigel, H., Plank, T., White, B., and Schmincke, H. U. (1996). “Geochemical fluxes during seafloor alteration of the basaltic upper oceanic crust: DSDP Sites 417 and 418,” in Subduction Top to Bottom, eds G. E. Bebout, D. W. Scholl, S. H. Kirby, and J. P. Platt (Washington, DC: American Geophysical Union), 19–38.
Teagle, D. A., Alt, J. C., Bach, W., and Erzinger, J. (1996). “Alteration of upper ocean crust in a ridge-flank hydrothermal upflow zone: mineral, chemical, and isotopic constraints from Hole 896-A,” in Proceedings of the Ocean Drilling Program, Scientific Results, Vol. 148 (Texas, TX: Ocean Drilling Program, College Station), 119–150.
Türke, A., Nakamura, K., and Bach, W. (2015). Palagonitization of basalt glass in the flanks of mid-ocean ridges: implications for the bioenergetics of oceanic intracrustal ecosystems. Astrobiology 15, 793–803. doi: 10.1089/ast.2014.1255
Wang, K. (2004). “Applying fundamental principles and mathematical models to understand processes and estimate parameters,” in Hydrogeology of the Oceanic Lithosphere, eds E. E. Davis and H. Elderfield (New York, NY: Cambridge University Press), 497–595.
Zhang, G.-L., Chen, L.-H., and Li, S.-Z. (2013). Mantle dynamics and generation of a geochemical mantle boundary along the east pacific rise – pacific/antarctic ridge. Earth Planet. Sci. Lett. 383, 153–163. doi: 10.1016/j.epsl.2013.09.045
Keywords: radiolysis, hydrogen, basalt, ocean crust, geochemistry, deep biosphere
Citation: Dzaugis ME, Spivack AJ, Dunlea AG, Murray RW and D’Hondt S (2016) Radiolytic Hydrogen Production in the Subseafloor Basaltic Aquifer. Front. Microbiol. 7:76. doi: 10.3389/fmicb.2016.00076
Received: 01 September 2015; Accepted: 15 January 2016;
Published: 04 February 2016.
Edited by:Jason B. Sylvan, Texas A&M University, USA
Reviewed by:Craig Lee Moyer, Western Washington University, USA
Wolfgang Bach, University of Bremen, Germany
Lisa Mayhew, University of Colorado Boulder, USA
Copyright © 2016 Dzaugis, Spivack, Dunlea, Murray and D’Hondt. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mary E. Dzaugis, firstname.lastname@example.org