Regional groundwater focusing of nitrogen and noble gases into the Hugoton-Panhandle giant gas field, USA

The Hugoton-Panhandle giant gas field, located across SW Kansas and the Texas and Oklahoma panhandles in the USA, is the case type example of high nitrogen concentrations in a natural gas being linked with high helium concentrations. We collected 31 samples from producing wells in a north-south traverse of the 350-km-long field. The samples reflect the previously observed north-south change in ⁴He/N₂, with values changing from 0.020 to 0.049 respectively. ³He/⁴He, ²¹Ne/²²Ne, and ⁴⁰Ar/³⁶Ar vary between 0.14-0.25 Ra, 0.0373-0.0508, and 818-1156 respectively, and are caused by quantifiable contributions from mantle, crustal, and atmosphere-derived sources. The atmosphere-derived ²⁰Ne/³⁶Ar ratios are indistinguishable from groundwater values. The crustal ⁴He/²¹Ne* and ⁴He/⁴⁰Ar* ratios show a 60% excess of ⁴He compared to predicted production ratios in the crust and are typical of noble gases released from the shallow crust. The mantle ³He/N₂ and groundwater-recharge ³⁶Ar/N₂ ratios enable us to rule out significant magmatic or atmosphere contributions to the gas field N₂, which is dominantly crustal in origin.Correlated ²⁰Ne/N₂ and ⁴He/N₂ shows mixing between two distinct crustal N₂ components. One N₂ component (N₂*) is associated with the crustal ⁴He and groundwater-derived ²⁰Ne, and the other with no resolvable noble gas contribution. Measured δ¹⁵N_N2 values vary from +2.7‰ to +9.4‰. The N₂* and non-He-associated N₂ endmembers are inferred to have δ¹⁵N_N2 = −3‰ and +13‰ and contribute from between 25-60% and 75-40% of the nitrogen respectively. The non-He-associated nitrogen is probably derived from relatively mature organic matter in the sedimentary column. The δ¹⁵N_N2* value is not compatible with a crystalline or high-grade metamorphic source and, similar to the ⁴He, is inferred to be from a shallow or low metamorphic-grade source rock. ⁴He mass balance requires a regional crustal source, its association with significant magmatic ³He pointing to a tectonically active source to the west of the Hugoton system. The volume of groundwater required to source the ²⁰Ne in the gas field demonstrates the viability of the groundwater system in providing the collection, transport, and focusing mechanism for the ⁴He and N₂*. The N₂*/²⁰Ne ratio shows that the N₂* transport must be in the aqueous phase, and that the degassing mechanism is probably contact between the regional groundwater system and the preexisting reservoir hydrocarbon gas phase.     

Large-scale moraine deformation at the Athabasca Glacier,Jasper National Park,Alberta, Canada

In this paper the development of a large-scale gravitational deformation involving the eastern lateral moraine of the Athabasca Glacier in Jasper National Park, Alberta, Canada, is described. Interpretation and analysis of sequential aerial photographs indicates that a 540-m-wide segment of the eastern lateral moraine began to deform in the early 1950s; however, significant movement only began in the late 1960s. Since then, the moraine has undergone progressive gravitational deformation leading to a network of fractures, bulging, and the development of a large gap in the moraine crest. Geographic information system analysis of topographic changes between 1967 and 2006 indicates that the displaced volume of the moraine is approximately 9.0 × 10⁵ m³. In the last 39 years, the moraine crest has displaced 55 m (1.4 m yr⁻¹) down towards the glacier. The development of slope instability is linked to a combination of debuttressing from recent glacier recession, deformation of the moraine, as well as the movement of a large, mobile, debris-mantled slope impinging the upslope margin of the lateral moraine. This case study illustrates the importance of glacial conditioning and local geomorphological factors in creating conditions for large-scale moraine instability in recently deglacierized alpine basins.     

Origins of Xenolithic Eclogites and Pyroxenites from the Central Slave Craton, Canada

Major- and trace-element and Sr–Nd–Hf isotopic compositionsof garnet and clinopyroxene in kimberlite-borne eclogite andpyroxenite xenoliths were used to establish their origins andevolution in the subcontinental lithospheric mantle beneaththe central Slave Craton, Canada. The majority of eclogitescan be assigned to three groups (high-Mg, high-Ca or low-Mgeclogites) that have distinct trace-element patterns. Althoughpost-formation metasomatism involving high field strength element(HFSE) and light rare earth element (LREE) addition has partiallyobscured the primary compositional features of the high-Mg andhigh-Ca eclogites, trace-element features, such as unfractionatedmiddle REE (MREE) to heavy REE (HREE) patterns suggestive ofgarnet-free residues and low Zr/Sm consistent with plagioclaseaccumulation, could indicate a subduction origin from a broadlygabbroic protolith. In this scenario, the low REE and smallpositive Eu anomalies of the high-Mg eclogites suggest moreprimitive, plagioclase-rich protoliths, whereas the high-Caeclogites are proposed to have more evolved protoliths withhigher (normative) clinopyroxene/plagioclase ratios plus trappedmelt, consistent with their lower Mg-numbers, higher REE andabsence of Eu anomalies. In contrast, the subchondritic Zr/Hfand positive slope in the HREE of the low-Mg eclogites are similarto Archaean second-stage melts and point to a previously depletedsource for their precursors. Low ratios of fluid-mobile to lessfluid-mobile elements and of LREE to HREE are consistent withdehydration and partial melt loss for some eclogites. The trace-elementcharacteristics of the different eclogite types translate intolower _Nd for high-Mg eclogites than for low-Mg eclogites. Withinthe low-Mg group, samples that show evidence for metasomaticenrichment in LREE and HFSE have lower _Nd and _Hf than a samplethat was apparently not enriched, pointing to long-term evolutionat their respective parent–daughter ratios. Garnet andclinopyroxene in pyroxenites show different major-element relationshipsfrom those in eclogites, such as an opposite CaO–Na₂Otrend and the presence of a CaO–Cr₂O₃ trend, independentof whether or not opx is part of the assemblage. Therefore,these two rock types are probably not related by fractionationprocesses. The presence of opx in about half of the samplesprecludes direct crystallization from eclogite-derived melts.They probably formed from hybridized melts that reacted withthe peridotitic mantle. KEY WORDS: eclogites; pyroxenite xenoliths; mantle xenoliths; eclogite trace elements; eclogite Sr isotopes; eclogite Hf isotopes; eclogite Nd isotopes     

Evaluation of the importance of clay confining units on groundwater flow in alluvial basins using solute and isotope tracers: the case of Middle San Pedro Basin in southeastern Arizona (USA)

As groundwater becomes an increasingly important water resource worldwide, it is essential to understand how local geology affects groundwater quality, flowpaths and residence times. This study utilized multiple tracers to improve conceptual and numerical models of groundwater flow in the Middle San Pedro Basin in southeastern Arizona (USA) by determining recharge areas, compartmentalization of water sources, flowpaths and residence times. Ninety-five groundwater and surface-water samples were analyzed for major ion chemistry (water type and Ca/Sr ratios) and stable (¹⁸O, ²H, ¹³C) and radiogenic (³H, ¹⁴C) isotopes, and resulting data were used in conjunction with hydrogeologic information (e.g. hydraulic head and hydrostratigraphy). Results show that recent recharge (<60 years) has occurred within mountain systems along the basin margins and in shallow floodplain aquifers adjacent to the San Pedro River. Groundwater in the lower basin fill aquifer (semi confined) was recharged at high elevation in the fractured bedrock and has been extensively modified by water-rock reactions (increasing F and Sr, decreasing ¹⁴C) over long timescales (up to 35,000 years BP). Distinct solute and isotope geochemistries between the lower and upper basin fill aquifers show the importance of a clay confining unit on groundwater flow in the basin, which minimizes vertical groundwater movement.     

Frontogenesis and Frontal Progression of a Trapping-Generated Estuarine Convergence Front and Its Influence on Mixing and Stratification

Estuarine fronts are well known to influence transport of waterborne constituents such as phytoplankton and sediment, yet due to their ephemeral nature, capturing the physical driving mechanisms and their influence on stratification and mixing is difficult. We investigate a repetitive estuarine frontal feature in the Snohomish River Estuary that results from complex bathymetric shoal/channel interactions. In particular, we highlight a trapping mechanism by which mid-density water trapped over intertidal mudflats converges with dense water in the main channel forming a sharp front. The frontal density interface is maintained via convergent transverse circulation driven by the competition of lateral baroclinic and centrifugal forcing. The frontal presence and propagation give rise to spatial and temporal variations in stratification and vertical mixing. Importantly, this front leads to enhanced stratification and suppressed vertical mixing at the end of the large flood tide, in contrast to what is found in many estuarine systems. The observed mechanism fits within the broader context of frontogenesis mechanisms in which varying bathymetry drives lateral convergence and baroclinic forcing. We expect similar trapping-generated fronts may occur in a wide variety of estuaries with shoal/channel morphology and/or braided channels and will similarly influence stratification, mixing, and transport.     

Exploring controls of the early and stepped deglaciation on the western margin of the British Irish Ice Sheet

New optically stimulated luminescence dating and Bayesian models integrating all legacy and BRITICE-CHRONO geochronology facilitated exploration of the controls on the deglaciation of two former sectors of the British–Irish Ice Sheet, the Donegal Bay (DBIS) and Malin Sea ice-streams (MSIS). Shelf-edge glaciation occurred ~27 ka, before the global Last Glacial Maximum, and shelf-wide retreat began 26–26.5 ka at a rate of ~18.7–20.7 m a^–1. MSIS grounding zone wedges and DBIS recessional moraines show episodic retreat punctuated by prolonged still-stands. By ~23–22 ka the outer shelf (~25 000 km²) was free of grounded ice. After this time, MSIS retreat was faster (~20 m a^–1 vs. ~2–6 m a^–1 of DBIS). Separation of Irish and Scottish ice sources occurred ~20–19.5 ka, leaving an autonomous Donegal ice dome. Inner Malin shelf deglaciation followed the submarine troughs reaching the Hebridean coast ~19 ka. DBIS retreat formed the extensive complex of moraines in outer Donegal Bay at 20.5–19 ka. DBIS retreated on land by ~17–16 ka. Isolated ice caps in Scotland and Ireland persisted until ~14.5 ka. Early retreat of this marine-terminating margin is best explained by local ice loading increasing water depths and promoting calving ice losses rather than by changes in global temperatures. Topographical controls governed the differences between the ice-stream retreat from mid-shelf to the coast.     

Mixed fluid sources involved in diamond growth constrained by Sr–Nd–Pb–C–N isotopes and trace elements

Sub-micrometer inclusions in diamonds carry high-density fluids (HDF) from which the host diamonds have precipitated. The chemistry of these fluids is our best opportunity of characterizing the diamond-forming environment. The trace element patterns of diamond fluids vary within a limited range and are similar to those of carbonatitic/kimberlitic melts that originate from beneath the lithospheric mantle. A convecting mantle origin for the fluid is also implied by C isotopic compositions and by a preliminary Sr isotopic study (Akagi, T., Masuda, A., 1988. Isotopic and elemental evidence for a relationship between kimberlite and Zaire cubic diamonds. Nature 336, 665–667.). Nevertheless, the major element chemistry of HDFs is very different from that of kimberlites and carbonatites, varying widely and being characterized by extreme K enrichment (up to ～ 39 wt.% on a water and carbonate free basis) and high volatile contents. The broad spectrum of major element compositions in diamond-forming fluids has been related to fluid–rock interaction and to immiscibility processes.Elemental signatures can be easily modified by a variety of mantle processes whereas radiogenic isotopes give a clear fingerprint of the time-integrated evolution of the fluid source region. Here we present the results of the first multi radiogenic-isotope (Sr, Nd, Pb) and trace element study on fluid-rich diamonds, implemented using a newly developed off-line laser sampling technique. The data are combined with N and C isotope analysis of the diamond matrix to better understand the possible sources of fluid involved in the formation of these diamonds. Sr isotope ratios vary significantly within single diamonds. The highly varied but unsupported Sr isotope ratios cannot be explained by immiscibility processes or fluid-mineral elemental fractionations occurring at the time of diamond growth. Our results demonstrate the clear involvement of a mixed fluid, with one component originating from ancient incompatible element-enriched parts of the lithospheric mantle while the trigger for releasing this fluid source was probably carbonatitic/kimberlitic melts derived from greater depths. We suggest that phlogopite mica was an integral part of the enriched lithospheric fluid source and that breakdown of this mica releases K and radiogenic Sr into a fluid phase. The resulting fluids operate as a major metasomatic agent in the sub-continental lithospheric mantle as reflected by the isotopic composition and trace element patterns of G10 garnets.     

Numerical computations of a theoretical model of ribbed moraine formation

We develop numerical solutions of a theoretical model which has been proposed to explain the formation of subglacial bedforms. The model has been shown to have the capability of producing bedforms in two dimensions, when they may be interpreted as ribbed moraine. However, these investigations have left unanswered the question of whether the theory is capable of producing fully three‐dimensional bedforms such as drumlins. We show that, while the three‐dimensional calculations show realistic quasi‐three‐dimensional features such as dislocations in the ribbing pattern, they do not produce genuine three‐dimensional drumlins. We suggest that this inadequacy is due to the treatment of subglacial drainage in the theory as a passive variable, and thus that the three‐dimensional forms may be associated with conditions of sufficient subglacial water flux. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling wetland connectivity during overbank flooding in a tropical floodplain in north Queensland,Australia

Hydrological connectivity between floodplain wetlands and rivers is one of the principal driving mechanisms for the diversity, productivity and interactions of the major biota in river–floodplain systems. This article describes a method of quantifying flood‐induced overbank connectivity using a hydrodynamic model (MIKE 21) to calculate the timing, the duration and the spatial extent of the connections between several floodplain wetlands and rivers in the Tully–Murray catchment, north Queensland, Australia. Areal photogrammetry and field surveyed stream cross data were used to reproduce floodplain topography and rivers in the model. Laser altimetry (LiDAR)–derived fine resolution elevation data, for the central floodplain, were added to the topography model to improve the resolution of key features including wetlands, flow pathways and natural and artificial flow barriers. The hydrodynamic model was calibrated using a combination of in‐stream and floodplain gauge records. A range of off‐stream wetlands including natural and artificial, small and large were investigated for their connectivity with two main rivers (Tully and Murray) flowing over the floodplain for flood events of 1‐, 20‐ and 50‐year recurrence intervals. The duration of the connection of individual wetlands varied from 1 to 12 days, depending on flood magnitude and location in the floodplain, with some wetlands only connected during large floods. All of the wetlands studied were connected to the Tully River for shorter periods than they were to the Murray River because of the higher bank heights and levees on the Tully River and wetland proximity to the Murray River. Other than hydrology, land relief, riverbank elevation and levee banks along the river were found key factors controlling the degree of connectivity. These variations in wetland connectivity could have important implications for aquatic biota that move between rivers and off‐stream habitats during floods. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial analysis of sand dunes with a new global topographic dataset: new approaches and opportunities

On June 29, 2009, version 1 of the ASTER GDEM (Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model) was made available at no cost or restrictions to users worldwide via electronic download. The ASTER GDEM provides expanded spatial coverage and better resolution than other global digital elevation models (DEMs). In this paper we demonstrate how the ASTER GDEM provides new opportunities for investigating large aeolian sand dunes in three‐dimensions. Two dune‐specific spatial analysis methods are presented to illustrate potential applications of these data for discriminating dune generations and quantifying spatial variations of sediment supply. Moreover, we review how existing and emerging fields of dune pattern analysis and simulation modeling will be able to make significant advances through application of these data, potentially leading to future progress in studies of dune morphodynamics, environmental controls, and paleoenvironmental reconstructions. Copyright © 2010 John Wiley & Sons, Ltd.