Energetic Charged Particles Above Thunderclouds

The French government has committed to launch the satellite TARANIS to study transient coupling processes between the Earth’s atmosphere and near-Earth space. The prime objective of TARANIS is to detect energetic charged particles and hard radiation emanating from thunderclouds. The British Nobel prize winner C.T.R. Wilson predicted lightning discharges from the top of thunderclouds into space almost a century ago. However, new experiments have only recently confirmed energetic discharge processes which transfer energy from the top of thunderclouds into the upper atmosphere and near-Earth space; they are now denoted as transient luminous events, terrestrial gamma-ray flashes and relativistic electron beams. This meeting report builds on the current state of scientific knowledge on the physics of plasmas in the laboratory and naturally occurring plasmas in the Earth’s atmosphere to propose areas of future research. The report specifically reflects presentations delivered by the members of a novel Franco-British collaboration during a meeting at the French Embassy in London held in November 2011. The scientific subjects of the report tackle ionization processes leading to electrical discharge processes, observations of transient luminous events, electromagnetic emissions, energetic charged particles and their impact on the Earth’s atmosphere. The importance of future research in this area for science and society, and towards spacecraft protection, is emphasized.     

An automated cross‐correlation based event detection technique and its application to a surface passive data set

In studies on heavy oil, shale reservoirs, tight gas and enhanced geothermal systems, the use of surface passive seismic data to monitor induced microseismicity due to the fluid flow in the subsurface is becoming more common. However, in most studies passive seismic records contain days and months of data and manually analysing the data can be expensive and inaccurate. Moreover, in the presence of noise, detecting the arrival of weak microseismic events becomes challenging. Hence, the use of an automated, accurate and computationally fast technique for event detection in passive seismic data is essential. The conventional automatic event identification algorithm computes a running‐window energy ratio of the short‐term average to the long‐term average of the passive seismic data for each trace. We show that for the common case of a low signal‐to‐noise ratio in surface passive records, the conventional method is not sufficiently effective at event identification. Here, we extend the conventional algorithm by introducing a technique that is based on the cross‐correlation of the energy ratios computed by the conventional method. With our technique we can measure the similarities amongst the computed energy ratios at different traces. Our approach is successful at improving the detectability of events with a low signal‐to‐noise ratio that are not detectable with the conventional algorithm. Also, our algorithm has the advantage to identify if an event is common to all stations (a regional event) or to a limited number of stations (a local event). We provide examples of applying our technique to synthetic data and a field surface passive data set recorded at a geothermal site.     

Airborne MSS for land cover classification II

A basic methodology for land cover classification using airborne multispectral scanner (MSS) imagery is outlined. This includes waveband selection and radiometric calibration; correction for scan angle and atmosphere; training and classification and accuracy assessment. Refinements to this basic methodology include per‐field sampling and the addition of low‐pass filtering, image texture, prior probabilities and two dates of imagery.

For a study area in upland England, eight land covers were classified with a mean accuracy of 52.6 percent using the basic methodology. This was increased to 79.0 percent by using a suitability refined methodology. Per‐field sampling accounted for the largest proportion of this increase.     

U–Pb and Th–Pb dating of apatite by LA-ICPMS

Apatite is a common U- and Th-bearing accessory mineral in igneous and metamorphic rocks, and a minor but widespread detrital component in clastic sedimentary rocks. U–Pb and Th–Pb dating of apatite has potential application in sedimentary provenance studies, as it likely represents first cycle detritus compared to the polycyclic behavior of zircon. However, low U, Th and radiogenic Pb concentrations, elevated common Pb and the lack of a U–Th–Pb apatite standard remain significant challenges in dating apatite by LA-ICPMS, and consequently in developing the chronometer as a provenance tool.This study has determined U–Pb and Th–Pb ages for seven well known apatite occurrences (Durango, Emerald Lake, Kovdor, Mineville, Mud Tank, Otter Lake and Slyudyanka) by LA-ICPMS. Analytical procedures involved rastering a 10 μm spot over a 40 × 40 μm square to a depth of 10 μm using a Geolas 193 nm ArF excimer laser coupled to a Thermo ElementXR single-collector ICPMS. These raster conditions minimized laser-induced inter-element fractionation, which was corrected for using the back-calculated intercept of the time-resolved signal. A Tl–U–Bi–Np tracer solution was aspirated with the sample into the plasma to correct for instrument mass bias. External standards (Ple?ovice and 91500 zircon, NIST SRM 610 and 612 silicate glasses and STDP5 phosphate glass) along with Kovdor apatite were analyzed to monitor U–Pb, Th–Pb, U–Th and Pb–Pb ratiosCommon Pb correction employed the ²⁰⁷Pb method, and also a ²⁰⁸Pb correction method for samples with low Th/U. The ²⁰⁷Pb and ²⁰⁸Pb corrections employed either the initial Pb isotopic composition or the Stacey and Kramers model and propagated conservative uncertainties in the initial Pb isotopic composition. Common Pb correction using the Stacey and Kramers (1975) model employed an initial Pb isotopic composition calculated from either the estimated U–Pb age of the sample or an iterative approach. The age difference between these two methods is typically less than 2%, suggesting that the iterative approach works well for samples where there are no constraints on the initial Pb composition, such as a detrital sample. No ²⁰⁴Pb correction was undertaken because of low ²⁰⁴Pb counts on single collector instruments and ²⁰⁴Pb interference by ²⁰⁴Hg in the argon gas supply.Age calculations employed between 11 and 33 analyses per sample and used a weighted average of the common Pb-corrected ages, a Tera–Wasserburg Concordia intercept age and a Tera–Wasserburg Concordia intercept age anchored through common Pb. The samples in general yield ages consistent (at the 2σ level) with independent estimates of the U–Pb apatite age, which demonstrates the suitability of the analytical protocol employed. Weighted mean age uncertainties are as low as 1–2% for U- and/or Th-rich Palaeozoic–Neoproterozoic samples; the uncertainty on the youngest sample, the Cenozoic (31.44 Ma) Durango apatite, ranges from 3.7–7.6% according to the common Pb correction method employed. The accurate and relatively precise common Pb-corrected ages demonstrate the U–Pb and Th–Pb apatite chronometers are suitable as sedimentary provenance tools. The Kovdor carbonatite apatite is recommended as a potential U–Pb and Th–Pb apatite standard as it yields precise and reproducible ²⁰⁷Pb-corrected, ²³²Th–²⁰⁸Pb, and common Pb-anchored Tera–Wasserburg Concordia intercept ages.     

Integration of the life of field seismic data with the reservoir model at the Valhall Field‡

The frequent time‐lapse observations from the life of field seismic system across the Valhall field provide a wealth of information. The responses from the production and injection wells can be observed through time‐shift and amplitude changes. These observations can be compared to modelled synthetic seismic responses from a reservoir simulation model of the Valhall Field. The observed differences between the observations and the modelling are used to update and improve the history match of the reservoir model. The uncertainty of the resulting model is reduced and a more confident prediction of future reservoir performance is provided. A workflow is presented to convert the reservoir model to a synthetic seismic response and compare the results to the observed time‐lapse responses for any time range and area of interest. Correlation based match quality factors are calculated to quantify the visual differences. This match quality factor allows us to quantitatively compare alternative reservoir models to help identify the parameters that best match the seismic observations. Three different case studies are shown where this workflow has helped to reduce the uncertainty range associated with specific reservoir parameters. By updating various reservoir model parameters we have been able to improve the match to the observations and thereby improve the overall reservoir model predictability. The examples show positive results in a range of different reservoir modelling issues, which indicates the flexibility of this workflow and the ability to have an impact in most reservoir modelling challenges.     

Three‐dimensional gravity‐current flow within a subaqueous bend: Spatial evolution and force balance variations

The nature of three‐dimensional flow in submarine channel bends is poorly understood, largely due to the absence of detailed data from natural channels. Herein, data from density‐driven flows in a large reservoir on the Huanghe (Yellow) River are presented showing the spatio‐temporal variation of flow around a subaqueous bend. The data demonstrate for the first time that reversed helical flow, relative to that found in river channel bends, can occur from the centrifugal forcing of flow, even when the Coriolis force acts in the opposite direction. The data also suggest that reversed helical flow fields in submarine channels may be more frequent than currently estimated, notably for bends where Coriolis and centrifugal forces combine in the same direction. In addition, this study provides the first field evidence suggesting that sinuous submarine channels can exhibit an asymmetry in helical flow orientation between left and right‐turning bends, which will have major implications for the morphodynamics of submarine channels, their resultant patterns of sedimentation and, ultimately, the distribution of depositional units across submarine fan systems.     

Neoarchean suprasubduction zone arc magmatism in southern India: Geochemistry,zircon U-Pb geochronology and Hf isotopes of the Sittampundi Anorthosite Complex

The Sittampundi Anorthosite Complex (SAC) in southern India is one of the well exposed Archean layered anorthosite-gabbro-ultramafic rock associations. Here we present high precision geochemical data for the various units of SAC, coupled with zircon U-Pb geochronology and Hf isotopic data for the anorthosite. The zircon ages define two populations, the older yield a concordia age of 2541 ± 13 Ma, which is interpreted as the best estimate of the magmatic crystallization age for the Sittampundi anorthosite. A high-grade metamorphic event at 2461 ± 15 Ma is suggested by the upper intercept age of the younger zircon population. A Neoproterozoic event at 715 ± 180 Ma resulted in Pb loss from some of the metamorphic zircons. The magmatic age of the anorthosite correlates well with the timing of crystallization of the arc-related ~ 2530 Ma magmatic charnockites in the adjacent Salem Block, while the metamorphic age is synchronous with the regional metamorphic event. The geochemical data suggest that the rocks were derived from a depleted mantle source. Sub-arc mantle metasomatism of slab derived fluids and subsequent partial melting produced hydrous, aluminous basalt magma. The magma fractionated at depth to produce a variety of high-alumina basalt compositions, from which the anorthositic complex with its chromite-rich and amphibole-rich layers formed as cumulates within the magma chamber of a supra-subduction zone arc. The coherent initial¹⁷⁶Hf/¹⁷⁷Hf ratios and positive εHf values (1.7 – 4.5) of the magmatic zircons in the anorthosite are consistent with derivation of a rather homogeneous juvenile parent magma from a depleted mantle source. Our study further confirms that the southern part of the Dharwar Craton was an active convergent margin during the Neoarchean with the generation and emplacement of suprasubduction zone arc magmas which played a significant role in continental growth.     

Geostatistical characterization of soil moisture and chloride distribution in deep vadose profiles of the Badain Jaran Desert,Northwestern China

Geostatistical analysis of soil moisture and ion accumulation can provide a deep understanding of aspects of a region’s hydrology such as recharge rates and flow paths. By combining this approach with detrended canonical correspondence analysis, the present study investigated the vertical and spatial variation of soil moisture and chloride (Cl^?) in vadose profiles of megadunes near lakes in the southeastern Badain Jaran Desert, northwestern China. Soil moisture varied vertically within the dune profiles following exponential, logarithmic, or log-normal distributions at different sites as a result of regulation by different local controls. The main local controls for Cl^? concentration were slope and orientation of the dune. Cl^? was not affected by local controls at depths >5 m. Sills of spatial variation in soil moisture and Cl^? in the variogram analysis were highly sensitive to local controls in the top 4 m of the profiles, and the values were three times those at depths >4 m. The lag for soil moisture was insensitive to local controls, and was 50–60 m in the top 4 m of the profiles and 20 m in deeper layers, whereas the lag for Cl^? was about 20 m at all depths. The results demonstrate that understanding recharge rates using factors such as Cl^? accumulation requires a careful survey to ensure that the local controls that regulate these factors are adequately accounted for.     

The Oligocene Lund Tuff, Great Basin, USA: a very large volume monotonous intermediate

Unusual monotonous intermediate ignimbrites consist of phenocryst-rich dacite that occurs as very large volume (>1000 km³) deposits that lack systematic compositional zonation, comagmatic rhyolite precursors, and underlying plinian beds. They are distinct from countless, usually smaller volume, zoned rhyolite–dacite–andesite deposits that are conventionally believed to have erupted from magma chambers in which thermal and compositional gradients were established because of sidewall crystallization and associated convective fractionation. Despite their great volume, or because of it, monotonous intermediates have received little attention. Documentation of the stratigraphy, composition, and geologic setting of the Lund Tuff – one of four monotonous intermediate tuffs in the middle-Tertiary Great Basin ignimbrite province – provides insight into its unusual origin and, by implication, the origin of other similar monotonous intermediates.The Lund Tuff is a single cooling unit with normal magnetic polarity whose volume likely exceeded 3000 km³. It was emplaced 29.02±0.04 Ma in and around the coeval White Rock caldera which has an unextended north–south diameter of about 50 km. The tuff is monotonous in that its phenocryst assemblage is virtually uniform throughout the deposit: plagioclase>quartz≈hornblende>biotite>Fe–Ti oxides≈sanidine>titanite, zircon, and apatite. However, ratios of phenocrysts vary by as much as an order of magnitude in a manner consistent with progressive crystallization in the pre-eruption chamber. A significant range in whole-rock chemical composition (e.g., 63–71 wt% SiO₂) is poorly correlated with phenocryst abundance.These compositional attributes cannot have been caused wholly by winnowing of glass from phenocrysts during eruption, as has been suggested for the monotonous intermediate Fish Canyon Tuff. Pumice fragments are also crystal-rich, and chemically and mineralogically indistinguishable from bulk tuff. We postulate that convective mixing in a sill-like magma chamber precluded development of a zoned chamber with a rhyolitic top or of a zoned pyroclastic deposit. Chemical variations in the Lund Tuff are consistent with equilibrium crystallization of a parental dacitic magma followed by eruptive mixing of compositionally diverse crystals and high-silica rhyolite vitroclasts during evacuation and emplacement. This model contrasts with the more systematic withdrawal from a bottle-shaped chamber in which sidewall crystallization creates a marked vertical compositional gradient and a substantial volume of capping-evolved rhyolite magma. Eruption at exceptionally high discharge rates precluded development of an underlying plinian deposit.The generation of the monotonous intermediate Lund magma and others like it in the middle Tertiary of the western USA reflects an unusually high flux of mantle-derived mafic magma into unusually thick and warm crust above a subducting slab of oceanic lithosphere.