Quantification of the relation between surface morphodynamics and subsurface sedimentological product in sandy braided rivers

This study uses digital elevation models and ground-penetrating radar to quantify the relation between the surface morphodynamics and subsurface sedimentology in the sandy braided South Saskatchewan River, Canada. A unique aspect of the methodology is that both digital elevation model and ground-penetrating radar data were collected from the same locations in 2004, 2005, 2006 and 2007, thus enabling the surface morphodynamics to be tied explicitly to the associated evolving depositional product. The occurrence of a large flood in 2005 also allowed the influence of discharge to be assessed with respect to the processproduct relationship. The data demonstrate that the morphology of the study reach evolved even during modest discharges, but more extensive erosion was caused by the large flood. In addition, the study reach was dominated by compound bars before the flood, but switched to being dominated by unit bars during and after the flood. The extent to which the subsurface deposits (the product') were modified by the surface morphodynamics (the process') was quantified using the changes in radar-facies recorded in sequential ground-penetrating radar surveys. These surveys reveal that during the large flood there was an increase in the proportion of facies associated with bar margin accretion and larger dunes. In subsequent years, these facies became truncated and replaced with facies associated with smaller dune sets. This analysis shows that unit bars generally become truncated more laterally than vertically and, thus, they lose the high-angle bar margin deposits and smaller scale bar-top deposits. In general, the only fragments that remain of the unit bars are dune sets, thus making identification of the original unit barform problematic. This novel data set has implications for what may ultimately become preserved in the rock record.

Interpretação de dados de GPR com base na hierarquização de superfícies limitantes e na adaptação de critérios sismoestratigráficos

Due to its high resolution, Ground Penetrating Radar (GPR) has been used to image subsurface sedimentary deposits. Because GPR and Seismic methods share some principles of image construction, the classic seismostratigraphic interpretation method has been also applied as an attempt to interpret GPR data. Nonetheless some advances in few particular contexts, the adaptations from seismic to GPR of seismostratigraphic tools and concepts unsuitable because the meaning given to the termination criteria in seismic stratigraphy do not represent the adequate geologic record in the GPR scale. Essentially, the open question relies in proposing a interpretation method for GPR data which allow not only relating product and sedimentary process in the GPR scale but also identifying or proposing depositional environments and correlating these results with the well known Sequence Stratigraphy cornerstones. The goal of this dissertation is to propose an interpretation methodology of GPR data able to perform this task at least for siliciclastic deposits. In order to do so, the proposed GPR interpretation method is based both on seismostratigraphic concepts and on the bounding surface hierarchy tool from Miall (1988). As consequence of this joint use, the results of GPR interpretation can be associated to the sedimentary facies in a genetic context, so that it is possible to: (i) individualize radar facies and correlate them to the sedimentary facies by using depositional models; (ii) characterize a given depositional system, and (iii) determine its stratigraphic framework highligthing how it evolved through geologic time. To illustrate its use the proposed methodology was applied in a GPR data set from Galos area which is part of the Galinhos spit, located in Rio Grande do Norte state, Northeastern Brazil. This spit presents high lateral sedimentary facies variation, containing in its sedimentary record from 4th to 6th cicles caused by high frequency sea level oscillation. The interpretation process was done throughout the following phases: (i) identification of a vertical facies succession, (ii) characterization of radar facies and its associated sedimentary products, (iii) recognition of the associated sedimentary process in a genetic context, and finally (iv) proposal of an evolutionay model for the Galinhos spit. This model proposes that the Galinhos spit is a barrier island constituted, from base to top, of the following sedimentary facies: tidal channel facies, tidal flat facies, shore facies, and aeolic facies (dunes). The tidal channel facies, in the base, is constituted of lateral accretion bars and filling deposits of the channels. The base facies is laterally truncated by the tidal flat facies. In the foreshore zone, the tidal flat facies is covered by the shore facies which is the register of a sea transgression. Finally, on the top of the stratigraphic column, aeolic dunes are deposited due to areal exposition caused by a sea regression

Conditioning of multiple-point statistics facies simulations to tomographic images

Geophysical tomography captures the spatial distribution of the underlying geophysical property at a relatively high resolution, but the tomographic images tend to be blurred representations of reality and generally fail to reproduce sharp interfaces. Such models may cause significant bias when taken as a basis for predictive flow and transport modeling and are unsuitable for uncertainty assessment. We present a methodology in which tomograms are used to condition multiple-point statistics (MPS) simulations. A large set of geologically reasonable facies realizations and their corresponding synthetically calculated cross-hole radar tomograms are used as a training image. The training image is scanned with a direct sampling algorithm for patterns in the conditioning tomogram, while accounting for the spatially varying resolution of the tomograms. In a post-processing step, only those conditional simulations that predicted the radar traveltimes within the expected data error levels are accepted. The methodology is demonstrated on a two-facies example featuring channels and an aquifer analog of alluvial sedimentary structures with five facies. For both cases, MPS simulations exhibit the sharp interfaces and the geological patterns found in the training image. Compared to unconditioned MPS simulations, the uncertainty in transport predictions is markedly decreased for simulations conditioned to tomograms. As an improvement to other approaches relying on classical smoothness-constrained geophysical tomography, the proposed method allows for: (1) reproduction of sharp interfaces, (2) incorporation of realistic geological constraints and (3) generation of multiple realizations that enables uncertainty assessment.

Architectural characterization of a delta-front reservoir analogue combining ground penetrating radar and electrical resistivity tomography: Roda Sandstone (Lower Eocene, Graus-Tremp basin, Spain)

Three-dimensional reconstruction of reservoir analogues can be improved combining data from different geophysical methods. Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT) data are valuable tools, since they provide subsurface information from internal architecture and facies distribution of sedimentary rock bodies, enabling the upgrading of depositional models and heterogeneity reconstruction. The Lower Eocene Roda Sandstone is a well-known deltaic complex widely studied as a reservoir analogue that displays a series of sandstone wedges with a general NE to SW progradational trend. To provide a better understanding of internal heterogeneity of a 10m-thick progradational delta-front sandstone unit, 3D GPR data were acquired. In addition, common midpoints (CMP) to measure the sandstone subsoil velocity, test profiles with different frequency antennas (25, 50 and 100MHz) and topographic data for subsequent correction in the geophysical data were also obtained. Three ERT profiles were also acquired to further constrain GPR analysis. These geophysical results illustrate the geometry of reservoir analogue heterogeneities both depositional and diagenetic in nature, improving and complementing previous outcrop-derived data. GPR interpretation using radar stratigraphy principles and attributes analysis provided: 1)tridimensional geometry of major stratigraphic surfaces that define four units in the GPR Prism, 2) image the internal architecture of the units and their statistical study of azimuth and dips, useful for a quick determination of paleocurrent directions. These results were used to define the depositional architecture of the progradational sandbody that shows an arrangement in very-high-frequency sequences characterized by clockwise paleocurrent variations and decrease of the sedimentary flow, similar to those observed at a greater scale in the same system. This high-frequency sequential arrangement has been attributed to the autocyclic dynamics of a supply-dominated delta- front where fluvial and tidal currents are in competition. The resistivity models enhanced the viewing of reservoir quality associated with cement distribution caused by depositional and early diagenetic processes related to the development of transgressive and regressive systems tracts in igh-frequency sequences.

Estratigrafia da sucessão sedimentar Pós-Barreiras (Zona Bragantina, Pará) com base em radar de penetração no solo

Investigação pioneira aplicando Radar de Penetração no Solo (GPR) na área da Praia do Atalaia, norte do Brasil, nos permitiu caracterizar, pela primeira vez, fácies e estratigrafia dos depósitos conhecidos informalmente como Sedimentos Pós-Barreiras (Plioceno e mais recente). Esta sucessão sedimentar recobre discordantemente o embasamento miocênico, representado pelas formações Pirabas/Barreiras. Três unidades estratigráficas foram reconhecidas. A Unidade 1, inferior, consiste em um intervalo com 6 m de espessura média dominado por reflexões pobremente definidas e de baixa amplitude, e que intergradam com reflexões de média escala dos tipos tangencial, obliquo e hummocky. A Unidade 2, intermediária, possui cerca de 9 m de espessura e inclui, principalmente, reflexões oblíquas de larga escala, cujas configurações variam de paralela, tangencial, sigmoidal a sigmoidal-complexa. A Unidade 3, superior, corresponde a um intervalo entre 3,5 e 9 m de espessura, sendo dominada por reflexões hummocky, seguidas por reflexões de média escala dos tipos oblíquo, paralelo a sub-paralelo, e em corte-e-preenchimento. A análise da configuração das reflexões internas e geometria das reflexões nos leva a propor que a unidade correspondente aos Sedimentos Pós-Barreiras é mais variável faciologicamente que inicialmente imaginado, incluindo depósitos eólicos (dunas costeiras), bem como depósitos de cordão litorâneo, planície de maré, canal e mangue. Além disto, o mapeamento das três unidades descritas acima é importante para desvendar a complexidade de sedimentação versus erosão durante o Neógeno tardio no norte do Brasil.

Classificação de sistemas meteorológicos e comparação da precipitação estimada pelo radar e medida pela rede telemétrica na bacia hidrográfica do alto Tietê

Foram analisadas características da precipitação estimada a partir de 145.194 campos de refletividade, de um total de 827 dias entre 1998 e 2003, obtidos do Radar Meteorológico de São Paulo (RSP). Os eventos foram classificados de acordo com intensidades de precipitação; em Convectivos (EC) e Estratiformes (EE). Quanto à morfologia, cinco tipos de sistemas foram identificados; Convecção Isolada (CI), Brisa Marítima (BM), Linhas de Instabilidade (LI), Bandas Dispersas (BD) e Frentes Frias (FF). Eventos convectivos dominam na primavera e verão e estratiformes no outono e inverno. A CI e a BM tiveram maiores picos de atuação entre outubro e março enquanto as FF de abril a setembro. BD atuam durante todo o ano e as LI só não foram observadas nos meses de junho e julho. Uma comparação pontual entre a precipitação medida pela telemetria e estimada com o radar foi realizada e, mostrou haver, na maioria dos casos, um viés positivo do RSP, para acumulações de 10, 30 e 60 minutos. Com o objetivo de integrar as estimativas de precipitação do radar com as medidas da rede telemétrica, por meio de uma análise objetiva estatística, foram obtidas dos campos de precipitação do radar as estruturas das correlações espaciais em função da distância para acumulações de chuva de 15, 30, 60 e 120 minutos para os cinco tipos de sistemas precipitantes que foram caracterizados. As curvas das correlações espaciais médias de todos os eventos de precipitação de cada sistema foram ajustadas por funções polinomiais de sexta ordem. Os resultados indicam diferenças significativas na estrutura espacial das correlações entre os sistemas precipitantes.

Facies, delta(13)C, delta(15)N and C/N analyses in a late Quaternary compound estuarine fill, northern Brazil and relation to sea level

The present work integrates sedimentary facies, (14)C dating, delta(13)C, delta(15)N, and C/N with geologic and geomorphologic data available from literature. The aim was to characterize the depositional settings of a late Quaternary estuary in northeastern Marajo Island and analyze its evolution within the context of relative sea level fluctuations. The data derive from four continuous cores along a proximal-to-distal transect of a paleoestuary, previously recognized using remote sensing information. Fifteen sediment samples recorded ages ranging from 42,580 +/- 1430 to 3184 +/- 37 (14)C yr B.P. Fades analysis indicated fine- to coarse-grained sands with parallel lamination or cross stratification, massive or laminated muds and heterolithic deposits. delta(13)C (-28.1 parts per thousand to -19.7 parts per thousand, mean = -23.0 parts per thousand), delta(15)N (+ 14.8 parts per thousand to + 4.7 parts per thousand, mean = + 9.2 parts per thousand) and C/N (14.5 to 1.5, mean = 7.9) indicate mostly marine and freshwater phytoplankton sources for the organic matter. The results confirm a large late Quaternary paleoestuary in northeastern Marajo Island. The distribution of delta(13)C, delta(15)N, and C/N, together with fades associations, led to identify depositional settings related to fluvial channel, floodplain, tidal channel/tidal flat, central basin, tidal delta, and tidal inlet/sand barrier. These deposits are consistent with a wave-dominated estuary. Variations in stratigraphy and geochemistry are controlled by changes in relative sea level, revealing a main transgression from an undetermined time around 42,000 (14)C yr B.P. and 29,340 (+/- 200) (14)C yr B.P., which is synchronous to the overall drop in sea level after the last interglacial. Following this period, and probably until 9110 +/- 37 (14)C yr B.P., i.e., during a time interval encompassing two glacial episodes including the Last Glacial and the Younger Dryas, there was a pronounced drop in sea level, recorded by the development of a major erosional discontinuity due to valley re-incision. Sea level rose again until 5464 +/- 40 (14)C yr B.P, just before the main worldwide mid-Holocene transgressive peak. Mid to late Holocene coastal progradation ended the Marajo paleoestuarine history, and promoted the establishment of continental conditions throughout the island. The divergence comparing the Marajo sea level behavior with the eustatic curve allows hypothesizing that post-rifting tectonics along the Brazilian Equatorial margin influenced the sedimentary evolution of the studied paleoestuary. Considering that sedimentary facies in estuarine settings are highly variable both laterally and vertically, the present integration of facies with isotope and elemental analyses was crucial to provide a more precise interpretation of the Late Pleistocene and Holocene Marajo paleoestuary, and analyze its sea level history within the eustatic and tectonic context. (C) 2010 Elsevier B.V. All rights reserved.

Quaternary paleoenvironments and relative sea-level changes in Marajo Island (Northern Brazil): Facies, delta(13)C, delta(15)N and C/N

The reconstruction of physical environments of Amazonian areas is of great interest to determine the dynamic evolution of the Amazon drainage basin. However. few studies have emphasized the Quaternary deposits in this region. which is mostly due to the lack of natural exposures imposed by the low topography. This work integrates facies analysis. radiocarbon dating, delta(13)C, delta(15)N, and C/N of an 124 m-thick core from an area located at the mouth of the Amazon River. northeastern Amazonia. The study records deposits up to 50.795 (14)C yr B P. in age. which formed in a variety of depositional environments including fluvial channel, tidal flat, outer estuarine basin to shallow marine. inner estuarine basin, estuarine channel and lagoon. Facies interpretation was significantly improved with the inclusion of delta(13)C, delta(15)N, and C/N analyses of organic matter extracted from the sediments The obtained values conform to a transitional. mostly estuarine paleosetting evolved during successive relative sea-level fluctuations. The results suggest fluvial deposition between 40,950 (+/- 590) and 50.795 (14)C yr B P, with a rise in relative sea level that commenced between 35,567 (+/- 649) and 39,079 (+/- 1114) (14)C yr B P. An overall transgression took place until 29,340 (+/- 340) (14)C yr B P., after which the relative sea level dropped, favoring valley rejuvenation and incision. Following this time up to 10,479 (+/- 34) (14)C yr B.P. a rise in relative sea level filled up the valley with estuarine deposits After 10.479(+/- 34) (14)C yr B.P., the estuary was replaced by a lagoon At the end of the Holocene, the coastline prograided approximately 45 km northward, replaci ng the lagoon by a lake system Despite the influence of eustatic fluctuations. regional tectonics played a significant role to create new space where these Late Pleistocene and Holocene sediments accumulated. (C) 2009 Elsevier B V All rights reserved.

Facies architecture of the CRP-3 drillhole, Victoria Land Basin, Antarctica

The Cenozoic Victoria Land Basin (VLB) stratigraphic section penetrated by CRP-3 is mostly of Early Oligocene age. It contains an array of lithofacies comprising fine-grained mudrocks, interlaminated and interbedded mudrocks/sandstones, mud-rich and mud-poor sandstones, conglomerates and diametites that are together interpreted as the products of shallow marine to possibly non-marine environments of deposition, affected by the periodic advance and retreat of tidewater glaciers. This lithofacies assemblage can be readily rationalised using the facies scheme designed originally for CRP-2/2A, and published previously. The uppermost 330 metres below sea floor (mbsf) shows a cyclical arrangement of lithofacies also similar to that recognised throughout CRP-2/2A, and interpreted to reflect cyclical variations in relative sea-level driven by ice volume fluctuations ("Motif A"). Between 330 and 480 mbsf, a series of less clearly cyclical units, generally fining-upward but nonetheless incorporating a significant subset of the facies assemblage, has been identified and noted in the Initial Report as "Motif B. Below 480 mbsf, the section is arranged into a repetitive succession of fining-upward units, each of which comprises dolerite clast conglomerate at the base passing upward into relatively thick intervals of sandstones. The cycles present down 480 mbsf are defined as sequences, each interpreted to record cyclical variation of relative sea-level. The thickness distribution of sequences in CRP-3 provides some insights into the geological variables controlling sediment accumulation in the Early Oligocene section. The uppermost part of the section in CRP-3 comprises two or three thick, complete sequences that show a broadly symmetrical arrangement of lithofacies (similar to Sequences 9-11 in CRP-2/2A). This suggests a period of relatively rapid tectonic subsidence, which allowed preservation of the complete facies cycle. Below Sequence 3, however, is a considerable interval of thin, incomplete and erosionally truncated sequences (4-23), which incorporates both the remainder of Motif A sequences and all Motif B sequences recognised. The thinner and more truncated sequences suggest sediment accumulation under conditions of reduced accommodation, and given the lack of evidence for glacial conditions (see Powell et al., this volume) tends to argue for a period of reduced tectonic subsidence. The section below 480 mbsf consists of a series of fining-upward, conglomerate to sandstone intervals which cannot be readily interpreted in terms of relative sea-level change. A relatively mudrock-rich interval above the basal conglomerate/breccia (782-762 mbsf) may record initial flooding of the basin during early rift subsidence. The lithostratigraphy summarised above has been linked to seismic reflection data using depth conversion techniques (Henrys et al., this volume). The three uppermost reflectors ("o", "p" and "q") correlate to the package of thick sequences 1-3, and several deeper reflectors can also be correlated to sequence boundaries. The package of thick Sequences 1-3 shows a sheet-like cross-sectional geometry on seismic reflection lines, unlike the similar package recognised in CRP-2/2A.

Delineation of sulphide ore-zones by borehole radar tomography at Hellyer Mine, Australia

Velocity and absorption tomograms are the two most common forms of presentation of radar tomographic data. However, mining personnel, geophysicists included, are often unfamiliar with radar velocity and absorption. In this paper, general formulae are introduced, relating velocity and attenuation coefficient to conductivity and dielectric constant. The formulae are valid for lossy media as well as high-resistivity materials. The transformation of velocity and absorption to conductivity and dielectric constant is illustrated via application of the formulae to radar tomograms from the Hellyer zinc-lead-silver mine, Tasmania, Australia. The resulting conductivity and dielectric constant tomograms constructed at Hellyer demonstrated the potential of radar tomography to delineate sulphide ore zones. (C) 2001 Elsevier Science B.V. All rights reserved.

Air target detection via bistatic radar based on LEOS communication signals

The paper discusses the bistatic radar parameters for the case when the transmitter is a satellite emitting communication signals. The model utilises signals from an Iridium-like low earth orbiting satellite system. The maximum detection range, when thermal noise-limited, is discussed at the theoretical level and these results are compared with experimentation. Satellite-radar signal levels and the power of ground reflections are evaluated.

Three-dimensional space-borne synthetic aperture radar (SAR) imaging with multiple pass processing

Three-dimensional (3D) synthetic aperture radar (SAR) imaging via multiple-pass processing is an extension of interferometric SAR imaging. It exploits more than two flight passes to achieve a desired resolution in elevation. In this paper, a novel approach is developed to reconstruct a 3D space-borne SAR image with multiple-pass processing. It involves image registration, phase correction and elevational imaging. An image model matching is developed for multiple image registration, an eigenvector method is proposed for the phase correction and the elevational imaging is conducted using a Fourier transform or a super-resolution method for enhancement of elevational resolution. 3D SAR images are obtained by processing simulated data and real data from the first European Remote Sensing satellite (ERS-1) with the proposed approaches.

Resistively loaded helical antennas for ground-penetrating radar

Resistively loaded helical antennas, used in the normal mode and horizontally polarised, are modelled using the moment method above typical lossy ground. The distributed resistive loading was adjusted to maintain a two octave bandwidth. The centre frequency of 1 m dipoles was reduced from 250 MHz for the straight resistive wire to 50 MHz for a helix of pitch 2.5 cm and diameter 5 cm. The reduction in efficiency required to maintain the bandwidth for this helix was 12 dB. This agrees reasonably with the theory for small antennas in free space. The results were also verified by comparing measurements performed on a monopole resistively loaded helical antenna in a watertank with the numerical model used elsewhere.