Mantle heat flow and geotherms for major tectonic units in central Europe

The results of seismic measurements along the deep seismic sounding profile VII and terrestrial heat flow measurements used for construction of heat generation models for the crust in the Paleozoic Platform region, the Sudetic Mountains (Variscan Internides) and the European Precambrian Platform show considerable differences in mantle heat flow and temperatures. At the base of the crust variations from 440–510°C in the models of Precambrian Platform to 700–820°C for the Paleozoic Platform and the Variscan Internides (Sudets) are found. These differences are associated with considerable mantle heat flow variations.The calculated models show mantle heat flow of about 8.4–12.6 mW m^–2 for the Precambrian Platform and 31 mW m^–2 to 40.2 mW m^–2 for Paleozoic orogenic areas. The heat flow contribution originating from crustal radioactivity is almost the same for the different tectonic units (from 33.5 mW m^–2 to 37.6 mW m^–2). Considerable physical differences in the lower crust and upper mantle between the Precambrian Platform and the adjacent areas, produced by lateral temperature variations, could be expected. On the basis of carbon ratio data it can be concluded that the Carboniferous paleogeothermal gradient was much lower in the Precambrian Platform area than in the Paleozoic Platform region.     

Hydraulic parameters and solute velocities in triple-porosity karstic-fissured-porous carbonate aquifers: case studies in southern Poland

 Rock and flow parameters of three karstic-fissured-porous aquifers in the Krakow-Silesian Triassic formations were measured using various methods and compared. Though cavern and fissure porosities are shown to be very low (cavern porosity below 0.5% and fracture porosity below 0.2%), they contribute dominantly to the hydraulic conductivity (from about 1.3×10^–6 to about 11×10^–6 m/s). Matrix porosity (2–11%) is shown to be the main water reservoir for solute transport and the main or significant contributor to the specific yield (<2%). Though the matrix porosity is shown to be much larger than the sum of the cavern and fissure porosities, its contribution to the total hydraulic conductivity is practically negligible (hydraulic conductivity of the matrix is from about 5×10^–11 m/s to about 2×10^–8 m/s). On the other hand, the matrix porosity (for neglected cavern and fissure porosities) when combined with tracer ages (or mean travel times) is shown to yield proper values of the hydraulic conductivity (K) by applying the following formula: K≅(matrix porosity×mean travel distance)/(mean hydraulic gradient×mean tracer age). Confirming earlier findings of the authors, this equation is shown to be of great practical importance because matrix porosity is easily measured in the laboratory on rock samples, whereas cavern and fracture porosities usually remain unmeasurable. Received: 21 February 1997 · Accepted: 13 May 1997     

Coupled HM analysis using zero‐thickness interface elements with double nodes—Part II: Verification and application

In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd.     

A complex variable solution for stress and displacement field around a lined circular tunnel at great depth

By applying the series expansion technique in the complex variable method established by Muskhelishvili, the plane elasticity problem for the stress and displacement field around a lined circular tunnel in conjunction with the consideration of misfit and interaction between the liner and the surrounding geomaterial is dealt with. The tunnel is assumed to be driven in a homogeneous and isotropic geomaterial. The coefficients in the Laurent series expansion of the stress functions are determined. The complex potentials in the liner and the surrounding geomaterial are explicitly derived, respectively. As an example, the case of a lined circular tunnel located in an isotropic initial stress field but subjected to uniform internal pressure is numerically considered. Numerical results indicate that the installation of tunnel liner can reduce the influences of the tunnel excavation on the in situ displacement and stress fields. However, the relative thickness and rigidity of the liner should be in an appropriate range. In addition, the effect of the tunnel excavation upon the displacement field is more significant than that upon the stress field. As far as the stress field in the surrounding geomaterial is concerned, when the ratio between the cover depth of tunnel and the tunnel radius is larger than 5, the results for the stress field in the paper are applicable. When the ratio between the tunnel depth and the tunnel radius is larger than 20, the results are applicable for the displacement field. Copyright © 2008 John Wiley & Sons, Ltd.     

Autochthonous block fields in southern Norway: Implications for the geometry,thickness, and isostatic loading of the Late Weichselian Scandinavian ice sheet

The consistent geographical and altitudinal distribution of autochthonous block fields (mantle of bedrock weathered in situ) and trimlines in southern Norway suggests a multi-domed and asymmetric Late Weichselian ice sheet. Low-gradient ice-sheet profiles in the southern Baltic region, in the North Sea, and along the outer fjord areas of southern Norway, are best explained by movement of ice on a bed of deforming sediment, although water lubricated sliding or a combination of the two, may not be excluded. The ice-thickness distribution of the Late Weichselian Scandinavian ice sheet is not in correspondence with the modern uplift pattern of Fennoscandia. Early Holocene crustal rebound was apparently determined by an exponential, glacio-isostatic rise. Later, however, crustal movements appear to have been dominated by large-scale tectonic uplift of the Fennoscandian Shield, centred on the Gulf of Bothnia, the region of maximum lithosphere thickness.     

Numerical estimation of effective diffusion coefficients for charged porous materials based on micro-scale analyses

In order to describe diffusive transport of solutes through a porous material, estimation of effective diffusion coefficients is required. It has been shown theoretically that in the case of uncharged porous materials the effective diffusion coefficient of solutes is a function of the pore morphology of the material and can be described by the tortuosity (tensor) (Bear, 1988 [1]). Given detailed information of the pore geometry at the micro-scale the tortuosity of different materials can be accurately estimated using homogenization procedures. However, many engineering materials (e.g., clays and shales) are characterized by electrical surface charges on particles of the porous material which strongly affect the (diffusive) transport properties of ions. For these type of materials, estimation of effective diffusion coefficients have been mostly based on phenomenological equations with no link to underlying micro-scale properties of these charged materials although a few recent studies have used alternative methods to obtain the diffusion parameters (Jougnot et al., 2009; Pivonka et al., 2009; Revil and Linde, 2006 2, 3 and 4). In this paper we employ a recently proposed up-scaled Poisson–Nernst–Planck type of equation (PNP) and its micro-scale counterpart to estimate effective ion diffusion coefficients in thin charged membranes. We investigate a variety of different pore geometries together with different surface charges on particles. Here, we show that independent of the charges on particles, a (generalized) tortuosity factor can be identified as function of the pore morphology only using the new PNP model. On the other hand, all electro-static interactions of ions and charges on particles can consistently be captured by the ratio of average concentration to effective intrinsic concentration in the macroscopic PNP equations. Using this formulation allows to consistently take into account electrochemical interactions of ions and charges on particles and so excludes any ambiguity generally encountered in phenomenological equations.     

滑动构造影响下煤层厚度与构造破碎带关系分析

徐庄滑动构造是河南省荥巩煤田主要构造之一,对区内二1煤层的赋存特征影响较大。经对大量计河井田勘查资料分析,发现滑动构造破碎带厚度与其下伏二1煤层厚度相关程度较高,表现为两者在中、浅部厚度变化大,并沿倾向呈薄、厚相间的条带状展布,大致在二1煤层-300m水平附近,滑动构造破碎带厚度变薄尖灭,煤层厚度亦相应变小且趋于稳定。通过对钻孔数据的回归分析,认为二1煤层厚度与其上滑动构造破碎带厚度呈正相关关系,这一结论对后期矿山开发中工作面的布置及做好安全防范措施有指导意义。     

综合物探技术在煤系上覆薄硬岩矿区的应用初探——以山西潞安集团司马矿区为例

山西潞安集团司马矿区东高西低,西部延至长治盆地,矿区水系发育。主采3#煤层埋藏较浅,层位稳定。新生界不整合于下伏基岩之上,厚度变化较大,属典型的薄硬岩煤矿区。为准确提供矿区煤炭开采上限的地质依据,采用了综合勘探技术,其中利用三维地震勘探技术追踪新生界底面反射波及第四系"天窗"位置,利用直流电法与瞬变电磁法提供电性分层,并对其解释结果进行相互印证。勘探实例表明,三维地震勘探、电法勘探与钻探的综合运用不仅可解决薄硬岩煤矿区新生界厚度、含水层富水性及第四系缺失等问题,还可提供新生界含水层与基岩含水层之间的水力联系、制约因素等资料。     

Seasonal changes in sea ice conditions along the Northeast Passage in 2007 and 2012

Remote sensing data from passive microwave and satellite-based altimeters, associated with the data measured underway, were used to characterize seasonal and spatial changes in sea ice conditions along...     

基于卫星高度计的北极海冰厚度变化研究

A modified algorithm taking into account the first year(FY) and multiyear(MY) ice densities is used to derive a sea ice thickness from freeboard measurements acquired by satellite altimetry ICESat(2003–2008). Estimates agree with various independent in situ measurements within 0.21 m. Both the fall and winter campaigns see a dramatic extent retreat of thicker MY ice that survives at least one summer melting season. There were strong seasonal and interannual variabilities with regard to the mean thickness. Seasonal increases of 0.53 m for FY the ice and 0.29 m for the MY ice between the autumn and the winter ICESat campaigns, roughly 4–5 month separation, were found. Interannually, the significant MY ice thickness declines over the consecutive four ICESat winter campaigns(2005–2008) leads to a pronounced thickness drop of 0.8 m in MY sea ice zones. No clear trend was identified from the averaged thickness of thinner, FY ice that emerges in autumn and winter and melts in summer. Uncertainty estimates for our calculated thickness, caused by the standard deviations of multiple input parameters including freeboard, ice density, snow density, snow depth, show large errors more than 0.5 m in thicker MY ice zones and relatively small standard deviations under 0.5 m elsewhere. Moreover, a sensitivity analysis is implemented to determine the separate impact on the thickness estimate in the dependence of an individual input variable as mentioned above. The results show systematic bias of the estimated ice thickness appears to be mainly caused by the variations of freeboard as well as the ice density whereas the snow density and depth brings about relatively insignificant errors.