Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).     

Mais comment s'écoule donc un glacier ? Aperçu historique

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).     

Up-temperature flow of surface-derived fluids in the mid-crust: the role of pre-orogenic burial of hydrated fault rocks

The Walter‐Outalpa shear zone in the southern Curnamona Province of NE South Australia is an example of a shear zone that has undergone intensely focused fluid flow and alteration at mid‐crustal depths. Results from this study have demonstrated that the intense deformation and ductile shear zone reactivation, at amphibolite facies conditions of 534 ± 20 °C and 500 ± 82 MPa, that overprint the Proterozoic Willyama Supergroup occurred during the Delamerian Orogeny (c. 500 Ma) (EPMA monazite ages of 501 ± 16 and 491 ± 19 Ma). This is in contrast to the general belief that the majority of basement deformation and alteration in the southern Curnamona Province occurred during the waning stages of the Olarian Orogeny (c. 1610–1580 Ma). These shear zones contain hydrous mineral assemblages that cut wall rocks that have experienced amphibolite facies metamorphism during the Olarian Orogeny. The shear zone rock volumes have much lower δ¹⁸O values (as low as 1‰) than their unsheared counterparts (7–9‰), and calculated fluid δ¹⁸O values (5–8‰) consistent with a surface‐derived fluid source. Hydrous minerals show a decrease in δD_(H2O) from ?14 to ?22‰, for minerals outside the shear zones, to ?28 to ?40‰, for minerals within the shear zones consistent with a contribution from a meteoric source. It is unclear how near‐surface fluids initially under hydrostatic pressure penetrate into the middle crust where fluid pressures approach lithostatic, and where fluid flow is expected to be dominantly upward because of pressure gradients. We propose a mechanism whereby faulting during basin formation associated with the Adelaidean Rift Complex (c. 700 Ma) created broad hydrous zones containing mineral assemblages in equilibrium with surface waters. These panels of fault rock were subsequently buried to depths where the onset of metamorphism begins to dehydrate the fault rock volumes evolving a low δ¹⁸O fluid that is channelled through shear zones related to Delamerian Orogenic activity.     

2D Self-Potential tomographies for studying groundwater flows in the Varco d''Izzo landslide (Basilicata, southern Italy)

A geoelectrical monitoring activity has been carried out to improve the geological and hydrogeological knowledge about the Varco d'Izzo landslide (Potenza, Basilicata, Southern Apennine, Italy), an active rotational–translational slide evolved in earth-flow. In this work we have focused on the Self-Potential (SP) method by applying three different SP measuring techniques and combining modern technologies for data acquisition and new methods for tomographic inversion. A SP map and three static SP tomographies have been carried out to better analyse the groundwater circulation system and to better reconstruct the geometry of the landslide body. In the accumulation zone, which is the area more exposed to the geomorphological activity, a new SP measuring strategy has been applied. This strategy, based on time-continuous 2D SP tomographies, helps identify water flow changes in subsurface by studying the time series of SP tomographic images. The analysis of time-dependent changes of water infiltration in near surface is the key to better understand the hydrogeological processes underlying the ground instability phenomena. The time-lapse analysis of tomographic images has allowed us not only to investigate the correlation between the temporal changes of SP signals and rainfall events, but also to quantify the range of these changes. The modification of the distribution of the SP source accumulation zones is associated with the dynamics of the groundwater flows. These preliminary results allow us to consider the SP tomographic method as a tool for geophysical monitoring of landslide areas and encourage to develop new measuring systems for near-real time applications.     

The distribution of radiogenic heat production as a function of depth in the Sierra Nevada Batholith, California

Geochemical analyses and geobarometric determinations have been combined to create a depth vs. radiogenic heat production database for the Sierra Nevada batholith, California. This database shows that mean heat production values first increase, then decrease, with increasing depth. Heat production is 2 μW/m³ within the 3-km-thick volcanic pile at the top of the batholith, below which it increases to an average value of 3.5 μW/m³ at 5.5 km depth, then decreases to 0.5–1 μW/m³ at 15 km depth and remains at these values through the entire crust below 15 km. Below the crust, from depths of 40–125 km, the batholith's root and mantle wedge that coevolved beneath the batholith appears to have an average radiogenic heat production rate of 0.14 μW/m³. This is higher than the rates from most published xenolith studies, but reasonable given the presence of crustal components in the arc root assemblages. The pattern of radiogenic heat production interpreted from the depth vs. heat production database is not consistent with the downward-decreasing exponential distribution predicted from modeling of surface heat flow data. The interpreted distribution predicts a reasonable range of geothermal gradients and shows that essentially all of the present day surface heat flow from the Sierra Nevada could be generated within the 35 km thick crust. This requires a very low heat flux from the mantle, which is consistent with a model of cessation of Sierran magmatism during Laramide flat-slab subduction, followed by conductive cooling of the upper mantle for 70 m.y. The heat production variation with depth is principally due to large variations in uranium and thorium concentration; potassium is less variable in concentration within the Sierran crust, and produces relatively little of the heat in high heat production rocks. Because silica content is relatively constant through the upper 30 km of the Sierran batholith, while U, Th, and K concentrations are highly variable, radiogenic heat production does not vary directly with silica content.     

Evidence for deep groundwater flow and convective heat transport in mountainous terrain, Delta County, Colorado, USA

The Tongue Creek watershed lies on the south flank of Grand Mesa in western Colorado, USA and is a site with 1.5 km of topographic relief, heat flow of 100 mW/m², thermal conductivity of 3.3 W m^–1 °C^–1, hydraulic conductivity of 10^-8 m/s, a water table that closely follows surface topography, and groundwater temperatures 3–15°C above mean surface temperatures. These data suggest that convective heat transport by groundwater flow has modified the thermal regime of the site. Steady state three-dimensional numerical simulations of heat flow, groundwater flow, and convective transport were used to model these thermal and hydrological data. The simulations provided estimates for the scale of hydraulic conductivity and bedrock base flow discharge within the watershed. The numerical models show that (1) complex three-dimensional flow systems develop with a range of scales from tens of meters to tens of kilometers; (2) mapped springs are frequently found at locations where contours of hydraulic head indicate strong vertical flow at the water table, and; (3) the distribution of groundwater temperatures in water wells as a function of surface elevation is predicted by the model.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Characterization of water bursting and discharge into underground mines with multilayered groundwater flow systems in the North China coal basin

A conceptual framework for characterization of water bursting and discharge into underground mines with multilayered groundwater flow systems is presented, based on the features of the site conditions and analyses of the water bursting and mine inundation events of the North China coal basin. Comprehensive analyses of the hydrogeological conditions for establishment of the three-dimensional groundwater flow systems of the North China coal basin revealed different vertical hydraulic connection paths or channels. These connections include karst collapse columns, fault or fracture zones, buried weathering zones, and fracture networks in aquitards, within the multilayered groundwater flow systems. Also examined, was the effect of the primary features of the flow system and those different connection paths on water bursting and discharge into the underground coal mines. It was demonstrated that appropriately identifying and adequately understanding the three-dimensional multilayered groundwater flow systems and those various vertical hydraulic connection mechanisms are critical for appropriately preventing water bursting hazards and predicting water discharge into underground mines. A valuable guideline is presented on characterization of water bursting and discharge into underground mines with multilayered groundwater flow systems.

---

Resumen Un cadre de travail conceptuel pour la caractérisation du jaillissement et de la vidange des eaux dans les mines à plusieurs niveaux d'écoulement est présenté, basé sur les conditions et l'analyse sur site du jaillissement des eaux et des inondations dans les bassins charbonneux du Nord de la Chine. L'analyse compréhensive des conditions hydrogéologiques des systèmes d'écoulement tridimensionnel de ces bassins, révèle des trajets ou des chenaux qui sont autant de connections hydrauliques verticales. Ces connections comprennent des effondrement karstiques en colonne, des zones de fractures et de failles, des zones d'altération, des réseaux de fractures dans les aquitards, dans le système hydrogéologique multicouche. L'influence des caractéristiques primaires du système d'écoulement et de ces différentes connections sur le jaillissement et la vidange dans les mines, est également examiné. Il a été démontré que l'identification appropriée et la compréhension adéquate du système d'écoulement souterrain multicouche et tridimensionnel, et des mécanismes variés de connection hydraulique verticale est critiquée pour prévenir de manière appropriée les risques de jaillissement et de vidange dans les mines. Des guidelines valorisables sont présentées, pour la caractérisation du jaillissement et de la vidange dans de tels systèmes.

---

Résumé Se presenta un marco conceptual para la de caracterización de la despresurización y descarga de agua a minas subterráneas con sistemas de acuíferos de múltiples niveles, el cual se basa en las condiciones del sitio, análisis de la despresurización de agua e inundaciones de minas en la cuenca de carbón del norte de la China.Un análisis comprehensivo de las condiciones hidrogeológicas para el establecimiento de los sistemas de flujo de agua de tri-dimensionales en la cuenca de carbón del norte de la China ha revelado diferentes canales hidráulicos verticales. Estas conexiones incluyen columnas de karst colapsadas, fallas o zonas de fracturas, zonas de regalita hundida y redes de fracturas en zonas de poca permeabilidad dentro de los sistemas de flujo de múltiples niveles.También se investigó el efecto de los rasgos primarios del sistema de flujo y los diferentes canales conectivos sobre la despresurización y descarga de agua a minas subterráneas. Se demostró que el entendimiento e identificación adecuada de los sistemas de acuíferos de múltiples niveles en tres dimensiones es crítico para prevenir apropiadamente la despresurización y descarga de agua a minas subterráneas. Se presentan lineamientos valiosos en relación con la despresurización y descarga de agua a minas subterráneas en sistemas de flujo de múltiple nivel.

     

A Connectivity Index for Discrete Fracture Networks

Connectivity is an important measure for assessing flow transport in rock, especially through fractures. In this paper, rock fracture systems are modelled by a discrete fracture model simulated by a marked point process. A connectivity index is then introduced to quantify the connectivity between any two points in space. Monte Carlo simulation is used to evaluate the connectivity index for stationary cases and relationships between the connectivity index and the parameters of the discrete fracture model are analysed. The average number of intersections per fracture, X_f, and the fracture intensity, P12 (P32), are calculated and the relationships between these parameters and the connectivity index are investigated, concluding that X_f is the more suitable parameter for the classification of rock mass flow properties. The relationships between the connectivity index and the percolation state of the fractured medium are also discussed. An edge correction is briefly discussed and a practical example is used to demonstrate the method of computing the connectivity index.     

西北地区泥石流沟形态因素的非线性分析——以兰州市为例

以兰州市区96条主要的泥石流沟为例,利用分形理论的统计自相似性对泥石流沟数的累计分布随沟谷形态要素沟道比降、汇水面积和沟道长度的变化进行了分析;并与不同流域形态要素的分维进行了比较,探讨了分维不同的原因。兰州市区泥石流沟三个形态要素均具有统计意义上的自相似性,其中以汇水面积显示出分维最小(0.46),沟道长度次之(0.71),沟道比降最大(1.09)。分维的大小描述了作为系统的沟谷及其泥石流的形成的空间结构,即分形结构,同时也反映了不同形态要素对系统所起作用的差异,分维越小,所起作用越显著。其结果对西北干旱半干旱地区石流灾害的进一步研究与防治提供了新的思路。