Implementation aspects of sequential Gaussian simulation on irregular points

The increasing use of unstructured grids for reservoir modeling motivates the development of geostatistical techniques to populate them with properties such as facies proportions, porosity and permeability. Unstructured grids are often populated by upscaling high-resolution regular grid models, but the size of the regular grid becomes unreasonably large to ensure that there is sufficient resolution for small unstructured grid elements. The properties could be modeled directly on the unstructured grid, which leads to an irregular configuration of points in the three-dimensional reservoir volume. Current implementations of Gaussian simulation for geostatistics are for regular grids. This paper addresses important implementation details involved in adapting sequential Gaussian simulation to populate irregular point configurations including general storage and computation issues, generating random paths for improved long range variogram reproduction, and search strategies including the superblock search and the k-dimensional tree. An efficient algorithm for computing the variogram of very large irregular point sets is developed for model checking.     

Influence of artificial channels on the source and extent of saline water intrusion in the wind tide dominated wetlands of the southern Albemarle estuarine system (USA)

Saline water intrusion is degrading water quality in the channelized coastal wetlands of the southern Albemarle estuarine system (AES). The source, transport and fate of the saline water in the southern AES was determined by monitoring specific conductivity and water levels in small artificial channels, the Alligator River, the Alligator–Pungo Canal and the groundwater system for ~12 months. Results indicate that water levels are affected by wind tides which trigger the movement of saline water into the interior of the wetlands via the small canals. The wind tides are mostly driven by episodic southerly winds pushing saline water into the canal network and the groundwater regime proximal to the Alligator River. The saline waters persist in the canals as long as the wind tide events last. Specific conductivities from canals and groundwater are shown to be unexpectedly higher closer to the source of the Alligator River than toward the Albemarle Sound, suggesting that the large Alligator–Pungo Canal facilitates the northward migration of saline water from the Pamlico Sound to the Alligator River. Overwash and reversals in the flow directions between groundwater and surface water bodies suggest that saline water that is present during wind tide events may migrate into the groundwater system from surface water bodies. The results of this study reveal that whereas the large Alligator–Pungo Canal channels saline water to the AES, small artificial channels may also play significant roles in degrading water quality in the interior of channelized coastal wetlands.     

A re-examination of petrogenesis and 40Ar/39Ar systematics in the Chain of Ponds K-feldspar: “diffusion domain” archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, upon which thermochronological modeling of Ar diffusion in discrete “domains” was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. ⁴⁰Ar/³⁹Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest “diffusion domains.” Heterochemical phases also play an important role in producing nonlinear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick’s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the “diffusion domain” model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution among different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for ⁴⁰Ar/³⁹Ar hygrochronology as it is for U–Pb geochronology.     

Nucleosynthetic zirconium isotope anomalies in acid leachates of carbonaceous chondrites

Stepwise dissolutions of the carbonaceous chondrites Orgueil (CI), Murchison (CM) and Allende (CV) reveal large nucleosynthetic anomalies for Zr isotopes that contrast with the uniform compositions found in bulk meteorites. Two complementary nucleosynthetic components are observed: one enriched and one depleted in s-process nuclides. The latter component, characterized by excess ⁹⁶Zr, is most distinctive in the acetic acid leachate (up to ε⁹⁶Zr ≈ 50). The excess decreases with increasing acid strength and the final leaching steps of the experiment are depleted in ⁹⁶Zr and thus enriched in s-process nuclides. Presolar silicon carbide grains are likely host phases for part of the anomalous Zr released during these later stages. However, by mass balance they cannot account for the ⁹⁶Zr excesses observed in the early leaching steps and this therefore hints at the presence of at least one additional carrier phase with significant amounts of anomalous Zr. The data provide evidence that average solar system material consists of a homogenized mixture of different nucleosynthetic components, which can be partly resolved by leaching experiments of carbonaceous chondrites.     

High-Resolution Records of the Holocene Paleoenvironmental Variation Reflected by Carbonate and Its Isotopic Compositions in Bosten Lake and Response to Glacial Activities

Abstract: The Early Holocene paleoclimate in Bosten Lake on the northern margin of the Tarim Basin, southern Xinjiang, is reconstructed through an analysis of a 953 cm long core (BSTC2000) taken from Bosten Lake. Multiple proxies of this core, including the mineral components of carbonate, carbonate content, stable isotopic compositions of carbonate, Ca/Sr, TOC and C/N and C/S of organic matter, are used to reconstruct the climatic change since 8500 a B.P. The chronology model is made by nine AMS 14C ages of leaves, seeds and organic matter contained in two parallel cores. The climate was cold and wet during 8500 to 8100 a B.P. Temperature increased from 8100 to 6400 a B.P., the climate was warm and humid, and the lake expanded. The lake level was highest during this stage. Then from 6400 to 5100 a B.P., the climate became cold and the lake level decreased slightly. During the late mid-Holocene, the climate was hot and dry from 5100 to 3100 a B.P., but there was a short cold period during 4400 to 3800 a B.P. At this temporal interval, a mass of ice and snow melting water supplied the lake at the early time and made the lake level rise. The second highest lake level stage occurred during 5200 to 3800 a B.P. The climate was cool and wet during 3100 to 2200 a B.P., when the lake expanded with decreasing evaporation. The lake had the last short-term high level during 3100 to 2800 a B.P. After this short high lake level period, the lake shrank because of the long-term lower temperature and reduced water supply. From 2200 to 1200 a B.P., the climate was hot and dry, and the lake shrank greatly. Although the temperature decreased somewhat from 1200 a B.P. to the present, the climate was warm and dry. The lake level began to rise a little again, but it did not reach the river bed altitude of the Konqi River, an outflow river of the Bosten Lake.     

Late Quaternary glaciation and equilibrium-line altitudes of the Mayan Ice Cap, Guatemala, Central America

The Sierra los Cuchumatanes (3837 m), Guatemala, supported a plateau ice cap and valley glaciers around Montaña San Juan (3784 m) that totaled ∼ 43 km² in area during the last local glacial maximum. Former ice limits are defined by sharp-crested lateral and terminal moraines that extend to elevations of ∼ 3450 m along the ice cap margin, and to ca. 3000-3300 m for the valley glaciers. Equilibrium-line altitudes (ELAs) estimated using the area-altitude balance ratio method for the maximum late Quaternary glaciation reached as low as 3470 m for the valley glaciers and 3670 m for the Mayan Ice Cap. Relative to the modern altitude of the 0°C isotherm of ∼ 4840 m, we determined ELA depressions of 1110-1436 m. If interpreted in terms of a depression of the freezing level during maximal glaciation along the modern lapse rate of − 5.3°C km^− 1, this ΔELA indicates tropical highland cooling of ∼ 5.9 to 7.6 ± 1.2°C. Our data support greater glacial highland cooling than at sea level, implying a high tropical sensitivity to global climate changes. The large magnitude of ELA depression in Guatemala may have been partially forced by enhanced wetness associated with southward excursions of the boreal winter polar air mass.     

Equilibrium H/H fractionations in organic molecules: I. Experimental calibration of ab initio calculations

Carbon-bound hydrogen in sedimentary organic matter can undergo exchange over geologic timescales, altering its isotopic composition. Studies investigating the natural abundance distribution of ¹H and ²H in such molecules must account for this exchange, which in turn requires quantitative knowledge regarding the endpoint of exchange, i.e., the equilibrium isotopic fractionation factor (α_eq). To date, relevant data have been lacking for molecules larger than methane. Here we describe an experimental method to measure α_eq for C-bound H positions adjacent to carbonyl group (H_α) in ketones. H at these positions equilibrates on a timescale of days as a result of keto-enol tautomerism, allowing equilibrium ²H/¹H distributions to be indirectly measured. Molecular vibrations for the same ketone molecules are then computed using Density Functional Theory at the B3LYP/6-311G** level and used to calculate α_eq values for H_α. Comparison of experimental and computational results for six different straight and branched ketones yields a temperature-dependent linear calibration curve with slope = 1.081−0.00376T and intercept = 8.404−0.387T, where T is temperature in degrees Celsius. Since the dominant systematic error in the calculation (omission of anharmonicity) is of the same size for ketones and C-bound H in most other linear compounds, we propose that this calibration can be applied to analogous calculations for a wide variety of organic molecules with linear carbon skeletons for temperatures below 100 °C. In a companion paper (Wang et al., 2009) we use this new calibration dataset to calculate the temperature-dependent equilibrium isotopic fractionation factors for a range of linear hydrocarbons, alcohols, ethers, ketones, esters and acids.     

Nature of magnetic grains in basalts and implications for palaeomagnetism

Journal of Earth System Science - Investigations involving temperature dependence of low-field susceptibility and of low- and high-field hysteresis have been carried out on more than one thousand...     

Robustness of the holistic seismic risk evaluation in urban centers using the USRi

The Urban Seismic Risk index (USRi) published in a previous article (Carreño et al., Nat Hazards 40:137–172, 2007) is a composite indicator that measures risk from an integrated perspective and guides decision-making for identifying the main interdisciplinary factors of vulnerability to be reduced or intervened. The first step of the method is the evaluation of the potential physical damage (hard approach) as a result of the convolution of the seismic hazard with the physical vulnerability of buildings and infrastructure. Subsequently, a set of social context conditions that aggravate the physical effects is also considered (soft approach). According to this procedure, the physical risk index is evaluated for each unit of analysis from existing loss scenarios, whereas the total risk index is obtained by multiplying the former index by an impact factor using an aggravating coefficient, based on variables associated with the socio-economic conditions of each unit of analysis. The USRi has been developed using the underlying holistic and multi-hazard approach of the Urban Risk Index framework proposed for the evaluation of disaster risk in different megacities worldwide. This article presents the sensitivity analysis of the index to different parameters such as input data, weights and transformation functions used for the scaling or normalization of variables. This analysis has been performed using the Monte Carlo simulation to validate the robustness of this composite indicator, understanding as robustness how the cities maintain the ranking as well as predefined risk level ranges, when compared with the deterministic results of risk. Results are shown for different cities of the world.