Geothermal evolution of an intruded dike in the rift zone of Kilauea volcano, Hawaii from VLF and self-potential measurements

NASA Astrophysics Data System (ADS)

Davis, Paul M.

2015-09-01

Self-potential (SP) and VLF measurements were made in 1973, 1975, 1995, 1997 and 2012 across a basaltic dike that intruded into the Koae fault zone of Kilauea volcano, Hawaii in May 1973. The SP anomaly remained strong throughout. In 2012 it was at about 60% of the strength it had in 1973. In contrast, the VLF anomaly, though diminished, was still observable in 1995/1997, but by 2012 it had disappeared. A hydrothermal dike model, with parameters calibrated by modeling the solidification of Kilauea Iki lava lake, is used to calculate temperatures and conductivity variation. Following Jaeger's (1957) method, we find that the time in years for a dike of width W (m) to solidify is 0.0075W2. Thus, a 1 m dike solidifies within the first few days, and after 39 years is only tens of degrees above ambient. Given the orders of magnitude difference between the conductivities of wet and dry basalt, we infer, that after solidification, the VLF anomalies were caused by induction in a localized veil of wet, hot basalt enveloping the dike, that was generated initially by condensation of steam, and subsequently by condensation of evaporated water as temperatures reduced. The conductivity anomaly persisted until the mid-nineties. By 2012, temperatures and condensation were too small for a VLF signal. The persistent SP anomaly is attributed to localized fluid disruption, with evaporation mainly at the water table and in the vadose zone. Streaming potentials are associated with evaporative circulation in the vadose zone. Next to the dike a positive potential is generated by upward flow of moisture-laden air, with a smaller negative potential on its flanks from downward infiltrating rainwater. The analysis indicates that the combination of SP and VLF measurements can characterize the evolving geothermal regime of intrusions above the water table.

Conduits and dike distribution analysis in San Rafael Swell, Utah

NASA Astrophysics Data System (ADS)

Kiyosugi, K.; Connor, C.; Wetmore, P. H.; Ferwerda, B. P.; Germa, A.

2011-12-01

Volcanic fields generally consist of scattered monogenetic volcanoes, such as cinder cones and maars. The temporal and spatial distribution of monogenetic volcanoes and probability of future activity within volcanic fields is studied with the goals of understanding the origins of these volcano groups, and forecasting potential future volcanic hazards. The subsurface magmatic plumbing systems associated with volcanic fields, however, are rarely observed or studied. Therefore, we investigated a highly eroded and exposed magmatic plumbing system on the San Rafael Swell (UT) that consists of dikes, volcano conduits and sills. San Rafael Swell is part of the Colorado Plateau and is located east of the Rocky Mountain seismic belt and the Basin and Range. The overburden thickness at the time of mafic magma intrusion (Pliocene; ca. 4 Ma) into Jurassic sandstone is estimated to be ~800 m based on paleotopographical reconstructions. Based on a geologic map by P. Delaney and colleagues, and new field research, a total of 63 conduits are mapped in this former volcanic field. The conduits each reveal features of root zone and / or lower diatremes, including rapid dike expansion, peperite and brecciated intrusive and host rocks. Recrystallized baked zone of host rock is also observed around many conduits. Most conduits are basaltic or shonkinitic with thickness of >10 m and associated with feeder dikes intruded along N-S trend joints in the host rock, whereas two conduits are syenitic and suggesting development from underlying cognate sills. Conduit distribution, which is analyzed by a kernel function method with elliptical bandwidth, illustrates a N-S elongate higher conduit density area regardless of the azimuth of closely distributed conduits alignment (nearest neighbor distance <200 m). In addition, dike density was calculated as total dike length in unit area (km/km^2). Conduit and sill distribution is concordant with the high dike density area. Especially, the distribution

Constraints on dike propagation from continuous GPS measurements

NASA Astrophysics Data System (ADS)

Segall, Paul; Cervelli, Peter; Owen, Susan; Lisowski, Mike; Miklius, Asta

2001-09-01

The January 1997 East Rift Zone eruption on Kilauea volcano, Hawaii, occurred within a network of continuous Global Positioning System (GPS) receivers. The GPS measurements reveal the temporal history of deformation during dike intrusion, beginning ˜8 hours prior to the onset of the eruption. The dike volume as a function of time, estimated from the GPS data using elastic Green's functions for a homogeneous half-space, shows that only two thirds of the final dike volume accumulated prior to the eruption and the rate of volume change decreased with time. These observations are inconsistent with simple models of dike propagation, which predict accelerating dike volume up to the time of the eruption and little or no change thereafter. Deflationary tilt changes at Kilauea summit mirror the inferred dike volume history, suggesting that the rate of dike propagation is limited by flow of magma into the dike. A simple, lumped parameter model of a coupled dike magma chamber system shows that the tendency for a dike to end in an eruption (rather than intrusion) is favored by high initial dike pressures, compressional stress states, large, compressible magma reservoirs, and highly conductive conduits linking the dike and source reservoirs. Comparison of model predictions to the observed dike volume history, the ratio of erupted to intruded magma, and the deflationary history of the summit magma chamber suggest that most of the magma supplied to the growing dike came from sources near to the eruption through highly conductive conduits. Interpretation is complicated by the presence of multiple source reservoirs, magma vesiculation and cooling, as well as spatial variations in dike-normal stress. Reinflation of the summit magma chamber following the eruption was measured by GPS and accompanied a rise in the level of the Pu'u O'o lava lake. For a spheroidal chamber these data imply a summit magma chamber volume of ˜20 km3, consistent with recent estimates from seismic

Constraints on dike propagation from continuous GPS measurements

USGS Publications Warehouse

Segall, P.; Cervelli, Peter; Owen, S.; Lisowski, M.; Miklius, Asta

2001-01-01

The January 1997 East Rift Zone eruption on Kilauea volcano, Hawaii, occurred within a network of continuous Global Positioning System (GPS) receivers. The GPS measurements reveal the temporal history of deformation during dike intrusion, beginning ??? 8 hours prior to the onset of the eruption. The dike volume as a function of time, estimated from the GPS data using elastic Green's functions for a homogeneous half-space, shows that only two thirds of the final dike volume accumulated prior to the eruption and the rate of volume change decreased with time. These observations are inconsistent with simple models of dike propagation, which predict accelerating dike volume up to the time of the eruption and little or no change thereafter. Deflationary tilt changes at Kilauea summit mirror the inferred dike volume history, suggesting that the rate of dike propagation is limited by flow of magma into the dike. A simple, lumped parameter model of a coupled dike magma chamber system shows that the tendency for a dike to end in an eruption (rather than intrusion) is favored by high initial dike pressures, compressional stress states, large, compressible magma reservoirs, and highly conductive conduits linking the dike and source reservoirs. Comparison of model predictions to the observed dike volume history, the ratio of erupted to intruded magma, and the deflationary history of the summit magma chamber suggest that most of the magma supplied to the growing dike came from sources near to the eruption through highly conductive conduits. Interpretation is complicated by the presence of multiple source reservoirs, magma vesiculation and cooling, as well as spatial variations in dike-normal stress. Reinflation of the summit magma chamber following the eruption was measured by GPS and accompanied a rise in the level of the Pu'u O'o lava lake. For a spheroidal chamber these data imply a summit magma chamber volume of ??? 20 km3, consistent with recent estimates from seismic

Volcanotectonic earthquakes induced by propagating dikes

NASA Astrophysics Data System (ADS)

Gudmundsson, Agust

2016-04-01

Volcanotectonic earthquakes are of high frequency and mostly generated by slip on faults. During chamber expansion/contraction earthquakes are distribution in the chamber roof. Following magma-chamber rupture and dike injection, however, earthquakes tend to concentrate around the dike and follow its propagation path, resulting in an earthquake swarm characterised by a number of earthquakes of similar magnitudes. I distinguish between two basic processes by which propagating dikes induce earthquakes. One is due to stress concentration in the process zone at the tip of the dike, the other relates to stresses induced in the walls and surrounding rocks on either side of the dike. As to the first process, some earthquakes generated at the dike tip are related to pure extension fracturing as the tip advances and the dike-path forms. Formation of pure extension fractures normally induces non-double couple earthquakes. There is also shear fracturing in the process zone, however, particularly normal faulting, which produces double-couple earthquakes. The second process relates primarily to slip on existing fractures in the host rock induced by the driving pressure of the propagating dike. Such pressures easily reach 5-20 MPa and induce compressive and shear stresses in the adjacent host rock, which already contains numerous fractures (mainly joints) of different attitudes. In piles of lava flows or sedimentary beds the original joints are primarily vertical and horizontal. Similarly, the contacts between the layers/beds are originally horizontal. As the layers/beds become buried, the joints and contacts become gradually tilted so that the joints and contacts become oblique to the horizontal compressive stress induced by a driving pressure of the (vertical) dike. Also, most of the hexagonal (or pentagonal) columnar joints in the lava flows are, from the beginning, oblique to an intrusive sheet of any attitude. Consequently, the joints and contacts function as potential shear

Modeling Thermal Pressurization Around Shallow Dikes Using Temperature-Dependent Hydraulic Properties: Implications for Deformation Around Intrusions

NASA Astrophysics Data System (ADS)

Townsend, Meredith R.

2018-01-01

Pressurization and flow of groundwater around igneous intrusions depend in part on the hydraulic diffusivity of the host rocks and processes that enhance diffusivity, such as fracturing, or decrease diffusivity, such as mineral precipitation during chemical alteration. Characterizing and quantifying the coupled effects of alteration, pore pressurization, and deformation have significant implications for deformation around intrusions, geothermal energy, contact metamorphism, and heat transfer at mid-ocean ridges. Fractures around dikes at Ship Rock, New Mexico, indicate that pore pressures in the host rocks exceeded hydrostatic conditions by at least 15 MPa following dike emplacement. Hydraulic measurements and petrographic analysis indicate that mineral precipitation clogged the pores of the host rock, reducing porosity from 0.25 to <0.10 and reducing permeability by 5 orders of magnitude. Field data from Ship Rock are used to motivate and constrain numerical models for thermal pore fluid pressurization adjacent to a meter-scale dike, using temperature-dependent hydraulic properties in the host rock as a proxy for porosity loss by mineral precipitation during chemical alteration. Reduction in permeability by chemical alteration has a negligible effect on pressurization. However, reduction in porosity by mineral precipitation increases fluid pressure by constricting pore volume and is identified as a potentially significant source of pressure. A scaling relationship is derived to determine when porosity loss becomes important; if permeability is low enough, pressurization by porosity loss outweighs pressurization by thermal expansion of fluids.

Quantification of the intrusion process at Kīlauea volcano, Hawai'i

NASA Astrophysics Data System (ADS)

Wright, Thomas L.; Marsh, Bruce

2016-12-01

The characteristic size of two types of intrusions identified beneath Kīlauea's East Rift zone are uniquely estimated by combining time constraints from fractional crystallization and the rates of magma solidification during cooling. Some intrusions were rapidly emplaced as dikes, but stalled before reaching the surface, and cooled and crystallized to feed later fractionated eruptions. More specifically, using the observed time interval between initial emplacement and eruption of fractionated lava, whose degree of fractionation is estimated from petrologic mixing calculations, the extent of solidification or cooling needed to produce this amount of fractionation can be directly inferred. And from the known erupted volumes the spatial extent or size of this fractionated volume can be analytically related to the full size of the source body itself. Two examples yield dike widths of 82 and 68 m. Other intrusions remain close to the east rift magma transport path and are observed to last for decades or longer as viable magma bodies that may participate in feeding later eruptions. The thickness of semi-permanent reservoirs near the East Rift Zone magma transport path can be estimated by assuming a resupply rate that is sufficiently frequent to restrict cooling to < 10 °C. It is inferred that both types of intrusions likely began as dike offshoots from the East Rift Zone magma transport path, but the frequently resupplied bodies may have later been converted to sills or laccoliths of heights estimated at 43-62 m. Our modeled intrusions contrast with models of rapidly emplaced thinner dikes feeding shallow intrusions, which are accompanied by intense rift earthquake swarms and are often associated with eruptions. These calculations show that long-term heating of the wallrock of the magma transport paths serves to slow conduit cooling, which may be partly responsible for sustaining long East Rift Zone eruptions. Adjacent to the vertical transport path beneath K�

Coupled heat and silica transport associated with dike intrusion into sedimentary rock: effects on isotherm location and permeability evolution

NASA Astrophysics Data System (ADS)

Dutrow, Barbara L.; Travis, Bryan J.; Gable, Carl W.; Henry, Darrell J.

2001-11-01

An 11-meter-wide alkalic monchiquite dike recovered from the subsurface of Louisiana has produced a metasomatic aureole in the adjacent interbedded carbonate mudstones and siltstones. The asymmetric contact aureole, which extends nearly 6 m above and 4 m below the intrusion, contains the metamorphic minerals, diopside, pectolite, fluor-apophyllite, fluorite, and garnet. A series of coupled heat and mass transport calculations was undertaken to provide thermal constraints for the aureole, in the absence of robust geothermometric assemblages, and insights into accompanying mass transport associated with the sedimentary rock-dike system. Calculations were completed for systems with homogeneous, anisotropic, and layered permeability, κ. Transport, dissolution, and precipitation of silica were also incorporated into calculations. All systems modeled indicate that the thermal pulse waned in ∼3 yr with a return to background temperatures in ∼10 yr. Heat and fluid transport produce maximum temperature isotherms that are distinctly different in spatial extent and lateral variability for each numerical system. The homogeneous κ case produced isotherms that pinch and swell vertically above the dike and have large lateral variations, in contrast to the anisotropic κ case that produced a single large plume above the dike. The layered system κ case produced the most spatially extensive thermal aureole, unlike that recorded in the rocks. Addition of dissolved silica to the flow system significantly impacts the calculated transport of heat and fluid, primarily due to density changes that affect upwelling dynamics. Although precipitation and dissolution of SiO2 can affect flow through the feedback to permeability, κ changes were found to be minor for these system conditions. Where κ decreased, flow was refocused into higher κ zones, thus mitigating the κ differences over time. This negative feedback tends to defocus flow and provides a mechanism for lateral migration of

Dike intrusions into bituminous coal, Illinois Basin: H, C, N, O isotopic responses to rapid and brief heating

USGS Publications Warehouse

Schimmelmann, A.; Mastalerz, Maria; Gao, L.; Sauer, P.E.; Topalov, K.

2009-01-01

Unlike long-term heating in subsiding sedimentary basins, the near-instantaneous thermal maturation of sedimentary organic matter near magmatic intrusions is comparable to artificial thermal maturation in the laboratory in terms of short duration and limited extent. This study investigates chemical and H, C, N, O isotopic changes in high volatile bituminous coal near two Illinois dike contacts and compares observed patterns and trends with data from other published studies and from artificial maturation experiments. Our study pioneers in quantifying isotopically exchangeable hydrogen and measuring the D/H (i.e., 2H/1H) ratio of isotopically non-exchangeable organic hydrogen in kerogen near magmatic contacts. Thermal stress in coal caused a reduction of isotopically exchangeable hydrogen in kerogen from 5% to 6% in unaltered coal to 2-3% at contacts, mostly due to elimination of functional groups (e.g., {single bond}OH, {single bond}COOH, {single bond}NH2). In contrast to all previously published data on D/H in thermally matured organic matter, the more mature kerogen near the two dike contacts is D-depleted, which is attributed to (i) thermal elimination of D-enriched functional groups, and (ii) thermal drying of hydrologically isolated coal prior to the onset of cracking reactions, thereby precluding D-transfer from relatively D-enriched water into kerogen. Maxima in organic nitrogen concentration and in the atomic N/C ratio of kerogen at a distance of ???2.5 to ???3.5 m from the thicker dike indicate that reactive N-compounds had been pyrolytically liberated at high temperature closer to the contact, migrated through the coal seam, and recombined with coal kerogen in a zone of lower temperature. The same principle extends to organic carbon, because a strong ??13Ckerogen vs. ??15Nkerogen correlation across 5.5 m of coal adjacent to the thicker dike indicates that coal was functioning as a flow-through reactor along a dynamic thermal gradient facilitating back

Dike propagation energy balance from deformation modeling and seismic release

NASA Astrophysics Data System (ADS)

Bonaccorso, Alessandro; Aoki, Yosuke; Rivalta, Eleonora

2017-06-01

Magma is transported in the crust mainly by dike intrusions. In volcanic areas, dikes can ascend toward the free surface and also move by lateral propagation, eventually feeding flank eruptions. Understanding dike mechanics is a key to forecasting the expected propagation and associated hazard. Several studies have been conducted on dike mechanisms and propagation; however, a less in-depth investigated aspect is the relation between measured dike-induced deformation and the seismicity released during its propagation. We individuated a simple x that can be used as a proxy of the expected mechanical energy released by a propagating dike and is related to its average thickness. For several intrusions around the world (Afar, Japan, and Mount Etna), we correlate such mechanical energy to the seismic moment released by the induced earthquakes. We obtain an empirical law that quantifies the expected seismic energy released before arrest. The proposed approach may be helpful to predict the total seismic moment that will be released by an intrusion and thus to control the energy status during its propagation and the time of dike arrest.