Stratigraphic framework of gas accumulations in the Ría de Pontevedra (NW Spain)

Knowledge of marine geological environments in which shallow gas is accumulating is becoming increasingly important in global studies of climate change because a measurable proportion of the total methane source comes from continental margins. Previous studies have revealed that coastal environments represent important geological environments where microbial methane is being generated, is accumulating, and is being released. In the Ría de Pontevedra, at least 4.5 km² of seafloor in the innermost part of the ría is underlain by sediments containing natural gas. Seismic interpretation contributes new findings for the definition of periods and geological environments in which the gas could have been generated, and is accumulating and released in the Ría de Pontevedra. Groundtruthing the seismic data (facies, environments) makes it possible to identify favourable geological environments for gas generation in the sedimentary infill of the Ría de Pontevedra. Sequence stratigraphy based on high-resolution seismic profiles and post-Last Glacial Maximum sea-level records makes it possible to establish the stratigraphic architecture of the ría and to define the periods in which gas could have been generated. The results of this study show that the sedimentary infill is composed of a fifth-order sequence developed since the Last Glacial Maximum. Within this sequence, gas appears to have accumulated in the Holocene deposits associated with the latest transgressive and highstand system tracts. Seismic analysis shows that gas could have been generated in different geological environments in the Ría de Pontevedra. If coastal environments at times of lower sea level were similar to those of the present, organic-rich mud deposits (deposited mainly in lakes, estuaries and floodplains) could have survived transgression and remained buried as potential gas sources in the inner part of the ría.     

Spin rate estimation of sounding rockets using GPS wind-up

Carrier phase wind-up is a well-known effect that arises from the relative rotation between a transmitting and receiving antenna. In GPS measurements at L1 frequency, this effect translates into an error of 19.029 cm per full relative rotation of antennas. Since this effect is independent of the satellite elevation for pure rotation about the antenna boresight axis, it is usually absorbed by the clock estimation in navigation algorithms. Therefore, the impact of wind-up is usually neglected for applications that do not require accuracies to the cm level like RTK. However, in receiving platforms with high rotation rate, the accumulated wind-up value can be important and actually be larger than receiver noise or even ionospheric variations. Therefore, in such scenarios, the wind-up contribution can be isolated and used as a source of information to compute the spin rate of such platforms using an appropriate combination of GPS observables. This work shows some results of a coarse, yet simple, approach to monitor the rotation angle and spin-rate of spin stabilized sounding rockets flown by DLR.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Water balance in the playa-lakes of an arid environment,Monegros, NE Spain

The playa-lakes of the Monegros desert in north-east Spain are saline wetlands in an arid environment, a rare phenomenon in Europe. These extremely valuable habitats are threatened by changes associated with agricultural expansion and incorporation of new irrigated areas. An understanding of the present hydrologic regime will enable changes to be identified, particularly those brought about by flooding and pollution caused by irrigation surplus. This study sets out to show the results of applying a daily water balance in three selected playa-lakes. The balance was in two parts and consisted of: (1) the average balance for all the endorheic basin using the BALAN_11 program, and (2) the water balance in some playa-lakes, applying discharge flows obtained from the previous balance. The resulting volumes of water were converted to water depths and contrasted with reference volumes taken from field and Landsat images. The model was calibrated by applying various hypotheses of function which enabled the results to be adjusted. The proposed balance is an acceptable reproduction of field water measurements during this period, and underlines the consistency of the conceptual model. The methodology used is appropriate for understanding the playa-lakes function and for monitoring them for conservation purposes.     

Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

Sediment Organic Carbon in Todos Santos Bay, Baja California, Mexico

Sediment grain size and organic carbon (OC) data collected over the past 50 years, together with δ¹³C values of OC in recently collected samples, were analyzed to improve understanding of sediment OC distribution and abundance in Todos Santos Bay. Sediments in the submarine canyon at the mouth of the bay and in quiet-water locations along the shore are fine grained, high in OC, and have generally low δ¹³C values; sediments in high-energy environments are low in OC and have high δ¹³C values. A bivariate isotopic mixing model indicates that none of the sediments contain >50% terrigenous OC (average ~30%), and that the terrigenous OC content of the sediments is a small proportion of the OC content of local soils. Sediment OC composition is apparently controlled by energy-related sorting and deposition, oxidation of much of the original terrigenous OC, and replacement of some terrigenous OC by marine OC.     

Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.