Late Cretaceous to Eocene geology of the South Central American forearc area (southern Costa Rica and western Panama) : initiation and evolution of an intra-oceanic convergent margin

Résumé : L'arc volcanique du sud de l'Amérique Centrale se situe sur la marge SW de la Plaque Caraïbe, au-dessus des plaques subduites de Cocos et Nazca. Il s'agit de l'un des arcs intra-océaniques les plus étudiés au monde, qui est généralement considéré comme s'étant développé à la fin du Crétacé le long d'un plateau océanique (le Plateau Caraïbe ou CLIP) et se trouvant actuellement dans un régime de subduction érosive. Au cours des dernières décennies, des efforts particuliers ont été faits pour comprendre les processus liés à la subduction sur la base d'études géophysiques et géochimiques. Au sud du Costa Rica et à l'ouest du Panama, des complexes d'accrétions et structures à la base de l'arc volcanique ont été exposés grâce à la subduction de rides asismiques et de failles transformantes. Des affleurements, situés jusqu'à seulement 15 km de la fosse, offrent une possibilité unique de mieux comprendre quelques uns des processus ayant lieu le long de la zone de subduction. Nous présentons de nouvelles contraintes sur l'origine de ces affleurements en alliant une étude de terrain poussée, de nouvelles données géochimiques, sédimentaires et paléontologiques, ainsi que des observations structurales effectuées en télédétection. Une nouvelle stratigraphie tectonique entre le Campanien et l'Éocène est définie pour la région d'avant-arc située entre la Péninsule d'Osa (Costa Rica) et la Péninsule d'Azuero (Panama). Nos résultats montrent que la partie externe de la marge est composée d'un arrangement complexe de roches ignées et de séquences sédimentaires de recouvrement qui comprennent principalement le socle de l'arc, des roches d'arc primitif, des fragments de monts sous-marins accrétés et des mélanges d'accrétion. Des preuves sont données pour le développement de l'arc volcanique du sud de l'Amérique Centrale sur un plateau océanique. Le début de la subduction le long de la marge SW de la Plaque Caraïbe a eu lieu au Campanien et a généré des roches d'arc primitif caractérisées par des affinités géochimiques particulières, globalement intermédiaires entre des affinités de plateau et d'arc insulaire. L'arc était mature au Maastrichtien et formait un isthme essentiellement continu entre l'Amérique du Nord et l'Amérique du Sud. Ceci a permis la migration de faunes terrestres entre les Amériques et pourrait avoir contribué à la crise fin Crétacé -Tertiaire en réduisant les courants océaniques subéquatoriaux entre le Pacifique et l'Atlantique. Plusieurs unités composées de fragments de monts sous-marins accrétés sont définies. La nature et l'arrangement structural de ces unités définissent de nouvelles contraintes sur les modes d'accrétion des monts sous-marins/îles océaniques et sur l'évolution de la marge depuis la formation de la zone de subduction. Entre la fin du Crétacé et l'Éocène moyen, la marge a enregistré plusieurs épisodes ponctuels d'accrétion de monts sous-marins alternant avec de la subduction érosive. A l'Éocène moyen, un événement tectonique régional pourrait avoir causé un fort couplage entre les plaques supérieure et inférieure, menant à des taux plus important d'accrétion de monts sous-marins. Durant cette période, la situation le long de la marge était très semblable à la situation actuelle et caractérisée par la présence de monts sous-marins subductants et l'absence d'accrétion de sédiments. L'enregistrement géologique montre qu'il n'est pas possible d'attribuer une nature érosive ou accrétionnaire à la marge dans le passé ou -par analogie- aujourd'hui, parce que (1) les processus d'accrétion et érosifs varient fortement spatialement et temporellement et (2) il est impossible d'évaluer la quantité exacte de matériel tectoniquement enlevé à la marge depuis le début de la subduction. Au sud du Costa Rica, certains fragments de monts sous-marins accrétés sont représentatifs d'une interaction entre une ride et un point chaud dans le Pacifique au Crétacé terminal/Paléocène. L'existence de ces fragments de monts sous-marins et la morphologie du fond de l'Océan Pacifique indiquent que la formation de la ride de Cocos-Nazca s'est formée au moins ~40 Ma avant l'âge proposé par les modèles tectoniques actuels. Au Panama, nous avons identifié une île océanique d'âge début Éocène qui a été accrétée à l'Éocène moyen. L'accrétion a eu lieu à très faible profondeur par détachement de l'île dans la fosse, et a mené à une exceptionnelle préservation des structures volcaniques. Des affleurement comprenant aussi bien des parties basses et hautes de l'édifice volcanique on été étudiées, depuis la phase sous-marine bouclier jusqu'à la phase subaérienne post-bouclier. La stratigraphie nous a permis de différencier les laves de la phase sous-marine de celles de la phase subaérienne. La composition des laves indique une diminution progressive de l'intensité de la fusion partielle de la source et une diminution de la température des laves produites durant les derniers stades de l'activité volcanique. Nous interprétons ces changements comme étant liés à l'éloignement progressif de l'île océanique de la zone de fusion ou point chaud. Abstract The southern Central American volcanic front lies on the SW edge of the Caribbean Plate, inboard of the subducting Cocos and Nazca Plates. It is one of the most studied intra-oceanic convergent margins around the world, which is generally interpreted to have developed in the late Cretaceous along an oceanic plateau (the Caribbean Large Igneous Province or CLIP) and to be currently undergoing a regime of subduction erosion. In the last decades a particular effort has been made to understand subduction-related processes on the basis of geophysical and geochemical studies. In southern Costa Rica and western Panama accretionary complexes and structures at the base of the volcanic front have been exposed in response to subduction of aseismic ridges and transforms. Onland exposures are located as close as to 15 km from the trench and provide a unique opportunity to better understand some of the processes occurring along the subduction zone. We provide new constraints on the origins of these exposures by integrating a comprehensive field work, new geochemical, sedimentary and paleontological data, as well as structural observations based on remote imaging. A new Campanian to Eocene tectonostratigraphy is defined for the forearc area located between the Osa Peninsula (Costa Rica) and the Azuero Peninsula (Panama). Our results show that the outer margin is composed of a complicated arrangement of igneous complexes and overlapping sedimentary sequences that essentially comprise an arc basement, primitive island-arc rocks, accreted seamount fragments and accretionary mélanges. Evidences are provided for the development of the southern Central American arc on the top an oceanic plateau. The subduction initiation along the SW edge of the Caribbean Plate occurred in the Campanian and led to formation of primitive island-arc rocks characterized by unusual geochemical affinities broadly intermediate between plateau and arc affinities. The arc was mature in the Maastrichtian and was forming a predominantly continuous landbridge between the North and South Americas. This allowed migration of terrestrial fauna between the Americas and may have contributed to the Cretaceous-Tertiary crisis by limiting trans-equatorial oceanic currents between the Pacific and the Atlantic. Several units composed of accreted seamount fragments are defined. The nature of the units and their structural arrangement provide new constraints on the modes of accretion of seamounts/oceanic islands and on the evolution of the margin since subduction initiation. Between the late Cretaceous and the middle Eocene, the margin recorded several local episodes of seamount accretion alternating with tectonic erosion. In the middle Eocene a regional tectonic event may have triggered strong coupling between the overriding and subducting plates, leading to higher rates of seamount accretion. During this period the situation along the margin was very similar to the present and characterized by subducting seamounts and absence of sediment accretion. The geological record shows that it is not possible to ascribe an overall erosive or accretionary nature to the margin in the past and, by analogy, today, because (1) accretionary and erosive processes exhibit significant lateral and temporal variations and (2) it is impossible to estimate the exact amount of material tectonically eroded from the margin since subduction initiation. In southern Costa Rica, accreted seamount fragments point toward a plume-ridge interaction in the Pacific in the late Cretaceous/Paleocene. This occurrence of accreted seamount fragments and morphology of the Pacific Ocean floor is indicative of the formation of the Cocos-Nazca spreading system at least ~40 Ma prior to the age proposed in current tectonic models. In Panama, we identified a remarkably-well preserved early Eocene oceanic island that accreted in the middle Eocene. The accretion probably occurred at very shallow depth by detachment of the island in the trench and led to an exceptional preservation of the volcanic structures. Exposures of both deep and superficial parts of the volcanic edifice have been studied, from the submarine-shield to subaerial-postshield stages. The stratigraphy allowed us to distinguish lavas produced during the submarine and subaerial stages. The lava compositions likely define a progressive diminution of source melting and a decrease in the temperature of erupted melts in the latest stages of volcanic activity. We interpret these changes to primarily reflect the progressive migration of the oceanic island out of the melting region or hotspot.

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

Detection and characterization of hydraulically active fractures in a carbonate aquifer: Results from self-potential, temperature and fluid conductivity logging along a 260-m-deep borehole in the Combioula hydrothermal system in the southwestern Swiss Alps

A geophysical and geochemical study has been conducted in a fractured carbonate aquifer located at Combioula in the southwestern Swiss Alps with the objective to detect and characterize hydraulically active fractures along a 260-m-deep borehole. Hydrochemical analyses, borehole diameter, temperature and fluid electrical conductivity logging data were integrated in order to relate electrokinetic self-potential signals to groundwater flow inside the fracture network. The results show a generally good, albeit locally variable correlation of variations of the self-potential signals with variations in temperature, fluid electrical conductivity and borehole diameter. Together with the hydrochemical evidence, which was found to be critical for the interpretation of the self-potential data, these measurements not only made it possible to detect the hydraulically active fractures but also to characterize them as zones of fluid gain or fluid loss. The results complement the available information from the corresponding litholog and illustrate the potential of electrokinetic self-potential signals in conjunction with temperature, fluid electrical conductivity and hydrochemical analyses for the characterization of fractured aquifers, and thus may offer a perspective for an effective quantitative characterization of this increasingly important class of aquifers and geothermal reservoirs.

Quantifying sediment storage in a high Alpine valley (Turtmanntal, Switzerland)

The determination of sediment storage is a critical parameter in sediment budget analyses. But, in many sediment budget studies the quantification of magnitude and time-scale of sediment storage is still the weakest part and often relies on crude estimations only, especially in large drainage basins (>100km2). We present a new approach to storage quantification in a meso-scale alpine catchment of the Swiss Alps (Turtmann Valley, 110km2). The quantification of depositional volumes was performed by combining geophysical surveys and geographic information system (GIS) modelling techniques. Mean thickness values of each landform type calculated from these data was used to estimate the sediment volume in the hanging valleys and the trough slopes. Sediment volume of the remaining subsystems was determined by modelling an assumed parabolic bedrock surface using digital elevation model (DEM) data. A total sediment volume of 781·3×106?1005·7×106m3 is deposited in the Turtmann Valley. Over 60% of this volume is stored in the 13 hanging valleys. Moraine landforms contain over 60% of the deposits in the hanging valleys followed by sediment stored on slopes (20%) and rock glaciers (15%). For the first time, a detailed quantification of different storage types was achieved in a catchment of this size. Sediment volumes have been used to calculate mean denudation rates for the different processes ranging from 0·1 to 2·6mm/a based on a time span of 10ka. As the quantification approach includes a number of assumptions and various sources of error the values given represent the order of magnitude of sediment storage that has to be expected in a catchment of this size.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

Textural controls on induced polarization and permeability in granular media

The spatial resolution visualized with hydrological models and the conceptualized images of subsurface hydrological processes often exceed resolution of the data collected with classical instrumentation at the field scale. In recent years it was possible to increasingly diminish the inherent gap to information from point like field data through the application of hydrogeophysical methods at field-scale. With regards to all common geophysical exploration techniques, electric and electromagnetic methods have arguably to greatest sensitivity to hydrologically relevant parameters. Of particular interest in this context are induced polarisation (IP) measurements, which essentially constrain the capacity of a probed subsurface region to store an electrical charge. In the absence of metallic conductors the IP- response is largely driven by current conduction along the grain surfaces. This offers the perspective to link such measurements to the characteristics of the solid-fluid-interface and thus, at least in unconsolidated sediments, should allow for first-order estimates of the permeability structure.¦While the IP-effect is well explored through laboratory experiments and in part verified through field data for clay-rich environments, the applicability of IP-based characterizations to clay-poor aquifers is not clear. For example, polarization mechanisms like membrane polarization are not applicable in the rather wide pore-systems of clay free sands, and the direct transposition of Schwarz' theory relating polarization of spheres to the relaxation mechanism of polarized cells to complex natural sediments yields ambiguous results.¦In order to improve our understanding of the structural origins of IP-signals in such environments as well as their correlation with pertinent hydrological parameters, various laboratory measurements have been conducted. We consider saturated quartz samples with a grain size spectrum varying from fine sand to fine gravel, that is grain diameters between 0,09 and 5,6 mm, as well as corresponding pertinent mixtures which can be regarded as proxies for widespread alluvial deposits. The pore space characteristics are altered by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples, while keeping any electrochemical variability during the measurements as small as possible. The results do not follow simple assumptions on relationships to single parameters such as grain size. It was found that the complexity of natural occurring media is not yet sufficiently represented when modelling IP. At the same time simple correlation to permeability was found to be strong and consistent. Hence, adaptations with the aim of better representing the geo-structure of natural porous media were applied to the simplified model space used in Schwarz' IP-effect-theory. The resulting semi- empiric relationship was found to more accurately predict the IP-effect and its relation to the parameters grain size and permeability. If combined with recent findings about the effect of pore fluid electrochemistry together with advanced complex resistivity tomography, these results will allow us to picture diverse aspects of the subsurface with relative certainty. Within the framework of single measurement campaigns, hydrologiste can than collect data with information about the geo-structure and geo-chemistry of the subsurface. However, additional research efforts will be necessary to further improve the understanding of the physical origins of IP-effect and minimize the potential for false interpretations.¦-¦Dans l'étude des processus et caractéristiques hydrologiques des subsurfaces, la résolution spatiale donnée par les modèles hydrologiques dépasse souvent la résolution des données du terrain récoltées avec des méthodes classiques d'hydrologie. Récemment il est possible de réduire de plus en plus cet divergence spatiale entre modèles numériques et données du terrain par l'utilisation de méthodes géophysiques, notamment celles géoélectriques. Parmi les méthodes électriques, la polarisation provoquée (PP) permet de représenter la capacité des roches poreuses et des sols à stocker une charge électrique. En l'absence des métaux dans le sous-sol, cet effet est largement influencé par des caractéristiques de surface des matériaux. En conséquence les mesures PP offrent une information des interfaces entre solides et fluides dans les matériaux poreux que nous pouvons lier à la perméabilité également dirigée par ces mêmes paramètres. L'effet de la polarisation provoquée à été étudié dans différentes études de laboratoire, ainsi que sur le terrain. A cause d'une faible capacité de polarisation des matériaux sableux, comparé aux argiles, leur caractérisation par l'effet-PP reste difficile a interpréter d'une manière cohérente pour les environnements hétérogènes.¦Pour améliorer les connaissances sur l'importance de la structure du sous-sol sableux envers l'effet PP et des paramètres hydrologiques, nous avons fait des mesures de laboratoire variées. En détail, nous avons considéré des échantillons sableux de quartz avec des distributions de taille de grain entre sables fins et graviers fins, en diamètre cela fait entre 0,09 et 5,6 mm. Les caractéristiques de l'espace poreux sont changées en modifiant (i) la distribution de taille des grains, (ii) le degré de compaction, et (iii) le niveau d'hétérogénéité dans la distribution de taille de grains. En suite nous étudions comment ces changements influencent l'effet-PP, la perméabilité et la surface spécifique des échantillons. Les paramètres électrochimiques sont gardés à un minimum pendant les mesures. Les résultats ne montrent pas de relation simple entre les paramètres pétro-physiques comme par exemples la taille des grains. La complexité des media naturels n'est pas encore suffisamment représenté par les modèles des processus PP. Néanmoins, la simple corrélation entre effet PP et perméabilité est fort et consistant. En conséquence la théorie de Schwarz sur l'effet-PP a été adapté de manière semi-empirique pour mieux pouvoir estimer la relation entre les résultats de l'effet-PP et les paramètres taille de graines et perméabilité. Nos résultats concernant l'influence de la texture des matériaux et celles de l'effet de l'électrochimie des fluides dans les pores, permettront de visualiser des divers aspects du sous-sol. Avec des telles mesures géo-électriques, les hydrologues peuvent collectionner des données contenant des informations sur la structure et la chimie des fluides des sous-sols. Néanmoins, plus de recherches sur les origines physiques de l'effet-PP sont nécessaires afin de minimiser le risque potentiel d'une mauvaise interprétation des données.

Soil diffuse degassing and thermal energy fluxes from the southern Lakki plain, Nisyros (Greece)

Two diffuse soil CO2 flux surveys from the southern Lakki plain show that CO2 is mainly released from the hydrothermal explosion craters. The correspondence between high CO2 fluxes and elevated soil temperatures suggests that a flux of hot hydrothermal fluids ascends towards the surface. Steam mostly condenses near the surface and the heat given off is conductively transferred to the atmosphere through the soil, accompanied by a large CO2 flux. Tt was calculated, that 68 t d(-1) of hydrothermal CO2 are released through the total surveyed area of similar to1.3 km(2) Admitting that a steam flux of 2200 t d(-1) accompanies this CO2 flux, the thermal energy released through steam condensation amounts to 58 MW.

European climate variations over the past half-millennium reconstructed from groundwater

Temperature reconstructions for recent centuries are the basis of estimations of the natural variability in the climate system before and during the onset of anthropogenic perturbation. Here we present, for the first time, an independent and physically based reconstruction of mean annual temperature over the past half millennium obtained from groundwater in France. The reconstructed noble gas temperature (NGT) record suggests cooler than present climate conditions throughout the 16th-19th centuries. Periods of warming occur in the 17th-18th and 20th century, while cooling is reconstructed in the 19th century. A noticeable coincidence with other temperature records is demonstrated. Deuterium excess varies in parallel with the NGT, and indicates variation in the seasonality of the aquifer recharge; whereas high excess air in groundwater indicates periods with high oscillations of the water table.

Evolution and sedimentology of a channel fill in the sandy braided South Saskatchewan River and its comparison to the deposits of an adjacent compound bar

The depositional stratigraphy of within-channel deposits in sandy braided rivers is dominated by a variety of barforms (both singular `unit' bars and complex `compound' bars), as well as the infill of individual channels (herein termed `channel fills'). The deposits of bars and channel fills define the key components of facies models for braided rivers and their within-channel heterogeneity, knowledge of which is important for reservoir characterization. However, few studies have sought to address the question of whether the deposits of bars and channel fills can be readily differentiated from each other. This paper presents the first quantitative study to achieve this aim, using aerial images of an evolving modern sandy braided river and geophysical imaging of its subsurface deposits. Aerial photographs taken between 2000 and 2004 document the abandonment and fill of a 1 3 km long, 80 m wide anabranch channel in the sandy braided South Saskatchewan River, Canada. Upstream river regulation traps the majority of very fine sediment and there is little clay (<1%) in the bed sediments. Channel abandonment was initiated by a series of unit bars that stalled and progressively blocked the anabranch entrance, together with dune deposition and stacking at the anabranch entrance and exit. Complete channel abandonment and subsequent fill of up to 3 m of sediment took approximately two years. Thirteen kilometres of ground-penetrating radar surveys, coupled with 18 cores, were obtained over the channel fill and an adjacent 750 m long, 400 m wide, compound bar, enabling a quantitative analysis of the channel and bar deposits. Results show that, in terms of grain-size trends, facies proportions and scale of deposits, there are only subtle differences between the channel fill and bar deposits which, therefore, renders them indistinguishable. Thus, it may be inappropriate to assign different geometric and sedimentological attributes to channel fill and bar facies in object-based models of sandy braided river alluvial architecture.