Application of grouped detrital zircon analyses to determine provenance and closely approximate true depositional age: Early Cretaceous McMurray-Clearwater succession,Canada

In the Lower CretaceousMcMurray-Clearwater succession of the intracontinental Alberta Foreland Basin,Canada,detrital zircon U-Pb geochronology samples(referred to herein as DZ samples)have been used to interpret the strata as representing a paleo-continental-scale drainage system.However,the majority of DZ samples are relatively small(n≈90–100),and syndepositional DZ(i.e.,crystallization age<5 Ma older than depositional age)are rare.This has forced a reliance on dinocysts with long stratigraphic ranges to chronostratigraphically subdivide the McMurray-Clearwater succession rather than employing maximum depositional ages(MDAs)derived from DZ samples.Herein,43 DZ samples(taken from20 subsurface cores)are assigned to 1 of 5 stratigraphic intervals,and in each stratigraphic interval all associated DZ samples are combined to produce a grouped DZ sample.Analysis and comparison of individual and grouped DZ samples are used to(1)assess variability in provenance through time and space,and(2)assess the accuracy of chronostratigraphically subdividing the succession using MDAs.Along the main paleo-drainage axis,a comparison of dissimilarity between DZ samples from the same stratigraphic interval,as well as between stratigraphic intervals,reveals increasing average dissimilarity between individual DZ samples and their respective grouped DZ samplewith increasing spatial separation of samples.These data indicate that in the McMurray Depocenter some sediment is sourced from local tributaries,leading to geographical provenance variability.CalculatedMDAs for all 43 DZ samples and groupedMDAs(gMDAs)for the 5 grouped DZ samples are compared to an ash-derived absolute age and existing biostratigraphy.In theMcMurray Formation,comparison of MDAs to gMDAs shows that in basins with rare syndepositional DZ,the gMDA method improved depositional age estimates by transforming low-confidence MDAs(e.g.,youngest single grains)into high-confidence(multi-grain)gMDAs.In the Clearwater Formationwhere syndepositional DZ are plentiful(i.e.,>5%of the total DZ population),calculating maximumlikelihood ages fromgrouped DZ samples avoids negatively biased(i.e.,too young)MDAs.We suggest grouped DZ samples and the gMDA method be used in systems with multiple DZ samples from a well-defined stratigraphic interval as a means of assessing variability in provenancewithin a depositional system and for improving estimates of depositional ages using DZ.     

Modeling the kinetics of silica nanocolloid formation and precipitation in geologically relevant aqueous solutions

The kinetics of the formation and precipitation of nanocolloidal silica from geologically relevant aqueous solutions is investigated. Changes in monomeric (SiO_2(mono)), nanocolloidal (SiO_2(nano)) and precipitated silica (SiO_2(ppt)) concentrations in aqueous solutions from pH 3 to 7, ionic strengths (IS) of 0.01 and 0.24 molal, and initial SiO₂ concentrations of 20.8, 12.5 and 4.2 mmolal (reported in [Icopini, G.A., Brantley, S.L., Heaney, P.J., 2005. Kinetics of silica oligomerization and nanocolloid formation as a function of pH and ionic strength at 25 °C. Geochim. Cosmochim. Acta69(2), 293-303.]) were fit using two kinetic models. The first model, termed the concentration model, is taken from Icopini et al. (2005) and assumes that the rate of change of SiO_2(mono) as a function of time has a fourth-order dependence on the concentration of SiO_2(mono) in solution. The second model, termed the supersaturation model, incorporates the equilibrium concentration of amorphous silica and predicts that polymerization will be a function of the degree of silica supersaturation in solution with respect to amorphous silica. While both models generally predicted similar rate constants for a given set of experimental conditions, the supersaturation model described the long-term equilibrium behavior of the SiO_2(mono) fraction more accurately, resulting in significantly better fits of the monomeric data. No difference was seen between the model fits of the nanocolloidal silica fraction. At lower pH values (3-4), a metastable equilibrium was observed between SiO_2(mono) and SiO_2(nano). This equilibrium SiO_2(mono) concentration was found to be 6 mmolal, or three times the reported solubility of bulk amorphous silica under the experimental conditions studied and corresponds to the predicted solubility of amorphous silica colloids approximately 3 nm in diameter. Atomic force microscopy was used to determine the average size of the primary nanocolloidal particles to be ∼3 nm, which is in direct agreement with the solubility calculations. Larger aggregates of the primary nanocolloids were also observed to range in size from 30 to 40 nm. This work provides the first kinetic models describing the formation and evolution of nanocolloidal silica in environmentally relevant aqueous solutions. Results indicate that nanocolloidal silica is an important species at low pH and neutral pH at low ionic strengths and may play a more important role in geochemical cycles in natural aqueous systems than previously considered.     

Surface-Based Geological Reservoir Modelling Using Grid-Free NURBS Curves and Surfaces

Mathematical Geosciences - Building geometrically realistic representations of geological heterogeneity in reservoir models is a challenging task that is limited by the inflexibility of pre-defined...     

Geomorphology and environmental management of the Yellow River source zone

Geologic,climatic and anthropogenic factors have fashioned the distinctive landscapes and ecosystems of the Yellow River source zone(Brierley and Huang2013;Figure1).Located at thenorth-east corner of the Qinghai-Tibetan Plateau(QTP),at elevations between 2500 and 4500m,     

Analysis of longitudinal profiles along the eastern margin of the Qinghai-Tibetan Plateau

Resulting from the collision of the Eurasian and Indian plates, the Qinghai-Tibetan Plateau is commonly known as the ’roof of the world’. Collectively the Yarlung Tsangpo, Nu, Lancang, Yangtze, Yalong, and Yellow River basins drain the eastern margin of the plateau. In this paper, we utilize Shuttle Radar Topography Mission elevation data to examine morphometric and relief attributes of these basins to reveal insights into tectonic activity and rates of incision. A robust technique using Matlab is proposed to alleviate errors associated with SRTM data in the derivation of river longitudinal profiles. Convex longitudinal profiles are interpreted to be a product of uplift rates that exceed rates of channel incision along the entire margin of the Qinghai-Tibetan Plateau. Highest relief towards the south reflects extensive fluvial incision. High relief is also prominent along major active faults. Erosion patterns are related to distance from knickpoints. Highest rates of erosion and incision are evident towards the south, with decreasing values towards the north, suggesting a link between tectonic activity and erosion.     

Estimating aquifer channel recharge using optical data interpretation

Recharge through intermittent and ephemeral stream channels is believed to be a primary aquifer recharge process in arid and semiarid environments. The intermittent nature of precipitation and flow events in these channels, and their often remote locations, makes direct flow and loss measurements difficult and expensive. Airborne and satellite optical images were interpreted to evaluate aquifer recharge due to stream losses on the Frio River in south-central Texas. Losses in the Frio River are believed to be a major contributor of recharge to the Edwards Aquifer. The results of this work indicate that interpretation of readily available remote sensing optical images can offer important insights into the spatial distribution of aquifer recharge from losing streams. In cases where upstream gauging data are available, simple visual analysis of the length of the flowing reach downstream from the gauging station can be used to estimate channel losses. In the case of the Frio River, the rate of channel loss estimated from the length of the flowing reach at low flows was about half of the loss rate calculated from in-stream gain-loss measurements. Analysis based on water-surface width and channel slope indicated that losses were mainly in a reach downstream of the mapped recharge zone. The analysis based on water-surface width, however, did not indicate that this method could yield accurate estimates of actual flow in pool and riffle streams, such as the Frio River and similar rivers draining the Edwards Plateau.