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).     

Radar determination of the spatial structure of hydraulic conductivity

Spatial variability of hydraulic conductivity exerts a predominant control on the flow of fluid through porous media. Heterogeneities influence advective pathways, hydrodynamic dispersion, and density-dependent dispersion; they are, therefore, a key concern for studies of ground water resource development, contaminant transport, and reservoir engineering. Ground-penetrating radar contributes to the remote, geophysical characterization of the macroscale variability of natural porous media. On a controlled excavation of a glacial-fluvial sand and gravel deposit in the Fanshawe Delta area (Ontario, Canada), the hydraulic conductivity field of a 45 x 3 m vertical exposure was characterized using constant-head permeameter measurements performed on undisturbed horizontal sediment cores. Ground-penetrating radar data were collected along the excavation face in the form of both reflection and common midpoint surveys. Comparison of geostatistical analyses of the permeameter measurements and the radar data suggests thatthe horizontal correlation structure of radar stack velocity can be used to directly infer the horizontal correlation structure of hydraulic conductivity. The averaging nature of the common midpoint survey is manifest in the vertical correlation structure of stack velocity, making it less useful. Radar reflection data do not exhibit a spatial structure similar to that of hydraulic conductivity possibly because reflections are a result of material property contrasts rather than the material properties themselves.     

Geomagnetic evolution: 400 years of change on planet Earth

Spherical harmonic coefficients of the geomagnetic field, calculated from historical observations of declination, inclination and intensity, and from archaeomagnetic inclination results, have been used to produce a film of geomagnetic change since 1600 A.D. The non-dipole geomagnetic field is found to be constantly changing: no fixed or standing non-dipole features are observed. Non-dipole foci are seen to have lifetimes of a few hundred years. The westward drift, which was an important feature of the 18th and early 19th century geomagnetic field, was less pronounced in the 17th century. The growth, evolution, decay and replacement of non-dipole foci, but not their movement are found to have been the major features producing century-long secular directional magnetic variation. Most of the low degree and order spherical harmonic coefficients have changed significantly over the last few hundred years. In particular the change in sign of the axisymmetric quadrupole around 1837 A.D. is noted. Sustained, century-long, intensity changes, however, appear to have been dominated by variations in the intensity of the centred dipole, rather than by non-dipole field fluctuations.     

Geochemistry of waters in the Valley of Ten Thousand Smokes region, Alaska

Meteoric waters from cold springs and streams outside of the 1912 eruptive deposits filling the Valley of Ten Thousand Smokes (VTTS) and in the upper parts of the two major rivers draining the 1912 deposits have similar chemical trends. Thermal springs issue in the mid-valley area along a 300-m lateral section of ash-flow tuff, and range in temperature from 21 to 29.8°C in early summer and from 15 to 17°C in mid-summer. Concentrations of major and minor chemical constituents in the thermal waters are nearly identical regardless of temperature. Waters in the downvalley parts of the rivers draining the 1912 deposits are mainly mixtures of cold meteoric waters and thermal waters of which the mid-valley thermal spring waters are representative. The weathering reactions of cold waters with the 1912 deposits appear to have stabilized and add only subordinate amounts of chemical constituents to the rivers relative to those contributed by the thermal waters. Isotopic data indicate that the mid-valley thermal spring waters are meteoric, but data is inconclusive regarding the heat source. The thermal waters could be either from a shallow part of a hydrothermal system beneath the 1912 vent region or from an incompletely cooled, welded tuff lens deep in the 1912 ash-flow sheet of the upper River Lethe area.Bicarbonate-sulfate waters resulting from interaction of near-surface waters and the cooling 1953–1968 southwest Trident plug issue from thermal springs south of Katmai Pass and near Mageik Creek, although the Mageik Creek spring waters are from a well-established, more deeply circulating hydrothermal system. Katmai caldera lake waters are a result of acid gases from vigorous drowned fumaroles dissolving in lake waters composed of snowmelt and precipitation.     

An evaluation of the health status of the Lavaca Bay, Texas ecosystem using Crassostrea virginica as the sentinel species

Locational risks for compromised ecosystem health for the eastern oysters (Crassostrea virginica) harvested from Lavaca Bay, Texas were estimated. Flow cytometric evaluation of variations in DNA content and the lysosomal destabilization assay were used for evaluation of genotoxicity and stress, respectively. Bayesian geo-statistical methods were utilized to estimate and evaluate spatial effects. For models with spatial risks, continuous surface maps of predicted parameter values were created to evaluate risk location. Lysosomal destabilization assay results were spatially oriented whereas flow cytometry results were fit best with the random effects model. While not spatially oriented, the highest levels of variations in DNA content were also present near industrial facilities. Locational risks of increased biomarkers of genotoxicity and stress in the eastern oyster (Crassostrea virginica) were increased with proximity to industrial facilities     

Evaluating an Analytical Model to Predict Subsurface LNAPL Distributions and Transmissivity from Current and Historic Fluid Levels in Groundwater Wells: Comparing Results to Numerical Simulations

A recent analytical model predicts free, entrapped, and residual LNAPL saturations and the LNAPL transmissivity in the subsurface from current and historic fluid levels in groundwater wells. As such, the model accounts for effects of fluid level fluctuations in a well. The model was developed to predict LNAPL specific volumes and transmissivities from current fluid level measurements in wells and either recorded historic fluid level fluctuations in wells or estimates. An assumption is made in the model that the predictions are not dependent on whether the historic highest or lowest fluid level elevations in a well occur first. To test the assumption, we conduct two simulations with a modified multiphase flow numerical code TMVOC that incorporates relative permeability‐saturation‐capillary head relations employed in the model. In one simulation, the initial condition is for fluid levels in a well at the historic highest elevations. In the other simulation, the initial condition is for fluid levels in a well at the historic lowest elevations. We change the boundary conditions so both historical conditions occur followed by generating the current condition. Results from the numerical simulations are compared to model predictions and show the assumption in the analytical model is reasonable. The analytical model can be used to develop/refine conceptual site models and for assessing potential LNAPL recovery endpoints, especially on sites with fluctuating fluid levels in wells.     

A reinterpretation of the entrainment process in some laboratory experiments

Pollard et al. (1973) proposed a closure hypothesis for the entrainment process in the stress-mixed layer of a stratified fluid. Some laboratory experiments suitable for testing the hypothesis were carried out by Kantha et al. (1977) and by Kato and Phillips (1969). In the early period of each experiment, before side-wall friction becomes important, these experiments support the closure hypothesis, with a critical Froude number close to 1, perhaps F = 1.2 ± 0.2.