Diesel transport monitoring in simulated unconfined aquifers using miniature resistivity arrays

The principal aim of this study is to assess the scope of monitoring diesel plume migration in a scaled aquifer model with a miniaturised electrical resistivity array. Respectively 1000 and 500 ml of diesel were injected in both the unsaturated and water-saturated zones of a sand body overlying a clay aquitard, and diesel migration was monitored with a miniature electrode array and an off-the-shelf resistivity meter. Inverted time-lapse electrical resistivity tomography (ERT) data reflect downward and lateral spreading of the diesel plume away from the injection point in the unsaturated zone. Diesel was also imaged to spread upwards and laterally away from the injection point in the saturated zone, as controlled by capillary rise. In both cases later-time ERT images reflected preferential pooling of diesel on the water table, as well as vertical smearing of pooled diesel in response to simulated water-table fluctuations. Repeat fluid electrical conductivity (EC) and dissolved oxygen (DO) measurements validate the observed changes in bulk resistivity caused by both diesel injections. Artefacts introduced by 2D inversion of 3D contaminant transport were abound. Time-lapse ERT imaging of diesel transport is therefore inferred to be feasible and well-suited to complementing conventional techniques of intrusive site investigation, although time-lapse 3D or 4D ERT imaging is strongly advocated.     

3D geological modeling for prediction of subsurface Mo targets in the Luanchuan district,China

Three-dimensional (3D) district-scale geoscience information for the Luanchuan Mo district was integrated for understanding the development of its regional geology and ore-forming processes and for decision-making about potential targets for mineral exploration. The methodology and datasets used were: (1) construction of an initial geological model (25 km × 20 km × 2.5 km) using 1:10,000 scale geological map, nine geological cross-sections and gravity and magnetic data; (2) construction of three large-scale Mo deposits model (5 km × 4 km × 2.5 km) using 1:2000 scale geological and topographic maps, 288 boreholes (total core length of 158,700 m), and 32 1:2000 scale cross-sections; (3) 3D inversion of 1:25,000 scale gravity and magnetic data for identification metallogenic anomaly zones which are associated with Jurassic intrusions; (4) extraction of ore-controlling formation and sequence of the Luanchuan Group using the large-scale 3D models of Mo deposits and results of analysis of lithogeochemical samples from outcrops and borehole cores; (5) identification of ore-forming and ore-controlling faults using the large-scale 3D model of Mo deposits and mineralized Jurassic granite porphyry stocks; (6) boost weights-of-evidence and concentration–volume (C–V) fractal analyses to integrate metallogenic information and to identify and classify potential Mo targets. Four classes of exploration targets were identified using C–V modeling and 3D known orebodies model: the first and second class targets are mainly located in three large magma-skarn type deposit camps, occupying ~ 1.4 km³ with total estimated reserve of ~ 2.3 Mt; the third class targets, which are mainly located in Huangbeiling and Yuku deposit camps comprising concealed magma-skarn type deposits, occupy ~ 2.8 km³ and represent a new target exploration zone in the Luanchuan district; the fourth class targets, which are located in the Huoshenmiao, Majuan, and Daping zones, occupy ~ 15 km³ and represent potential mineral resources with likely similar orebody features as the Yuku deposit.     

Influence of shallow infiltration on time-lapse ERT: Experience of advanced interpretation

Previous time-lapse Electrical Resistivity Tomography (ERT) studies have experienced difficulties in reconstructing reliable calculated resistivity changes in the subsurface. Increases or decreases of resistivity appear in the calculated ERT image where no changes were noted in the subsurface, leading to erroneous hydrological interpretations of the geophysical results. In this article, we investigate how a variation of actual resistivity with time and at shallow depth can influence time-lapse ERT results and produce resistivity artefacts at depth. We use 1 and 2-D numerical modelling to simulate infiltration scenarios. Using a standard time-lapse inversion, we demonstrate the resistivity artefact production according to the electrode spacing parameter. We used an advanced inversion methodology with a decoupling line at shallow depth to attenuate or remove resistivity artefacts. We also applied this methodology to a field data set obtained in a semi-arid environment in Burkina Faso, West Africa. Here, time-lapse ERT shows several resistivity artefacts of calculated resistivity if a standard inversion is used. We demonstrate the importance of a dense sampling of shallow resistivity variations at shallow depth. Advanced interpretation allows us to significantly attenuate or remove the resistivity artefact production at intermediate depth and produce reliable interpretation of hydrological processes.     

Joint use of TEM and MRS methods in a complex geological setting

Transient Electromagnetic (TEM), known also as Time Domain Electromagnetic (TDEM) and Magnetic Resonance Sounding (MRS) methods were applied jointly to investigate variations in lithology and groundwater salinity in the Nahal Hever South area (Dead Sea coast of Israel). The subsurface in this area is highly heterogeneous and composed of intercalated sand and clay layers over a salt formation, which is partly karstified. Groundwater is very saline, with a chloride concentration of 100–225 g/l. TEM is known as an efficient tool for investigating electrically conductive targets like saline water, but it is sensitive to both the salinity of groundwater and the porosity of rocks. MRS, however, is sensitive primarily to groundwater volume, but it also allows delineating of lithological variations in water-saturated formations. MRS is much less sensitive to variations in groundwater salinity in comparison with TEM. We show that MRS enables us to resolve the fundamental uncertainty in TEM interpretation caused by the equivalence between groundwater resistivity and lithology. Combining TEM and MRS, we found that the sandy Dead Sea aquifer filled with Dead Sea brine is characterized by a bulk resistivity of ρ_x > 0.4 Ωm, whereas zones with silt and clay in the subsurface are characterized by a bulk resistivity of ρ_x < 0.4 Ωm. These observations are confirmed by calibration of the TEM method performed near 18 boreholes.     

Late Cenozoic basalt and gabbro in the subsurface in the Phetchabun Basin,Thailand: Implications for the Southeast Asian Volcanic Province

Fragments of basaltic and gabbroic rocks were obtained in cuttings from 15 exploration wells in the Na Sanun area of the Wichian Buri Sub-basin of the Phetchabun Basin in central Thailand. The samples represent flows and sills in lacustrine and fluvial sedimentary rocks of the Lower to mid-Miocene Wichian Buri Group. Mafic igneous units were identified in the sections based on their typically high-amplitude seismic reflections, confirmed by the examination of several hundred well cuttings and magnetic susceptibility measurements. Cross-sections of the sub-basin were constructed on the basis of previously published subsurface interpretations, seismic and well data, and petrological observations. Basaltic flows A, B, E, and F have ages of ca. 2 Ma, 16 Ma, 24 Ma and 18 Ma, based on inferred stratigraphic position. Gabbroic sill C and dioritic sill G are inferred to be correlative at ca. 11.6 Ma, and differ petrologically from ca. 12.8 Ma gabbroic sill D. Major minerals in both basaltic and gabbroic samples are plagioclase (ca. An₅₀), anorthoclase, and augite, with pervasive alteration to Na- and Ca- zeolite minerals and analcime. Leucodioritic sill G also contains amphibole and high Ti-phlogopite. Overall, the rocks show within-plate tholeiitic to alkalic characteristics, and show similarities to basaltic surface outcrops of similar ages in the Wichian Buri-Lop Buri area. No evidence was seen in the subsurface for the andesitic to rhyolitic rocks of similar ages that occur at surface, but their presence cannot be precluded based on our limited data.     

Human Health Risk Assessment of a landfill based on volatile organic compounds emission,immission and soil gas concentration measurements

A Human Health Risk Assessment (HHRA) was required for a closed landfill located in Cerdanyola del Vallès (Barcelona, Spain). The HHRA had two objectives, to evaluate the present risk of the identified receptors in the area and to safely develop the future urban planning of the area, therefore 3 scenarios for the current situation and 4 for the future situation were developed.After reviewing the existing data and exploring the needs of information, the assessment in this study was focused on the measurement of volatile organic compounds (VOCs) fluxes from the subsoil (emission from the landfill at 5 points), concentrations of VOCs in the air (immission in 4 urban sites) and concentration of VOCs in soil–gas (measurements at 5 m below ground surface outside the landfill at 8 sites). Around 70 VOCs were analyzed by using multi-sorbent tubes and Thermal Desorption Gas Chromatography (TD–GC–MS). The VOCs that were detected and quantified include alkanes, aromatic hydrocarbons, alcohols, ketones, halocarbons, aldehydes, esters, terpenoids, ethers and some nitrogenated and sulfur compounds, furans and carboxylic acids. Specific mercury flux measurements were performed in a hot spot by using carulite tubes, that were also analyzed by using Thermal Decomposition, Amalgamation, and Atomic Absorption Spectrophotometry.Results showed average values of volatile emission fluxes ranging from non-detected to 331 μg m⁻² day⁻¹ (dichlorodifluoromethane). In the case of immission, the concentration of VOCs measured in the air of populated area surrounding the landfill ranged values from non-detected to 42.0 μg m⁻³ (acetic acid). The soil–gas measurements in piezometers around the landfill showed individual VOC values with a maximum 830 μg m⁻³ for dichlorodifluoromethane.With the obtained fluxes and concentrations in air and soil–gas, USEPA methodology and modeling was used to evaluate equivalent concentration in the scenarios considered. Toxicity values from IRIS database were used to finally obtain chemical risk indicators. Admissible risk indicators were obtained in all scenarios. The VOCs that contributed more to risk indexes in RH2 were trichloroethylene, trimethylbenzene, chloroform, 1,2-dichloroethane and carbon tetrachloride. The carcinogenic risk in RH7 was linked to the presence of benzene and chloroform. The comparison of the measurements of the present work with other landfills evidence that HHRA in ambient air would be needed in order to perform a correct landfill management.     

Heat flow in the rifted continental margin of the South China Sea near Taiwan and its tectonic implications

Temperature measurements carried out on 9 hydrocarbon exploration boreholes together with Bottom Simulating Reflectors (BSRs) from reflection seismic images are used in this study to derive geothermal gradients and heat flows in the northern margin of the South China Sea near Taiwan. The method of Horner plot is applied to obtain true formation temperatures from measured borehole temperatures, which are disturbed by drilling processes. Sub-seafloor depths of BSRs are used to calculate sub-bottom temperatures using theoretical pressure/temperature phase boundary that marks the base of gas hydrate stability zone. Our results show that the geothermal gradients and heat flows in the study area range from 28 to 128 °C/km and 40 to 159 mW/m², respectively. There is a marked difference in geothermal gradients and heat flow beneath the shelf and slope regions. It is cooler beneath the shelf with an average geothermal gradient of 34.5 °C/km, and 62.7 mW/m² heat flow. The continental slope shows a higher average geothermal gradient of 56.4 °C/km, and 70.9 mW/m² heat flow. Lower heat flow on the shelf is most likely caused by thicker sediments that have accumulated there compared to the sediment thickness beneath the slope. In addition, the continental crust is highly extended beneath the continental slope, yielding higher heat flow in this region. A half graben exists beneath the continental slope with a north-dipping graben-bounding fault. A high heat-flow anomaly coincides at the location of this graben-bounding fault at the Jiulong Ridge, indicating vigorous vertical fluid convection which may take place along this fault.     

Naturally occurring arsenic: Mobilization at a landfill in Maine and implications for remediation

Elevated levels of dissolved arsenic (∼300 μg L⁻¹) have been detected beneath and in groundwater plumes extending away from a closed landfill in southern Maine. This study sought to determine the source of arsenic to the aquifer, the processes responsible for arsenic mobilization, and to evaluate the effectiveness of remediation efforts that have occurred at this site. The As appears to originate in the natural (glacial) aquifer solids, which contain ∼5 mg kg⁻¹ As on a dry weight basis. This conclusion is supported by the relatively uniform distribution of As in sediment samples, results of laboratory batch incubation experiments, and comparisons with groundwaters in nearby wetlands, which also have high levels of dissolved As that do not appear to originate within the landfill. The As is mobilized in the subsurface by strongly reducing conditions beneath the landfill and in nearby wetlands. In the aquifer beneath the landfill, the average oxidation–reduction potential (ORP) is −95 mV (Eh + 105 mV), and these reducing conditions were primarily induced by landfill leachate. Remediation efforts at this site have included installation of a low permeability clay cap; groundwater extraction, oxidation, and re-injection; and subsurface oxidation by injection of magnesium peroxide. The natural source of arsenic within the aquifer solids, coupled with widespread reducing conditions, has severely limited the effectiveness of these interventions on groundwater arsenic concentrations.     

Highstands during Marine Isotope Stage 5: evidence from the Ironshore Formation of Grand Cayman,British West Indies

The Ironshore Formation on Grand Cayman is formed of six unconformity-bounded packages (units A–F). Units A, B, C, and D, known from the subsurface in the northeastern part of Grand Cayman, formed during Marine Isotope Stages (MIS) 11(?), 9, 7, and 5e, respectively. Unconformities at the tops of units A, B, and C are highlighted by terra rossa and/or calcrete layers. Strata in core obtained from wells drilled in George Town Harbour and exposed on the west part of Grand Cayman belong to unit D, and the newly defined units E and F. Corals from unit E yielded Th/U ages of ∼104 ka whereas conch shells from unit F gave ages of ∼84 ka. Unit E equates to MIS 5c whereas unit F developed during MIS 5a.Th/U dating of corals and conchs from the Ironshore Formation on the western part of Grand Cayman shows that unit D formed during the MIS 5e highstand whereas units E and F developed in association with highstands at 95–110 ka (MIS 5c) and 73–87 ka (MIS 5a). Unit E, ∼5 m thick in the offshore cores, is poorly represented in onshore exposures. Unit F, which unconformably overlies unit D at most localities, is formed largely of fossil-poor, cross-bedded ooid grainstones. The unconformity at the top of unit D, a marine erosional surface with up to 2.5 m relief, is not characterized by terra rossa or calcrete in the offshore cores or onshore exposures. Unit D formed with a highstand of +6 m asl, whereas units E and F developed when sea level was +2 to +5 asl and +3 to +6 m asl, respectively. Thus, the highstands associated with MIS 5e, 5c, and 5a were at similar elevations.     

On the origin of saline fluids in the KTB (Continental Deep Drilling Project of Germany)

Highly saline fluids were encountered during the German Continental Deep Drilling Project (KTB) from depths ranging between 2 and 3 km to about 9 km. The most reliable data were obtained from samples extracted during a long-term pumping test in the 4000-m deep KTB pilot hole. Some 460 m³ Ca–Na–Cl brines with about 68 g l⁻¹ total dissolved solids (TDS) and some 270 m³ associated gases, mainly N₂ and CH₄ were pumped to the surface from the main fracture system situated near the bottom of the pilot hole. Geochemical and isotopic data support the hydraulic tests which suggest the presence of an open and large fluid reservoir at depth. The pumped fluids from this main fracture system were released from a deep reservoir situated at more than 5500 m depth which is hydraulically connected with the 9101 m deep KTB main hole, drilled some 250 m to the northeast of the pilot hole.While Ca and Sr contents of the extracted brines may be the result of water–rock interaction, Cl is most likely of external origin. The Cl is hypothesized to derive from geotectonic processes rather than to descending infiltration of paleo-seawater (evaporitic brines). The sampled fluids have probably migrated from a deeper reservoir to their present position since the Cretaceous–Tertiary period due to tectonic activity. However, several isotopic studies have identified an admixture of descending paleowaters down to more than 4000 m depth. The high ³⁶Cl/Cl ratio of the fluids sampled during the long-term pumping test point to a host rock highly enriched in U–Th, unlike the sampled KTB country rocks. The fluid reservoir is believed to be in contact with the Falkenberg granite massif situated about 2 km to the E of the KTB holes, capable of supplying sufficient neutron flux for considerable subsurface production of ³⁶Cl. The Na–Cl–(K-, SO₄) precursor fluids of the Ca–Na–Cl brines were produced in the course of extensive tectonic processes since the Late Caledonian within the Bohemian Massif.     

Depositional sequences in a foreland basin (north-western domain of the continental Duero basin,Spain)

The Cenozoic record of the north-western domain of the Duero basin is articulated at the surface through a set of continental depositional sequences called, from base to top, the Vegaquemada sequence, the Candanedo sequence, and the Barrillos sequence. These depositional sequences were deposited in continental sedimentary environments. The deposition of the first sequence occurred through a fluvial system with floodplains cut by low-sinuosity channels. The Vegaquemada sequence was developed between the Middle Eocene and the Early Agenian. The second sequence was formed by a set of highly efficient transport alluvial fans that evolved laterally towards fluvial systems with low-sinuosity fluvial channels and an extensive floodplain, where several types of palaeosols were formed. This sequence developed between the Early Agenian and the Late Vallesian. The third unit–the Barrillos sequence (between the Late Vallesian and the Turolian/Ruscinian transition), was generated by a set of highly efficient transport alluvial fans dominated by low-sinuosity fluvial channels.In subsurface geology, seismic and well data are used to rebuild the stratigraphic architecture. The two basal depositional sequences can be identified with two seismic units: the Palaeogene Seismic Unit (PgSU) and the Neogene Seismic Unit (NgSU), respectively. In the present work, we obtained the isovelocity, isochron, and isobath maps for the top and base of the two Cenozoic units. The Palaeozoic (PzSU) and Mesozoic (MzSU) seismic units are found under these two units. Through study of the logs of the various boreholes, it was only possible to analyse the upper 700 m of the Candanedo Sequence (NgSU), without encompassing the total thickness of the unit. Several middle-order sequences were differentiated, in general showing a sequential fining-upwards evolutionary character. Additionally, for the boreholes analysed two main types of electrofacies were identified, both representing fluvial channels and floodplain deposits.The north-western domain of the Duero basin is interpreted to have been formed in response to the tectonic uplifting of the Cantabrian Mountains since Middle-Eocene times. Integration of the data concerning the surface and subsurface geology in this domain reveals that this basin edge behaved as a foreland basin during Cenozoic stages. The foredeep, with a depth of 2800 m, is oriented east–west and has a sediment thickness of up to 3500 m. The forebulge is located in the southwestern zone and represents an area of basement uplifting in which a minimum thickness of materials from the Cenozoic depositional sequences has accumulated.     

Characterizing the vadose zone and a perched aquifer near the Vosges ridge at the La Soutte experimental site,Obernai, France

We consider a series of hydrogeophysical techniques that provide a multiscale investigation of the water content in the vadose zone and of the perched aquifer at the experimental site of “La Soutte” in the Vosges Mountains (France). It is located in a catchment area where several springs and streams occur along fractured volcanic and weathered plutonic rocks. The site is the object of a long-term study that uses both continuous and repeated measurements to monitor hydrogeological processes. The main results from AMT and DC resistivity techniques allow the determination of a high-resolution 3D resistivity model over a large range of depths (from 10⁰ to 10³ m). We discuss their use and propose a hydrogeological model (porosity, water conductivity and water content). We also use MRS and GPR for a detailed investigation of the shallow part of the catchment that consists of soil and weathered rocks of highly varying thickness (0 to 15 m). MRS is used to map the thickness and total water volume content by unit surface of the saturated weathered zone. It also yields estimates of the vadose zone thickness through the depth to the top of the saturated zone. Moreover, we show results from GPR CMP measurements that yield estimates of the water content and porosity in the shallowest layer (0–30 cm) by simple interpretation of the ground direct wave.     

Seismogenesis of clustered seismicity beneath the Kangra–Chamba sector of northwest Himalaya: Constraints from 3D local earthquake tomography

To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km² centered NE to the epicenter of Kangra earthquake.     

The insoluble organic matter in carbonaceous chondrites: Chemical structure,isotopic composition and origin

Carbonaceous chondrites are characterized by their enrichment in organic matter, mainly represented by insoluble organic matter (IOM), which consists of small aromatic units linked by short-branched aliphatic chains. Furthermore, IOM contains organic radicals heterogeneously distributed along with diradicaloids. These chemical features discriminate IOM from terrestrial counterparts. Isotopic compositions, especially the D/H isotopic ratio, are also distinct. IOM is highly enriched in D (D/H > 350 × 10⁻⁶), and the D/H isotopic ratio is heterogeneous. The isotopic composition is the result of interstellar-like processes that could have taken place during the first ages of the protosolar nebula. Chemical structure and isotopic composition clearly show that IOM is synthesized by an abiotic process and is subsequently affected by aqueous alteration or high-temperature metamorphism on the parent body.     

The relationship between fluid flow and mineral weathering in heterogeneous unsaturated porous media: A physical and geochemical characterization of a waste-rock pile

The relationships between factors that control subsurface flow and the timing, duration, and intensity of acidity generation and leaching of metals from waste-rock dumps are investigated. A 12 m high waste-rock pile that had been constructed in 1994 at Key Lake, Saskatchewan, Canada was disassembled, sampled and characterized in 2000. Physical properties that control water flow were characterized by measuring soil–water suction, volumetric water content, and the grain-size distribution at 60 randomized sites within the pile. Grain-size distribution was also measured at an additional 20 grid locations within the pile. Paste pH, pore-water geochemistry, mineralogy, and water-soluble extractions were used to investigate geochemical processes and sulfide oxidation at each of the 20 grid locations. A field-based soil–water characteristic curve could not be developed from the spatially variable and hysteretic field data; consequently, the grain-size distribution was used as a relative measure of subsurface flow and of the tendency to contain water under unsaturated conditions. The geochemical characterization demonstrated that marcasite underwent preferential weathering relative to pyrite and chalcopyrite, that dolomite was the main buffering carbonate mineral, and that gypsum, jarosite, and Fe oxyhydroxides were the main secondary (supergene) minerals. The pore waters contained up to 78,000 mg L⁻¹ SO₄, 690 mg L⁻¹ Ni and 1400 mg L⁻¹ U (800, 11.7 and 6 mM, respectively), suggesting that significant weathering had occurred. The pore water chemistry varied considerably between sampling sites. However, neither a correlation of pore-water chemistry with grain-size distribution nor a spatial relationship within the sampled grid was discernible.     

Groundwater fluxes into a submerged sinkhole area,Central Italy,using radon and water chemistry

The groundwater contribution into Green Lake and Black Lake (Vescovo Lakes Group), two cover collapse sinkholes in Pontina Plain (Central Italy), was estimated using water chemistry and a ²²²Rn budget. These data can constrain the interactions between sinkholes and deep seated fluid circulation, with a special focus on the possibility of the bedrock karst aquifer feeding the lake. The Rn budget accounted for all quantifiable surface and subsurface input and output fluxes including the flux across the sediment–water interface. The total value of groundwater discharge into Green Lake and Black Lake (∼540 ± 160 L s⁻¹) obtained from the Rn budget is lower than, but comparable with historical data on the springs group discharge estimated in the same period of the year (800 ± 90 L s⁻¹). Besides being an indirect test for the reliability of the Rn-budget “tool”, it confirms that both Green and Black Lake are effectively springs and not simply “water filled” sinkholes. New data on the water chemistry and the groundwater fluxes into the sinkhole area of Vescovo Lakes allows the assessment of the mechanism responsible for sinkhole formation in Pontina Plain and suggests the necessity of monitoring the changes of physical and chemical parameters of groundwater below the plain in order to mitigate the associated risk.     

Kinematics of the crust around the Tanggula Shan in North–Central Tibet: Constraints from paleomagnetic data

We have conducted a paleomagnetic investigation on the Middle–Upper Jurassic marine strata exposed in the hanging wall of the Tanggula Thrust system near the Yanshiping area, northern Tibet. Progressive demagnetization experiments successfully isolated stable magnetization over a broad spectrum of demagnetization temperatures. The mean direction of the characteristic remanent magnetizations for the Middle–Late Jurassic Yanshiping Group in stratigraphic coordinates (D/I (Declination/Inclination) = 5.6°/60.3°, k = 22.9, α95 = 12.9°, N = 7 s) is much more clustered than the mean direction in geographic coordinates (D/I = 345.5°/37.2°, k = 2.5, α95 = 48.4°), indicating magnetization was not acquired after folding. Although the conventional fold test is positive, incremental untilting test on the characteristic remanent magnetization reveals that a maximum value of precision parameter k occurs at 82.1 ± 4.6% untilting (D/I = 3.3°/57.8°, k = 43.9, α95 = 9.2°), which indicates the ChRMs are probably acquired during Late Cretaceous folding. This synfolding magnetization component is therefore secondary. The corresponding pole position (84.4°N, 119.4°E with dp/dm = 13.5/9.9°) is inconsistent with Jurassic–Early Cretaceous paleopoles of the region, but the paleolatitude is consistent with the Late Cretaceous paleolatitude observed in the Qiangtang terrane and its periphery. The synfolding component is carried by both magnetite and hematite, which were identified by isothermal remnant magnetization acquisition experiments, unblocking temperatures of stable magnetic components, and Curie temperature determination and correlated with observed hydrothermal veins. Available geological evidences indicate that the synfolding magnetization is probably the result of chemical remagnetization caused by orogenic fluids or hydrothermal sources during the early uplift of the Tibetan Plateau.     

Dissolution kinetics of selected natural minerals relevant to potential CO2-injection sites – Part 2: Dissolution and alteration of carbonates and feldspars in CO2-bearing brines

Chemical interaction processes among injected CO₂, saline fluids and potential reservoir materials are experimentally simulated to derive dissolution rates of natural materials (minerals) that can be used as input parameters for modeling of CO₂ storage in deep saline formations and risk analyses. In order to study dissolution processes, mineral aliquots were exposed to CO₂-bearing brines at elevated temperature (60, 100, 150 °C) and pressure (85 bar) and at various run durations. Several potential reservoir rocks include carbonates as cement. Calcite and dolomite grains were therefore mainly used as solid starting material. Experiments with the two feldspar varieties alkali feldspar and almost pure anorthite were performed in addition. Grain sizes of the mineral starting materials varied between <63 μm and 500 μm with most experiments performed at grain size fractions of 160 – 250 μm and 250 – 500 μm. All experiments run with a complex synthetic brine (total dissolved solids: ∼156 g/l) according to a natural upper cretaceous formation water. Dry ice was used as CO₂-source. All experiments were done in closed batch reactors. These reactors allow mimicking reservoir conditions far from the injection site as well as reservoir conditions after finishing the CO₂ injection. The concentration changes during the experiment were monitored by ICP-OES measurements of the initial and the post-run fluids. Dissolution rates were derived based on the concentration changes of the brine.Most of the studied experimental variables and parameters (temperature, run duration, grain size, brine composition – expressed as pH-value and ionic strength) impact alteration of the reacting agents, i.e. they change the chemical composition of the brine, change the surfaces of the mineral aliquots exposed to the CO₂-bearing brine, and induce formation of secondary minerals. Hence, all influencing parameters on dissolution processes have to be considered and time-resolved changes of the dissolution behavior have to be implemented in numerical simulations of processes at CO₂ injection sites and CO₂ storage reservoirs.     

Method for deep temperature estimation with regard to the paleoclimate influence on heat flow

New results of the study of the distorting effect of paleoclimate on the formation of the thermal field of the Earth’s uppermost crust inthe South Urals are presented. We discuss some consequences of taking into account the paleoclimate influence on estimates of deep heat flow and possible applications of these data. The dependence of the heat flow density on the studied depth of boreholes is considered, and deep heat loss and deep temperatures are estimated. A new method for deep temperature estimation is proposed, which takes into account the paleoclimate influence on the heat flow. The method is tested on the available high-quality temperature data for deep boreholes. Deep temperatures to a depth of –10,000 m are estimated for the platform part of the Republic of Bashkortostan using the proposed method. Isotherm schemes are constructed for elevations of –5000 and –10,000 m below sea level.The necessity of using heat flow values corrected for the paleoclimate influence to estimate the Earth’s thermal state is justified. Some examples illustrate that underestimation of heat flow values measured in shallow boreholes might lead to underestimation of deep temperatures and global heat losses.     

Mobilisation of Iron from rocks in a fractured aquifer: Lithological and geochemical controls

Iron mobilisation from aquifer rocks in an important fractured aquifer system in South Africa is resulting in clogging of boreholes by Fe oxide minerals. Leach experiments using natural waters were conducted to determine the effects of redox conditions, pH lithology and presence of organic acids on the rate and extent of Fe dissolution from aquifer rocks, with the aim of clarifying the association of Fe clogging with geological formations that show Fe staining on weathering. The results indicate that the greatest amount of Fe (>30 mmol/kg rock) is leached from arenaceous rocks with low total Fe contents (49.0–75.0 mmol/kg) under anoxic conditions. Rocks with the highest Fe contents (>800 mmol/kg) generated low concentrations of Fe (<10 mmol/kg) even under favourable conditions of 0 mg/L DO and pH 3. The extent of Fe dissolution from the rocks was found to be most strongly dependent on the redox conditions, and the form of Fe present in the rock, with ascorbate-extracted amorphous Fe being the most mobile. The rate of dissolution is affected by pH and the presence of natural organic acids in the leachate. However, the effect of organic acids was only noticeable on arenaceous rocks.