Potential of future seismogenesis in Hebei Province (NE China) due to stress interactions between strong earthquakes

Northeast China, a densely populated area, is affected by intense seismic activity, which includes large events that caused extensive disaster and tremendous loss of life. For contributing to the continuous efforts for seismic hazard assessment, the earthquake potential from the active faults near the cities of Zhangjiakou and Langfang in Hebei Province is examined. We estimate the effect of the coseismic stress changes of strong (M ⩾ 5.0) earthquakes on the major regional active faults, and mapped Coulomb stress change onto these target faults. More importantly our calculations reveal that positive stress changes caused by the largest events of the 1976 Tangshan sequence make the Xiadian and part of Daxing fault, thus considered the most likely sites of the next strong earthquake in the study area. The accumulated static stress changes that reached a value of up to 0.4 bar onto these faults, were subsequently incorporated in earthquake probability estimates for the next 30 years.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Detection of large,Holocene earthquakes using diatom analysis of coastal sedimentary sequences,Wellington, New Zealand

Paleoenvironmental reconstructions of three coastal waterbodies in Wellington, New Zealand, reveal that sites were isolated from the sea within the last 7500 years through coseismic uplift and barrier growth. Evidence for coseismic uplift consists of distinct transitions in diatom assemblages representing large changes in relative sea-level or water-table level, commonly in association with sedimentological evidence for catchment disturbance or marine influx. Transitions are abrupt, laterally extensive and synchronous within each waterbody. Amount of change across transition horizons is assessed using quantitative estimates of paleosalinity and waterbody type as proxies for relative sea-level change. Seven transitions involve large paleoenvironmental changes and provide evidence for earthquakes occurring at approximately 5100, 3200 (recorded at two sites), 2300 (recorded at two sites), 1000 cal years BP and 1855 AD. Five other transitions involve smaller paleoenvironmental changes and are considered to be consistent with effects of earthquakes but do not provide independent evidence for earthquake occurrence. These smaller transitions occur at approximately 6800, 3600, 2200, 1000 (coincident with a large transition) and 500 cal years BP. The data refine ages and provide information about the extent and effects of past large earthquakes in the region. These are the first paleoecologically derived earthquake signatures for Wellington and they contribute to the sparse collection worldwide of off-fault sedimentary earthquake records for predominantly strike-slip faults.     

Influence of the 2008 Wenchuan earthquake (Mw 7.9) on the occurrence probability of future earthquakes on neighboring faults

The Longquan–Shan fault and the Huya fault are two major neighboring faults of the Longmen–Shan fault zone where the 12 May 2008 Wenchuan earthquake (M_w 7.9) occurred. To study the influence of the Wenchuan event on these two active faults, we calculate changes of Coulomb stress on the Longquan–Shan fault and the Huya fault caused by the Wenchuan mainshock. Our results indicate that the Coulomb stress in the northern section (Zone A) of the Longquan–Shan fault is increased by 0.07–0.10 bars, that in the middle section (Zone B) by 0.04–0.11 bars, and that in the southern section (Zone C) shows almost no change. For the Huya fault, the Coulomb stress is decreased by 0.01–0.03 bars in the northern section (Zone A), 0.10–0.35 bars in the middle section (Zone B), and nearly 0.5 bars in the southern section (Zone C). The epicenter distribution of small earthquakes (M_L ⩾ 1.5) on the Longquan–Shan fault and the Huya fault after the Wenchuan earthquake is consistent with the distribution of the Coulomb stress change. This implies that the Wenchuan earthquake may have triggered small events on the Longquan–Shan fault, but inhibited those on the Huya fault. We then use the rate/state friction law to calculate the occurrence probability of future earthquakes in the study region for the next decade. They include the distribution of b-values, magnitude of completeness (M_c), the background seismicity rate, a value of Aσ_n and the duration for the transient effect (t_a) in the study region. We also estimate the earthquake occurrence probabilities on the neighboring faults after the Wenchuan earthquake. Our results show that, the occurrence probability of future earthquakes in the Longquan–Shan has a slight increase, being 7% for M ⩾ 5.0 shocks during the next decade, but the earthquake probability in the Huya region is reduced obviously, being 5–20%, 7–26% and 3–9% for M ⩾ 5.0 shocks during the next decade in sections A, B and C of the Huya fault, respectively.     

Late Holocene relative sea-level changes and the earthquake deformation cycle around upper Cook Inlet,Alaska

Multiple peat-silt couplets preserved in tidal marsh sediment sequences suggest that numerous great plate boundary earthquakes caused the coast around Cook Inlet, Alaska, to subside over the past 3500 years. Field and laboratory analyses of the two youngest couplets record the well-documented earthquake of AD 1964 and the penultimate one, approximately 850 cal yr BP. Diatom assemblages from a range of modern day estuarine environments from tidal flat through salt marsh to acidic bog produce quantitative diatom transfer function models for elevation reconstructions based on fossil samples. Only nine out of 124 fossil assemblages analysed, including previously published data for the AD 1964 earthquake, have a poor modern analogue. Calibration of fossil samples indicate co-seismic subsidence of 1.50±0.32 m for AD 1964, similar to measurements taken after the earthquake, and 1.45±0.34 m for the ∼850 cal yr BP earthquake. Elevation standard errors for individual fossil samples range from ∼0.08 m in peat layers to ∼0.35 m in silt units. Lack of a chronology within fossil silt units prevents identification of changes in the rate of recovery and land uplift between the post-seismic and inter-seismic periods. However, preservation of multiple peat-silt couplets indicates no net emergence over multiple earthquake cycles. Glacio-isostatic movements from Little Ice Age glacier advance and retreat explains a ∼0.15 m relative sea-level oscillation recorded within the peat layer subsequently submerged as a result of the AD 1964 earthquake. Before both this and the ∼850 cal yr BP earthquake, diatom assemblages suggest pre-seismic relative sea-level rise of ∼0.12±0.13 m, representing possible precursors to great earthquakes.     

Arsenic and other drinking water quality issues,Muzaffargarh District,Pakistan

In 49 samples of groundwater, sampled in Muzaffargarh District of south-western Punjab, central Pakistan, concentrations of As exceeded the World Health Organisation provisional guideline value, and United States Environmental Protection Agency (USEPA) Maximum Contaminant Level (MCL), of 10 μg L⁻¹ in 58% of samples and reached up to 906 μg L⁻¹. In this semi-arid region canal irrigation has lead to widespread water-logging, and evaporative concentration of salts has the potential to raise As concentrations in shallow groundwater well above 10 μg L⁻¹. In fact, in rural areas, concentrations stay below 25 μg L⁻¹ because As in the oxic shallow groundwater, and in recharging water, is sorbed to aquifer sediments. In some urban areas, however, shallow groundwater is found to contain elevated levels of As. The spatial distribution of As-rich shallow groundwater indicates either direct contamination with industrial or agricultural chemicals, or some other anthropogenic influence. Geochemical evidence suggests that pollutant organics from unconfined sewage and other sources drives reduction of hydrous ferric oxide (HFO) releasing sorbed As to shallow groundwater. The situation is slightly less clear for seven wells sampled which tap deeper groundwater, all of which were found with >50 μg L⁻¹ As. Here As concentrations seem to increase with depth and differing geochemical signatures are seen, suggesting that As concentrations in older groundwater may be governed by different processes. Other data on parameters of potential concern in drinking water are discussed briefly at the end of the paper.     

Isotopic and geochemical identifications of groundwater salinisation processes in Salalah coastal plain,Sultanate of Oman

A detail investigation was carried out to improve the current knowledge of groundwater salinisation processes in coastal aquifers using hydrochemical and isotopic parameters. Data of major ions for 40 wells located in the Salalah plain aquifer, Sultanate of Oman, were collected during pre-monsoon 2004 and analysed. The groundwater changes along the general flow path towards the coast from fresh (EC < 1500 μS/cm), brackish (EC: 1500–3000 μS/cm) and saline (EC > 3000 μS/cm). Results of inverse modeling simulations using PHREEQC show that dissolution of halite may be the main source of Cl and Na in the study area. Ionic delta calculation indicates that the depletion of Na and K and enrichment of Ca and Mg in groundwater were probably attributed to reverse ion exchange reactions. During a sampling campaign conducted in October 2015, 11 groundwater samples were collected for Cl, Br and isotopic analysis (²H/¹⁸O). Molar Cl/Br ratios in fresh groundwater were higher than those of seawater, indicating the impact of halite dissolution on the groundwater quality. For saline groundwater, these ratios were less than those of seawater, showing the influence of anthropogenic input from agriculture on the same. Relatively depleted isotopic signature of all groundwater samples show that the monsoon precipitation is the main source of groundwater recharge in the study area.     

A conceptual framework of groundwater flow in some crystalline aquifers in Southeastern Ghana

A conceptual groundwater flow model was developed for the crystalline aquifers in southeastern part of the Eastern region, Ghana. The objective was to determine approximate levels of groundwater recharge, estimate aquifer hydraulic parameters, and then test various scenarios of groundwater extraction under the current conditions of recharge. A steady state groundwater flow model has been calibrated against measured water levels of 19 wells in the area. The resulting recharge is estimated to range from 8.97 × 10^?5 m/d to 7.14 × 10^?4 m/d resulting in a basin wide average recharge of about 9.6% of total annual precipitation, which results in a basin wide quantitative recharge of about 2.4 million m³/d in the area. This compares to recharge estimated from the chloride mass balance of 7.6% of precipitation determined in this study. The general groundwater flow in the area has also been determined to conform to the general northeast–southwest structural grain of the country. The implication is that the general hydrogeology is controlled by post genetic structural entities imposed on the rocks to create ingresses for sufficient groundwater storage and transport. Calibrated aquifer hydraulic conductivities range between 0.99 m/d and over 19.4 m/d. There is a significant contribution of groundwater discharge to stream flow in the study area. Increasing groundwater extraction will have an effect on stream flow. This study finds that the current groundwater extraction levels represent only 0.17% of the annual recharge from precipitation, and that groundwater can sustain future increased groundwater demands from population growth and industrialization.     

Source complexity of the 4 March 2010 Jiashian,Taiwan, Earthquake determined by joint inversion of teleseismic and near field data

The Jiashian earthquake (M_L 6.4) occurred on 4 March 2010. It was the largest inland event in southern Taiwan of 2010. The mainshock location was unexpected since it occurred in an area with relatively low background seismicity. In addition, reports of earthquake focal mechanisms do not fit with any known active fault geometry. In order to understand the origin of this earthquake, especially its rupture process, we perform a joint source inversion by using teleseismic body wave, GPS coseismic displacements and near field ground motion data. In this study, we considered a northwest–southeast trending fault with a northeast dip retrieved from GPS coseismic data and aftershocks distribution. To analyze the detailed slip distribution in space and time, we used near field 3D Green’s functions provided by spectral-element method and a full time–space inversion technique. We find a complex rupture process with several slip patches distributed inside two main asperities. The slip map reveals a mean slip of 12.9 cm for a maximum slip of 27.3 cm leading to a M_w 6.47 for this event. The rupture initiates in the deepest portion of the fault at 20 km depth, and propagated upward up to 2 km depth to form the two asperities. The source time function of this event revealed two pulses corresponding to the two asperities, for a total duration time of about 16 s. Most aftershocks occurred near the upper boundary of the deepest asperity while no aftershocks are located close to the shallowest one. We infer that the locations of these slip patches are related to the surrounding fault systems that may have restricted the rupture propagation during the earthquake.     

Seismic behavior in central Taiwan: Response to stress evolution following the 1999 Mw 7.6 Chi-Chi earthquake

Following the 1999 M_w 7.6 Chi-Chi earthquake, a large amount of seismicity occurred in the Nantou region of central Taiwan. Among the seismic activities, eight M_w ⩾ 5.8 earthquakes took place following the Chi-Chi earthquake, whereas only four earthquakes with comparable magnitudes took place from 1900 to 1998. Since the seismicity rate during the Chi-Chi postseismic period has never returned to the background level, such seismicity activation cannot simply be attributed to modified Omori’s Law decay. In this work, we attempted to associate seismic activities with stress evolution. Based on our work, it appears that the spatial distribution of the consequent seismicity can be associated with increasing coseismic stress. On the contrary, the stress changes imparted by the afterslip; lower crust–upper mantle viscoelastic relaxation; and sequent events resulted in a stress drop in most of the study region. Understanding seismogenic mechanisms in terms of stress evolution would be beneficial to seismic hazard mitigation.     

Seismic anisotropy in the wedge above the Philippine Sea slab beneath Kanto and southwest Japan derived from shear wave splitting

We conduct shear wave splitting measurements on waveform data from the Hi-net and the broadband F-net seismic stations in Kanto and SW Japan generated by shallow and intermediate-depth earthquakes occurring in the subducting Philippine Sea and Pacific slabs. We obtain 1115 shear wave splitting parameter pairs. The results are divided into those from the shallow (depth < 50 km) and the deep (depth > 50 km) events. The deep events beneath Kanto are further divided into PHS1 and PHS2 (upper and lower planes of the double seismic zone in the Philippine Sea slab, respectively), PAC1 and PAC2 (western and eastern Pacific slab, respectively) events. The results from the shallow events represent the crustal anisotropy, and their fast directions are more or less aligned in the σ_Hmax directions, implying that the anisotropy is produced by the alignment of the vertical cracks in the crust induced by the compressive stresses. In Kanto, Kii Peninsula and Kyushu regions, the results from the deep events suggest a contribution from the mantle wedge anisotropy. Events from all groups beneath Kanto show NW, NE and EW fast directions. This complex pattern seems to be produced by the corner flows induced by both the WNW PAC plate subduction and the oblique NNW PHS slab subduction with the associated olivine lattice-preferred orientations (LPOs), and the anisotropy frozen in the PHS slab. The deep events beneath Kii Peninsula show NE and NW fast directions and may be produced by the corner flow produced by the NNW PHS slab subduction with the associated olivine LPOs. The NE directions might also be produced by the segregated melts in the thin layers parallel to the PHS slab subduction. The deep events beneath N Kyushu show NNW fast directions, which may result from the southeastward flow in the upper mantle inferred from the stresses in the upper plate. Results from the deep events beneath middle-south Kyushu show dominantly E–W fast directions, in both the fore- and back-arcs. They may be produced by the corner flow of the westward PHS slab subduction with the olivine LPOs. Because the source regions with multiple fast directions are not resolved in this study, further detailed analyses of shear wave splitting are necessary for a better understanding of the stress state, the induced mantle flow, and the melt-segregation processes.     

Bearing of plate geometry and rheology on shallow-focus mega-thrust seismicity with special reference to 26 December 2004 Sumatra event

We propose a model pertaining to the generation of 26th December 2004 off Sumatra mega-event in the backdrop of other similar type earthquakes along subduction zones around the world. Reconstructions of Benioff trajectories through the hypocenters of historical earthquakes including six mega-earthquakes indicate (i) confinement of hypocenters right within the descending lithosphere, and (ii) natural coincidence of foci of the mega-events around the zones of plate flexing. These observations are discussed in detail with special emphasis on the Sumatra margin considering the role of rheological anomaly across the cross-section of the descending lithosphere; yield strength envelope and residual stress accumulation through time. The intraplate origin of shallow mega-thrust earthquakes allowed us to advocate the ‘zone of flexing’ along the profiles of the subducting plates as nodal area for stress concentration. We propose here that at elevated confining pressure and temperature, loading of unidirectional cyclic stress on time-average bending stress enhanced the material yield strength (i.e., strain-hardening), and leads the semi-brittle portion of the lithosphere into near-brittle condition through rheological transformation. Under subsequent rise in neutral surface and increase in compressive stress field, non-coaxial deformation triggered shear failure on 26th December 2004 preferably at the rheological interface between strain-hardened near-brittle layer and deformed ductile layer within the sub-oceanic mantle.A two-stage fracture mechanism viz. a slow (～1.1 km/s) bilateral initiation in an essentially strain-hardened near-brittle domain and a follow-up very rapid progression (3.3 km/s) in the brittle, crustal domain was mainly involved in the generation of 2004 off Sumatra mega-event. Estimation shows an amount of 3.38 × 10²² to 4.50 × 10²² N m seismic moment (M_o) and 8.95–9.03 moment magnitude (M_w) for the southern part of the 1300 km extended rupture i.e. between the North Andaman to the north and the Sumatra at its south. The study necessitates the reassessment of other shallow-focus mega-thrust earthquakes along the subduction margins around the globe.     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

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.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Geochemical and isotopic variations in shallow groundwater in areas of the Fayetteville Shale development,north-central Arkansas

Exploration of unconventional natural gas reservoirs such as impermeable shale basins through the use of horizontal drilling and hydraulic fracturing has changed the energy landscape in the USA providing a vast new energy source. The accelerated production of natural gas has triggered a debate concerning the safety and possible environmental impacts of these operations. This study investigates one of the critical aspects of the environmental effects; the possible degradation of water quality in shallow aquifers overlying producing shale formations. The geochemistry of domestic groundwater wells was investigated in aquifers overlying the Fayetteville Shale in north-central Arkansas, where approximately 4000 wells have been drilled since 2004 to extract unconventional natural gas. Monitoring was performed on 127 drinking water wells and the geochemistry of major ions, trace metals, CH₄ gas content and its C isotopes (δ¹³C_CH4), and select isotope tracers (δ¹¹B, ⁸⁷Sr/⁸⁶Sr, δ²H, δ¹⁸O, δ¹³C_DIC) compared to the composition of flowback-water samples directly from Fayetteville Shale gas wells. Dissolved CH₄ was detected in 63% of the drinking-water wells (32 of 51 samples), but only six wells exceeded concentrations of 0.5 mg CH₄/L. The δ¹³C_CH4 of dissolved CH₄ ranged from −42.3‰ to −74.7‰, with the most negative values characteristic of a biogenic source also associated with the highest observed CH₄ concentrations, with a possible minor contribution of trace amounts of thermogenic CH₄. The majority of these values are distinct from the reported thermogenic composition of the Fayetteville Shale gas (δ¹³C_CH4 = −35.4‰ to −41.9‰). Based on major element chemistry, four shallow groundwater types were identified: (1) low (<100 mg/L) total dissolved solids (TDS), (2) TDS > 100 mg/L and Ca–HCO₃ dominated, (3) TDS > 100 mg/L and Na–HCO₃ dominated, and (4) slightly saline groundwater with TDS > 100 mg/L and Cl > 20 mg/L with elevated Br/Cl ratios (>0.001). The Sr (⁸⁷Sr/⁸⁶Sr = 0.7097–0.7166), C (δ¹³C_DIC = −21.3‰ to −4.7‰), and B (δ¹¹B = 3.9–32.9‰) isotopes clearly reflect water–rock interactions within the aquifer rocks, while the stable O and H isotopic composition mimics the local meteoric water composition. Overall, there was a geochemical gradient from low-mineralized recharge water to more evolved Ca–HCO₃, and higher-mineralized Na–HCO₃ composition generated by a combination of carbonate dissolution, silicate weathering, and reverse base-exchange reactions. The chemical and isotopic compositions of the bulk shallow groundwater samples were distinct from the Na–Cl type Fayetteville flowback/produced waters (TDS ∼10,000–20,000 mg/L). Yet, the high Br/Cl variations in a small subset of saline shallow groundwater suggest that they were derived from dilution of saline water similar to the brine in the Fayetteville Shale. Nonetheless, no spatial relationship was found between CH₄ and salinity occurrences in shallow drinking water wells with proximity to shale-gas drilling sites. The integration of multiple geochemical and isotopic proxies shows no direct evidence of contamination in shallow drinking-water aquifers associated with natural gas extraction from the Fayetteville Shale.     

Strontium isotope quantification of siderite,brine and acid mine drainage contributions to abandoned gas well discharges in the Appalachian Plateau

Unplugged abandoned oil and gas wells in the Appalachian region can serve as conduits for the movement of waters impacted by fossil fuel extraction. Strontium isotope and geochemical analysis indicate that artesian discharges of water with high total dissolved solids (TDS) from a series of gas wells in western Pennsylvania result from the infiltration of acidic, low Fe (Fe < 10 mg/L) coal mine drainage (AMD) into shallow, siderite (iron carbonate)-cemented sandstone aquifers. The acidity from the AMD promotes dissolution of the carbonate, and metal- and sulfate-contaminated waters rise to the surface through compromised abandoned gas well casings. Strontium isotope mixing models suggest that neither upward migration of oil and gas brines from Devonian reservoirs associated with the wells nor dissolution of abundant nodular siderite present in the mine spoil through which recharge water percolates contribute significantly to the artesian gas well discharges. Natural Sr isotope composition can be a sensitive tool in the characterization of complex groundwater interactions and can be used to distinguish between inputs from deep and shallow contamination sources, as well as between groundwater and mineralogically similar but stratigraphically distinct rock units. This is of particular relevance to regions such as the Appalachian Basin, where a legacy of coal, oil and gas exploration is coupled with ongoing and future natural gas drilling into deep reservoirs.     

Late Paleozoic peperites in West Junggar,China, and how they constrain regional tectonic and palaeoenvironmental setting

Late Paleozoic peperites have been identified for the first time at the bottom of Tailegula Formation in West Junggar, China. This finding is significant for the reconstruction of Late Paleozoic evolution in the Junggar region. The peperites form successions up to 500 m thick interbedded with basaltic lava and sedimentary rocks. Four types of peperites are described and interpreted as resulting from basaltic lava bulldozed into wet, unconsolidated sediments at their basal contacts. Zircon LA-ICP-MS U–Pb dating of a tuff lens enclosed by lava showed that the peperites formed in the Late Devonian (ca. 364 Ma). The peperite-bearing units probably formed at a water depth of less than 3 km and are generally undeformed, occurring in continuous stratigraphic sections distributed regionally over a distance of 100 km on either side of the Darbut and Baijiantan ophiolitic belts, in contrast to the highly deformed slices of ophiolite. They demonstrate that the Darbut and Baijiantan ophiolitic belts should not be interpreted as significant plate boundaries and represent the underlying ocean crust uplifted along tectonic lineaments within a continuous shallow remnant ocean basin. The peperites formed during the spreading phase of the remnant ocean basin and represent the final stages of creation of oceanic crust.     

Hydrogeochemical evolution of confined groundwater in northeastern Osaka Basin,Japan: estimation of confined groundwater flux based on a cation exchange mass balance method

A confined aquifer system has developed in argillaceous marine and freshwater sediments of Pliocene–Holocene age in the northeastern Osaka Basin (NEOB) in central Japan. The shallow groundwater (<100 m) in the system is recharged in a northern hilly to mountainous area with dominantly Ca-HCO₃ type water, which changes as it flows toward the SW to Mg-HCO₃ type and then to Na-HCO₃ type water. Comparison of the chemical and Sr isotopic compositions of the groundwater with those of the bulk and exchangeable components of the underground sediments indicates that elements leached from the sediments contribute negligibly to the NEOB aquifer system. Moreover, model calculations show that contributions of paleo-seawater in the deep horizon and of river water at the surface are not major factors of chemical change of the groundwater. Instead, the zonal pattern of the HCO₃-dominant groundwater is caused by the loss of Ca²⁺ from the water as it is exchanged for Mg²⁺ in clays, followed by loss of Mg + Ca as they are exchanged for Na + K in clays between the Ca-HCO₃ type recharge water and the exchangeable cations in the clay layers, which were initially enriched in Na⁺. Part of this process was reproduced in a chromatographic experiment in which Na type water with high ⁸⁷Sr/⁸⁶Sr was obtained from Mg type water with low ⁸⁷Sr/⁸⁶Sr by passing it through marine clay packed in a column. The flux of recharge water into the confined aquifer system according to this chromatographic model is estimated to be 0.99 mm/day, which is compatible with the average recharge flux to unconfined groundwater in Japan (1 mm/day).     

Origin,distribution and hydrogeochemical controls on methane occurrences in shallow aquifers,southwestern Ontario,Canada

Natural gas reservoirs in organic-rich shales in the Appalachian and Michigan basins in the United States are currently being produced via hydraulic fracturing. Stratigraphically-equivalent shales occur in the Canadian portion of the basins in southwestern Ontario with anecdotal evidence of gas shows, yet there has been no commercial shale gas production to date. To provide baseline data in the case of future environmental issues related to hydraulic fracturing and shale gas production, such as leakage of natural gas, saline water, and/or hydraulic fracturing fluids, and to evaluate hydrogeochemical controls on natural gas accumulations in shallow groundwater in general, this study investigates the origin and distribution of natural gas and brine in shallow aquifers across southwestern Ontario. An extensive geochemical database of major ion and trace metal chemistry and methane concentrations of 1010 groundwater samples from shallow, domestic wells in bedrock and overburden aquifers throughout southwestern Ontario was utilized. In addition, select wells (n = 36) were resampled for detailed dissolved gas composition, δ¹³C of CH₄, C₂, C₃, and CO₂, and δD of CH₄. Dissolved gases in groundwater from bedrock and overburden wells were composed primarily of CH₄ (29.7–98.6 mol% of total gas volume), N₂ (0.8–66.2 mol%), Ar + O₂ (0.2–3.4 mol%), and CO₂ (0–1.2 mol%). Ethane was detected, but only in low concentrations (<0.041 mol%), and no other higher chain hydrocarbons were present, except for one well in overburden overlying the Dundee Formation, which contained 0.81 mol% ethane and 0.21 mol% propane. The highest methane concentrations (30 to >100 in situ % saturation) were found in bedrock wells completed in the Upper Devonian Kettle Point Formation, Middle Devonian Hamilton Group and Dundee Formation, and in surficial aquifers overlying these organic-rich shale-bearing formations, indicating that bedrock geology is the primary control on methane occurrences. A few (n = 40) samples showed Na–Cl–Br evidence of brine mixing with dilute groundwater, however only one of these samples contained high (>60 in situ % saturation) CH₄. The relatively low δ¹³C values of CH₄ (−89.9‰ to −57.3‰), covariance of δD values of CH₄ and H₂O, positive correlation between δ¹³C values of CH₄ and CO₂, and lack of higher chain hydrocarbons (C₃₊) in all but one dissolved gas sample indicates that the methane in groundwater throughout the study area is primarily microbial in origin. The presence or absence of alternative electron acceptors (e.g. dissolved oxygen, Fe, NO₃, SO₄), in addition to organic substrates, controls the occurrence of microbial CH₄ in shallow aquifers. Microbial methane has likely been accumulating in the study area, since at least the Late Pleistocene to the present, as indicated by the co-variance and range of δD values of CH₄ (−314‰ to −263‰) and associated groundwater (−19‰ to −6‰ δD-H₂O).