Abstraction and use of groundwater for heating buildings in a hard rock environment

In order to turn more efficient the heating of class rooms in the lower floor of the old building of the University of Évora (a XVI century building), five drillings were organised inside the area of the university (Figure 1). The purpose was to use the temperature differential of groundwater in relation to air, by means of a heat exchanger, and use this process to heat the rooms using less energy, turning the heating process less expensive. The wells were drilled in fractured rocks (gneisses), and the purpose was to locate them at least around 100 m one from each other, whilst trying to have a hydraulic connection in-between. From the five initial wells, four were successful in terms of productivity, but just two of them (RA1 and RA2) proved to be hydraulically connected. The wells were equipped with screens for all their drilled depth (100 m), except for the first six meters and some two or three pipes of six meters each, to allow space for the installation for submersible pumps. The length of the installed screens guarantees a good efficiency of the system. In the wells with no connection, the heating system can work using each single well for abstraction and injection, but the process is much less efficient than in the cases where interaction between wells is possible through the rock’s fracture network.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

复杂裂隙岩体统计断裂本构关系与变形参数研究

Rock mass is widely recognized as a kind of geologic body which consists of rock blocks and discontinuities. The deformation and failure of rock mass is not only determined by rock block，but also by discontinuity which is virtually more important. Mutual cutting and combination of discontinuities controlled mechanical property of rock mass. The complex cutting of discontinuities determine the intense anisotropy on mechanical property of rock mass，especially under the effect of ground stress. Engineering practice has show that the brittle failure of hard rock always occurs when its working stress is far lower than the yield strength and compressive strength，the failure always directly related to the fracture propagation of discontinuities. Fracture propagation of discontinuities is the virtue of hard rock’s failure. We can research the rock mass discontinuous mechanical properties precisely by the methods of statistical analysis of discontinuities and Fracture Mechanics. According to Superposition Principle in Fracture Mechanics，A Problem or C Problem could be chosen to research. Problem A mainly calculates the crack-tip stress field and displacement field on internal discontinuities by numerical method. Problem C calculate the crack-tip stress field and displacement field under the assumption of that the mainly rock mass stress field has been known. So the Problem C avoid the complex mutual interference of stress fields of discontinuities，which is called crack system problem in Fracture Mechanics. To solve Problem C, field test on stress field in the rock mass is needed. The linear Superposition of discontinuities strain energies are Scientific and Rational. The difference of Fracture Mechanics between rock mass and other materials can mostly expression as：other materials Fracture Mechanics mostly face the problem A，and can’t avoid multi-crack puzzle, while the Rock mass Fracture Mechanics answer to the Problem C. Problem C can avoid multi-discontinuities mutual interference puzzle via the ground stress test. On the basis of Problem C, Fracture Mechanics could be used conveniently in rock mass. The rock mass statistics fracture constitutive relations, which introduced in this article, are based on the Problem C and the Discontinuity Strain Energy linear superposition. This constitutive relation has several merits: first, it is physical constitutive relation rather than empirical; second, it is very fit to describe the rock mass anisotropy properties; third, it elaborates the exogenous factors such as ground stress. The rock mass statistics fracture constitutive relation is the available approach to answer to the physical, anisotropic and ground stress impacted rock mass problems. This article stand on the foundation of predecessor’s statistics fractures constitutive relation, and improved the discontinuity distributive function. This article had derived the limitation of negative exponential distribution in the course of regression analysis, and advocated to using the two parameter negative exponential distribution for instead. In order to solve the problems of two-dimension stability on engineering key cross-sectional view in rock mass, this article derived the rock mass planar flexibility tensor, and established rock mass two-dimension penetrate statistics fracture constitutive relation on the basis of penetrate fracture mechanics. Based on the crack tip plasticity research production of penetrate fracture, for example the Irwin plasticity equifinality crack, this article established the way to deal with the discontinuity stress singularity and plastic yielding problem at discontinuity tip. The research on deformation parameters is always the high light region of rock mass mechanics field. After the dam foundation excavation of XiaoWan hydroelectric power station, dam foundation rock mass upgrowthed a great deal of unload cracks, rock mass mechanical property gotten intricacy and strong anisotropy. The dam foundation rock mass mostly upgrowthed three group discontinuities： the decantation discontinuity, the steep pitch discontinuity, and the schistosity plane. Most of the discontinuities have got partial unload looseness. In accordance with ground stress field data, the dam foundation stress field greatly non-uniform, which felled under the great impaction of tectonic stress field, self-weight stress field, excavation geometric boundary condition, and excavation, unload. The discontinuity complexity and stress field heterogeneity, created the rock mass mechanical property of dam foundation intricacy and levity. The research on the rock mass mechanics, if not take every respected influencing factor into consideration as best as we can, major errors likely to be created. This article calculated the rock mass elastic modulus that after Xiao Wan hydroelectric power station dam foundation gutter excavation finished. The calculation region covered possession monolith of Xiao Wan concrete double-curvature arch dam. Different monolith were adopted the penetrate fracture statistics constitutive relation or bury fracture statistics constitutive relation selectively. Statistics fracture constitutive relation is fit for the intensity anisotropy and heterogeneity rock mass of Xiao Wan hydroelectric power station dam foundation. This article had contrastive analysis the statistics fracture constitutive relation result with the inclined plane load test actual measurement elastic modulus and RMR method estimated elastic modulus, and find that the three methods elastic modulus have got greatly comparability. So, the statistics fracture constitutive relations are qualified for trust. Generally speaking，this article had finished following works based on predecessors job: “Argumentation the C Problems of superposition principle in Fracture Mechanics, establish two-dimension penetrate statistics fracture constitutive relation of rock mass, argue the negative exponential distribution limitation and improve it, improve of the three-dimension berry statistics fracture constitutive relation of rock mass, discontinuity-tip plastic zone isoeffect calculation, calculate the rock mass elastic modulus on two-dimension cross-sectional view”. The whole research clue of this article inherited from the “statistics rock mass mechanics” of Wu Faquan（1992）.

Hidrogeomorfologia e sustentabilidade de recursos hídricos subterrâneos

A água subterrânea de rochas duras é uma fonte importante para fins domésticos, industriais e agrícolas e mesmo para o consumo humano. A geologia, a tectónica, a geomorfologia e as características hidrológicas controlam o fluxo, ocorrência e armazenamento das águas subterrâneas. A disponibilidade da água subterrânea no meio geológico está totalmente dependente das áreas de recarga e de descarga numa determinada bacia. A precipitação é a principal fonte de recarga em aquíferos descontínuos, enquanto que a descarga depende dos declives do terreno e dos gradientes do nível hidrostático e ainda das condições hidrogeológicas do solo. A hidrogeomorfologia é um domínio interdisciplinar emergente, que estuda as relações entre as unidades geomorfológicas e o regime das águas superficiais e subterrâneas de uma determinada área. A compreensão do papel da geomorfologia é essencial para avaliar de forma rigorosa os sistemas hidrogeológicos e os recursos hídricos. Os dados de detecção remota providenciam uma informação espacial valiosa e actualizada da superfície terrestre e dos recursos naturais. Os recentes avanços tecnológicos colocaram as técnicas de detecção remota e os sistemas de informação geográfica (SIG) numa posição cimeira como ferramentas de gestão metodológica. Foi criada, em ambiente SIG, uma base de geo-dados, essencialmente derivada da detecção remota, da cartografia e do trabalho de campo. Esta base de dados, organizada em diferentes níveis de informação, inclui uma avaliação principalmente focalizada no uso do solo, climatologia, declives, geologia, geomorfologia e hidrogeologia. No presente estudo foram cruzados diversos níveis de informação, com a geração de múltiplos mapas temáticos para atingir um quadro integrado dos diversos sectores no Norte e Centro de Portugal. Os sectores em estudo (Caldas da Cavaca, Termas de Entre-os-Rios, Águas de Arouca e Águas do Alardo) estão localizados em sistemas hidrogeológicos predominantemente constituídos por rochas graníticas, por vezes intersectadas por filões de quartzo, aplito-pegmatíticos e doleríticos. Para apoiar a elaboração dos mapas hidrogeomorfológicos foi criada uma base SIG, contendo diversa informação, nomeadamente topografia, hidrografia, litologia, tectónica, morfoestrutura, hidrogeologia, geofísica e uso do solo. Além disso, foram realizadas várias campanhas de campo, as quais permitiram: o estabelecimento dum mapeamento geológico, geomorfológico e hidrogeológico; a caracterização in situ do grau de alteração, resistência e grau de fracturação dos maciços rochosos; o desenvolvimento de um inventário hidrogeológico em conjunto com alguns ensaios expeditos in situ. A interligação entre os parâmetros geomorfológicos, hidrológicos e hidrogeológicos dos sistemas de água subterrânea “normal” e hidromineral destaca a importância de uma cartografia e duma modelação conceptual hidrogeomorfológica. Além disso, contribuirá para um melhor apoio à decisão na gestão sustentável dos recursos hídricos.

Caracterização geológico-geofísica do meio aqüífero fissural : uma contribuição aos modelos de fluxo e armazenamento de água subterrânea

Northeastern Brazil is mainly formed by crystalline terrains (around 60% in area). Moreover, this region presents a semi-arid climate so that it is periodically subject to drought seasons. Furthermore, ground water quality extracted fromwells usually presents poor quality because of their high salinity contents. Nevertheless, ground water is still a very important source of water for human and animal consumption in this region. Well sitting in hard rocks terrains in Northeastern Brazil offers a mean success index of aboul 60%, given that a successful siting is defined by a well producing at least 0.5 m³/h. This low index reveals lack of knowledga about the true conditions of storage and percolation of ground water in crystalline rocks. Two models for structures storing and producing ground water in crystalline rocks in Northeastem Brazil have been proposed in the literature. The first model,tradnionally used for well sitting since the sixties are controlled by faults or fractures zones. This model is commonly referred, in Brazilian hydrogeological literature, as the "creek-crack" model (riacho-fenda in Portuguese). Sites appearing to present dense drainage network are preferred for water well siting - particularly at points where the drainages cross-cul each other. Field follow up work is usually based only on geological criteria. The second model is the "eluvio-alluvial through" (calha eluvio-aluvionar in Portuguese); it is also described in the literature but it is not yet incorporated in well sitting practice. This model is based on the hypothesis that reclilinear drainages can also be controlled by the folietion of the rock. Eventually, depending upon the degree of weathering, a through-shaped structure filled with sediments (alluvium and regolith) can be developed which can store and water can be produced from. Using severalfield case studies, this Thesis presents a thorough analysis ofthe two above cited models and proposes a new model. The analysis is based on an integrated methodological approach using geophysics and structural geology. Both land (Resitiviy and Ground Penetrating Radar- GPR) and aerogeophysical (magnetics and frequency domain eletromagnetics) surveys were used. Slructural analysis emphasized neolectonic aspects; in general, itwas found that fractures in the E-W direction are relatively open, as compared to fracturas inthe N-S direction, probably because E-W fractures were opened by the neotectonic stress regime in Northeastern Brazil, which is controlled by E-W compression and N-S extension. The riacho-fenda model is valid where drainages are controlled by fractures. The degree of fracturing and associated weathering dictale the hydrogeological potential of the structure. Field work in structural analogues reveals that subvertical fractures show consistent directions both in outcrop and aerophotograph scales. Geophysical surveys reveal subvertical conductive anomalies associated to the fracture network controlling the drainage; one of the borders of the conductive anomaly usually coincide wih the drainage. An aspect of particular importance to the validation of fracture control are the possible presence of relalively deep conductive anomalies wihoul continuation or propagalion to the surface. The conductive nature of lhe anomaly is due to the presence of wealhered rock and sedirnenls (alluvium and/or regolilh) storing ground waler which occur associated to the fracture network. Magnetic surveys are not very sensisnive to these structures.lf soil or covering sedirnents are resislive (> 100 Ohm.m), GPR can ba used to image precisely lhe fracture network. A major limialion of riacho-fenda model, revealed by GPR images, is associated to the fact thal subhorizontal fractures do play a very important role in connecting the fracture network, besides connect shallow recharge zones to relalively deep subvertical frecture zones. Iffractures play just a secondary control on the drainage, however, r/acho-fenda model may have a very limiled validny; in these cases, large portions oflhe drainage do nol coincide wilh frectures and mosl oflhewells localed in lhe drainage surrounding would resull dry. Usually, a secondary conlrol on lhe drainage by Ihefraclure networkcan be revealed only wilh detailed geophysical survey. The calha elClv1o-aluvlonarmodel is valid where drainages are conlrolled by folialion. The degree 01 wealhering 01 lhe lolialion planes dictales lhe hydrogeological polenlial 01 lhe slruclure. Outcrop analysis reveals Ihal lolialion and drainage direclions are parallel and Ihal no Iraclures, orfraclures wilh diflerent directions 01 lhe drainage direclion occur. Geophysical surveys reveal conduclive anomalies in a slab lorm associaled 10 lhe Ihrough 01 lhe wealhered rock and sedimenls (alluvium and/or regolith). Magnelic surveys can ofler a very good conlrol on lolialion direclion. An importanl aspect 10 validale lolialion conlrol are lhe presence 01 conductive anomalies showing shallow and deep portions area which are linked. Illhere is an exlensive soil cover, r/acho-fenda and calha eIClv1o-aluv/onar conlrols can be easily misinlerpreled in lhe absence 01 geophysical conlrol. Certainly, Ihis lacl could explain at leasl a part of lhe failure index in well sitting. The model wealhering sack (bolsllo de Intempertsmo in Portuguese) is proposed to explain cases where a very inlensive wealhering occur over lhe crystalline rock so Ihal a secondary inlerslilial porosity is crealed. The waler is Ihen stored in lhe porous of lhe regolilh in a similar mannerlo sedimentary rocks. A possible example ofthis model was delecled by using land geophysical survey where a relalivelyvery deep isolaled conduclive anomaly, in a slab form, was delected. Iflhis structure does store ground waler, certainly Ihere must be a link 01 lhe deep slructure wilh lhe surface in orderlo provide walerfeeding. This model mighl explain anomalous waler yields as greal as 50 m³/h Ihalsomelimescan occur in crystalline rocks in Northeaslern Brazil

Beurteilung der Eignung der NPEMFE-Methode (Natural Pulsed Electromagnetic Field of Earth) mit dem "Cereskop" in Rutschungen und in Locker- und Festgesteinen mit Spannungsänderungen im Mittel- und Hochgebirge

Spannungsumlagerungen in Mineralen und Gesteinen induzieren in geologisch aktiven Bereichen mikromechanische und seismische Prozesse, wodurch eine schwache natürliche elektromagnetische Strahlung im Niederfrequenzbereich emittiert wird. Die elektromagnetischen Emissionen von nichtleitenden Mineralen sind auf dielektrische Polarisation durch mehrere physikalische Effekte zurückzuführen. Eine gerichtete mechanische Spannung führt zu einer ebenso gerichteten elektromagnetischen Emission. Die Quellen der elektromagnetischen Emissionen sind bekannt, jedoch können sie noch nicht eindeutig den verschiedenen Prozessen in der Natur zugeordnet werden, weshalb im Folgenden von einem seismo-elektromagnetischen Phänomen (SEM) gesprochen wird. Mit der neuentwickelten NPEMFE-Methode (Natural Pulsed Electromagnetic Field of Earth) können die elektromagnetischen Impulse ohne Bodenkontakt registriert werden. Bereiche der Erdkruste mit Spannungsumlagerungen (z.B. tektonisch aktive Störungen, potenzielle Hangrutschungen, Erdfälle, Bergsenkungen, Firstschläge) können als Anomalie erkannt und abgegrenzt werden. Basierend auf dem heutigen Kenntnisstand dieser Prozesse wurden Hangrutschungen und Locker- und Festgesteine, in denen Spannungsumlagerungen stattfinden, mit einem neuentwickelten Messgerät, dem "Cereskop", im Mittelgebirgsraum (Rheinland-Pfalz, Deutschland) und im alpinen Raum (Vorarlberg, Österreich, und Fürstentum Liechtenstein) erkundet und die gewonnenen Messergebnisse mit klassischen Verfahren aus Ingenieurgeologie, Geotechnik und Geophysik in Bezug gesetzt. Unter Feldbedingungen zeigte sich großenteils eine gute Übereinstimmung zwischen den mit dem "Cereskop" erkundeten Anomalien und den mit den konventionellen Verfahren erkundeten Spannungszonen. Auf Grundlage der bisherigen Kenntnis und unter Einbeziehung von Mehrdeutigkeiten werden die Messergebnisse analysiert und kritisch beurteilt.

(Table 1) List of samples with carbonate mineralogy, stable isotope, and accessory-mineral information

Petrographic and stable-isotope (d13C, d18O) patterns of carbonates from the Logatchev Hydrothermal Field (LHF), the Gakkel Ridge (GR), and a Late Devonian outcrop from the Frankenwald (Germany) were compared in an attempt to understand the genesis of carbonate minerals in marine volcanic rocks. Specifically, were the carbonate samples from modern sea floor settings and the Devonian analog of hydrothermal origin, low-temperature abiogenic origin (as inferred for aragonite in serpentinites from elsewhere on the Mid-Atlantic Ridge), or biogenic origin? Aragonite is the most abundant carbonate mineral in serpentinites from the two modern spreading ridges and occurs within massive sulfides of the LHF. The precipitation and preservation of aragonite suggests high Mg2+ and sulfate concentrations in fluids. Values of d18OPDB as high as +5.3 per mill for serpentinite-hosted aragonite and as high as +4.2 per mill for sulfide-hosted aragonite are consistent with precipitation from cold seawater. Most of the corresponding d13C values indicate a marine carbon source, whereas d13C values for sulfide-hosted aragonite as high as +3.6 per mill may reflect residual carbon dioxide in the zone of methanogenesis. Calcite veins from the LHF, by contrast, have low d18OPDB (-20.0 per mill to -16.1 per mill) and d13C values (-5.8 per mill to -4.5 per mill), indicative of precipitation from hydrothermal solutions (~129°-186°C) dominated by magmatic CO2. Calcite formation was probably favored by fluid rock interactions at elevated temperatures, which tend to remove solutes that inhibit calcite precipitation in seawater (Mg2+ and sulfate). Devonian Frankenwald calcites show low d18O values, reflecting diagenetic and metamorphic overprinting. Values of d13C around 0 per mill for basalt-hosted calcite indicate seawater-derived inorganic carbon, whereas d13C values for serpentinite-hosted calcite agree with mantle-derived CO2 (for values as low as -6 per mill) with a contribution of amagmatic carbon (for values as low as -8.6 per mill), presumably methane. Secondary mineral phases from the LHF for which a biogenic origin appears feasible include dolomite dumbbells, clotted carbonate, and a network of iron- and silica-rich filaments.

Packaging and testing of fiber Bragg gratings for use as strain sensor for rock specimens

This paper reports a packaging and calibration procedure for surface mounting of fiber Bragg grating (FBG) sensors to measure strain in rocks. The packaging of FBG sensors is performed with glass fiber and polyester resin, and then subjected to tensile loads in order to obtain strength and deformability parameters, necessaries to assess the mechanical performance of the sensor packaging. For a specific package, an optimal curing condition has been found, showing good repeatability and adaptability for non-planar surfaces, such as occurs in rock engineering. The successfully packaged sensors and electrical strain gages were attached to standard rock specimens of gabbro. Longitudinal and transversal strains under compression loads were measured with both techniques, showing that response of FBG sensors is linear and reliable. An analytical model is used to characterize the influences of rock substrate and FBG packaging in strain transmission. As a result, we obtained a sensor packaging for non-planar and complex natural material under acceptable sensitivity suitable for very small strains as occurs in hard rocks.

Cálculo de la recarga natural en grandes áreas, en función de la litología y la precipitación = Calculation of natural recharge in large areas, as a function of lithology and rainfall

En una región amplia como España se demuestra —mediante inferencias estadísticas sobre una muestra completa de 875 manantiales en los que se conoce su caudal medio y la litología de su área de alimentación y que han sido agrupados en regiones de distinta pluviometría— que la recarga media anual es una fracción fija de la precipitación media para cada litología. Se han establecido así unas tasas de recarga respecto de la precipitación para seis grupos litológicos de diferente permeabilidad: arenas, gravas y formaciones aluviales en general, 8.3%; conglomerados, 5.6%; areniscas, 7.3%; calizas y dolomías, 34.3%; margas, margocalizas, limos y arcillas, 3.3%; otras rocas, 1.3%. Teniendo en cuenta la representatividad de España, la cual tiene una gran variabilidad de litología, pluviometría, topografía, etcétera, estas tasas de recarga respecto de la precipita-ción son probablemente valores cuasi universales que pueden ser utilizados para estimar la recarga media o los recursos hídricos subterráneos medios de regiones amplias en cualquier parte del mundo, salvo en regiones especiales, como las que tienen permafrost, por ejemplo. En todo caso, estas tasas de recarga podrían ser retocadas para cada región según sus particulares características. Los datos de precipitación y litología son muy corrientes, por lo que el método puede ser ampliamente utilizado para completar balances hidráulicos.In a region as large as Spain, annual mean recharge is shown to be a fixed proportion of the mean rainfall for each lithology. This determination is based on statistical inferences from a complete sample of 875 springs for which mean flow and catchment areas are known and which have been grouped into distinct rainfall regions. Recharge rates have thus been established with respect to rainfall for six lithological groups with different permeability: sands, gravels and generally alluvial formations, 8.3%; conglomerates, 5.6%; sandstones, 7.3%; limestone and dolomite 34.3%; marls, marly limestones, silts and clays, 3.3%; and hard rocks, 1.3%. Considering the representativeness of Spain, which is large in size and has a highly varied lithology, topography and rainfall, these recharge rates for rainfall are probably quasi-universal values that can be used to estimate average recharge or average groundwater resources of large regions in any part of the world (except in special cases such as areas with permafrost, for example). For any case, these recharge rates can be adapted to each region according to its particular characteristics. Rainfall and lithology data are very common, and so the method can be widely used to calculate hydraulic balances.

Between rocks and hard places

A travel article about the Aland Islands, Finland, that discusses the mix of Swedish, Finnish, and Russian cultural influences in the area. On the map, Finland seems to end in fragments. The gods have stomped their heels on the southwestern corner, and between the cities of Helsinki and Turku it is jagged, rocky islands that form the final landmarks...

Stress-Strain Prediction of Jointed Rocks using Artificial Neural Networks

The applicability of Artificial Neural Networks for predicting the stress-strain response of jointed rocks at varied confining pressures, strength properties and joint properties (frequency, orientation and strength of joints) has been studied in the present paper. The database is formed from the triaxial compression tests on different jointed rocks with different confining pressures and different joint properties reported by various researchers. This input data covers a wide range of rock strengths, varying from very soft to very hard. The network was trained using a 3 layered network with feed forward back propagation algorithm. About 85% of the data was used for training and remaining15% for testing the predicting capabilities of the network. Results from the analyses were very encouraging and demonstrated that the neural network approach is efficient in capturing the complex stress-strain behaviour of jointed rocks. A single neural network is demonstrated to be capable of predicting the stress-strain response of different rocks, whose intact strength vary from 11.32 MPa to 123 MPa and spacing of joints vary from 10 cm to 100 cm for confining pressures ranging from 0 to 13.8 MPa.

The typology of Irish hard rock aquifers based on an integrated hydrogeological and geophysical approach

Groundwater flow in hard-rock aquifers is strongly controlled by the characteristics and distribution of structural heterogeneity. A methodology for catchment-scale characterisation is presented, based on the integration of complementary, multi-scale hydrogeological, geophysical and geological approaches. This was applied to three contrasting catchments underlain by metamorphic rocks in the northern parts of Ireland (Republic of Ireland and Northern Ireland, UK). Cross-validated surface and borehole geophysical investigations confirm the discontinuous overburden, lithological compartmentalisation of the bedrock and important spatial variations of the weathered bedrock profiles at macro-scale. Fracture analysis suggests that the recent (Alpine) tectonic fabric exerts strong control on the internal aquifer structure at meso-scale, which is likely to impact on the anisotropy of aquifer properties. The combination of the interpretation of depth-specific hydraulic-test data with the structural information provided by geophysical tests allows characterisation of the hydrodynamic properties of the identified aquifer units. Regionally, the distribution of hydraulic conductivities can be described by inverse power laws specific to the aquifer litho-type. Observed groundwater flow directions reflect this multi-scale structure. The proposed integrated approach applies widely available investigative tools to identify key dominant structures controlling groundwater flow, characterising the aquifer type for each catchment and resolving the spatial distribution of relevant aquifer units and associated hydrodynamic parameters.