Geodetic intersection on the ellipsoid

Through each of two known points on the ellipsoid a geodesic is passing in a known azimuth. We solve the problem of intersection of the two geodesics. The solution for the latitude is obtained as a closed formula for the sphere plus a small correction, of the order of the eccentricity of the ellipsoid, which is determined by numerical integration. The solution is iterative. Once the latitude is obtained, the longitude is determined without iteration.     

Fluid inclusion and stable isotope characteristics of the Karalar (Gazipaşa,Antalya) barite-galena deposits,Southern Turkey

The Karalar galena-barite deposits are typical examples of the carbonate-hosted barite-galena deposits widely occurring in the Central Tauride Belt. These deposits are located in Permian limestones as ore veins along fault zones that are discordant to the bedding planes and as thin veinlets within hardly brecciated bottom zones of Permian limestones. The ore deposits contain mainly barite and galena and small amounts of sphalerite, pyrite, fahlore, limonite, quartz, and calcite. Barites occur during the earlier episode of mineralization and were mylonitized before the formation of galena. Galena and other minerals occur epigenetically with barite along porous zones between brecciated barite crystals and are especially abundant in the hardly mylonitized zones. Fluid inclusion studies indicate that the ore deposits of the area were developed by hydrothermal fluids with following characteristics: they contained NaCl, CaCl₂, and MgCl₂; the salinity of the fluids was relatively high and their temperature was low during the crystallization of barites in the earlier episode of mineralization; and the salinity of the fluids decreased and their temperature increased during the crystallization of sulfide minerals through the later episode of mineralization. δ₁₈O and the δD results indicate that the water in the mineralizing fluid was of meteoric origin. The deposits were formed by deep meteoric water circulating through marine sediments. Shallow circulating and slightly warmed fluids dissolved Ba and sulfate from marine sediments of the basement and transported them to the mineralization environment during the early episode, while deep circulating and more heated fluids dissolved Pb, Zn, and other elements from the basement and reached the environment during the later episode of mineralization. The text was submitted by the authors in English.     

Comparison of dissolved and particulate arsenic distributions in shallow aquifers of Chakdaha, India, and Araihazar, Bangladesh

Background  

The origin of the spatial variability of dissolved As concentrations in shallow aquifers of the Bengal Basin remains poorly understood. To address this, we compare here transects of simultaneously-collected groundwater and aquifer solids perpendicular to the banks of the Hooghly River in Chakdaha, India, and the Old Brahmaputra River in Araihazar, Bangladesh.     

Seasonal variation of carbon fluxes in a sparse savanna in semi arid Sudan

Background

Large spatial, seasonal and annual variability of major drivers of the carbon cycle (precipitation, temperature, fire regime and nutrient availability) are common in the Sahel region. This causes large variability in net ecosystem exchange and in vegetation productivity, the subsistence basis for a major part of the rural population in Sahel. This study compares the 2005 dry and wet season fluxes of CO₂ for a grass land/sparse savanna site in semi arid Sudan and relates these fluxes to water availability and incoming photosynthetic photon flux density (PPFD). Data from this site could complement the current sparse observation network in Africa, a continent where climatic change could significantly impact the future and which constitute a weak link in our understanding of the global carbon cycle.

Results

The dry season (represented by Julian day 35–46, February 2005) was characterized by low soil moisture availability, low evapotranspiration and a high vapor pressure deficit. The mean daily NEE (net ecosystem exchange, Eq. 1) was -14.7 mmol d^-1 for the 12 day period (negative numbers denote sinks, i.e. flux from the atmosphere to the biosphere). The water use efficiency (WUE) was 1.6 mmol CO₂ mol H₂O^-1 and the light use efficiency (LUE) was 0.95 mmol CO₂ mol PPFD^-1. Photosynthesis is a weak, but linear function of PPFD. The wet season (represented by Julian day 266–273, September 2005) was, compared to the dry season, characterized by slightly higher soil moisture availability, higher evapotranspiration and a slightly lower vapor pressure deficit. The mean daily NEE was -152 mmol d^-1 for the 8 day period. The WUE was lower, 0.97 mmol CO₂ mol H₂O^-1 and the LUE was higher, 7.2 μmol CO₂ mmol PPFD^-1 during the wet season compared to the dry season. During the wet season photosynthesis increases with PPFD to about 1600 μmol m^-2s^-1 and then levels off.

Conclusion

Based on data collected during two short periods, the studied ecosystem was a sink of carbon both during the dry and wet season 2005. The small sink during the dry season is surprising and similar dry season sinks have not to our knowledge been reported from other similar savanna ecosystems and could have potential management implications for agroforestry. A strong response of NEE versus small changes in plant available soil water content was found. Collection and analysis of flux data for several consecutive years including variations in precipitation, available soil moisture and labile soil carbon are needed for understanding the year to year variation of the carbon budget of this grass land/sparse savanna site in semi arid Sudan.     

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.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Contributions of boron isotopes to understanding the hydrogeochemistry of the coastal detritic aquifer of Castellón Plain,Spain

The Castellón Plain alluvial aquifer, Spain, is intensively exploited to meet the demand for agricultural irrigation and industrial water supply. The geochemistry of its groundwater shows complex salinization in the northern and southern parts of the aquifer, with significant pollution from human origin in the central portion. Boron content and B isotope geochemistry are useful for distinguishing between various sources of pollution and their relative importance in different parts of this aquifer. Boron concentrations in the groundwater vary between 0.01 and 0.85 mg/L. In the more saline groundwaters, found at the northern and southern ends of the study area, the presence of B is linked to inputs from seawater and water with a calcium-magnesium sulphate facies, which feed the aquifer and clearly influence the chemistry of its waters. Evidence of B adsorption processes in some samples is shown by the low B/Cl ratios and the high values of δ¹¹B. In the central portion of the aquifer, the high B/Cl ratios and the strongly negative δ¹¹B are related to pollution of human origin.     

New age constraints on emplacement of the Cévenol granitoids, South French Massif Central

During the development of the Variscan orogeny, large amounts of granitic melt were produced, giving rise to the intrusion of granitoids at different structural levels. Despite numerous studies, ages available from previous work on the Cévennes granites remain largely imprecise. In order to better constrain the age and emplacement mode of these granites, we have combined U–Pb dating on monazites and zircons and ⁴⁰Ar/³⁹Ar dating on biotites with petrological observations, major element chemical analysis and SEM zircon imaging on five samples from the Aigoual–St Guiral–Liron and Mont Lozère granitic massifs. The results revealed that granitic intrusions and cooling in Southern Cévennes occurred in a short time span at ∼306 Ma after the main episode of regional metamorphism. Petrological and chemical data suggest that they result from a mixing between mantle-derived basic magmas (lamprophyres) and lower crust acid magmas. At a regional scale the production of these melts occurred at the end of crustal thickening induced by nappe stacking, at the same time as the late anatectic events recorded further north in the Velay dome and the granulite facies metamorphism recorded in metasedimentary granulite enclaves brought up by Tertiary volcanoes of the Velay area (Bournac).