Petrography and uplift history of the Quaternary Takidani Granodiorite: could it have hosted a supercritical (HDR) geothermal reservoir?

The Quaternary Takidani Granodiorite (Japan Alps) is analogous to the type of deep-seated (3–5 km deep) intrusive-hosted fracture network system that might support (supercritical) hot dry/wet rock (HDR/HWR) energy extraction. The I-type Takidani Granodiorite comprises: porphyritic granodiorite, porphyritic granite, biotite-hornblende granodiorite, hornblende-biotite granodiorite, biotite-hornblende granite and biotite granite facies; the intrusion has a reverse chemical zonation, characterized by >70 wt% SiO₂ at its inferred margin and <67 wt% SiO₂ at the core. Fluid inclusion evidence indicates that fractured Takidani Granodiorite at one time hosted a liquid-dominated, convective hydrothermal system, with <380°C, low-salinity reservoir fluids at hydrostatic (mesothermal) pressure conditions. ‘Healed’ microfractures also trapped >600°C, hypersaline (35 wt% NaCl_eq) fluids of magmatic origin, with inferred minimum pressures of formation being 600–750 bar, which corresponds to fluid entrapment at 2.4–3.0 km depth. Al-in-hornblende geobarometry indicates that hornblende crystallization occurred at about 1.45 Ma (7.7–9.4 km depth) in the (marginal) eastern Takidani Granodiorite, but later (at 1.25 Ma) and shallower (6.5–7.0 km) near the core of the intrusion. The average rate of uplift across the Takidani Granodiorite from the time of hornblende crystallization has been 5.1–5.9 mm/yr (although uplift was about 7.5 mm/yr prior to 1.2 Ma), which is faster than average uplift rates in the Japan Alps (3 mm/yr during the last 2 million years). A temperature–depth–time window, when the Takidani Granodiorite had potential to host an HDR system, would have been when the internal temperature of the intrusive was cooling from 500°C to 400°C. Taking into account the initial (7.5 mm/yr) rate of uplift and effects of erosion, an optimal temperature–time–depth window is proposed: for 500°C at 1.54–1.57 Ma and 5.2±0.9 km (drilling) depth; and 400°C at 1.36–1.38 Ma and 3.3±0.8 km (drilling) depth, which is within the capabilities of modern drilling technologies, and similar to measured temperature–depth profiles in other active hydrothermal systems (e.g. at Kakkonda, Japan).     

Hydrothermal Quartz Vein Formation, Revealed by Coupled SEM-CL Imaging and Fluid Inclusion Microthermometry: Shuteen Complex, South Gobi, Mongolia

Abstract. Scanning electron microscopy-cathodoluminescence (SEM-CL) imaging of vein quartz in the Cu-mineralised, Shuteen Complex (South Gobi, Mongolia) has revealed a complex history of crystal growth, dissolution and microfracture healing, associated with several hydrothermal events that could not be detected using other observational techniques (e.g. transmitted/reflected light microscopy, back-scattered electron imaging, or secondary electron imaging).
The quartz initially grew as CL-bright/grey crystals in a 345±30C liquid reservoir, as inferred by the analysis of primary liquid fluid inclusions (average Th of 343C; 6.6∼7.7 wt% NaCl_eq). Quartz precipitation occurred at the edge of the crystals as reservoir fluids cooled to 260±25C, as indicated by micron-scale CL-dark/CL-bright quartz growth bands containing abundant fluid inclusions (with an average T_h values of 261C). Pressure fluctuations were the likely cause of dissolution, as SEM-CL imaging reveals the quartz have corroded or rounded crystal edges, and precipitation of later quartz into open space. SEM-CL imaging shows the quartz contains healed microfractures that trapped low salinity fluids (3.9 wt% NaC1_eq) with Th values of 173±15C.
SEM-CL imaging provides a means of deciphering the thermal and chemical evolution of the fossil Shuteen hydrothermal system, and the nature of hydrothermal quartz vein-forming processes, by facilitating the correlation of distinct fluid inclusion populations and their relative chronology, with specific hydrothermal events.     

Mineralogical and textural variation of silica minerals in hydrothermal flow-through experiments: Implications for quartz vein formation

We conducted hydrothermal flow-through experiments at 430 °C and 31 MPa to investigate the mechanism of silica precipitation on granite under crustal conditions. Two experiments were performed using different input solutions: a single-component Si solution, and a multi-component solution with minor Al, Na, and K. The degree of supersaturation with respect to quartz, Ω = C_Si/C_Si,Qtz,eq, where C_Si and C_Si,Qtz,eq indicate Si concentration in solutions and the solubility of quartz within water, respectively, decreased from 3-3.5 to <1.1 along the flow path. A variety of silica minerals formed during the experiments (opal-A, opal-C, chalcedony, and quartz), and their occurrences and modal abundances changed in response to Ω and the presence of additives in the solution.For near-equilibrium solutions (Ω < ∼1.2), silica precipitation occurred in a simple way in both experiments, being restricted to overgrowths on pre-existing quartz surfaces in the granite. At higher saturation levels (Ω > ∼1.2), silica minerals were deposited on other surfaces in addition to quartz. In the single-component experiment, the dominant silica minerals changed in the order of opal-A → opal-C → quartz with decreasing Si concentration along the flow path. In contrast, in the multi-component experiment, quartz and minor chalcedony formed throughout the entire reaction vessel. This finding indicates that impurities (Na, K, and Al) in the solutions inhibited the precipitation of opal and enhanced the direct nucleation of quartz. The systematic appearance of metastable silica minerals were examined by nucleation processes and macroscopic precipitation kinetics. Our experimental results indicate that different precipitation mechanisms yield contrasting textures, which in turn suggests that vein textures can be used as indicators of solution chemistry within the fracture.     

The significance of silica precipitation on the formation of the permeable–impermeable boundary within Earth's crust

The hydrological system within Earth's crust is divided into the permeable zone and the underlying, much less permeable zone. In this study, we investigated the solubility and precipitation kinetics of silica in water under the conditions of the Kakkonda geothermal field, Japan, where well WD‐1a penetrates the boundary between the hydrothermal convection zone and the heat conduction zone. We found that quartz solubility has a local minimum at the depth of the hydrological boundary, at ~3100‐m depth (380 °C and 24 MPa), in the case of either hydrostatic conditions or fluid pressure increase above hydrostatic at deeper levels. The hydrothermal experiments reveal that rapid quartz precipitation occurs via nucleation when fluids are brought from the liquid region to the supercritical region. The preferential precipitation of quartz at a specific depth plays a significant role in forming and sustaining the permeable–impermeable boundary in the crust.     

Estimation and exclusion of multipath range error for robust positioning

In integrated systems for accurate positioning, which consist of GNSS, INS, and other sensors, the GNSS positioning accuracy has a decisive influence on the performance of the entire system and thus is very important. However, GNSS usually exhibits poor positioning results in urban canyon environments due to pseudorange measurement errors caused by multipath creation, which leads to performance degradation of the entire positioning system. For this reason, in order to maintain the accuracy of an integrated positioning system, it is necessary to determine when the GNSS positioning is accurate and which satellites can have their pseudorange measured accurately without multipath errors. Thus, the objective of our work is to detect the multipath errors in the satellite signals and exclude these signals to improve the positioning accuracy of GNSS, especially in an urban canyon environment. One of the previous technologies for tackling this problem is RAIM, which checks the residual of the least square and identifies the suspicious satellites. However, it presumes a Gaussian measurement error that is more common in an open-sky environment than in the urban canyon environment. On the other hand, our proposed method can estimate the size of the pseudorange error directly from the information of altitude positioning error, which is available with an altitude map. This method can estimate even the size of non-Gaussian error due to multipath in the urban canyon environment. Then, the estimated pseudorange error is utilized to weight satellite signals and improve the positioning accuracy. The proposed method was tested with a low-cost GNSS receiver mounted on a test vehicle in a test drive in Nagoya, Japan, which is a typical urban canyon environment. The experimental result shows that the estimated pseudorange error is accurate enough to exclude erroneous satellites and improve the GNSS positioning accuracy.     

Mineral distribution within polymineralic veins in the Sanbagawa belt,Japan: implications for mass transfer during vein formation

Pelitic schists of the Sanbagawa metamorphic belt contain several types of polymineralic veins that formed during the late stages of exhumation. The vein mineral assemblages are quartz + albite + K-feldspar + chlorite ± calcite (Type I, II) and quartz + albite + calcite (Type III). Type I and II veins contain quartz and albite with stretched-crystal and elongate-blocky textures, respectively. The mineral species within Type I veins vary with compositional bands within the host rocks. Type III veins are characterized by euhedral to subhedral quartz grains with concentric zoning and a homogeneous distribution along the vein length. The vein textures vary depending on the crack aperture during multiple crack-seal events: <0.08 mm for Type I, and 0.5–10 mm for Type III. Type II veins show intermediate features between Type I and III veins in terms of mineral distribution (weak dependence on the host rock composition) and apparent crack aperture (less than 1–15 mm). These observations suggest a transition in the dominant transport mechanism of vein components with increasing crack aperture, from diffusion from host rocks to fluid advection along cracks.     

Factors Controlling the Fractionation and Seasonal Mobility Variations of Ga and In in Systems Impacted by Acidic Thermal Waters: Effects of Thermodynamics and Bacterial Activity

This work assessed both the fractionation and the seasonal mobility variations of Ga and In in systems impacted by acidic thermal waters. This was accomplished by performing thermodynamic calculations using the PHREEQC algorithm and by assessing the activity of acidophilic iron-oxidizing bacteria. The pH of the Kusatsu thermal waters in Gunma Prefecture, central Japan, is rapidly increased following the addition of a lime suspension. After an abrupt pH increase, under which conditions free ions of Ga and In and their complexes with Cl^? and SO₄^2? exist only in negligible quantities, the majority of dissolved Ga and In is removed by sorption onto suspended hydrous ferric oxides (HFOs). These HFOs are then transported to an artificial lake without significant sedimentation along the river. Subsequently, the suspended HFOs settle out and are added to sediments without significant fractionation between Ga and In. The Tamagawa thermal waters in Akita Prefecture, northeast Japan, are also treated with lime. However, complete neutralization requires mixing with some tributary streams, leading to a gradual downstream increase in pH. Dissolved Ga is, in general, sorbed by HFOs in upstream areas, leading to wide dispersal of Ga across the entire watershed. In comparison, In is transported to the lake inlet predominantly as a Cl^? complex species without significant removal along the river, with the majority being precipitated in an artificial lake, where Cl^? concentrations are too low to form stable complex species with In, and thus, dissolved In is sorbed by HFOs. As a result, In is effectively concentrated within downstream lakebed sediments, whereas Ga is dispersed along the river. Seasonal variations in Ga mobility within the Tamagawa field between snowmelt and low-flow seasons are primarily controlled by pH, because hydrolysis reactions of these metals, which are related to sorption reactions, tend to occur in the upstream regions in the snowmelt season. However, under warmer conditions, HFO formation preferably occurs due to the activity of acidophilic iron-oxidizing bacteria. Thus, under similar pH variations, dissolved Ga is more effectively removed by HFOs during warmer seasons. On the contrary, because HFOs are abundantly formed in low-flow season, even under colder conditions, before In hydrolysis reaction starts to occur, In mobility is less affected by water temperature and then bacterial activity.     

Style of fluid flow and deformation in and around an ancient out-of-sequence thrust: An example from the Nobeoka Tectonic Line in the Shimanto accretionary complex, Southwest Japan

To understand the characteristics of deformation of an out-of-sequence thrust (OST) and the style of fluid flow along it, we investigated the Nobeoka Tectonic Line, which has been interpreted as a deep OST (7–9 km), in the Shimanto accretionary complex, Southwest Japan. The shear zone in the footwall differs significantly in the along-strike direction not only in thickness, which varied from 100 to 300 m, but also in lithology and mineral vein development. These variations might reflect primarily differences in lithology; that is, the sandstone-dominant shear zone with a large amount of mineral veins precipitated in microcracks is relatively thick, whereas the shale-dominant shear zone with a small amount of veins and with textures indicating highly pressurized pore fluid, is thinner. By comparison with characteristics of a shallow OST (3–5 km), we conclude that the shallow OST has experienced repeated brittle failure with rapid slip and focused fluid flow whereas the deep OST has experienced both brittle and ductile deformation, followed by fluid flow of various styles, depending on the lithology.     

X-ray CT based numerical analysis of fracture flow for core samples under various confining pressures

The X-ray CT based numerical analysis of fracture flow for core samples, recently developed by the authors, was applied to two granite core samples having either a mated artificial or a mated natural fracture at confining pressures of 5 to 50 MPa. A third-generation medical X-ray CT scanner was used to image the samples within a core holder consisting of an aluminum liner and a carbon fiber overwrap. Fracture models (i.e., aperture distributions) were obtained by the CT images, the resolution of which was coarser than the apertures, and a single-phase flow simulation was performed using a local cubic law-based fracture flow model. Numerical results were evaluated by a fracture porosity measurement and a solution displacement experiment using NaCl and NaI aqueous solutions. These numerical results coincided only qualitatively with the experimental results, primarily due to image noise from the aluminum liner of the core holder. Nevertheless, the numerical results revealed flow paths within the fractures and their changes with confining pressure, whereas the experimental results did not provide such results. Different stress-dependencies in the flow paths were observed between the two samples despite the similar stress-dependency in fracture porosity and permeability. The changes in total area of the flow paths with confining pressure coincided qualitatively with changes in breakthrough points in the solution displacement experiment. Although the data is limited, the results of the present study suggest the importance of analyzing fluid flows within naturally fractured core samples under in situ conditions in order to better understand the fracture flow characteristics in a specific field. As demonstrated herein, X-ray CT-based numerical analysis is effective for addressing this concern. Using a multi-phase flow model, as well as a core holder constructed of an engineered plastic, should provide a useful, non-destructive, and non-contaminative X-ray CT-based fracture flow analysis for core samples under in situ conditions in future studies.