The bursting behavior of gas slugs: laboratory and analytical insights into Strombolian volcanic eruptions.

The aim of this thesis is to study how explosive behavior and geophysical signals in a volcanic conduit are related to the development of overpressure in slug-driven eruptions. A first suite of laboratory experiments of gas slugs ascending in analogue conduits was performed. Slugs ascended into a range of analogue liquids and conduit diameters to allow proper scaling to the natural volcanoes. The geometrical variation of the slug in response to the explored variables was parameterised. Volume of gas slug and rheology of the liquid phase revealed the key parameters in controlling slug overpressure at bursting. Founded on these results, a theoretical model to calculate burst overpressure for slug-driven eruptions was developed. The dimensionless approach adopted allowed to apply the model to predict bursting pressure of slugs at Stromboli. Comparison of predicted values with measured data from Stromboli volcano showed that the model can explain the entire spectrum of observed eruptive styles at Stromboli – from low-energy puffing, through normal Strombolian eruptions, up to paroxysmal explosions – as manifestations of a single underlying physical process. Finally, another suite of laboratory experiments was performed to observe oscillatory pressure and forces variations generated during the expansion and bursting of gas slugs ascending in a conduit. Two end-member boundary conditions were imposed at the base of the pipe, simulating slug ascent in closed base (zero magma flux) and open base (constant flux) conduit. At the top of the pipe, a range of boundary conditions that are relevant at a volcanic vent were imposed, going from open to plugged vent. The results obtained illustrate that a change in boundary conditions in the conduit concur to affect the dynamic of slug expansion and burst: an upward flux at the base of the conduit attenuates the magnitude of the pressure transients, while a rheological stiffening in the top-most region of conduit changes dramatically the magnitude of the observed pressure transients, favoring a sudden, and more energetic pressure release into the overlying atmosphere. Finally, a discussion on the implication of changing boundary on the oscillatory processes generated at the volcanic scale is also given.

Studio dei Vulcani di Fango per la definizione della Migrazione dei Fluidi Profondi

L’unione di approcci differenti durante lo studio dei sistemi petroliferi, come l’accoppiamento dello studio delle emissioni in superficie con le analisi geochimiche e strutturali, è un aspetto principale nelle strategie di sviluppo per la ricerca degli idrocarburi. La presenza di acqua connata nelle sequenze sedimentarie profonde e la sua sovrappressione che viene generata dalle spesse coperture sedimentarie, incrementata inoltre dalla generazione di idrocarburi in profondità, sono fattori di controllo primari per la migrazione e l’emissione di fluidi in superficie. I risultati ottenuti da questo studio forniscono nuovi elementi per la comprensione del ruolo dello studio dei vulcani di fango nell’esplorazione petrolifera, e nuove importanti prove per la caratterizzazione dei sistemi petroliferi nelle aree considerate.

Polygonal cracks in the seasonal semi-translucent CO2 ice layer in Martian polar areas

In this paper, we use morphological and numerical methods to test the hypothesis that seasonally formed fracture patterns in the Martian polar regions result from the brittle failure of seasonal CO2 slab ice. The observations by the High Resolution Imaging Science Experiment (HiRISE) of polar regions of Mars show very narrow dark elongated linear patterns that are observed during some periods of time in spring, disappear in summer and re-appear again in the following spring. They are repeatedly formed in the same areas but they do not repeat the exact pattern from year to year. This leads to the conclusion that they are cracks formed in the seasonal ice layer. Some of models of seasonal surface processes rely on the existence of a transparent form of CO2 ice, so-called slab ice. For the creation of the observed cracks the ice is required to be a continuous media, not an agglomeration of relatively separate particles like a firn. The best explanation for our observations is a slab ice with relatively high transparency in the visible wavelength range. This transparency allows a solid state green-house effect to act underneath the ice sheet raising the pressure by sublimation from below. The trapped gas creates overpressure and the ice sheet breaks at some point creating the observed cracks. We show that the times when the cracks appear are in agreement with the model calculation, providing one more piece of evidence that CO2 slab ice covers polar areas in spring.

Dynamics and kinematics of eruptive activity at Fuego volcano, Guatemala 2005-2009

Volcanoes are the surficial expressions of complex pathways that vent magma and gasses generated deep in the Earth. Geophysical data record at least the partial history of magma and gas movement in the conduit and venting to the atmosphere. This work focuses on developing a more comprehensive understanding of explosive degassing at Fuego volcano, Guatemala through observations and analysis of geophysical data collected in 2005 – 2009. A pattern of eruptive activity was observed during 2005 – 2007 and quantified with seismic and infrasound, satellite thermal and gas measurements, and lava flow lengths. Eruptive styles are related to variable magma flux and accumulation of gas. Explosive degassing was recorded on broadband seismic and infrasound sensors in 2008 and 2009. Explosion energy partitioning between the ground and the atmosphere shows an increase in acoustic energy from 2008 to 2009, indicating a shift toward increased gas pressure in the conduit. Very-long-period (VLP) seismic signals are associated with the strongest explosions recorded in 2009 and waveform modeling in the 10 – 30 s band produces a best-fit source location 300 m west and 300 m below the summit crater. The calculated moment tensor indicates a volumetric source, which is modeled as a dike feeding a SW-dipping (35°) sill. The sill is the dominant component and its projection to the surface nearly intersects the summit crater. The deformation history of the sill is interpreted as: 1) an initial inflation due to pressurization, followed by 2) a rapid deflation as overpressure is explosively release, and finally 3) a reinflation as fresh magma flows into the sill and degasses. Tilt signals are derived from the horizontal components of the seismometer and show repetitive inflation deflation cycles with a 20 minute period coincident with strong explosions. These cycles represent the pressurization of the shallow conduit and explosive venting of overpressure that develops beneath a partially crystallized plug of magma. The energy released during the strong explosions has allowed for imaging of Fuego’s shallow conduit, which appears to have migrated west of the summit crater. In summary, Fuego is becoming more gas charged and its summit centered vent is shifting to the west - serious hazard consequences are likely.

Geophysical and Hydrological Evaluation of Two Bog Complexes in a Northern Peatland: Implications for the Distribution of Biogenic Gases at the Basin Scale

Ground penetrating radar (GPR) was used to determine peat basin geometry and the spatial distribution of free-phase biogenic gasses in two separate units of a northern peatland (Central and Southern Unit of Caribou Bog, Maine). The Central Unit is characterized by a deep basin structure (15 m maximum depth) and a raised (eccentric) bog topographic profile (up to 2 m topographic variation). Here numerous regions of electromagnetic (EM) wave scattering are considered diagnostic of the presence of extensive free-phase biogenic gas. In contrast, the Southern Unit is shallower (8 m maximum depth) and has a slightly convex upwards bog profile (less than 1 m topographic variation), and areas of EM wave scattering are notably absent. The biogenic gas zones interpreted from GPR in the Central Unit are associated with: (1) wooded heath vegetation at the surface, (2) open pools at the surface, (3) high water table elevations near the center of the basin, and (4) a region of overpressure (at approximately 5 m depth) immediately below the zone of free-phase gas accumulation. The latter suggests (1) a transient pressure head associated with low hydraulic conductivity resulting from the biogenic gasses themselves or confining layers in the peat that restrict both gas release and groundwater flow and/or (2) overpressure in the peat column as a result of the gas buildup itself. In contrast, the Southern Unit, where zones of EM scattering are absent, is characterized by: (1) predominantly shrub vegetation, (2) a lack of open pools, (3) only minor variations (less than 1 m) in water table elevation throughout the entire unit; and (4) generally upward groundwater flow throughout the basin. The results illustrate the nonuniformity of free-phase biogenic gas distribution at the peat basin scale and provide insights into the processes and controls associated with CH4 and CO2 accumulation in peatlands.

Diagenesis of a light, tight-oil chert reservoir at the Ediacaran/Cambrian boundary, Sultanate of Oman

The Al Shomou Silicilyte Member (Athel Formation) in the South Oman Salt Basin shares many of the characteristics of a light, tight-oil (LTO) reservoir: it is a prolifi c source rock mature for light oil, it produces light oil from a very tight matrix and reservoir, and hydraulic fracking technology is required to produce the oil. What is intriguing about the Al Shomou Silicilyte, and different from other LTO reservoirs, is its position related to the Precambrian/Cambrian Boundary (PCB) and the fact that it is a ‘laminated chert‘ rather than a shale. In an integrated diagenetic study we applied microstructural analyses (SEM, BSE) combined with state-of-the-art stable isotope and trace element analysis of the silicilyte matrix and fractures. Fluid inclusion microthermometry was applied to record the salinity and minimum trapping temperatures. The microstructural investigations reveal a fi ne lamination of the silicilyte matrix with a mean lamina thickness of ca. 20 μm consisting of predominantly organic matter-rich and fi nely crystalline quartz-rich layers, respectively. Authigenic, micron-sized idiomorphic quartz crystals are the main matrix components of the silicilyte. Other diagenetic phases are pyrite, apatite, dolomite, magnesite and barite cements. Porosity values based on neutron density logs and core plug data indicate porosity in the silicilyte ranges from less than 2% to almost to 40%. The majority of the pore space in the silicilyte is related to (primary) inter-crystalline pores, with locally important oversized secondary pores. Pore casts of the silica matrix show that pores are extremely irregular in three dimensions, and are generally interconnected by a complex web or meshwork of fi ne elongate pore throats. Mercury injection capillary data are in line with the microstructural observations suggesting two populations of pore throats, with an effective average modal diameter of 0.4 μm. The acquired geochemical data support the interpretation that the primary source of the silica is the ambient seawater rather than hydrothermal or biogenic. A maximum temperature of ca. 45°C for the formation of microcrystalline quartz in the silicilyte is good evidence that the lithifi cation and crystallization of quartz occurred in the fi rst 5 Ma after deposition. Several phases of brittle fracturing and mineralization occurred in response to salt tectonics during burial. The sequences of fracture-fi lling mineral phases (dolomite - layered chalcedony – quartz – apatite - magnesite I+II - barite – halite) indicates a complex fl uid evolution after silicilyte lithifi cation. Primary, all-liquid fl uid inclusions in the fracturefi lling quartz are good evidence of growth beginning at low temperatures, i.e. ≤ 50ºC. Continuous precipitation during increasing temperature and burial is documented by primary two-phase fl uid inclusions in quartz cements that show brines at 50°C and fi rst hydrocarbons at ca. 70°C. The absolute timing of each mineral phase can be constrained based on U-Pb geochronometry, and basin modelling. Secondary fl uid inclusions in quartz, magnesite and barite indicate reactivation of the fracture system after peak burial temperature during the major cooling event, i.e. uplift, between 450 and 310 Ma. A number of fi rst-order trends in porosity and reservoir-quality distribution are observed which are strongly related to the diagenetic and fl uid history of the reservoir: the early in-situ generation of hydrocarbons and overpressure development arrests diagenesis and preserves matrix porosity. Chemical compaction by pressure dissolution in the fl ank areas could be a valid hypothesis to explain the porosity variations in the silicilitye slabs resulting in lower porosity and poorer connectivity on the fl anks of the reservoir. Most of the hydrocarbon storage and production comes from intervals characterized by Amthor et al. 114488 preserved micropores, not hydrocarbon storage in a fracture system. The absence of oil expulsion results in present-day high oil saturations. The main diagenetic modifi cations of the silicilyte occurred and were completed relatively early in its history, i.e. before 300 Ma. An instrumental factor for preserving matrix porosity is the diffi culty for a given slab to evacuate all the fl uids (water and hydrocarbons), or in other words, the very good sealing capacity of the salt embedding the slab.

Stable isotope and geochemical record of sediments from the Bering Sea

here is controversy over the role of marine methane hydrates in atmospheric methane concentrations and climate change during the last glacial period. In this study of two sediment cores from the southeast Bering Sea (700 m and 1467 m water depth), we identify multiple episodes during the last glacial period of intense methane flux reaching the seafloor. Within the uncertainty of the radiocarbon age model, the episodes are contemporaneous in the two cores and have similar timing and duration as Dansgaard-Oeschger events. The episodes are marked by horizons of sediment containing 13C-depleted authigenic carbonate minerals; 13C-depleted archaeal and bacterial lipids, which resemble those found in ANME-1 type anaerobic methane oxidizing microbial consortia; and changes in the abundance and species distribution of benthic foraminifera. The similar timing and isotopic composition of the authigenic carbonates in the two cores is consistent with a region-wide increase in the upward flux of methane bearing fluids. This study is the first observation outside Santa Barbara Basin of pervasive, repeated methane flux in glacial sediments. However, contrary to the "Clathrate Gun Hypothesis" (Kennett et al., 2003), these coring sites are too deep for methane hydrate destabilization to be the cause, implying that a much larger part of the ocean's sedimentary methane may participate in climate or carbon cycle feedback at millennial timescales. We speculate that pulses of methane in these opal-rich sediments could be caused by the sudden release of overpressure in pore fluids that builds up gradually with silica diagenesis. The release could be triggered by seismic shaking on the Aleutian subduction zone caused by hydrostatic pressure increase associated with sea level rise at the start of interstadials.

Consolidation properties and hydraulic conductivities of sediments from ODP Leg 207 sites

The stress history, permeability, and compressibility of sediments from Demerara Rise recovered during Ocean Drilling Program Leg 207 were determined using one-dimensional incremental load consolidation and low-gradient flow pump permeability tests. Relationships among void ratio, effective stress, and hydraulic conductivity are presented for sampled lithologic units and used to reconstruct effective stress, permeability, and in situ void ratio profiles for a transect of three sites across Demerara Rise. Results confirm that a significant erosional event occurred on the northeastern flank of the rise during the late Miocene, resulting in the removal of ~220 m of upper Oligocene-Miocene deposits. Although Neogene and Paleogene sediments tend to be overconsolidated, Cretaceous sediments are normally consolidated to underconsolidated, suggesting the presence of overpressure. A pronounced drop in permeability occurs at the transition from the Cretaceous black shales into the overlying Maastrichtian-upper Paleocene chalks and clays. The development of a hydraulic seal at this boundary may be responsible for overpressure in the Cretaceous deposits, leading to the lower overconsolidation ratios of these sediments. Coupled with large regional variations in sediment thickness (overburden stresses), the higher permeability overpressured Cretaceous sediments represent a regional lateral fluid conduit on Demerara Rise, possibly venting methane-rich fluids where it outcrops on the margin's northeastern flank.

(Table 1) Geochemical and geothermometry data from fibrous calcite vein fills in the Dabashan Fold Belt, China

Fibrous calcite veins with organic inclusions have been widely considered as indicators of oil and gas generation and migration under overpressure. Abundant fibrous calcite veins containing organic-bearing inclusions occur in faulted Lower Paleozoic through Triassic hydrocarbon source rocks in the Dabashan Foreland Belt (DBF). d13CPDB and d18OPDB values of the fibrous calcite range from - 4.8 to -1.9 to per mil and - 12.8 to - 8.4 per mil respectively, which is lighter than that of associated carbonate host rocks ranging from - 1.7 to + 3.1 per mil and - 8.7 to - 4.5 per mil. A linear relationship between d13CPDB and d18OPDB indicates that the calcite veins were precipitated from a mixture of basinal and surface fluids. The fibrous calcite contains a variety of inclusions, such as solid bitumen, methane bearing all-liquid inclusions, and vapor-liquid aqueous inclusions. Homogenization temperatures of aqueous inclusions range from 140 to 196° with an average of 179°. Salinities of aqueous inclusions average 9.7 wt% NaCl. Independent temperatures from bitumen reflectance and inclusion phase relationships of aqueous and methane inclusions were used to determine fluid pressures. Results indicate high pressures, elevated above typical lithostatic confining pressure, from 150 to 200 MPa. The elevated salinity and high temperature and pressure conditions of the fibrous calcite veins argue against an origin solely from burial overpressure resulting from clay transformation and dehydration reactions. Instead fluid inclusion P-T data and geochemistry results and regional geology indicate abnormally high pressures during fluid migration. These findings indicate that tectonic stress generated fracture and fault fluid pathways and caused migration of organic bearing fluids from the DBF during the Yanshan orogeny.

Physical properties and bulk minerals from IODP Holes 308-U1322B, U1322D and U1324B

How the micro-scale fabric of clay-rich mudstone evolves during consolidation in early burial is critical to how they are interpreted in the deeper portions of sedimentary basins. Core samples from the Integrated Ocean Drilling Program Expedition 308, Ursa Basin, Gulf of Mexico, covering seafloor to 600 meters below sea floor (mbsf) are ideal for studying the micro-scale fabric of mudstones. Mudstones of consistent composition and grain size decrease in porosity from 80% at the seafloor to 37% at 600 mbsf. Argon-ion milling produces flat surfaces to image this pore evolution over a vertical effective stress range of 0.25 (71 mbsf) to 4.05 MPa (597 mbsf). With increasing burial, pores become elongated, mean pore size decreases, and there is preferential loss of the largest pores. There is a small increase in clay mineral preferred orientation as recorded by high resolution X-ray goniometry with burial.

Shear experiments from IODP Expedition 322 and 333 samples

The mechanical behavior of the plate boundary fault zone is of paramount importance in subduction zones, because it controls megathrust earthquake nucleation and propagation as well as the structural style of the forearc. In the Nankai area along the NanTroSEIZE (Kumano) drilling transect offshore SW Japan, a heterogeneous sedimentary sequence overlying the oceanic crust enters the subduction zone. In order to predict how variations in lithology, and thus mechanical properties, affect the formation and evolution of the plate boundary fault, we conducted laboratory tests measuring the shear strengths of sediments approaching the trench covering each major lithological sedimentary unit. We observe that shear strength increases nonlinearly with depth, such that the (apparent) coefficient of friction decreases. In combination with a critical taper analysis, the results imply that the plate boundary position is located on the main frontal thrust. Further landward, the plate boundary is expected to step down into progressively lower stratigraphic units, assisted by moderately elevated pore pressures. As seismogenic depths are approached, the décollement may further step down to lower volcaniclastic or pelagic strata but this requires specific overpressure conditions. High-taper angle and elevated strengths in the toe region may be local features restricted to the Kumano transect.

Pore pressure measurements on IODP Exp308 sites

Overpressures measured with pore pressure penetrometers during Integrated Ocean Drilling Program (IODP) Expedition 308 reach 70% and 60% of the hydrostatic effective stress (View the MathML source) in the first 200 meters below sea floor (mbsf) at Sites U1322 and U1324, respectively, in the deepwater Gulf of Mexico, offshore Louisiana. High overpressures are present within low permeability mudstones where there have been multiple, very large, submarine landslides during the Pleistocene. Beneath 200 mbsf at Site U1324, pore pressures drop significantly: there are no submarine landslides in this mixture of mudstone, siltstone, and sandstone. The penetrometer measurements did not reach the in situ pressure at the end of the deployment. We used a soil model to determine that an extrapolation approach based on the inverse of square route of time (View the MathML source) requires much less decay time to achieve a desirable accuracy than an inverse time (1/t) extrapolation. Expedition 308 examined how rapid and asymmetric sedimentation above a permeable aquifer drives lateral fluid flow, extreme pore pressures, and submarine landslides. We interpret that the high overpressures observed are driven by rapid sedimentation of low permeability material from the ancestral Mississippi River. Reduced overpressure at depth at Site U1324 suggests lateral flow (drainage) whereas high overpressure at Site U1322 requires inflow from below: lateral flow in the underlying permeable aquifer provides one mechanism for these observations. High overpressure near the seafloor reduces slope stability and provides a mechanism for the large submarine landslides and low regional gradient (2°) offshore from the Mississippi delta.

Consolidation tests, elastic properties and grain size measurements at the consolidations samples from different Holes of IODP Expedition 302

The Integrated OceanDrilling Program's Expedition 302, the Arctic Coring Expedition (ACEX), recovered the first Cenozoic sedimentary sequence from the central Arctic Ocean. ACEX provided ground truth for basin scale geophysical interpretations and for guiding future exploration targets in this largely unexplored ocean basin. Here, we present results from a series of consolidation tests used to characterize sediment compressibility and permeability and integrate these with high-resolution measurements of bulk density, porosity and shear strength to investigate the stress history and the nature of prominent lithostratigraphic and seismostratigraphic boundaries in the ACEX record. Despite moderate sedimentation rates (10-30 m/Myr) and high permeability values (10**-15 -10**-18 m**2), consolidation and shear strength measurements both suggest an overall state of underconsolidation or overpressure. One-dimensional compaction modelling shows that to maintain such excess pore pressures, an in situ fluid source is required that exceeds the rate of fluid expulsion generated by mechanical compaction alone. Geochemical and sedimentological evidence is presented that identifes the Opal A-C/T transformation of biosiliceous rich sediments as a potential additional in situ fluid source.However, the combined rat of chemical and mechanical compaction remain too low to fully account for the observed pore pressure gradients, implying an additional diagenetic fluid source from within or below the recovered Cenozoic sediments from ACEX. Recognition of the Opal A-C/T reaction front in the ACEX record has broad reaching regional implications on slope stability and subsurface pressure evolution, and provides an important consideration for interpreting and correlating the spatially limited seismic data from the Arctic Ocean.

Common depth points and associated geographic coordinates

Fluids in subduction zones can influence seismogenic behaviour and prism morphology. The Eastern Makran subduction zone, offshore Pakistan, has a very thick incoming sediment section of up to 7.5 km, providing a large potential fluid source to the accretionary prism. A hydrate-related bottom simulating reflector (BSR), zones of high amplitude reflectivity, seafloor seep sites and reflective thrust faults are present across the accretionary prism, indicating the presence of fluids and suggesting active fluid migration. High amplitude free gas zones and seep sites are primarily associated with anticlinal hinge traps, and fluids here appear to be sourced from shallow biogenic sources and migrate to the seafloor along minor normal faults. There are no observed seep sites associated with the surface expression of the wedge thrust faults, potentially due to burial of the surface trace by failure of the steep thrust ridge slopes. Thrust fault reflectivity is restricted to the upper 3 km of sediment and the deeper décollement is non-reflective. We interpret that fluids and overpressure are not common in the deeper stratigraphic section. Thermal modelling of sediments at the deformation front suggests that the deeper sediment section is relatively dewatered and not currently contributing to fluid expulsion in the Makran accretionary prism.

Molecular and isotopic properties of gases from ODP Holes of ODP Leg 204

We report and discuss molecular and isotopic properties of hydrate-bound gases from 55 samples and void gases from 494 samples collected during Ocean Drilling Program (ODP) Leg 204 at Hydrate Ridge offshore Oregon. Gas hydrates appear to crystallize in sediments from two end-member gas sources (deep allochthonous and in situ) as mixtures of different proportions. In an area of high gas flux at the Southern Summit of the ridge (Sites 1248-1250), shallow (0-40 m below the seafloor [mbsf]) gas hydrates are composed of mainly allochthonous mixed microbial and thermogenic methane and a small portion of thermogenic C2+ gases, which migrated vertically and laterally from as deep as 2- to 2.5-km depths. In contrast, deep (50-105 mbsf) gas hydrates at the Southern Summit (Sites 1248 and 1250) and on the flanks of the ridge (Sites 1244-1247) crystallize mainly from microbial methane and ethane generated dominantly in situ. A small contribution of allochthonous gas may also be present at sites where geologic and tectonic settings favor focused vertical gas migration from greater depth (e.g., Sites 1244 and 1245). Non-hydrocarbon gases such as CO2 and H2S are not abundant in sampled hydrates. The new gas geochemical data are inconsistent with earlier models suggesting that seafloor gas hydrates at Hydrate Ridge formed from gas derived from decomposition of deeper and older gas hydrates. Gas hydrate formation at the Southern Summit is explained by a model in which gas migrated from deep sediments, and perhaps was trapped by a gas hydrate seal at the base of the gas hydrate stability zone (GHSZ). Free gas migrated into the GHSZ when the overpressure in gas column exceeded sealing capacity of overlaying sediments, and precipitated as gas hydrate mainly within shallow sediments. The mushroom-like 3D shape of gas hydrate accumulation at the summit is possibly defined by the gas diffusion aureole surrounding the main migration conduit, the decrease of gas solubility in shallow sediment, and refocusing of gas by carbonate and gas hydrate seals near the seafloor to the crest of the local anticline structure.