Recognition of down-valley translation in tidally influenced meandering fluvial deposits,Athabasca Oil Sands (Cretaceous), Alberta,Canada

The vast Athabasca Oil Sands of Alberta, Canada has an estimated resource of more than 1.7 trillion barrels of bitumen in-place, the majority of which is hosted in the Lower Cretaceous McMurray Formation. Despite its economical significance the depositional environment of the formation, and particularly the middle part that is the primary reservoir in most areas, is still a matter of considerable debate. These strata of interest are widely known to comprise fluvial and estuarine point bar deposits that were subject to varying degrees of marine influence. The orientation of point bar strata from the formation is tabulated and the majority is observed to dip paleo-basinward, to the north. This observation has proven difficult to explain by geologists attempting to build predictive models for the formation over the last two decades. However, the basinward-dipping point bar strata can be attributed to widespread down-valley translation of point bars in confined north-south oriented valleys, which have previously been delineated in the region. Differentiating the deposits of lateral point bar migration and down-valley translation in the rock record is not possible at the facies scale, and thus translated point bars have not been previously recognized or interpreted in the rock record, despite their prevalence in numerous modern fluvial valley systems. Their identification in the McMurray Formation has important implications for the delineation and development of Canada’s economically significant oil sands resources.     

Numerical simulation of surf–swash zone motions and turbulent flow

A two-dimensional numerical model was presented for the simulation of wave breaking, runup and turbulence in the surf and swash zones. The main components of the model are the Reynolds-Averaged Navier–Stokes equations describing the average motion of a turbulent flow, a k–ε turbulence closure model describing the transformation and dissipation processes of turbulence and a volume of fluid technique for tracking the free surface motion. Nearshore wave evolution on a sloping bed, the velocity field and other wave characteristics were investigated. First, the results of the model were compared with experimental results for different surf zone hydrodynamic conditions. Spilling and plunging breakers were simulated and the numerical model investigated for different wave parameters. The turbulence field was also considered and the spatial and time-dependent variations of turbulence parameters were discussed. In the next stage of the study, numerical results were compared with two sets of experimental data in the swash zone. Generally, there is good agreement except for turbulence predictions near the breaking point where the model does not represent well the physical processes. On the other hand, turbulence predictions were found to be excellent for the swash zone. The model provides a precise and efficient tool for the simulation of the flow field and wave transformations in the nearshore, especially in the swash zone. The numerical model can simulate the surface elevation of the vertical shoreline excursion on sloping beaches, while swash–swash interactions within the swash zone are accounted for.     

An evaluation of Radarsat-2 individual and combined image dates for land use/cover mapping

Various land use/cover types exhibit seasonal characteristics which can be captured in remotely sensed imagery. This study examined how different seasons of Radarsat-2 data influence land use/cover classification accuracies for two study sites. Two dates of Radarsat-2 C-band quad-polarised images were obtained for Washington, DC, USA and Wad Madani, Sudan. Spectral signatures were extracted and used with a maximum likelihood decision rule for classification and thematic accuracies were then determined. Both despeckled radar and derived texture measures were examined. Thematic accuracies for the two despeckled image dates were similar with a difference of 3% for Washington and 6% for Sudan. Merging the despeckled images for both seasons increased overall accuracy by 2% for Washington and 9% for Sudan. Further combining the original radar for both seasons with derived texture measures increased overall accuracies by 9% for Washington and 16% for Sudan for final overall accuracy values of 73 and 82%.     

Imputing forest carbon stock estimates from inventory plots to a nationally continuous coverage

The U.S. has been providing national-scale estimates of forest carbon (C) stocks and stock change to meet United Nations Framework Convention on Climate Change (UNFCCC) reporting requirements for years. Although these currently are provided as national estimates by pool and year to meet greenhouse gas monitoring requirements, there is growing need to disaggregate these estimates to finer scales to enable strategic forest management and monitoring activities focused on various ecosystem services such as C storage enhancement. Through application of a nearest-neighbor imputation approach, spatially extant estimates of forest C density were developed for the conterminous U.S. using the U.S.’s annual forest inventory. Results suggest that an existing forest inventory plot imputation approach can be readily modified to provide raster maps of C density across a range of pools (e.g., live tree to soil organic carbon) and spatial scales (e.g., sub-county to biome). Comparisons among imputed maps indicate strong regional differences across C pools. The C density of pools closely related to detrital input (e.g., dead wood) is often highest in forests suffering from recent mortality events such as those in the northern Rocky Mountains (e.g., beetle infestations). In contrast, live tree carbon density is often highest on the highest quality forest sites such as those found in the Pacific Northwest. Validation results suggest strong agreement between the estimates produced from the forest inventory plots and those from the imputed maps, particularly when the C pool is closely associated with the imputation model (e.g., aboveground live biomass and live tree basal area), with weaker agreement for detrital pools (e.g., standing dead trees). Forest inventory imputed plot maps provide an efficient and flexible approach to monitoring diverse C pools at national (e.g., UNFCCC) and regional scales (e.g., Reducing Emissions from Deforestation and Forest Degradation projects) while allowing timely incorporation of empirical data (e.g., annual forest inventory).     

Geochemistry of highly-undersaturated ocean island basalt suites from the South Atlantic Ocean: Fernando de Noronha and Trindade islands

The volcanic rocks of the South Atlantic Ocean islands of Fernando de Noronha and Trindade comprise a diverse magmatic series ranging from nephelinites and basanites to phonolites and, on Fernando de Noronha, trachytes. All rock types are highly silica undersaturated with the exception of Fernando de Noronha trachytes_, and have high abundances of incompatible trace elements and strongly LREE (light rare earth element)-enriched REE patterns. Crystal fractionation of parental basanitic magmas produced evolved phonolites and trachytes which display severe trace-element fractionation, even among trace elements (Nb, Ta, Zr, Hf) which normally behave highly incompatibly during crystallisation of alkaline magmas. Moderately to highly evolved compositions develop strongly MREE (middle REE)-depleted REE patterns, and become increasingly depleted in elements such as Nb and, in particular, Ta. Ratios of Nb/Ta and Zr/Hf are highly fractionated in phonolites (60–65, 64–77 respectively in Fernando de Noronha phonolites) compared to ratios in basanites (14, 45 respectively). The compatibility of Nb, Ta, and the REE, and the strong fractionation of Nb/Ta and Zr/Hf ratios and the MREE, during crystallisation from basanite to phonolite are attributable to the crystallisation of small amounts (<5%) of sphene. Trace-element behaviour is relatively insensitive to the major phenocryst phases, and is controlled by minor phases in highly undersaturated alkaline suites. Incompatible trace-element ratios (e.g. La/Nb, Th/Ta) in nephelinites and basanites from Fernando de Noronha and Trindade are generally comparable with those in basaltic and hawaiitic OIB (ocean island basalt) lavas from other South Atlantic islands, but are distinct from those of Gough and Tristan da Cunha OIB. The mantle source for the highly undersaturated volcanism on Fernando de Noronha and Trindade is similar in trace-element characteristics to the typical OIB source which produces alkaline lavas with significant relative enrichment in Nb and Ta compared to other trace elements (as expressed by low La/Nb, Ba/Nb and Th/Ta ratios). The highly undersaturated nature of the magmas and the slight fractionation of some incompatible-element ratios (elevated Ba/Nb, Ba/Rb, Ba/Th etc.) is consistent with a smaller degree of melting of a typical OIB source, but with residual phlogopite in the source to account for significant K depletion and LIL-element fractionation.     

Fracture analysis near the mid-ocean plate boundary, Reykjavik-Hvalfjördur area, Iceland

The geometry and thermal history of fractures have been determined at 59 stations from Reykjavik to Hvalfjördur in southwestern Iceland. The data provide information on crustal stress regimes in the vicinity of mid-ocean ridges.Two major, generalized fracture orientations are present

1. (1) a northeast system, trend 010°–030°, except on Akranes where the orientation is 040°–060°

2. (2) a broad east—west system containing one or more sets with strike between 070°–130°.

Thermal history of the host rock and fractures was determined from secondary minerals in vugs and fractures. The thermal history indicates that the northeast fracture set opened while the area was within the relatively hot axial zone of active volcanism and rifting. Some of the east—west trending fractures also opened at this time but many formed later, after the area had begun to cool and drift from the active zone.The northeast fracture set is essentially parallel to the trend of dikes and normal faults in southwestern Iceland. They have been interpreted as extension fractures (resulting in about 0.4% maximum extension) forming generally from the same stress field associated with normal faulting and dike injection in the active zone. Fracturing in an east-west direction (estimated 0.1% maximum extension), mainly near the edge and outside the active zone, indicates a reorientation of this stress field. The dominant mechanism related to the origin of the east—west fractures may be thermoelastic stresses arising from axial and basal accretion and cooling of lithospheric plates.Both fracture systems are inferred to have formed, in the Griffiths idealization, under nearly biaxial effective compressive loading on the order of 200 bar. The discrepancy between this value and the kilobar-order strengths of short-time laboratory tests reflects such factors as high temperature stress corrosion and fatigue. Fracture propagation is assumed to have been stable, but governed primarily by lateral load-diminishing mechanisms rather than by progressive loading. These relaxation mechanisms may have been episodic (northeast-system fissure swarm activity) or steady-state (thermoelastic contraction) in time.     

The MARID (mica-amphibole-rutile-ilmenite-diopside) suite of xenoliths in kimberlite

Within the ‘glimmerite’ nodules occurring within kimberlite pipes we recognize the MARID suite consisting of varying proportions of mica, amphibole, rutile, ilmenite and diopside. Banding of some specimens is interpreted as cumulate layering. All specimens were deformed either before incorporation into the host kimberlite or during intrusion. Compared with minerals in peridotite xenoliths, the MARID ones are lower in Al₂O₃ and Cr₂O₃, but richer in total iron. The MARID micas, amphiboles, diopsides, ilmenites and probably rutiles contain substantial Fe₂O₃ indicative of oxidizing conditions. The amphibole is potassic richterite. Micas of the megacryst suite in kimberlite have less total iron and Fe₂O₃ than micas of the MARID suite. We suggest that the rocks of the MARID suite crystallized under oxidizing conditions from a magma, chemically similar to kimberlite, within the higher parts of the upper mantle: the presence of amphibole restricts the depth to less than ~ 100 km. A xenolith containing olivine and orthopyroxene as well as minerals similar to but not the same compositionally as MARID-types is interpreted as a metasomite, possibly representing wall-rock of a magma body from which MARID-suite rocks crystallized.