Heat flow and crustal thermal structure in the Late Archaean Closepet Granite batholith,south India

The Late Archaean Closepet Granite batholith in south India is exposed at different crustal levels grading from greenschist facies in the north through amphibolite and granulite facies in the south along a ∼400 km long segment in the Dharwar craton. Two areas, Pavagada and Magadi, located in the Main Mass of the batholith, best represent the granitoid of the greenschist and amphibolite facies crustal levels respectively. Heat flow estimates of 38 mW m⁻² from Pavagada and 25 mW m⁻² from Magadi have been obtained through measurements in deep (430 and 445 m) and carefully sited boreholes. Measurements made in four boreholes of opportunity in Pavagada area yield a mean heat flow of 39 ± 4 (s.d.) mW m⁻², which is in good agreement with the estimate from deep borehole. The study, therefore, demonstrates a clear-cut heat flow variation concomitant with the crustal levels exposed in the two areas. The mean heat production estimates for the greenschist facies and amphibolite facies layers constituting the Main Mass of the batholith are 2.9 and 1.8 μW m⁻³, respectively. The enhanced heat flow in the Pavagada area is consistent with the occurrence of a radioelement-enriched 2-km-thick greenschist facies layer granitoid overlying the granitoid of the amphibolite facies layer which is twice as thick as represented in the Magadi area. The crustal heat production models indicate similar mantle heat flow estimates in the range 12–14 mW m⁻², consistent with the other parts of the greenstone-granite-gneiss terrain of the Dharwar craton.     

Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

High heat flow anomalies on an old oceanic plate observed seaward of the Japan Trench

Thirty-three new measurements on the seaward slope and outer rise of the Japan Trench along a parallel of 38°45′N revealed the existence of high heat flow anomalies on the subducting Pacific plate, where the seafloor age is about 135 m.y.. The most prominent anomaly with the highest value of 114 mW/m² is associated with a small mound on the outer rise, which was reported to be a kind of mud volcano. On the seaward slope of the trench, heat flow is variable: high (70–90 mW/m²) at some locations and normal for the seafloor age (about 50 mW/m²) at others. The spatial variation of heat flow may be related to development of normal faults and horst/graben structures due to bending of the Pacific plate before subduction, with fluid flow along the fault zones enhancing the vertical heat transfer. Possible heat sources of the high heat flow anomalies are intra-plate volcanism in the last several million years like that discovered recently on the Pacific plate east of the Japan Trench.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Evaluation of change in subsurface thermal environment due to groundwater flow in the Tokyo Lowland,Japan

Information on the distribution of subsurface temperature and hydraulic heads at 24 observation wells in and around the Tokyo Lowland, the eastern part of the Tokyo Metropolitan area, were examined to make clear the relationship between groundwater and the subsurface thermal environment in the urban area. Minimums in temperature–depth logs due to subsurface temperature increasing at shallow parts were recognized in 21 wells. This fact shows subsurface temperature is affected by ground surface warming in almost all of this area. Deeper than minimums, where the effects of surface warming became relatively small, regional variation is observed as follows: high temperatures are shown in the central part to the southern part, and low temperatures shown in the inland to eastern part. The high temperature area corresponds to an area where the lower boundary of groundwater flow is relatively shallow. This area corresponds also to an area with severe land subsidence resulting from excessive groundwater pumping. It is considered that this high temperature area is formed by the effects of upward groundwater flow affected by hydrogeological conditions and pumping. On the other hand, a comparison between past data (1956–1967) and present data (2001–2003) revealed widespread decreasing temperature in the inland area. This is explained by downward groundwater flow based on an analysis of temperature–depth logs. This fact suggests that subsurface temperature is not only increasing from the effects of surface warming but also decreasing from the effects of groundwater environment change due to pumping.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

鄂尔多斯高原地下水流系统的多层结构循环模式—来自深孔中PACKER系统分层水头测定的证据

地下水流系统和循环模式分析是研究地下水形成机理的基础,对正确认识和评价地下水资源具有重要意义。鄂尔多斯高原在含水系统和众多地下水排泄区的控制下,形成了多个不同规模、不同循环深度、相互独立的地下水流系统。PACKER系统分层试验测定的不同深度水头的数据证明,鄂尔多斯高原地下水流系统存在托斯多层水流模式,区域性水流系统一般包含浅循环、中间循环和深循环3个循环系统。浅循环系统的发育深度在200m以内,深循环系统的发育深度大于400m。     

幔源CO_2释出机理、脱气模式及成藏机制研究进展

针对幔源CO2如何从地幔岩浆中脱出并进入沉积地层中形成CO2气藏聚集这一关键问题,总结了国内外研究进展和前缘方向。研究表明,地幔深部的碱性玄武岩浆和碱性岩浆才是深部流体和CO2等挥发份大量赋存、渗滤和释出的场所。浅成侵入岩、次火山岩和火山通道等是CO2释放和聚集的有利位置,岩浆期后和岩浆衰弱期的热液活动阶段是CO2大量释放和聚集的有利时期。幔源CO2进入沉积盆地中具有3种脱气模式,即沿岩石圈断裂直接脱气模式、热流底辟体脱气模式和壳内岩浆房-基底断裂组合脱气模式。CO2的固有物化性质决定其运移相态多样,具有运移和聚集过程同步的特征。只有在满足大量的化学消耗及地层水或原油的溶解和耗散之后才能形成CO2有效聚集。幔源CO2成藏和分布主要受岩浆气源体和气源断裂体系的控制。今后,在超临界CO2及其对油气运移聚集的作用、CO2与深大断裂及火山岩的关系、CO2脱气运移机制、CO与常规烃类油气的耦合差异成藏机制等方面仍需要进一步的研究和探索。     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Micromechanical approach to unsaturated water flow in structured geomaterials by two-scale computations

This article presents a micromechanical approach to the problem of unsaturated water flow in heterogeneous porous media in transient conditions. The numerical formulation is based on the two-scale model obtained previously by periodic homogenization. It allows for a coupled solution of the non-linear flow equations at macroscopic and microscopic scales and takes into account the macroscopic anisotropy of the medium and the local non-equilibrium of the capillary pressure. The model was applied to simulate two-dimensional water infiltration at constant flux into an initially dry medium containing inclusions of square and rectangular shapes. The numerical results showed the influence of the inclusion–matrix conductivity ratio and the local geometry on the macroscopic behavior. The influence of the conductivity ratio manifested itself by the acceleration or retardation of the onset of the macroscopic water flux at the outlet, while the local geometry (anisotropy) significantly affected the macroscopic spatial distribution of the water flux. Such type of approach can be extended to simulate coupled phenomena (for example hydro-mechanical problems) with evolving local geometry.