The geology of kimberlite pipes of the Ekati property, Northwest Territories, Canada

This paper reviews key characteristics of kimberlites on the Ekati property, NWT, Canada. To date 150 kimberlites have been discovered on the property, five of which are mined for diamonds. The kimberlites intrude Archean basement of the central Slave craton. Numerous Proterozoic diabase dykes intrude the area. The Precambrian rocks are overlain by Quaternary glacial sediments. No Phanerozoic rocks are present. However, mudstone xenoliths and disaggregated sediment within the kimberlites indicate that late-Cretaceous and Tertiary cover (likely <200 m) was present at the time of emplacement. The Ekati kimberlites range in age from 45 to 75 Ma. They are mostly small pipe-like bodies (surface area mostly <3 ha but up to 20 ha) that typically extend to projected depths of 400–600 m below current surface. Pipe morphologies are strongly controlled by joints and faults. The kimberlites consist primarily of variably bedded volcaniclastic kimberlite (VK). This is dominated by juvenile constituents (olivine and lesser kimberlitic ash) and variable amounts of exotic sediment (primarily mud), with minor amounts of xenolithic wall-rock material (generally <5%). Kimberlite types include: mud-rich resedimented VK (mRVK); olivine-rich VK (oVK); sedimentary kimberlite; primary VK (PVK); tuffisitic kimberlite (TK) and magmatic kimberlite (MK). The presence and arrangement of these rock types varies widely. The majority of bodies are dominated by oVK and mRVK, but PVK is prominent in the lower portions of certain kimberlites. TK is rare. MK occurs primarily as precursor dykes but, in a few cases, forms pipe-filling intrusions. The internal geology of the kimberlites ranges from simple single-phase pipes (RVK or MK), to complex bodies with multiple, distinct units of VK. The latter include pipes infilled with steep, irregular VK blocks/wedges and at least one case in which the pipe is occupied by well-defined sub-horizontal VK phases, including a unique, 100-m-thick graded sequence. The whole-rock compositions of VK samples suggest significant loss of kimberlitic fines during eruption followed by variable dilution by surface sediment and concurrent incorporation of kimberlitic ash. Diamond distribution within the kimberlites reflects the amount and nature of mantle material sampled by individual kimberlite phases, but is modified considerably by eruption and depositional processes. The characteristics of the Ekati kimberlites are consistent with a two-stage emplacement process: (1) explosive eruption/s causing vent clearing followed by formation of a significant tephra rim/cone of highly fragmented, olivine-enriched juvenile material with varying amounts of kimberlitic ash and surface sediments (predominantly mud); and (2) infilling of the vent by direct deposition from the eruption column and/or resedimentation of crater rim materials. The presence of less fragmented, juvenile-rich PVK in the lower portions of certain pipes and the intrusion of large volumes of MK to shallow levels in some bodies suggest emplacement of relatively volatile-depleted, less explosive kimberlite in the later stages of pipe formation and/or filling. Explosive devolatilisation of CO₂-rich kimberlite magma is interpreted to have been the dominant eruption mechanism, but phreatomagmatism is thought to have played a role and, in certain cases, may have been dominant.     

Assessing the effectiveness of non-structural flood management measures in the Thames Estuary under conditions of socio-economic and environmental change

Risk analysis and appraisal of the benefits of structural flood risk management measures such as embankments is well established. Here, a method to quantify, over extended timescales, the effectiveness of non-structural measures such as land use spatial planning, insurance and flood resilient construction is presented. The integrated approach couples socio-economic and climate change scenarios with long term land use modelling and flood risk analysis to generate maps and time series of expected annual damages. The analysis has been applied on a case study in the Thames Estuary in the UK. Stakeholders helped develop a number of scenarios that might lead to substantial changes in existing planning and insurance policies in the UK. The effectiveness of these changes was analysed and showed the substantial benefits in terms of reduction of future flood risks that are achievable with changes in planning policy, financial incentives and resilient property construction in the floodplain. Moreover, the reward can be increased through earlier action. Subsequently, the benefits of a range of policies are explored under the UK Foresight socio-economic scenarios. Different structural and non structural flood management interventions are tested and the results demonstrate that despite the potential for large increases in flood risk in the Thames Estuary, in all scenarios substantial flood risk reductions are possible. The effectiveness of non-structural measures is however sensitive to socio-economic changes and governance arrangements. The analysis described here will help to identify portfolios of non-structural and structural options that are robust to uncertainties.     

Changing patterns of Pennsylvanian coal-swamp vegetation and implications of climatic control on coal occurrence

Improved regional and interregional stratigraphic correlations of Pennsylvanian strata permit comparisons of vegetational changes in Euramerican coal swamps. The coal-swamp vegetation is known directly from in situ coal-ball peat deposits from more than 65 coals in the United States and Europe. Interpretations of coal-swamp floras on the basis of coal-ball peat studies are extended to broader regional and stratigraphic patterns by use of coal palynology. Objectives of the quantitative analyses of the vegetation in relation to coal are to determine the botanical constituents at the peat stage and their environmental implications for plant growth and peat accumulation. Morphological and paleoecological analyses provide a basis for deducing freshwater regimes of coal swamps.Changes in composition of Pennsylvanian coal-swamp vegetation are quire similar from one paralic coal region to another and show synchrony that is attributable to climate. Paleobotany and paleogeography of the Euramerican province indicate a moist tropical paleoclimate. Rainfall, runoff and evapotranspiration were the variable climatic controls in the distribution of coal-swamp vegetation, peat accumulation and coal resources. In relative terms of climatic wetness the Pennsylvanian Period is divisible into five intervals, which include two relatively drier intervals that developed during the Lower-Middle and Middle-Upper Pennsylvanian transitions. The climate during Early Pennsylvanian time was moderately wet and the median in moisture availability. Early Middle Pennsylvanian was drier, probably seasonally dry-wet; late Middle Pennsylvanian was the wettest in the Midcontinent; early Late Pennsylvanian was the driest; and late Late Pennsylvanian was probably the wettest in the Dunkard Basin. The five climatic intervals represent a general means of dividing coal resources within each region into groups with similar botanical constituents and environments of peat accumulation. Regional differences in basinal geology and climate were significant variables, but the synchronous control of paleoclimate was of primary importance.     

Clinoform geometry,geomorphology, facies character and stratigraphic architecture of a sand‐rich subaqueous delta: Jurassic Sognefjord Formation,offshore Norway

The integration of core sedimentology, seismic stratigraphy and seismic geomorphology has enabled interpretation of delta‐scale (i.e. tens of metres high) subaqueous clinoforms in the upper Jurassic Sognefjord Formation of the Troll Field. Mud‐prone subaqueous deltas characterized by a compound clinoform morphology and sandy delta‐scale subaqueous clinoforms are common in recent tide‐influenced, wave‐influenced and current‐influenced settings, but ancient examples are virtually unknown. The data presented help to fully comprehend the criteria for the recognition of other ancient delta‐scale subaqueous clinoforms, as well as refining the depositional model of the reservoir in the super‐giant Troll hydrocarbon field. Two 10 to 60 m thick, overall coarsening‐upward packages are distinguished in the lower Sognefjord Formation. Progressively higher energy, wave‐dominated or current‐dominated facies occur from the base to the top of each package. Each package corresponds to a set of seismically resolved, westerly dipping clinoforms, the bounding surfaces of which form the seismic ‘envelope’ of a clinoform set and the major marine flooding surfaces recognized in cores. The packages thicken westwards, until they reach a maximum where the clinoform ‘envelope’ rolls over to define a topset–foreset–toeset geometry. All clinoforms are consistently oriented sub‐parallel to the edge of the Horda Platform (N005–N030). In the eastern half of the field, individual foresets are relatively gently dipping (1° to 6°) and bound thin (10 to 30 m) clinothems. Core data indicate that these proximal clinothems are dominated by fine‐grained, hummocky cross‐stratified sandstones. Towards the west, clinoforms gradually become steeper (5° to 14°) and bound thicker (15 to 60 m) clinothems that comprise medium‐grained, cross‐bedded sandstones. Topsets are consistently well‐developed, except in the westernmost area. No seismic or sedimentological evidence of subaerial exposure is observed. Deposition created fully subaqueous, near‐linear clinoforms that prograded westwards across the Horda Platform. Subaqueous clinoforms were probably fed by a river outlet in the north‐east and sculpted by the action of currents sub‐parallel to the clinoform strike.     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

Lu–Hf and U–Pb isotope systematics of zircons from the Storgangen intrusion, Rogaland Intrusive Complex, SW Norway: implications for the composition and evolution of Precambrian lower crust in the Baltic Shield

The Storgangen orebody is a concordantly layered, sill-like body of ilmenite-rich norite, intruding anorthosites of the Rogaland Intrusive Complex (RIC), SW Norway. 17 zircon grains were separated from ca. 5 kg of sand-size flotation waste collected from the on-site repository from ilmenite mining. These zircons were analysed for major and trace elements by electron microprobe, and for U–Pb and Lu–Hf isotopes by laser ablation microprobe plasma source mass spectrometry. Eight of the zircons define a well-constrained (MSWD=0.37) concordant population with an age of 949±7 Ma, which is significantly older than the 920–930 Ma ages previously reported for zircon inclusions in orthopyroxene megacrysts from the RIC. The remaining zircons, interpreted as inherited grains, show a range of ²⁰⁷Pb/²⁰⁶Pb ages up to 1407±14 Ma, with an upper intercept age at ca. 1520 Ma. The concordant zircons have similar trace element patterns, and a mean initial Hf isotope composition of ¹⁷⁶Hf/¹⁷⁷Hf_{949 Ma}=0.28223±5 (_Hf=+2±2). This is similar to the Hf-isotope composition of zircons in a range of post-tectonic Sveconorwegian granites from South Norway, and slightly more radiogenic than expected for mid-Proterozoic juvenile crust. The older, inherited zircons show Lu–Hf crustal residence ages in the range 1.85–2.04 Ga. One (undated) zircon plots well within the field of Hf isotope evolution of Paleoproterozoic rocks of the Baltic Shield. These findings indicate the presence of Paleoproterozoic components in the deep crust of the Rogaland area, but do not demonstrate that such rocks, or a Sveconorwegian mantle-derived component, contributed significantly to the petrogenesis of the RIC. If the parent magma was derived from a homogeneous, lower crustal mafic granulite source, the lower crustal protolith must be at least 1.5 Ga old, and it must have an elevated Rb/Sr ratio. This component would be indistinguishable in Sr, Nd and Hf isotopes from some intermediate mixtures between Sveconorwegian mantle and Paleoprotoerzoic felsic crust, but it cannot account for the initial ¹⁴³Nd/¹⁴⁴Nd of the most primitive, late Sveconorwegian granite in the region, without the addition of mantle-derived material.     

Vertical Stacking Statistics of Multi-facies Object-Based Models

Mathematical Geosciences - Equations describing facies proportions and amalgamation ratios are derived for randomly placed objects belonging to two or three foreground facies embedded in a...