Effect of Multipoint Heterogeneity on Nonlinear Transformations for Geological Modeling： Porosity-Permeability Relations Revisited

An analysis of statistical expected values for transformations is performed in this study to quantify the effect of heterogeneity on spatial geological modeling and evaluations. Algebraic transformations are frequently applied to data from logging to allow for the modeling of geological properties. Transformations may be powers, products, and exponential operations which are commonly used in well-known relations (e.g., porosity-permeability transforms). The results of this study show that correct computations must account for residual transformation terms which arise due to lack of independence among heterogeneous geological properties. In the case of an exponential porosity-permeability transform, the values may be positive. This proves that a simple exponential model back-transformed from linear regression underestimates permeability. In the case of transformations involving two or more properties, residual terms may represent the contribution of heterogeneous components which occur when properties vary together, regardless of a pair-wise linear independence. A consequence of power- and product-transform models is that regression equationswithin those transformations need corrections via residual cumulants. A generalization of this result isthat transformations of multivariate spatial attributes require multiple-point random variable relations. This analysis provides practical solutions leading to a methodology for nonlinear modeling using correct back transformations in geology.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

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.     

High-resolution seismic analysis of the coastal Mecklenburg Bay (North German Basin): the pre-Alpine evolution

The pre-Alpine structural and geological evolution in the northern part of the North German Basin have been revealed on the basis of a very dense reflection seismic profile grid. The study area is situated in the coastal Mecklenburg Bay (Germany), part of the southwestern Baltic Sea. From the central part of the North German Basin to the northern basin margin in the Grimmen High area a series of high-resolution maps show the evolution from the base Zechstein to the Lower Jurassic. We present a map of basement faults affecting the pre-Zechstein. The pre-Alpine structural evolution of the region has been determined from digital mapping of post-Permian key horizons traced on the processed seismic time sections. The geological evolution of the North German Basin can be separated into four distinct periods in the Rerik study area. During Late Permian and Early Triassic evaporites and clastics were deposited. Salt movement was initiated after the deposition of the Middle Triassic Muschelkalk. Salt pillows, which were previously unmapped in the study area, are responsible for the creation of smaller subsidence centers and angular unconformities in the Late Triassic Keuper, especially in the vicinity of the fault-bounded Grimmen High. In this area, partly Lower Jurassic sediments overlie the Keuper unconformably. The change from extension to compression in the regional stress field remobilized the salt, leading to a major unconformity marked at the base of the Late Cretaceous.     

Giant ancient landslide in the Alma water gap (Crimean Mountains,Ukraine): notes to the predisposition,structure, and chronology

Large-scale ancient landslides of the area of more than 5 km² and volume exceeding 200 × 10⁶ m³ are characteristic features of the valleys incised in the northern periphery of the Crimean Mountains (Ukraine). The largely affected area is located in the outermost cuesta range of the Crimean Mountains which consists of rigid Sarmatian limestones overlying weak Middle Miocene and Upper Palaeogene deposits. A giant landslide arose in the Alma water gap as a reflection of several coincident preparatory factors such as suitable bedrock stratification, smectite-rich bedrock exposed to swelling activity, presence of faults parallel to the valley trend, and river capture event which preceded the landslide event. The occurrence of such ancient megaslides is particularly interesting in the area which is characterized by low precipitation (<500 mm/year) and weak contemporary seismicity. It probably reflects a more dynamic environment in humid phases of the Holocene; however, seismic triggering along the Mesozoic suture zone cannot be rejected. Compressional features such as gravitational folds in the central and distal parts of the landslide, which probably correlate with the whole landslide genesis or its significant reactivation, arose, according to the radiocarbon dating, during the Holocene climatic optimum in the Atlantic period. The slope deformation has been relatively quiescent since that time, except minor historic reactivization which took place in the frontal part of the landslide. We suppose that the studied landslide could be classified as a transitional type of slope deformation with some signs of spreading and translational block slides.