Numerical model of pore‐pressure diffusion associated with the initiation of the 2010–2011 Guy–Greenbrier,Arkansas earthquakes

We model pore‐pressure diffusion caused by pressurized waste‐fluid injection at two nearby wells and then compare the buildup of pressure with the observed initiation and migration of earthquakes during the early part of the 2010–2011 Guy–Greenbrier earthquake swarm. Pore‐pressure diffusion is calculated using MODFLOW 2005 that allows the actual injection histories (volume/day) at the two wells to diffuse through a fractured and faulted 3D aquifer system representing the eastern Arkoma basin. The aquifer system is calibrated using the observed water‐level recovery following well shut‐in at three wells. We estimate that the hydraulic conductivities of the Boone Formation and Arbuckle Group are 2.2 × 10^?2 and 2.03 × 10^?3 m day^?1, respectively, with a hydraulic conductivity of 1.92 × 10^?2 m day^?1 in the Hunton Group when considering 1.72 × 10^?3 m day^?1 in the Chattanooga Shale. Based on the simulated pressure field, injection near the relatively conductive Enders and Guy–Greenbrier faults (that hydraulically connect the Arbuckle Group with the underlying basement) permits pressure diffusion into the crystalline basement, but the effective radius of influence is limited in depth by the vertical anisotropy of the hydraulic diffusivity. Comparing spatial/temporal changes in the simulated pore‐pressure field to the observed seismicity suggests that minimum pore‐pressure changes of approximately 0.009 and 0.035 MPa are sufficient to initiate seismic activity within the basement and sedimentary sections of the Guy–Greenbrier fault, respectively. Further, the migration of a second front of seismicity appears to follow the approximately 0.012 MPa and 0.055 MPa pore‐pressure fronts within the basement and sedimentary sections, respectively.     

Effective gas permeability of Tight Gas Sandstones as a function of capillary pressure – a non‐steady‐state approach

Single‐ and two‐phase (gas/water) fluid transport in tight sandstones has been studied in a series of permeability tests on core plugs of nine tight sandstones of the southern North Sea. Absolute (Klinkenberg‐corrected) gas permeability coefficients (k_gas__inf) ranged between 3.8 × 10^?16 and 6.2 × 10^?19 m² and decreased with increasing confining pressure (10–30 MPa) by a factor 3–5. Klinkenberg‐corrected (intrinsic) gas permeability coefficients were consistently higher by factors from 1.4 to 10 than permeability coefficients determined with water. Non‐steady‐state two‐phase (He/water) flow experiments conducted up to differential pressures of 10 MPa document the dynamically changing conductivity for the gas phase, which is primarily capillary‐controlled (drainage and imbibition). Effective gas permeability coefficients in the two‐phase flow tests ranged between 1.1 × 10^?17 and 2.5 × 10^?22 m², corresponding to relative gas permeabilities of 0.03% and 10%. In the early phase of the nonstationary flow regime (before establishment of steady‐state conditions), they may be substantially (>50%) lower. Effective gas permeability measurements are affected by the following factors: (i) Capillary‐controlled drainage/imbibition, (ii) viscous–dynamic effects (iii) and slip flow.     

Regional groundwater flow and interactions with deep fluids in western Apennine: the case of Narni‐Amelia chain (Central Italy)

The elemental fluxes and heat flow associated with large aquifer systems can be significant both at local and at regional scales. In fact, large amounts of heat transported by regional groundwater flow can affect the subsurface thermal regime, and the amount of matter discharged towards the surface by large spring systems can be significant relative to the elemental fluxes of surface waters. The Narni‐Amelia regional aquifer system (Central Italy) discharges more than 13 m³ sec^?1 of groundwater characterised by a slight thermal anomaly, high salinity and high pCO₂. During circulation in the regional aquifer, groundwater reacts with the host rocks (dolostones, limestones and evaporites) and mixes with deep CO₂‐rich fluids of mantle origin. These processes transfer large amounts of dissolved substances, in particular carbon dioxide, and a considerable amount of heat towards the surface. Because practically all the water circulating in the Narni‐Amelia system is discharged by few large springs (Stifone‐Montoro), the mass and energy balance of these springs can give a good estimation of the mass and heat transported from the entire system towards the surface. By means of a detailed mass and balance of the aquifer and considering the soil CO₂ fluxes measured from the main gas emission of the region, we computed a total CO₂ discharge of about 7.8 × 10⁹ mol a^?1 for the whole Narni‐Amelia system. Finally, considering the enthalpy difference between infiltrating water and water discharged by the springs, we computed an advective heat transfer related to groundwater flow of 410 ± 50 MW.     

Residence times and source ages of deep crustal fluids: interpretation of 129I and 36Cl results from the KTB‐VB drill site,Germany

We present here ¹²⁹I/I and ³⁶Cl/Cl ratios, together with halogen concentrations in crustal fluids from the continental deep drill site (KTB‐VB) in Germany, where fluids were collected from 4000 m depth during a pump test carried out in 2002 and 2003. Compared with seawater, the fluids are enriched by factors of 2, 8 and 40 for Cl^?, Br and I, respectively, and show little variation over the test period. The ¹²⁹I/I ratios are between 1700 and 4100 × 10^?15; the ³⁶Cl/Cl ratios are below 10 × 10^?15. Co‐variation between ¹²⁹I and ³⁶Cl concentrations in the fluids indicates that anthropogenic components are absent and that the ratios reflect an addition from crustal sources. The results suggest residence times of 10 Ma or more for the fluids in formations with uranium concentrations of 1 ppm. A minimum age of 30 Ma for the iodine source was derived from the correlation between ¹²⁹I and ³⁶Cl concentrations in the fluids. The results demonstrate that the halogen characteristics of the KTB fluids are very similar to those of other deep crustal fluids and that the combination of ¹²⁹I and ³⁶Cl systematics allows determination of residence times and source ages of such fluids.     

The upper continental crust, an aquifer and its fluid: hydaulic and chemical data from 4 km depth in fractured crystalline basement rocks at the KTB test site

Detailed information on the hydrogeologic and hydraulic properties of the deeper parts of the upper continental crust is scarce. The pilot hole of the deep research drillhole (KTB) in crystalline basement of central Germany provided access to the crust for an exceptional pumping experiment of 1‐year duration. The hydraulic properties of fractured crystalline rocks at 4 km depth were derived from the well test and a total of 23100 m³ of saline fluid was pumped from the crustal reservoir. The experiment shows that the water‐saturated fracture pore space of the brittle upper crust is highly connected, hence, the continental upper crust is an aquifer. The pressure–time data from the well tests showed three distinct flow periods: the first period relates to wellbore storage and skin effects, the second flow period shows the typical characteristics of the homogeneous isotropic basement rock aquifer and the third flow period relates to the influence of a distant hydraulic border, probably an effect of the Franconian lineament, a steep dipping major thrust fault known from surface geology. The data analysis provided a transmissivity of the pumped aquifer T = 6.1 × 10^?6 m² sec^?1, the corresponding hydraulic conductivity (permeability) is K = 4.07 × 10^?8 m sec^?1 and the computed storage coefficient (storativity) of the aquifer of about S = 5 × 10^?6. This unexpected high permeability of the continental upper crust is well within the conditions of possible advective flow. The average flow porosity of the fractured basement aquifer is 0.6–0.7% and this range can be taken as a representative and characteristic values for the continental upper crust in general. The chemical composition of the pumped fluid was nearly constant during the 1‐year test. The total of dissolved solids amounts to 62 g l^?1 and comprise mainly a mixture of CaCl₂ and NaCl; all other dissolved components amount to about 2 g l^?1. The cation proportions of the fluid (X_Ca approximately 0.6) reflects the mineralogical composition of the reservoir rock and the high salinity results from desiccation (H₂O‐loss) due to the formation of abundant hydrate minerals during water–rock interaction. The constant fluid composition suggests that the fluid has been pumped from a rather homogeneous reservoir lithology dominated by metagabbros and amphibolites containing abundant Ca‐rich plagioclase.     

Impacts of Pleistocene glaciation on large‐scale groundwater flow and salinity in the Michigan Basin

Pleistocene melting of kilometer‐thick continental ice sheets significantly impacted regional‐scale groundwater flow in the low‐lying stable interiors of the North American and Eurasian cratons. Glacial meltwaters penetrated hundreds of meters into the underlying sedimentary basins and fractured crystalline bedrock, disrupting relatively stagnant saline fluids and creating a strong disequilibrium pattern in fluid salinity. To constrain the impact of continental glaciation on variable density fluid flow, heat and solute transport in the Michigan Basin, we constructed a transient two‐dimensional finite‐element model of the northern half of the basin and imposed modern versus Pleistocene recharge conditions. The sag‐type basin contains up to approximately 5 km of Paleozoic strata (carbonates, siliciclastics, and bedded evaporites) overlain by a thick veneer (up to 300 m) of glacial deposits. Formation water salinity increases exponentially from <0.5 g l^?1 total dissolved solids (TDS) near the surface to >350 g l^?1 TDS at over 800 m depth. Model simulations show that modern groundwater flow is primarily restricted to shallow glacial drift aquifers with discharge to the Great Lakes. During the Pleistocene, however, high hydraulic heads from melting of the Laurentide Ice Sheet reversed regional flow patterns and focused recharge into Paleozoic carbonate and siliciclastic aquifers. Dilute waters (<20 g l^?1 TDS) migrated approximately 110 km laterally into the Devonian carbonate aquifers, significantly depressing the freshwater‐saline water mixing zones. These results are consistent with ¹⁴C ages and oxygen isotope values of confined groundwaters in Devonian carbonates along the basin margin, which reflect past recharge beneath the Laurentide Ice Sheet (14–50 ka). Constraining the paleohydrology of glaciated sedimentary basins, such as the Michigan Basin, is important for determining the source and residence times of groundwater resources, in addition to resolving geologic forces responsible for basinal‐scale fluid and solute migration.     

Analytical and numerical models of hydrothermal fluid flow at fault intersections

Fault intersections are the locus of hot spring activity and Carlin‐type gold mineralization within the Basin and Range, USA. Analytical and numerical solutions to Stokes equation suggest that peak fluid velocities at fault intersections increase between 20% and 47% when fracture apertures have identical widths but increase by only about 1% and 8% when aperture widths vary by a factor of 2. This suggests that fault zone intersections must have enlarged apertures. Three‐dimensional finite element models that consider intersecting 10‐ to 20‐m wide fault planes resulted in hot spring activity being preferentially located at fault zone intersections when fault zones were assigned identical permeabilities. We found that the onset of convection at the intersections of the fault zones occurred in our hydrothermal model over a narrow permeability range between 5 × 10^?13 and 7 × 10^?13 m². Relatively high vertical fluid velocities (0.3–3 m year^?1) extended away from the fault intersections for about 0.5–1.5 km. For the boundary conditions and fault plane dimensions used, peak discharge temperatures of 112°C at the water table occurred with an intermediate fault zone permeability of 5 × 10^?13 m². When fault plane permeability differed by a factor of 2 or more, the locus of hot spring activity shifted away from the intersections. However, increasing the permeability at the core of the fault plane intersection by 40% shifted the discharge back to the intersections. When aquifer units were assigned a permeability value equal to those of the fault planes, convective rolls developed that extend about 3 km laterally along the fault plane and into the adjacent aquifer.     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

Molecular transport of methane, ethane and nitrogen and the influence of diffusion on the chemical and isotopic composition of natural gas accumulations

Laboratory experiments have been performed to determine diffusion coefficients of natural gas components (methane, ethane and nitrogen) and isotope fractionation effects under simulated in situ pressure (up to 45 MPa effective stress) and temperature conditions (50–200°C) in water‐saturated pelitic and coarse‐grained rocks. Effective diffusion coefficients of molecular nitrogen (0.39 × 10^?11 to 21.6 × 10^?11 m² sec^?1 at 90°C) are higher than those for methane (0.18 × 10^?11 to 18.2 × 10^?11 m² sec^?1 at 90°C). Diffusive flux rates expressed in mass units are generally higher for N₂ than for CH₄. Both methane and (to a lesser extent) nitrogen diffusion coefficients decrease with increasing total organic carbon (TOC) content of the rock samples because of sorption processes on the organic matter. This effect decreases with increasing temperature. Effective diffusion coefficients increase upon a temperature increase from 50 to 200°C by a factor of four. Effective diffusion coefficients and steady‐state diffusive flux decrease with effective stress. Stationary diffusive fluxes drop by 50–70% for methane and 45–62% for nitrogen while effective diffusion coefficients are reduced by 38% (CH₄) and 32–48% (N₂), respectively. Isotope fractionation coefficients of diffusive transport are higher for methane (?1.56 and ?2.77‰) than for ethane (?0.84 and ?1.62‰). Application of the experimental results to geological systems show that diffusive transport has only a low transport efficiency. Significant depletion of natural gas reservoirs by molecular diffusion is only expected in cases of very poor caprock qualities (in terms of thickness and/or porosity) and over extended periods of geological time. Under these circumstances, the chemical and isotopic composition of a gas reservoir will change and maturity estimates based on these parameters may be deceptive. To account for these potential effects, nomograms have been developed to estimate diffusive losses and apply maturity corrections.     

Estimation of Soil Erosion Using Field and Modelling Approaches in an Undisturbed Arnhem Land Catchment,Northern Territory,Australia

An assessment of slope erosion at Tin Camp Creek catchment, Arnhem Land, Northern Territory, Australia, was carried out using the fallout environmental radioisotope caesium‐137 (¹³⁷Cs) as an indicator of soil erosion status, two numerical models (SIBERIA and the Revised Universal Soil Loss Equation (RUSLE)) and erosion pins. This undisturbed drainage basin is situated in the seasonally wet‐dry tropics, with high energy storms and a mean annual rainfall of 1389 mm. Tin Camp Creek catchment is unaffected by European agriculture or pastoral activities, but often experiences fire during the dry season. Two transects were sampled for ¹³⁷Cs in 2002 and 2004, and two models were used to convert ¹³⁷Cs measurements into soil loss estimates. Two methods using the theoretical Profile Distribution Model (PDM) gave net soil redistribution rates between +2.72 and –22.19 t ha^?1 yr^?1 and +2.95 and –24.06 t ha^?1 yr^?1, respectively, while an Australian empirical model (AEM) for uncultivated soils produced estimates between +1.84 and –7.00 t ha^?1 yr^?1 (negative values indicate soil erosion, positive, deposition). The RUSLE gave estimated soil losses for the two transects of approximately 10 t ha^?1 yr^?1, while the SIBERIA model produced values between 0.5 and 2 t ha^?1 yr^?1 for the transects and between 3.5 and 11 t ha^?1 yr^?1 for the total catchment. Average net soil losses of 14 and 15 t ha^?1 yr^?1 for the total catchment and slopes, respectively, were measured by erosion pins. The soil losses in the catchment are similar to those for some other transects in the Northern Territory and the Kimberley region of Western Australia (measured by the ¹³⁷Cs AEM), even though these areas are affected by pastoral activities. This may be at least partly explained by erosion in Tin Camp Creek catchment during high intensity rainstorms at the commencement of the wet season, especially if the slopes have been affected by fire during the previous dry season.     

An Examination of Vertical Accretion of Floodplain Sediments in the Labasa River Sugarcane Belt of Northern Fiji: Rates,Influences and Contributing Processes

The ¹³⁷Cs (caesium‐137) method was used to investigate medium‐term rates of sediment deposition on the floodplain of the Labasa River, on Vanua Levu island in northern Fiji. The Labasa basin is commercially important for sugarcane farming, which provides much revenue and sustains the economy in the greater Labasa area. Alluvium was sampled at three riverbank sites in vertical increments of 30 mm. Measured net vertical accretion rates, based on analysis of depth‐profiles of ¹³⁷Cs activity, ranged from 10 mm yr^?1 at a low‐lying site near a tributary confluence with the main river, to 60 mm yr^?1 on a levée that was elevated slightly above the adjacent floodplain. These rates of accretion are high, but in general agreement with rates recorded using similar methods in other tropical Pacific island river systems. Soil erosion under sugarcane on hilly parts of the lower Labasa basin is an important factor in rapid floodplain development. Observations made during Cyclone Ami, which traversed Vanua Levu island in early 2003, illustrate the major contribution of severe storm events to floodplain sediment supply, by triggering numerous landslides on catchment slopes and in‐channel debris floods, and by generating large‐magnitude overbank floods.     

Geochemical fingerprinting of hydraulic fracturing fluids from Qusaiba Hot Shale and formation water from Paleozoic petroleum systems,Saudi Arabia

Provenance studies of produced water are essential to trace flow dynamics and reservoir compartmentalization in petroleum systems and to quantify fluid recovery rates from unconventional fracturing. Produced water from a hydraulically fractured well in the Qusaiba Hot Shale in the Northern Exploration Area, Saudi Arabia, was daily monitored and analyzed for water chemistry, and environmental (δ²H, δ¹³C, δ¹⁸O_H2O, δ¹⁸O_SO4, δ³⁴S_SO4, δ³⁷Cl, ⁸⁷Sr/⁸⁶Sr) and cosmogenic isotopes (³H, ¹⁴C, ³⁶Cl), to differentiate from reference fluids of supply water, fracturing fluids, and formation water from adjacent Paleozoic units. Initially, recovered water is composed of fracturing fluids and subsequently replaced by a homogeneous cut of pristine formation water. Formation water is composed of dominant meteoric water (approximately 84 vol%) and minor fossil evaporated seawater. The young ¹⁴C‐apparent age between 6000 and 6700 years BP and depleted δ¹⁸O/δ²H values for the meteoric component confirm the infiltration of surface water into the Qusaiba Hot Shale interval or adjacent units during the Early Holocene Pluvial Period under cooler and wetter climatic conditions than present, which suggest the presence of a very recent, dynamic hydraulic flow system. ³⁶Cl/Cl ratios between 102 × 10^?15 and 31 × 10^?15 are ambiguous and can be attributed to atmospheric recharge close to the coast, mixing of ³⁶Cl‐enriched Quaternary meteoric recharge with ³⁶Cl‐depleted fossil seawater, and/or hypogene production by U‐Th‐enriched host rock. Produced waters from Qusaiba Hot Shale are within the compositional range of Na‐Cl‐type formation water from Paleozoic reservoir units in northern Saudi Arabia with salinities from 30 000 to 130 000 mg l^?1. As a novel technological approach for exploration wells in Northern Saudi Arabia, multi‐isotopic methods were successfully implemented to quantify flowback volumes from hydraulic fracturing, and to fingerprint pristine formation water or pore water in Paleozoic systems on their provenance, residence time, migration pathways, and secondary alteration processes.     

Hydraulic observations from a 1 year fluid production test in the 4000 m deep KTB pilot borehole

A long‐term pump test was conducted in the KTB pilot borehole (KTB‐VB), located in the Oberpfalz area, Germany. It produced 22 300 m³ of formation fluid. Initially, fluid production rate was 29 l min^?1 for 4 months, but was then raised to an average of 57 l min^?1 for eight more months. The aim of this study was to examine the fluid parameters and hydraulic properties of fractured, crystalline crusts as part of the new KTB programme ‘Energy and Fluid Transport in Continental Fault Systems’. KTB‐VB has an open‐hole section from 3850 to 4000 m depth that is in hydraulic contact with a prominent continental fault system in the area, called SE2. Salinity and temperature of the fluid inside the borehole, and consequently hydrostatic pressure, changed significantly throughout the test. Influence of these quantities on variations in fluid density had to be taken into account for interpretation of the pump test. Modelling of the pressure response related to the pumping was achieved assuming the validity of linear Darcy flow and permeability to be independent of the flow rate. Following the principle ‘minimum in model dimension’, we first examined whether the pressure response can be explained by an equivalent model where rock properties around the borehole are axially symmetric. Calculations show that the observed pressure data in KTB‐VB can in fact be reproduced through such a configuration. For the period of high pumping rate (57 l min^?1) and the following recovery phase, the resulting parameters are 2.4 × 10^?13 m³ in hydraulic transmissivity and 3.7 × 10^?9 m Pa^?1 in storativity for radial distances up to 187 m, and 4.7 × 10^?14 m³ and 6.0 × 10^?9 m Pa^?1, respectively, for radial distances between 187 and 1200 m. The former pair of values mainly reflect the hydraulic properties of the fault zone SE2. For a more realistic hydraulic study on a greater scale, program FEFLOW was used. Parameter values were obtained by matching the calculated induced pressure signal to fluid‐level variations observed in the KTB main hole (KTB‐HB) located at 200 m radial distance from KTB‐VB. KTB‐HB is uncased from 9031 to 9100 m and shows indications of leakage in the casing at depths 5200–5600 m. Analysis of the pressure record and hydraulic modelling suggest the existence of a weak hydraulic communication between the two boreholes, probably at depths around the leakage. Hydraulic modelling of a major slug‐test in KTB‐HB that was run during the pumping in KTB‐VB reveals the effective transmissivity of the connected formation to be 1 to 2 orders of magnitude lower than the one determined for the SE2 fault zone.     

Modelling the Lost City hydrothermal field: influence of topography and permeability structure

The Lost City hydrothermal field (LCHF) is hosted in serpentinite at the crest of the Atlantis Massif, an oceanic core complex close to the mid‐Atlantic Ridge. It is remarkable for its longevity and for venting low‐temperature (40–91°C) alkaline fluids rich in hydrogen and methane. IODP Hole U1309D, 5 km north of the LCHF, penetrated 1415 m of gabbroic rocks and contains a near‐conductive thermal gradient close to 100°C km^?1. This is remarkable so close to an active hydrothermal field. We present hydrothermal modelling using a topographic profile through the vent field and IODP site U1309. Long‐lived circulation with vent temperatures similar to the LCHF can be sustained at moderate permeabilities of 10^?14 to 10^?15 m² with a basal heatflow of 0.22 W m^?2. Seafloor topography is an important control, with vents tending to form and remain in higher topography. Models with a uniform permeability throughout the Massif cannot simultaneously maintain circulation at the LCHF and the near‐conductive gradient in the borehole, where permeabilities <10^?16 m² are required. A steeply dipping permeability discontinuity between the LCHF and the drill hole is required to stabilize venting at the summit of the massif by creating a lateral conductive boundary layer. The discontinuity needs to be close to the vent site, supporting previous inferences that high permeability is most likely produced by faulting related to the transform fault. Rapid increases in modelled fluid temperatures with depth beneath the vent agree with previous estimates of reaction temperature based on geochemical modelling.     

OSL and TL characteristics of fine grain quartz from Mongolian prehistoric pottery used for dating

The OSL, post-IR OSL and pulsed post-IR OSL applied to polymineral grains and calculated by fitting to the data the contributions from fast, medium and slow components revealed that the polymineral samples under study are dominated by the medium component. An increase in D _e ’s with increasing integration intervals was observed, which is considered as an indication of increasing medium and decayed fast component; and the equivalent doses obtained using different components or minerals reflect also the shape of the dose distributions. The identified fast component in polymineral sample has photoionization cross section of 1.2±0.02 × 10^?17 cm². The present study shows the usefulness of the application of different luminescence techniques combined with fitting procedures as a check which should be adopted in dating protocols. Based on luminescence ages obtained on poly-mineral grains from prehistoric pottery samples from the Boroo settlement, Mongolia, which are in agreement with independent age control by ¹⁴C on charcoal material, it is argued that the manufacturing of Xiongnu — pottery at this site lasted until ca. 130±75 AD.     

Modeling of hydrogen production by serpentinization in ultramafic‐hosted hydrothermal systems: application to the Rainbow field

The production of hydrogen by serpentinization in ultramafic‐hosted hydrothermal systems is simulated by coupling thermodynamic and dynamic modeling in the framework of a thermo‐hydraulic single‐pass model where a high‐temperature hydrothermal fluid moves preferentially through a main canal of high permeability. The alteration of ultramafic rocks is modeled with a first‐order kinetic formulation, wherein the serpentinization rate coefficient, K_r, takes the form: K_r = A exp(?α(T ? T₀)²). In this formulation, α determines the temperature range of the reaction and T₀ is the temperature at which the serpentinization rate reaches its maximum. This model is applied to the Rainbow hydrothermal system, which is situated on the Mid‐Atlantic Ridge, and characterized by a high temperature, a high mass flux, and a very high hydrogen concentration. The results show that a first‐order kinetic law gives a useful representation of the kinetics of serpentinization. The estimated value for the parameter A in the temperature‐dependent formulation of the serpentinization rate coefficient lies in the range (1–5) × 10^?11 s^?1. This effective parameter is several orders of magnitude lower than the values obtained from small grain‐size experiments, but in agreement with other published modeling studies of natural systems. Numerical simulations show that the venting site is able to produce the observed high concentration of hydrogen during the whole continuous lifetime of the Rainbow site.     

Hydrogeochemical dynamics affecting steam‐heated pools at El Chichón Crater (Chiapas – Mexico)

El Chichón is an active volcano located in the north‐western Chiapas, southern Mexico. The crater hosts a lake, a spring, named Soap Pool, emerging from the underlying volcanic aquifer and several mud pools/hot springs on the internal flanks of the crater which strongly interact with the current fumarolic system (steam‐heated pools). Some of these pools, the crater lake and a cold spring emerging from the 1982 pumice deposits, have been sampled and analysed. Water–volcanic gas interactions determine the heating (43–99°C) and acidification (pH 2–4) of the springs, mainly by H₂S oxidation. Significantly, in the study area, a significant NH₃ partial pressure has been also detected. Such a geochemically aggressive environment enhances alteration of the rock in situ and strongly increases the mineralization of the waters (and therefore their electrical conductivity). Two different mineralization systems were detected for the crater waters: the soap pool‐lake (Na⁺/Cl^? = 0.4, Na/Mg>10) and the crater mud pools (Na⁺/Cl^? > 10, Na/Mg < 4). A deep boiling, Na⁺‐K⁺‐Cl^?‐rich water reservoir generally influences the Soap Pool‐lake, while the mud pool is mainly dominated by water‐gas–rock interactions. In the latter case, conductivity of sampled water is directly proportional to the presence of reactive gases in solution. Therefore, chemical evolution proceeds through neutralization due to both rock alteration and bacterial oxidation of ammonium to nitrate. The chemical compositions show that El Chichón aqueous fluids, within the crater, interact with gases fed by a geothermal reservoir, without clear additions of deep magmatic fluids. This new geochemical dataset, together with previously published data, can be used as a base line with which to follow‐up the activity of this deadly volcano.     

A simulation of the hydrothermal response to the Chesapeake Bay bolide impact

Groundwater more saline than seawater has been discovered in the tsunami breccia of the Chesapeake Bay Impact Crater. One hypothesis for the origin of this brine is that it may be a liquid residual following steam separation in a hydrothermal system that evolved following the impact. Initial scoping calculations have demonstrated that it is feasible such a residual brine could have remained in the crater for the 35 million years since impact. Numerical simulations have been conducted using the code HYDROTHERM to test whether or not conditions were suitable in the millennia following the impact for the development of a steam phase in the hydrothermal system. Hydraulic and thermal parameters were estimated for the bedrock underlying the crater and the tsunami breccia that fills the crater. Simulations at three different breccia permeabilities suggest that the type of hydrothermal system that might have developed would have been very sensitive to the permeability. A relatively low breccia permeability (1 × 10^?16 m²) results in a system partitioned into a shallow water phase and a deeper superheated steam phase. A moderate breccia permeability (1 × 10^?15 m²) results in a system with regionally extensive multiphase conditions. A relatively high breccia permeability (1 × 10^?14 m²) results in a system dominated by warm‐water convection cells. The permeability of the crater breccia could have had any of these values at given depths and times during the hydrothermal system evolution as the sediments compacted. The simulations were not able to take into account transient permeability conditions, or equations of state that account for the salt content of seawater. Results suggest, however, that it is likely that steam conditions existed at some time in the system following impact, providing additional evidence that is consistent with a hydrothermal origin for the crater brine.     

A Survey of Soil Erosion in Australia using Caesium‐137

A survey of soil erosion was conducted in Australia using the fallout radioisotope caesium‐137 as an indicator of topsoil redistribution. Two hundred and six sites were sampled, 100 within rotational cropping and horticultural use, 52 within uncultivated permanent pasture and forest, and 54 in rangelands. Average net soil losses were approximately equal for cultivated cropping lands and rangelands (ca. 5.5 t ha^?1 yr^?1), and just over 1 t ha^?1 yr^?1 for pasture and forest. The Mann Whitney U Test revealed that losses under cropping and rangeland conditions were significantly higher (p < 0.05) than under uncultivated pasture and forest. Soil loss was negatively correlated with mean annual rainfall and slope gradient, and positively correlated with slope length (Spearman's rank correlation). There was no correlation between rates of soil loss and a rainfall erosi‐vity index. An assessment of erosional events was provided by landholders for 104 sites, with their ranking being weakly but significantly correlated with soil loss estimates (r =+0.35). Sixty percent of sites had net soil losses greater than 1 t ha^?1 yr^?1, and 74% of sites had losses of more than 0.5 t ha^?1 yr^?1. This latter rate may be regarded as a limit for a tolerable level of soil loss. These high rates of soil loss have occurred since the mid‐1950s despite there being significant landholder awareness of the soil erosion hazard.     

A simple model of gravity-wave momentum and energy fluxes transferred through the middle atmosphere to the upper atmosphere

Vertical fluxes of momentum and energy through the middle atmosphere are calculated by using a simple semi-empirical model of quasi-monochromatic internal gravity waves with dominant vertical wavenumbers. In this model those dominant gravity waves are assumed to saturate and break at each observational altitude by an effective critical-layer mechanism. The dominant value of the vertical wave-number is expressed by an exponential function of altitude, decreasing upward with a scale height of 34 km. This expression gives the momentum and energy flux densities decreasing upward with scale heights of 12 and 18 km, respectively, and typical values at 100 km altitude are estimated as 4 × 10⁻⁵ Pa and 4 × 10⁻³ W/m². A heat flux induced by wavebreaking turbulence also has an order of magnitude similar to that of the wave energy flux. Variabilities around these values and comparisons with other momentum and heat inputs to the upper atmosphere are only briefly discussed.