Deep geoelectrical models: geological and electromagnetic principles

An important result of recent years is the normal resistivity profile. It was obtained by interpretation of the combined apparent resistivity curve (magnetotelluric sounding and geomagnetic deep sounding) for the East European platform. This profile has no highly conducting layer and resistivity is greater than 100 ohm-m at asthenospheric depths. It corresponds well with geothermal indications of the absence of partial melting beneath the Precambrian plates. Nearly the same profiles have been obtained for the Canadian shield, and the Siberian and Australian platforms. Investigations carried out in many “hot” regions confirm the existence of a well-developed low-resistivity asthenosphere. Partially molten zones have conductances of about several thousand Siemens in the Eastern Pacific, Iceland and in the North American rift zone. Magnetotelluric soundings show that in many continental areas the lower part of the crust has low resistivity, in the range 10–20 ohm-m. Usually this crustal conductive layer is observed in regions of recent activity. Its total conductivity changes from several hundred to several thousand S. Many investigators propose that the most natural explanation of this conductivity is water solutions.It is necessary to note the distorting role of near-surface inhomogeneities. Local distortions can be eliminated by simple averaging of the experimental data. These average apparent resistivity curves are the starting point for the construction of deep geoelectrical models.     

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Réunion)

 A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Réunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60–100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos, the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise. Received: 3 November 1999 / Accepted: 15 September 1999     

The main structure and tectonic features of the Chernorud-Mukhor site on the western shore of Lake Baikal from TEM and SP measurements

Transient electromagnetic (TEM), self-potential (SP) and geoelectrical mapping measurements were carried out at the Chernorud-Mukhor site in the Priolkhonje area on the western shore of the Lake Baikal. All measurements were made along several profiles across the main strike of the regional Primorsky fault. TEM measurements were carried out in a time range from a few tens of microseconds to several tens of milliseconds. The most important result of the 1D modelling of TEM soundings is the discovery of nearly horizontal boundaries that divide high resistive overlying and well conducting underlying rocks. The resistivity of the former is in the range from 100 Ωm to 1000 Ωm, while the resistivity of the latter varies from less than 1 Ωm to several tens of Ωm. This good conductive zone could also be verified by geoelectrical mapping using Schlumberger array (AB/2=100 m). Due to high conductivity of the underlying rocks only the upper boundary of the conductive layer could be determined by TEM soundings. A regional SP anomaly with amplitude of about −450 mV has also been observed above the low resistivity zone indicating the electron nature of its conductance. Geologically, the conductive zone is represented by a graphite-bearing layer within the region of archean rocks. Since that layer extends over a large area, it may be used as a key in studying structures and tectonics of the Priolkhonje area. A 1D TEM geoelectric section shows a wide, gently sloping syncline as a probable base structure of the Chernorud-Mukhor site. Neotectonic faults divide the syncline into vertically displaced blocks that form a wide complicated graben with a total amplitude of about 250 m.     

Crustal electrical conductivity of the Indian continental subduction zone: New data from the profile in the Garhwal Himalaya

We present the results of studying the geoelectrical structure of the zone of continental subduction of the Indian lithospheric plate within the Gahrwal Himalaya. In the framework of the Russian–Indian project, the data of the broadband magnetotelluric soundings conducted by the Indian Institute of Technology Roorkee on the regional profile across the structures of the orogen were expanded, processed, and interpreted by the new program tools adapted for the measurements in the mountain conditions and for the presence of industrial noise. The constructed model of the deep electrical conductivity cross section for Garhwal revealed its two-dimensional (2D) features and more accurately delineated the location of the midcrustal conductor associated with the ramp structure of the detachment plane. The correlations with the regional distribution of the earthquake hypocenters and the seismotomographic images suggest a common, fluid-related nature of the seismic and geoelectrical anomalies in the crust of the Garhwal Tectonic Corridor and enabled the identification of the seismogenerating zones. Among the data of the expanded profile set of magnetotelluric and magnetovariational transfer functions, the response of a poorly explored deep conductive body is revealed. This object is located east of the profile and is probably associated with the activation of the ancient trans-Himalayan cratonic structures which prepares the segmentation of the Himalayan arc.     

An oxygen isotope study of hydrothermal alteration in the Lake City caldera, San Juan Mountains, Colorado

A 23-m.y.-old, fossil meteoric-hydrothermal system in the Lake City caldera (11 × 14 km) has been mapped out by measuring δ ¹⁸O values of 300 rock and mineral samples. δ ¹⁸O varies systematically throughout the caldera, reaching values as low as −2. Great topographic relief, regional tilting, and variable degrees of erosion within the caldera all combine to give us a very complete section through the hydrothermal system, from the surface down to a depth of more than 2000 m. The initial δ ¹⁸O value of the caldera-fill Sunshine Peak Tuff was very uniform (+7.2 ± 0.1), making it easy to determine the exact amount of ¹⁸O depletion experienced by each sample during hydrothermal alteration. Also, we have excellent stratigraphic control on depths beneath the mid-Tertiary surface, quantitative information on mineralogical alteration products, and accurate data on the shape of the central resurgent intrusion, which was the principal ‘heat engine’ that drove the hydrothermal circulation. Major conclusions are: (1) Although pristine mid-Tertiary meteoric waters in this area had δ ¹⁸O −14, these fluids were ¹⁸O-shifted upward to about δ¹⁸O = −8 to −5 prior to entering the shallow convective system associated with the resurgent intrusive rocks. Although there was undoubtedly radial inflow toward the caldera from all directions, the highly fractured Eureka Graben, southwest of the caldera, was probably the principal source of recharge groundwater for the Lake City system. (2) Fluid flow within the caldera was dominated by three major categories of permeable zones: the porous megabreccia units (which dip outward from the resurgent dome), vertical fractures and faults related to resurgence, and the caldera ring fault itself. All of these zones exhibit marked ¹⁸O depletions, and they are also typically intensely mineralogically altered. (3) The resurgent intrusive stock and its contact metamorphic aureole of hornfels both experienced water/rock ratios lower than the permeable zones; however, they have similarly low δ ¹⁸O values because they were altered at higher temperatures. (4) Throughout the caldera, the δ ¹⁸O of Sunshine Peak Tuff decreases with increasing depth (about 6 per mil/km), indicative of a shallow thermal gradient, typical of a convective hydrothermal system. The near-surface portion of this gradient was controlled by the temperature drop associated with boiling in the uprising fluid. (5) Deeply circulating meteoric water rose along permeable ring fractures 3 to 5 km beneath the mid-Tertiary surface. These fluids were drawn into the shallow convective system through the lower, porous, megabreccia units. Near the resurgent intrusions, fluid flow was again directed upward where resurgence-related, near-vertical fractures intersect the megabreccia units.     

Recent activity of Nisyros volcano (Greece) inferred from structural, geochemical and seismological data

This study summarizes the results of structural, geochemical and seismological surveys carried out at Nisyros volcano (Aegean Sea, Greece) during 1999–2001. Field mapping and mesostructural measurements at the summit caldera (Lakki plain) indicate that faults follow two main strikes: NE-SW and N-S. The N-S striking fault depicts extensional features accommodating the left-lateral component of motion of the NE-SW- striking main faults. The NE-SW preferred strike of the Lakki faults and of the mineral-filled veins as well as the distribution and NE-SW elongation of the hydrothermal craters indicate that tectonics plays a major role in controlling the fluid pathway in the Nisyros caldera. The same NE-SW trend is depicted by CO₂ anomalies revealed through detailed soil CO₂ flux surveys, thus indicating a structural control on the pattern of the hydrothermal degassing. Degassing processes account for a thermal energy release of about 43 MW, most of which occurs at Lofos dome, an area that was affected by hydrothermal eruptions in historical times. The seismic study was conducted in June 2001, using a deployment specifically aimed at detecting signals of magmatic-hydrothermal origin. Our instruments recorded local and regional earthquakes, a few local long-period events (LP), and bursts of monochromatic tremor. Local earthquake activity is concentrated beneath the caldera, at depths generally shallower than 6 km. Plane-wave decomposition of tremor signal indicates a shallow (<200 m) source located in the eastern part of the caldera. Conversely, LP events depict a source located beneath the central part of the caldera, in the area of Lofos dome, at depths in the 1–2-km range. In agreement with geochemical and structural measurements, these data suggest that both the deeper and shallower part of the hydrothermal system are subjected to instability in the fluid flow regimes, probably consequent to transient pressurization of the reservoir. These instabilities may be related to input of hot fluids from the deeper magmatic system, as suggested by the variations in geochemical parameters observed after the 1997–1999 unrest episode. The significance of seismological and geochemical indicators as precursors of hydrothermal explosive activity at Nisyros is discussed.Editorial responsibility: H. Shinohara     

Multiple caldera collapses inferred from the shallow electrical resistivity signature of the Las Cañadas caldera, Tenerife, Canary Islands

The Las Cañadas caldera of Tenerife (LCC) is a well exposed caldera depression filled with pyroclastic deposits and lava flows from the active Teide–Pico Viejo complex (TPVC). The caldera's origin is controversial as both the formation by huge lateral flank collapse(s) and multiple vertical collapses have been proposed. Although vertical collapses may have facilitated lateral slope failures and thus jointly contribute to the exposed morphology, their joint contribution has not been clearly demonstrated. Using results from 185 audiomagnetotelluric (AMT) soundings carried out between 2004 and 2006 inside the LCC, our study provides consistent geophysical constraints in favour of multiple vertical caldera collapse. One-dimensional modelling reveals a conductive layer at shallow depth (30–1000 m), presumably resulting from hydrothermal alteration and weathering, underlying the infilling resistive top layer. We present the resistivity distribution of both layers (resistivity images), the topography of the conductive layer across the LCC, as well as a cross-section in order to highlight the caldera's evolution, including the distribution of earlier volcanic edifices. The AMT phase anisotropy reveals the structural and radial characteristics of the LCC.     

Structure of the Nemrut caldera (Eastern Anatolia, Turkey) and associated hydrothermal fluid circulation

Plio-Quaternary volcanism played an important role in the present physical state of Eastern Anatolia. Mount Nemrut, situated to the west of Lake Van is one of the main volcanic centers in the region, with a spectacular summit caldera 8.5 × 7 km in diameter. The most recent eruptions of the volcano were in 1441, 1597 and 1692. Nemrut Lake covers the western half of the caldera; it is a deep, half-bowl-shaped lake with a maximum depth of 176 m. Numerous eruption centers are exposed within the caldera as a consequence of magma–water interaction. Current activity of Nemrut caldera is revealed as hot springs, fumaroles and a small, hot lake.Self-potential and bathymetric surveys carried out in the caldera were used to characterize the structure of the caldera and the associated hydrothermal fluid circulation. In addition, analyses based on digital elevation models and satellite imagery were used to improve our knowledge about the structure of the caldera. According to SP results, the flanks of the volcano represent “the hydrogeologic zone”, whereas the intra-caldera region is an “active hydrothermal area” where the fluid circulation is controlled by structural discontinuities. There is also a northern fissure zone which exhibits hydrothermal signatures. Nemrut caldera collapsed piecemeal, with three main blocks. Stress controlling the collapse mechanism seems to be highly affected by the regional neotectonic regime. In addition to the historical activity, current hydrothermal and hydrogeologic conditions in the caldera, in which there is a large lake and shallow water table, increase the risk of the quiescent volcano.     

Fluid circulation and structural discontinuities inside Misti volcano (Peru) inferred from self-potential measurements

One of the seven potentially active andesite stratovolcanoes in southern Peru, Misti (5822 m), located 17 km northeast and 3.5 km above Arequipa, represents a major threat to the population (900,000 inhabitants). Our recent geophysical and geochemical research comprises an extensive self-potential (SP) data set, an audio–magnetotelluric (AMT) profile across the volcano and CO₂ concentrations in the soil along a radial profile. The SP survey is the first of its kind in providing a complete mapping of a large andesitic stratovolcano 20 km in diameter. The SP mapping enables us to analyze the SP signature associated with a subduction-related active volcano.The general SP pattern of Misti is similar to that of most volcanoes with a hydrogeologic zone in the lower flanks and a hydrothermal zone in the upper central area. A quasi-systematic relationship exists between SP and elevation. Zones with constant SP/altitude gradients (Ce) are observed in both hydrogeologic (negative Ce) and hydrothermal (positive Ce) zones. Transition zones between the different Ce zones, which form a concentric pattern around the summit, have been interpreted in terms of lateral heterogeneities in the lithology. The highest amplitudes of SP anomalies seem to coincide with highly resistive zones. The hydrothermal system 6 km in diameter, which extends over an area much larger than the summit caldera, may be constrained by an older, concealed collapse caldera. A sealed zone has apparently developed through alteration in the hydrothermal system, blocking the migration of CO₂ upward. Significant CO₂ emanations are thus observed on the lower flanks but are absent above the hydrothermal zone.     

On the Locality of Transient Electromagnetic Soundings with a Single-Loop Configuration

The possibilities of reconstructing two-dimensional (2D) cross sections based on the data of the profile soundings by the transient electromagnetic method (TEM) with a single ungrounded loop are illustrated on three-dimensional (3D) models. The process of reconstruction includes three main steps: transformation of the responses in the depth dependence of resistivity ρ(h) measured along the profile, with their subsequent stitching into the 2D pseudo section; point-by-point one-dimensional (1D) inversion of the responses with the starting model constructed based on the transformations; and correction of the 2D cross section with the use of 2.5-dimensional (2.5D) block inversion. It is shown that single-loop TEM soundings allow studying the geological media within a local domain the lateral dimensions of which are commensurate with the depth of the investigation. The structure of the medium beyond this domain insignificantly affects the sounding results. This locality enables the TEM to reconstruct the geoelectrical structure of the medium from the 2D cross sections with the minimal distortions caused by the lack of information beyond the profile of the transient response measurements.     

Fluid circulation in a complex volcano-tectonic setting,inferred from self-potential and soil CO2 flux surveys: The Santa María–Cerro Quemado–Zunil volcanoes and Xela caldera (Northwestern Guatemala)

The region encompassing Santa María, Cerro Quemado, and Zunil volcanoes, close to Quetzaltenango, the second largest city of Guatemala, is volcanically and tectonically complex. In addition, the huge Xela caldera, about 20 km in diameter, crosses this area and links up to the important Zunil fault zone located between the three volcanoes. Two highly active geothermal sites, named Zunil-I and Zunil-II, are also located between these three volcanic edifices at the southeastern boundary of Xela caldera. In order to determine the permeability variations and the main structural discontinuities within this complex volcano-tectonic setting, self-potential and soil CO₂ flux measurements have been coupled, with a step of 20 m, along a 16.880 km-long profile crossing the entire area. Two shallow hydrothermal systems, with maximum lateral extensions of 1.5 km in diameter, are indicated by positive self-potential/elevation gradients below Santa María and Cerro Quemado volcanoes. Such small hydrothermal systems cannot explain the intense geothermal manifestations at Zunil-I and Zunil-II. Another minor hydrothermal system is indicated by self-potential measurements on the flank of Santa María along the edge of the Xela caldera. CO₂ flux measurements display slight variations inside the caldera and decreasing values crossing outside the caldera boundary. We hypothesize the presence of a magmatic body, inside the southeastern border of Xela caldera, to explain the deeper and more intense hydrothermal system manifested by the Zunil-I and the Zunil-II geothermal fields. This magmatic system may be independent from Santa María and Cerro Quemado volcanoes. Alternatively, the hypothesized Xela magmatic system could have a common magmatic origin with the Cerro Quemado dome complex, consistent with previous findings on regional gas emissions. Sectors bordering the Cerro Quemado dome complex also have high amplitude minima-short wavelength anomalies in self-potential, interpreted as preferential rain water infiltration along faults of major permeability, probably related with the most recent stages of Cerro Quemado dome growth.     

Fluid inclusion evidence for rapid formation of the vapor-dominated zone at Sulphur Springs, Valles caldera, New Mexico, USA

Microthermometric measurements were obtained for 618 fluid inclusions in hydrothermal quartz, fluorite and calcite and magmatic quartz phenocrysts in intracaldera tuffs from the VC-2A core hole in order to study evolutionary processes of the Sulphur Springs hydrothermal system in the Valles caldera. Relatively high T_h values in samples from shallow depths indicate erosion of about 200 m of caldera fill since deposition of hydrothermal minerals at shallow depths in the Sulphur Springs hydrothermal system, accompanied by a descent in the water table of the liquid-dominated reservoir. For samples collected below the current water level of the well, the minimum values of homogenization temperature (T_h) fit the present thermal profile, whereas minimum T_h values of samples from above the water level are several tens of degrees higher than the present thermal profile and fit a paleo-thermal profile following the boiling point curve for pure water, as adjusted to 92 °C at 20 m below the present land surface. This is attributed to development of an evolving vapor zone that formed subsequent to a sudden drop in the water table of the liquid-dominated reservoir. We suggest that these events were caused by the drainage of an intracaldera lake when the southwestern wall of the caldera was breached about 0.5 Ma. This model indicates that vapor zones above major liquid-dominated geothermal reservoirs can be formed due to dramatic changes in geohydrology and not just from simple boiling.     

The role of regional-scale faults in controlling a trapdoor caldera,Coromandel Peninsula,New Zealand

Regional gravity data from an eroded Miocene to Pliocene volcanic arc exposed in the Coromandel Peninsula, New Zealand, delineate a circular − 26-mGal, 15-km-diameter gravity anomaly. This anomaly, which has steep gradients on its northern and western margins but shallow gradients elsewhere, correlates with relatively young volcanic and volcaniclastic rocks within a broad topographic depression. Gravity modelling, using an exponentially decreasing density contrast with depth profile, requires very low-density rocks (ca. 2280 kg m^− 3) in the near-surface to account for the observed anomaly, giving a total depth of ca. 2.8 km for these rocks. The northern and western margins of this body dip steeply inward at 70°, whereas the southern and eastern margins have shallow inward dips (20–30°). The western margin coincides with the regional-scale Mangakino Fault, but the northern margin, recognizable only in the geophysical data (and named here the Ohinemuri Fault), is partially buried under younger volcanic rocks. We interpret these deep and steeply bounded, low-density volcanics in terms of a trapdoor caldera, faulted on its northern and western margins, with its hinge on the southern and eastern margins. Epithermal deposits are spatially associated with the Mangakino and Ohinemuri Faults, suggesting that both structures may have influenced hydrothermal fluid flow. These deposits pre-date caldera fill, indicating that caldera development followed pre-existing regional faults. These results delineate the subsurface geometry of a trapdoor caldera and highlight the role of pre-existing, regional-scale faults in controlling such caldera location and collapse.     

Resistivity measurements on the sea bottom to map fracture zones in the bedrock underneath sediments1

Refraction seismics with the shotpoints and the hydrophone cable on the sea-bottom, have become the standard geophysical method for investigating rock quality before constructing offshore tunnels in Norway. In connection with the construction of a sub-sea tunnel by the Norwegian Public Roads Administration, research work was carried out to compare two low-velocity zones, indicated by refraction seismics with other methods. A special resistivity cable for pole-dipole measurements on the sea-floor, with 10 m between the electrodes, was constructed. A 200 m long profile, crossing the two low-velocity zones, was measured with all combinations of electrode distances. The two zones were detected as low-resistivity zones. A special data processing technique to enhance the anomalies is described. Resistivity soundings in a seawater environment to detect fracture zones in the bedrock underneath the bottom sediments, are discussed. It is concluded that severely fractured zones, which may cause difficulties for the tunnel construction, can be detected both with sea-surface and sea-floor arrays using long electrode spacings.     

New evidence on the hydrothermal system in Long Valley caldera, California, from wells, fluid sampling, electrical geophysics, and age determinations of hot-spring deposits

Data collected since 1985 from test drilling, fluid sampling, and geologic and geophysical investigations provide a clearer definition of the hydrothermal system in Long Valley caldera than was previously available. This information confirms the existence of high-temperature (> 200°C) reservoirs within the volcanic fill in parts of the west moat. These reservoirs contain fluids which are chemically similar to thermal fluids encountered in the central and eastern parts of the caldera. The roots of the present-day hydrothermal system (the source reservoir, principal zones of upflow, and the magmatic heat source) most likely occur within metamorphic basement rocks beneath the western part of the caldera. Geothermometer-temperature estimates for the source reservoir range from 214 to 248°C. Zones of upflow of hot water could exist beneath the plateau of moat rhyolite located west of the resurgent dome or beneath Mammoth Mountain. Lateral flow of thermal water away from such upflow zones through reservoirs in the Bishop Tuff and early rhyolite accounts for temperature reversals encountered in most existing wells. Dating of hot-spring deposits from active and inactive thermal areas confirms previous interpretations of the evolution of hydrothermal activity that suggest two periods of extensive hot-spring discharge, one peaking about 300 ka and another extending from about 40 ka to the present. The onset of hydrothermal activity around 40 ka coincides with the initiation of rhyolitic volcanism along the Mono-Inyo Craters volcanic chain that extends beneath the caldera's west moat.     

Evaluating fracture patterns within a resurgent caldera: Campi Flegrei, Italy

Understanding deformation of active calderas allows their dynamics to be defined and their hazard mitigated. The Campi Flegrei resurgent caldera (Italy) is one of the most active and hazardous volcanoes in the world, characterized by post-collapse resurgence, eruptions, ground deformation, and seismicity. An original structural analysis provides an overview of the main fracture zones. NW-SE and NE-SW fractures (normal or transtensive faults and extensional fractures) predominate along the rim and within the caldera, suggesting a regional control, both during and after the collapses. While the NE-SW fractures are ubiquitous in the deposits of the last ∼37 ka, NW-SE fractures predominate in the last 4.5 ka, during resurgence. The most recently (<4.5 ka) strained area lies in the caldera center (Solfatara area), where the faults, with an overall ∼ENE-WSW extension direction, appear to be associated with the bending due to resurgence. Solfatara lies immediately to the east of the most uplifted part of the caldera (Pozzuoli area), where domes form and culminate both on the long-term (resurgence, accompanied by volcanic activity) and short-term deformation (1982–1984 bradyseism, accompanied by seismic and hydrothermal activity). Similar volcano-tectonic behavior characterizes the short- and long-term uplifts, and only the intensity of the tectonic and volcanic activity varies, being related to varying amounts of uplift. Seismicity and hydrothermal manifestations occur during the bradyseisms, with moderate uplift, while surface faulting and eruptions occur during resurgence, with higher uplift. The features observed at Campi Flegrei are found at other major calderas, suggesting consistent behavior of large magmatic systems.     

Minerals produced during cooling and hydrothermal alteration of ash flow tuff from Yellowstone drill hole Y-5

A rhyolitic ash-flow tuff in a hydrothermally active area within the Yellowstone caldera was drilled in 1967, and cores were studied to determine the nature and distribution of primary and secondary mineral phases. The rocks have undergone a complex history of crystallization and hydrothermal alteration since their emplacement 600,000 years ago. During cooling from magmatic temperatures, the glassy groundmass underwent either devitrification to alkali feldspar + α-cristobalite ± tridymite or granophyric crystallization to alkali feldspar + quartz. Associated with the zones of granophyric crystallization are prismatic quartz crystals in cavities similar to those termed miarolitic in plutonic rocks. Vapor-phase alkali feldspar, tridymite, magnetite, and sporadic α-cristobalite were deposited in cavities and in void spaces of pumice fragments. Subsequently, some of the vapor-phase alkali feldspar crystals were replaced by microcrystalline quartz, and the vapor-phase minerals were frosted by a coating of saccharoidal quartz.Hydrothermal minerals occur primarily as linings and fillings of cavities and fractures and as altered mafic phenocrysts. Chalcedony is the dominant mineral related to the present hydrothermal regime and occurs as microcrystalline material mixed with various amounts of hematite and goethite. The chalcedony displays intricate layering and was apparently deposited as opal from silica-rich water. Hematite and goethite also replace both mafic phenocrysts and vapor-phase magnetite. Other conspicuous hydrothermal minerals include montmorillonite, pyrite, mordenite, calcite, and fluorite. Clinoptilolite, erionite, illite, kaolinite, and manganese oxides are sporadic. The hydrothermal minerals show little correlation with temperature, but bladed calcite is restricted to a zone of boiling in the tuff and clearly was deposited when CO₂ was lost during boiling.Fractures and breccias filled with chalcedony are common throughout Y-5 and may have been produced by rapid disruption of rock caused by sudden decrease of fluid pressure in fractures, most likely a result of fracturing during resurgent doming in this part of the Yellowstone caldera. The chalcedony probably was deposited as opal or β-cristobalite from a pre-existing silica floc that moved rapidly into the fractures and breccias immediately after the sudden pressure drop.     

Results of Crust and Mantle Soundings in Central and Northern Europe in the 21st Century (Review)

The first decade of 21st century is characterized by the appearance of new approaches to deep induction soundings. The theory of magnetovariation and magnetotelluric soundings was generalised or corrected. Spatial derivatives of response functions (induction arrows) were obtained for the ultra-long periods. New phenomena have been detected by this method: secular variations of the Earth’s apparent resistivity and the rapid changes of induction arrows over the last 50 years. The first one can be correlated with the number of earthquakes, and the second one–with geomagnetic jerks in Central Europe. The extensive studies of geoelectrical structure of the crust and mantle were realized in the frame of a series of international projects. New information about geoelectrical structures of the crust in Northern Europe and Ukraine was obtained by deep electromagnetic soundings involving controlled powerful sources. An influence of the crust magnetic permeability on the deep sounding results was confirmed.     

Mobility and depositional controls of radioelements in hydrothermal systems at the Long Valley and Valles calderas

The loci and abundance of U and Th were examined in tuffaceous rocks encompassing hydrothermal systems at the Long Valley caldera, California and the Valles caldera, New Mexico. Aspects of these systems may be analogous to conditions expected in a potential site for a high-level waste repository in welded tuff. Examination of radioelements in core from scientific drill holes at these sites was accomplished by gamma-ray spectrometry and fission-track radiography. In the lateral-flowing hydrothermal system at the Long Valley caldera, where temperatures range from 140 to 200 °C, U is concentrated to 20 ppm in Fe-rich zones of varved tuff and to 50 ppm with Fe-rich mineral phases in tuff fragments of a calcite-cemented breccia. U-series disequilibrium in some of these samples suggests mobilization/deposition of parent U and/or its daughters. In the vapor zone of the Valles caldera's hydrothermal system (temperature ˜ 100 °C), the concordance of high U, low Th/U and decreasing whole-rock O-isotope ratios suggests that U was concentrated in response to hydrothermal circulation when the system was formerly liquid-dominated. In the underlying present-day liquid-dominated zone (temperature to 210 °C), U, up to several tens of parts per million, occurs with pyrite and Fe-oxide minerals, and in concentrations to several percents with a Ti-Nb-Y-rare earth mineral. In the Valles system's outflow zone, U is also concentrated in Fe-rich zones as well as in carbonaceous-rich zones in the Paleozoic sedimentary rocks that underlie the Quaternary tuff. Th, associated with accessory minerals, predominates in breccia zones and in a mineralized fault zone near the base of the Paleozoic sedimentary sequence. Relatively high concentrations of U occur in springs representative of water recharging the Valles caldera's hydrothermal system. In contrast, considerably lower U concentrations occur in hot waters (> 220 °C) and in the system's outflow plume, suggesting that U is concentrating in the hotter part of the system. The Long Valley and Valles observations indicate that U and Ra are locally mobile under hydrothermal conditions, and that reducing conditions associated with Fe-rich minerals and carbonaceous material are important factors in the adsorption of U, and thus can retard its transport in water at elevated temperature.     

The subsurface structure of the AIRA caldera and its vicinity in Southern Kyushu, Japan

Southern Kyushu, Japan, includes a chain of large and small calderas and active volcanoes, and the greatest part of it is covered with thick pyroclastic ejecta. The regional and local structures of this area are discussed from the standpoint of physical volcanology, with consideration of all available data.The regional structure of this area is examined in the light of gravity and geomagnetic anomalies. Two layers of the earth's uppermost crust are defined by spectrum analysis of the gravity anomalies. These two layers are identical with the two identified by seismicwave velocities. The Bouguer gravity anomalies are relatively high and rather monotonous over outcrops of the Mesozoic basement and the granite, but are relatively low and perturbed over calderas and caldera-like structures. Two low-gravity anomalies in Kagoshima Bay are remarkable. One is circular, with its center on the Aira caldera. The other is elongated between the Satsuma and Oosumi peninsulas. The southern end of the latter anomaly is occupied by the Ata caldera. Discussion of the gravity anomalies of the Aira caldera suggests that the subsurface basement has a funnel shape and is overlain by ‘fallback’. The sub bottom geology of the caldera suggests that it is formed by a few smaller depressions, though the distribution of the overall gravity anomalies is parallel with its shape.The southern part of Kagoshima Bay is characterized by a graben-like topography and low-gravity anomalies and, moreover, by several calderas. The middle part, between the Aira and Ata calderas, may have a graben-like structure. A profile crossing the bay through Sakurajima volcano is modeled on the basis of results from drilling and gravity surveys. The basement has a graben-like structure and is filled with coarse and low-density deposits, and the structure continues northwards to the Aira caldera with a funnel shape.A comparison of this area with the Taupo-Rotorua depression in New Zealand and Lake Toba in Indonesia, leads the authors to the conclusion that such major volcanic depressions may have been formed by amalgamation of a series of caldera-like structures which were formed by multiple violent explosions accompanied by ejection of a tremendous amount of pyroclastic material.