Thermal structure and paleo-heat flow in the Shimanto accretionary prism, Southwest Japan

Abstract Thermal structural analysis and paleo-heat flow estimation provide clues to understanding the thermal evolution of the accretionary complex. The thermal structure and heat flow in the Jurassic Chichibu and Cretaceous to Tertiary Shimanto accretionary complex, Southwest Japan, have been investigated by vitrinite reflectance measurement and fluid inclusion analysis. As a result, the local and multistage metamorphisms were recognized as follows. First, the Tertiary complex around the Miocene Ashizuri granite underwent exposure to extra-high temperatures. Second, the Okitsu Melange underwent exposure to higher temperatures than the surrounding strata and was formed concurrently with the Kula-Pacific ridge subduction beneath the Japanese Islands in the Eocene. Finally, the thermal structure of most of the Cretaceous and southern Jurassic complexes is independent of the geologic structure, indicating that these areas suffered thermal overprint. Regional radiometric dating studies show that most of the Cretaceous Shimanto complex was heated in the Eocene; the thermal overprint might have occurred as a result of ridge subduction. The heat flow during peak heating was estimated to be 95–120 mW/m² except for the Cretaceous Okitsu melange and the Cretaceous Nonokawa formation, north of the Okitsu Melange; a much higher value of heat flow of ~200 mW/m² was estimated in the Okitsu Melange. An estimation of heat flow failed for the non-okawa formation because thermal equilibrium between the fluid and rocks has not yet been reached. It is probable that the southern strata underwent a higher heat flow. Such a trenchward increase in heat flow resembles the present situation of the Nankai Trough, although the heat flow in the Eocene was much higher.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

复合填料热传导过程数值模拟

复合填料由废铸砂、粉煤灰、聚苯乙烯（EPS）颗粒、水泥和水按一定质量比例混合制成,具有低导热性、抗冻胀和轻质特性。根据热阻力法则和比等效导热系数相等法则,将材料视为由大量正方形单元体组成,其中心为一个球形EPS颗粒,这种单元体与总体的导热系数相等。运用AN SY S软件对复合填料单元体的传热过程进行了有限元模拟,根据瞬态法导热系数测试原理,推求不同EPS掺入比情况下的复合填料导热系数,分析EPS颗粒与导热系数的关系,发现材料导热系数随着EPS掺入比的增加而快速降低。将材料导热系数数值模拟与实测结果作比较,发现有限元模拟值与实测值接近,说明采用有限元方法可以实现对材料传热过程的模拟。     

Thermal calculation for high-pressure contact metamorphism: application to eclogite formation in the Sebadani area, the Sambagawa belt, SW Japan

Eclogite, a high-pressure–temperature metamorphic rock characterized by garnet + omphacite, is usually considered to be a product of regional metamorphism under a low geothermal gradient. However, in the Sebadani area of the Sambagawa metamorphic belt most petrologists agree that the eclogite formed by localized contact metamorphism due to intrusion of a body in the solid-state (the Sebadani mass). This process is termed ‘high-pressure contact metamorphism'. However, geological considerations suggest that the effect of such a process would be limited, firstly because the speed of emplacement for solid-state material will generally be much lower than that for magma and secondly because in the solid-state there is no heat of fusion in the body available for thermal effects. Thermal modelling of a solid-state intrusion, based on the heat conduction equation, allows the relationship between size of intrusion, velocity of emplacement and thermal effects to be calculated. Two cases have been considered: (1) a hot model, where none of the heat conducted into the surroundings is lost during the rise of the body; and (2) a cold model where all the heat conducted into the surroundings is lost. These models bracket possible thermal histories of the body. Calculations suggest that in the Sebadani region, production of the observed metamorphic features requires unrealistically high velocity and a much larger intruded body than is observed. These conclusions suggest that it is unlikely that eclogite in the Sebadani area was formed by high-pressure contact metamorphism, but rather that it represents the highest-grade part of the regional Sambagawa metamorphism.     

Surficial fractures in the Navajo sandstone,south-western USA: the roles of thermal cycles,rainstorms, granular disintegration,and iterative cracking

Deep (> 5 m) sheeting fractures in the Navajo sandstone are evident at numerous sites in southern Utah and derive from tectonic stresses. Strong diurnal thermal cycles are, however, the likely triggers for shallow (< 0.3 m) sheeting fractures. Data from subsurface thermal sensors reveal that large temperature differences between sensors at 2 and 15 cm depth on clear summer afternoons are as great as those that trigger sheeting fractures in exposed California granite. Extensive polygonal patterns in the Navajo sandstone are composed of surface-perpendicular fractures and were produced by contractile stresses. Numerous studies have shown that porewater diminishes the tensile strength of sandstone. Based on our thermal records, we propose that cooling during monsoonal rainstorms triggers polygonal fracturing of temporarily weakened rock. On steep outcrops, polygonal patterns are rectilinear and orthogonal, with T-vertices. Lower-angle slopes host hexagonal patterns (defined by the dominance of Y-vertices). Intermediate patterns with rectangles and hexagons of similar scale are common. We posit that outcropping fractures are advancing downward by iterative steps, and that hexagons on sandstone surfaces (like prismatic columns of basalt) have evolved from ancestral orthogonal polygons of similar scale. In lava flows, fractures elongate intermittently as they follow a steep thermal gradient (the source of stress) as it rapidly moves through the rock mass. In our model, a steep, surficial thermal gradient descends through unfractured sandstone, but at the slow pace of granular disintegration. Through time, as the friable rock on stable slopes erodes, iterative cracking advances into new space. Hexagonal patterns form as new fractures, imperfectly guided by the older ones, propagate in new directions, and vertices drift into a configuration that minimizes the ratio of fracture length to polygon area. © 2020 John Wiley & Sons, Ltd.     

Climate change and resilience of tributary thermal refugia for salmonids in eastern Canadian rivers

Abstract

River water temperature regimes are expected to change along with climate over the next decades. This work focuses on three important salmon rivers of eastern Canada, two of which warm up most summers to temperatures higher than the Atlantic salmon lethal limit (>28°C). Water temperature was monitored at 53 sites on the three basins during 2–18 summers, with about half of these sites either known or potential thermal refugia for salmon. Site-specific statistical models predicting water temperature, based on 10 different climate scenarios, were developed in order to assess how many of these sites will remain cool enough to serve as refugia in the future (2046–2065). The results indicate that, while 19 of the 23 identified refugia will persist, important increases in the occurrence and duration of temperature events in excess of 24°C and 28°C, respectively, in the mainstems of the rivers, will lead to higher demands for thermal refugia in the salmonid populations.

Editor Z.W. Kundzewicz; Associate editor T. Okruszko     

Macrodynamics of α2 dynamos

Abstract

Two distributions of the α-effect in a sphere are considered. The inviscid limit is approached both by direct numerical solution and by solution of a simpler nonlinear eigenvalue problem deriving from asymptotic boundary layer analysis for the case of stress-free boundaries. The inviscid limit in both cases is dominated by the need to satisfy the Taylor constraint which states that the integral of the Lorentz force over cylindrical (geostrophic) contours in a homogeneous fluid must tend to zero. For a small supercritical range in α, this condition can only be met by magnetic fields which vanish as the viscosity goes to zero. In this range, the agreement of the two approaches is excellent. In a portion of this range, the method of finite amplitude perturbation expansion is useful, and serves as a guide for understanding the numerical results. For larger α, evidence from the nonlinear eigenvalue problem suggests both that the Taylor state exists, and that the transition from small to large amplitude can require a finite amplitude (oscillatory) instability in accord with the findings of Soward and Jones (1983). However, solutions of the full equations have not been found which are independent of viscosity at larger values of α.     

Thermal evolution of the Tertiary Shimanto Belt, Muroto Peninsula, Shikoku, Japan

Abstract The Shimanto accretionary complex on the Muroto Peninsula of Shikoku comprises two major units of Tertiary strata: the Murotohanto Sub-belt (Eocene-Oligocene) and the Nabae Sub-belt (Oligocene-Miocene). Both sub-belts have been affected by thermal overprints following the peak of accretion-related deformation. Palaeotemperatures for the entire Tertiary section range from ～ 140 to 315°C, based upon mean vitrinite reflectance values of 0.9–5.0%R_m. Values of illite crystallinity index are consistent with conditions of advanced diagenesis and anchimetamorphism. Illite/mica b₀ lattice dimensions indicate that burial pressures were probably no greater than 2.5kbar. In general, levels of thermal maturity are higher for the Murotohanto Sub-belt than for the Nabae Sub-belt. The Eocene-Oligocene strata also display a spatial decrease in thermal maturity from south to north and this pattern probably was caused by regional-scale differential uplift following peak heating. Conversely, the palaeothermal structure within the Nabae Sub-belt is fairly uniform, except for the local effects of mafic intrusions at the tip of Cape Muroto. There is a paleotemperature difference of ～ 90°C across the boundary between the Murotohanto and Nabae Sub-belts (Shiina-Narashi fault), and this contrast is consistent with approximately 1200 m of post-metamorphic vertical offset. Subduction prior to Middle Miocene probably involved the Kula or fused Kula-Pacific plate and the background geothermal gradient during the Eocene-Oligocene phase of accretion was ～ 30–35°C/km. Rapid heating of the Shimanto Belt evidently occurred immediately after a Middle Miocene reorganization of the subduction boundary. Hot oceanic lithosphere from the Shikoku Basin first entered the subduction zone at ～ 15 Ma; this event also coincided with the opening of the Sea of Japan and the rapid clockwise rotation of southwest Japan. The background geothermal gradient at that time was ～ 70°C/km. Whether or not all portions of the inherited (Eocene-Oligocene) palaeothermal structure were overprinted during the Middle Miocene remains controversial.     

Hydrogeochemistry of the Simav geothermal field, western Anatolia, Turkey

Thermal waters hosted by Menderes metamorphic rocks emerge along fault lineaments in the Simav geothermal area. Thermal springs and drilled wells are located in the Eynal, Çitgöl and Na a locations, which are part of the Simav geothermal field. Studies were carried out to obtain the main chemical and physical characteristics of thermal waters. These waters are used for heating of residences and greenhouses and for balneological purposes. Bottom temperatures of the drilled wells reach 163°C with total dissolved solids around 2225 mg/kg. Surface temperatures of thermal springs vary between 51°C and 90°C. All the thermal waters belong to Na–HCO₃–SO₄ facies. The cold groundwaters are Ca–Mg–HCO₃ type. Dissolution of host rock and ion-exchange reactions in the reservoir of the geothermal system shift the Ca–Mg–HCO₃ type cold groundwaters to the Na–HCO₃–SO₄ type thermal waters. Thermal waters are oversaturated at discharge temperatures for aragonite, calcite, quartz, chalcedony, magnesite and dolomite minerals giving rise to a carbonate-rich scale. Gypsum and anhydrite minerals are undersaturated with all of the thermal waters. Boiling during ascent of the thermal fluids produces steam and liquid waters resulting in an increase of the concentrations of the constituents in discharge waters. Steam fraction, y, of the thermal waters of which temperatures are above 100°C is between 0.075 and 0.119. Reservoir pH is much lower than pH measured in the liquid phase separated at atmospheric conditions, since the latter experienced heavy loss of acid gases, mainly CO₂. Assessment of the various empirical chemical geothermometers and geochemical modelling suggest that reservoir temperatures vary between 175°C and 200°C.