Modelling biofilm growth in the presence of carbon dioxide and water flow in the subsurface

The concentration of greenhouse gases – particularly carbon dioxide (CO₂) – in the atmosphere has been on the rise in the past decades. One of the methods which have been proposed to help reduce anthropogenic CO₂ emissions is the capture of CO₂from large, stationary point sources and storage in deep geological formations. The caprock is an impermeable geological layer which prevents the leakage of stored CO₂, and its integrity is of utmost importance for storage security. Due to the high pressure build-up during injection, the caprock in the vicinity of the well is particularly at risk of fracturing. Biofilms could be used as biobarriers which help prevent the leakage of CO₂ through the caprock in injection well vicinity by blocking leakage pathways. The biofilm could also protect well cement from corrosion by CO₂-rich brine.     

Petrology of Medicine Lake Highland volcanics: Characterization of endmembers of magma mixing

Lavas from Medicine Lake volcano, Northern California have been examined for evidence of magma mixing. Mixing of magmas has produced basaltic andesite, andesite, dacite and rhyolite lavas at the volcano. We are able to identify the compositional characteristics of the components that were mixed and to estimate the time lag between the mixing event and eruption of the mixed magma. Compositional data from pairs of phenocrysts identify a high alumina basalt (HAB) and a silicic rhyolite as endmembers of mixing. Mg-rich olivine or augite and Ca-rich plagioclase are associated with the HAB component, and Fe-rich orthopyroxene and Na-rich plagioclase are associated with the rhyolitic component. Some lavas contain multiple phenocryst assemblages suggesting the incorporation of several magmas intermediate between the HAB and silicic components. Glass inclusions trapped in Mg-rich olivine and Na-rich plagioclase are similar in composition to the proposed HAB and rhyolite end members and provide supportive evidence for mixing. Textural criteria are also consistent with magma mixing. Thermal curvature of the liquidus surfaces in the basalt-andesite-rhyolite system allows magmas produced by mixing to be either supercooled or superheated. Intergranular textures of basaltic andesites and andesites result from cooling initiated below the liquidus. The trachytic textures of silicic andesites form from cooling initiated above the liquidus. Reversed compositional zoning profiles in olivine crystals were produced by the mixing event, and the homogenization of the compositional zoning has been used to estimate the time interval between magma mixing and eruption. Time estimates are on the order of 80 to 90 h, suggesting that the mixing event triggered eruption.     

Evaluation of gases,condensates, and SO2 emissions from Augustine volcano,Alaska: the degassing of a Cl-rich volcanic system

After the March–April 1986 explosive eruption a comprehensive gas study at Augustine was undertaken in the summers of 1986 and 1987. Airborne COSPEC measurements indicate that passive SO₂ emission rates declined exponentially during this period from 380±45 metric tons/day (T/D) on 7/24/86 to 27±6 T/D on 8/24/87. These data are consistent with the hypothesis that the Augustine magma reservoir has become more degassed as volcanic activity decreased after the spring 1986 eruption. Gas samples collected in 1987 from an 870°C fumarole on the andesitic lava dome show various degrees of disequilibrium due to oxidation of reduced gas species and condensation (and loss) of H₂O in the intake tube of the sampling apparatus. Thermochemical restoration of the data permits removal of these effects to infer an equilibrium composition of the gases. Although not conclusive, this restoration is consistent with the idea that the gases were in equilibrium at 870°C with an oxygen fugacity near the Ni–NiO buffer. These restored gas compositions show that, relative to other convergent plate volcanoes, the Augustine gases are very HCl rich (5.3–6.0 mol% HCl), S rich (7.1 mol% total S), and H₂O poor (83.9–84.8 mol% H₂O). Values of D and ¹⁸O suggest that the H₂O in the dome gases is a mixture of primary magmatic water (PMW) and local seawater. Part of the Cl in the Augustine volcanic gases probably comes from this shallow seawater source. Additional Cl may come from subducted oceanic crust because data by Johnston (1978) show that Cl-rich glass inclusions in olivine crystals contain hornblende, which is evidence for a deep source (>25km) for part of the Cl. Gas samples collected in 1986 from 390°–642°C fumaroles on a ramp surrounding the inner summit crater have been oxidized so severely that restoration to an equilibrium composition is not possible. H and O isotope data suggest that these gases are variable mixtures of seawater, FMW, and meteoric steam. These samples are much more H₂O-rich (92%–97% H₂O) than the dome gases, possibly due to a larger meteoric steam component. The 1986 samples also have higher Cl/S, S/C, and F/Cl ratios, which imply that the magmatic component in these gases is from the more degassed 1976 magma. Thus, the 1987 samples from the lava dome are better indicators than the 1986 samples of degassing within the Augustine magma reservoir, even though they were collected a year later and contain a significant seawater component. Future gas studies at Augustine should emphasize fumaroles on active lava domes. Condensates collected from the same lava-dome fumarole have enrichments ot 10⁷–10² in Cl, Br, F, B, Cd, As, S, Bi, Pb, Sb, Mo, Zn, Cu, K, Li, Na, Si, and Ni. Lower-temperature (200°–650°C) fumaroles around the volcano are generally less enriched in highly volatile elements. However, these lower-termperature fumaroles have higher concentration of rock-forming elements, probably derived from the wall rock.     

Geochemical variations in lavas from Kahoolawe volcano,Hawaii: evidence for open system evolution of plume-derived magmas

The Kahoolawe shield volcano produced precaldera and caldera-filling tholeiites and mildly alkalic post-caldera lavas that petrographically and compositionally resemble such lavas from other Hawaiian shield volcanoes. However, Kahoolawe tholeiites display wide ranges in incompatible trace element ratios (e.g., Nb/Th=9–24, Th/Ta=0.6–1.3), ⁸⁷Sr/⁸⁶Sr (0.70379–0.70440), ¹⁴³Nd/¹⁴⁴Nd (0.51273–0.51298), and ²⁰⁶Pb/²⁰⁴Pb (17.92–18.37). The isotopic variation exceeds that at any other Hawaiian shield volcano, and spans about half the range for all Hawaiian tholeiites. Quasi-cyclic temporal evolution of Kahoolawe tholeiites is consistent with combined fractional crystallization and periodic recharge by primitive magmas. Ratios of highly incompatible trace elements and Sr, Nd, and Pb isotopic ratios from coherent sub-trends that reflect recurrent interactions between variably evolved magmas and two other mantle components whose compositions are constrained by intersections between these trends. The most MgO-rich Kahoolawe tholeiites are partial melts of a high Nb/Th (23.5) ascending plume, possibly comprising ancient subducted oceanic lithosphere. Slightly evolved tholeiites experienced combined crystal fractionation and assimilation (AFC) of material derived from a distinct reservoir (Nb/Th 9) of asthenospheric derivation. The most evolved tholeiites display compositional shifts toward a third component, having mid ocean ridge basalt-like isotopic ratios but enriched OIB-like trace element ratios, representing part of the lithospheric mantle (or melts thereof). Periodic recurrence of all three magma variants suggests that eruptions may have tapped coeval reservoirs distributed over a large depth range. Kahoolawe provides new evidence concerning the nature of the Hawaiian plume, the distribution of compositional heterogeneities in the suboeanic mantle, and the processes by which Hawaiian tholeiites form and evolve.     

Reconnaissance dating: A new radiocarbon method applied to assessing the temporal distribution of Southern Ocean deep-sea corals

We have developed a rapid ‘reconnaissance’ method of preparing graphite for ¹⁴C/¹²C analysis. Carbonate (～15 mg) is combusted using an elemental analyzer and the resulting CO₂ is converted to graphite using a sealed tube zinc reduction method. Over 85% (n=45 replicates on twenty-one individual corals) of reconnaissance ages measured on corals ranging in age from 500 to 33,000 radiocarbon years (Ryr) are within two standard deviations of ages generated using standard hydrolysis methods on the same corals, and all reconnaissance ages are within 300 Ryr of the standard hydrolysis ages. Replicate measurements on three individual aragonitic corals yielded ages of 1076±35 Ryr (standard deviation; n=5), 10,739±47 Ryr (n=8), and 40,146±3500 Ryr (n=9). No systematic biases were found using different cleaning methods or variable sample sizes. Analysis of ¹³C/¹²C was made concurrently with the ¹⁴C/¹²C measurement to correct for natural fractionation and for fractionation during sample processing and analysis. This technique provides a new, rapid method for making accurate, percent-level ¹⁴C/¹²C analyses that may be used to establish the rates and chronology of earth system processes where survey-type modes of age estimation are desirable. For example, applications may include creation of sediment core-top maps, preliminary age models for sediment cores, and growth rate studies of marine organisms such as corals or mollusks. We applied the reconnaissance method to more than 100 solitary deep-sea corals collected in the Drake Passage in the Southern Ocean to investigate their temporal and spatial distribution. The corals used in this study are part of a larger sample set, and the subset that was dated was chosen based on species as opposed to preservation state, so as to exclude obvious temporal biases. Similar to studies in other regions, the distribution of deep-sea corals is not constant through time across the Drake Passage. Most of the corals from the Burdwood Bank (continental shelf of Argentina) have ages ranging between 0 and 2500 calendar years, whereas most of the corals from the Sars Seamount in the Drake Passage have ages between 10,000 and 12,500 calendar years. Such differences may be caused in part by sampling biases, but may also be caused by changes in larval transport, nutrient supply, or other environmental pressures.     

Puhimau Thermal Area: A Window into the Upper East Rift Zone of Volcano, Hawaii?

We report the results of two soil CO₂ efflux surveys by the closed chamber circulation method at the Puhimau thermal area in the upper East Rift Zone (ERZ) of volcano, Hawaii. The surveys were undertaken in 1996 and 1998 to constrain how much CO₂ might be reaching the ERZ after degassing beneath the summit caldera and whether the Puhimau thermal area might be a significant contributor to the overall CO₂ budget of . The area was revisited in 2001 to determine the effects of surface disturbance on efflux values by the collar emplacement technique utilized in the earlier surveys. Utilizing a cutoff value of 50 g m⁻² d⁻¹ for the surrounding forest background efflux, the CO₂ emission rates for the anomaly at Puhimau thermal area were 27 t d⁻¹ in 1996 and 17 t d⁻¹ in 1998. Water vapor was removed before analysis in all cases in order to obtain CO₂ values on a dry air basis and mitigate the effect of water vapor dilution on the measurements. It is clear that Puhimau thermal area is not a significant contributor to CO₂ output and that most of CO₂ (8500 t d⁻¹) is degassed at the summit, leaving only magma with its remaining stored volatiles, such as SO₂, for injection down the ERZ. Because of the low CO₂ emission rate and the presence of a shallow water table in the upper ERZ that effectively scrubs SO₂ and other acid gases, Puhimau thermal area currently does not appear to be generally well suited for observing temporal changes in degassing at .     

Global height system unification with GOCE: a simulation study on the indirect bias term in the GBVP approach

One of the main objectives of ESA’s Gravity Field and Steady-State Ocean Circulation mission GOCE (Gravity field and steady-state ocean circulation mission, 1999) is to allow global unification of height systems by directly providing potential differences between benchmarks in different height datum zones. In other words, GOCE provides a globally consistent and unbiased geoid. If this information is combined with ellipsoidal (derived from geodetic space techniques) and physical heights (derived from leveling/gravimetry) at the same benchmarks, datum offsets between the datum zones can be determined and all zones unified. The expected accuracy of GOCE is around 2–3 cm up to spherical harmonic degree n _max ≈ 200. The omission error above this degree amounts to about 30 cm which cannot be neglected. Therefore, terrestrial residual gravity anomalies are necessary to evaluate the medium and short wavelengths of the geoid, i.e. one has to solve the Geodetic Boundary Value Problem (GBVP). The theory of height unification by the GBVP approach is well developed, see e.g. Colombo (A World Vertical Network. Report 296, Department of Geodetic Science and Surveying, 1980) or Rummel and Teunissen (Bull Geod 62:477–498, 1988). Thereby, it must be considered that terrestrial gravity anomalies referring to different datum zones are biased due to the respective datum offsets. Consequently, the height reference surface of a specific datum zone deviates from the unbiased geoid not only due to its own datum offset (direct bias term) but is also indirectly affected by the integration of biased gravity anomalies. The latter effect is called the indirect bias term and it considerably complicates the adjustment model for global height unification. If no satellite based gravity model is employed, this error amounts to about the same size as the datum offsets, i.e. 1–2 m globally. We show that this value decreases if a satellite-only gravity model is used. Specifically for GOCE with n _max ≈ 200, the error can be expected not to exceed the level of 1 cm, allowing the effect to be neglected in practical height unification. The results are supported by recent findings by Gatti et al. (J Geod, 2012).     

Implications for eruptive processes as indicated by sulfur dioxide emissions from K lauea Volcano, Hawai‘i, 1979–1997

K lauea Volcano, Hawai‘i, currently hosts the longest running SO₂ emission-rate data set on the planet, starting with initial surveys done in 1975 by Stoiber and his colleagues. The 17.5-year record of summit emissions, starting in 1979, shows the effects of summit and east rift eruptive processes, which define seven distinctly different periods of SO₂ release. Summit emissions jumped nearly 40% with the onset (3 January 1983) of the Pu‘u ‘ ‘ -K paianaha eruption on the east rift zone (ERZ). Summit SO₂ emissions from K lauea showed a strong positive correlation with short-period, shallow, caldera events, rather than with long-period seismicity as in more silicious systems. This correlation suggests a maturation process in the summit magma-transport system from 1986 through 1993. During a steady-state throughput-equilibrium interval of the summit magma reservoir, integration of summit-caldera and ERZ SO₂ emissions reveals an undegassed volume rate of effusion of 2.1×10⁵ m³/d. This value corroborates the volume-rate determined by geophysical methods, demonstrating that, for K lauea, SO₂ emission rates can be used to monitor effusion rate, supporting and supplementing other, more established geophysical methods. For the 17.5 years of continuous emission rate records at K lauea, the volcano has released 9.7×10⁶ t (metric tonnes) of SO₂, 1.7×10⁶ t from the summit and 8.0×10⁶ t from the east rift zone. On an annual basis, the average SO₂ release from K lauea is 4.6×10⁵ t/y, compared to the global annual volcanic emission rate of 1.2×10⁷ t/y.