Mapping Curvilinear Structures with Local Anisotropy Kriging

Interpolating geo-data with curvilinear structures using geostatistics is often disappointing. Channels, for example, become disconnected sets of lakes when interpolated from point data. In order to improve the interpolation of geological structures (e.g., curvilinear structures), we present a new form of kriging, local anisotropy kriging (LAK). Local anisotropy kriging combines a gradient algorithm from image analysis with kriging in an iterative way. After an initial standard kriging interpolation, the gradient algorithm determines the local anisotropy for each cell in the grid using a search area around the cell. Subsequently, kriging is carried out with the spatially varying anisotropy. The anisotropy calculation and subsequent kriging steps will then succeed until the result is satisfactory in the way of reproducing the curvilinear structures. Depending on the size of the search area more or less detail in the geological structures can be reproduced with LAK. Using test examples we show that LAK interpolates data with curvilinear structures more realistically than standard kriging. In a real world case, using bathymetric data of the Oosterschelde estuary, LAK also proves to be quantitatively superior to standard kriging. Absolute interpolation errors are decreased by 23%. Local anisotropy kriging only uses information from point data, which makes the method very objective, it only presents “what the data can tell.”     

Simulated impacts of artificial groundwater recharge and discharge on the source area and source volume of an Atlantic Coastal Plain stream, Delaware, USA

A numerical groundwater-flow model was used to characterize the source area and volume of Phillips Branch, a baseflow-dominated stream incising a highly permeable unconfined aquifer on the low relief Delmarva Peninsula, USA. Particle-tracking analyses indicate that the source area (5.51 km²) is ～20% smaller than the topographically defined watershed (6.85 km²), and recharge entering ～37% of the surface watershed does not discharge to Phillips Branch. Groundwater residence time within the source volume ranges from a few days to almost 100 years, with 95% of the volume “flushing” within 50 years. Artificial discharge from groundwater pumping alters the shape of the source area and reduces baseflow due to the interception of stream flow paths, but has limited impacts on the residence time of groundwater discharged as baseflow. In contrast, artificial recharge from land-based wastewater disposal substantially reduces the source area, lowers the range in residence time due to the elimination of older flow paths to the stream, and leads to increased discharge to adjacent surface-water bodies. This research suggests that, in this and similar hydrogeologic settings, the “watershed” approach to water-resource management may be limited, particularly where anthropogenic stresses alter the transport of soluble contaminants through highly permeable unconfined aquifers.     

3D model of fracture zones at Soultz-sous-Forêts based on geological data,image logs,induced microseismicity and vertical seismic profiles

This study presents a new deterministic 3D model of the fracture zones observed in the granitic reservoir of the Soultz European geothermal project. The major fracture zones encountered around 6 wells (4550, EPS1, GPK1, GPK2, GPK3 and GPK4) consist in 53 main structures that are located and characterized in terms of size and orientation: 39 fracture zones, 8 microseismic structures and 6 structures derived from vertical seismic profiles are represented in the 3D model using Discrete Fracture Network tools of the gOcad modelling platform (Paradigm?, Earth Decision?). This work illustrates the complexity of 3D fracture zone correlation and interpretation in crystalline rock masses characterized at meter scale (borehole) and at the reservoir scale (kilometer) thanks to geophysical imaging techniques.     

Stable and clumped isotopes in desert carbonate spring and lake deposits reveal palaeohydrology: A case study of the Lower Jurassic Navajo Sandstone,south‐western USA

Carbon, oxygen and clumped isotope (Δ₄₇) values were measured from lacustrine and tufa (spring)‐mound carbonate deposits in the Lower Jurassic Navajo Sandstone of southern Utah and northern Arizona in order to understand the palaeohydrology. These carbonate deposits are enriched in both ¹⁸O and ¹³C across the basin from east to west; neither isotope is strongly sensitive to the carbonate facies. However, ¹⁸O is enriched in lake carbonate deposits compared to the associated spring mounds. This is consistent with evaporation of the spring waters as they exited the mounds and were retained in interdune lakes. Clumped isotopes (Δ₄₇) exhibit minor systematic differences between lake and tufa‐mound temperatures, suggesting that the rate of carbonate formation under ambient conditions was moderate. These clumped isotope values imply palaeotemperature elevated beyond reasonable surface temperatures (54 to 86°C), which indicates limited bond reordering at estimated burial depths of ca 4 to 5 km, consistent with independent estimates of sediment thickness and burial depth gradients across the basin. Although clumped isotopes do not provide surface temperature information in this case, they still provide useful burial information and support interpretations of the evolution of groundwater locally. The findings of this study significantly extend the utility of combining stable isotope and clumped isotope methods into aeolian environments.     

Abundance and size distribution of transparent exopolymer particles (TEP) in a coccolithophorid bloom in the northern Bay of Biscay

The distribution of transparent exopolymer particles (TEP) was investigated during a coccolithophorid bloom in the northern Bay of Biscay (North Atlantic Ocean) in early June 2006. MODIS chlorophyll-a (Chl-a) and reflectance images before and during the cruise were used to localize areas of important biological activity and high reflectance (HR). TEP profiles along the continental margin, determined using microscopic (TEP_micro) and colorimetric (TEP_color) methods, showed abundant (6.1×10⁶–4.4×10⁷ L^?1) and relatively small (0.5–20 μm) particles, leading to a low total volume fraction (0.05–2.2 ppm) of TEP_micro and similar vertical profiles of TEP_color. Estimates of carbon content in TEP (TEP-C) derived from the microscopic approach yielded surface concentration of 1.50 μmol C L^?1. The contribution of TEP-C to particulate organic carbon (POC) was estimated to be 12% (molar C ratio) during this survey. Our results suggest that TEP formation is a probable first step to rapid and efficient export of C during declining coccolithophorid blooms.     

Stable carbon and oxygen isotopic evidence for late Pleistocene to middle Holocene climatic fluctuations in the interior of southern Africa

Stable carbon and oxygen isotope analyses of ungulate grazers from four archaeological sites located in different environs within the Caledon River Valley have provided a relatively well‐dated proxy palaeoenvironmental and palaeoclimatic sequence for the period between 16 000 and 6000 calendar (cal.) yr BP. Within the overall trend towards hot mid‐Holocene temperatures and a summer rainfall pattern, stable carbon isotope results show that there were three periods when growth season temperatures were cool enough for C₃ grasses to be present: 16 000–14 000; 10 200–9600, and 8400–8000 cal. yr BP. Similar trends were recorded in stable oxygen isotope values, reflecting shifts in either temperature or available moisture. Although having a similar pattern to that of the lower altitude site, sites situated in foothills and montane portions of the valley consistently maintained lower temperatures until the mid‐Holocene altithermal. At this time growth season temperatures warmed sufficiently for a 100% C₄ grassland to expand in altitude from the warmer low lying localities. In relation to present understanding of synoptic and global climatic patterning, these findings suggest that the early to middle Holocene transition was not a gradual warming trend, but rather it was marked by a series of climatic fluctuations. Of particular note is the possible global, rather than regional, occurrence of the 8200 cal. yr BP ‘event’. Copyright © 2002 John Wiley & Sons, Ltd.     

Deep vs. shallow controlling factors of the crustal thermal field – insights from 3D modelling of the Beaufort‐Mackenzie Basin (Arctic Canada)

Significant lateral variations in observed temperatures in the Beaufort‐Mackenzie Basin raise the question on the temperature‐controlling factors. Based on the structural configuration of the sediments and underlying crust in the area, we calculate the steady‐state 3D conductive thermal field. Integrated data include the base of the relic permafrost layer representing the 0 °C‐isotherm, public‐domain temperature data (from 227 wells) and thermal conductivity data. For >75% of the wells the predicted temperatures deviate by <10 K from the observed temperatures, which validates the overall model setup and adopted thermal properties. One important trend reproduced by the model is a decrease in temperatures from the western to the eastern basin. While in the west, a maximum temperature of 185 °C is reached at 5000 m below sea level, in the east the maximum temperature is 138 °C. The main cause for this pattern lies in lateral variations in thermal conductivity indicating differences in the shale and sand contents of the different juxtaposed sedimentary units. North‐to‐south temperature trends reveal the superposition of deep and shallow effects. At the southern margin, where the insulating effect of the low‐conductive sediments is missing, temperatures are lowest. Farther north, where the sub‐sedimentary continental crust is thick enough to produce considerable heat and a thick pile of sediments efficiently stores heat, temperatures tend to be highest. Temperatures decrease again towards the northernmost distal parts of the basin, where thinned continental and oceanic crust produce less radiogenic heat. Wells with larger deviations of the purely conductive model from the temperature observations (>15 K at 10% of the wells) and their basin‐wide pattern of misfit tendency (too cold vs. too warm temperature predictions) point to a locally restricted coupling of heat transport to groundwater flow.     

Vertical mixing of gravel over a long flood series

Vertical sediment exchange is a fundamental component of bedload transport in gravel‐bed channels. This paper describes the characteristic depth of exchange achieved over a long flood series. Analysis is based on 11 recoveries of magnetically tagged gravels deployed in Carnation Creek, Canada, completed between 1990 and 2008. Vertical grain exchange mixes gravels throughout the streambed relatively rapidly. Within one to eight floods the mean burial depth approaches two times the surface layer thickness, quantified by the 90th percentile of the size distribution. Finer gravels are mixed more rapidly into the bed than coarser gravels. Both active and passive grain exchanges throughout most of the bed produce the overall vertical distribution of marked grains. Gravel exchanges exhibit fairly consistent patterns once tracers are well mixed by large floods. Results highlight the role of flood sequence in determining exchange depths, support the notion of an upper limit to exchange, and underscore the importance of passive grain exchange. Copyright © 2011 John Wiley & Sons, Ltd.     

Why Web GIS May Not Be Enough: A Case Study with the Virtual Research Vessel

During several decades of investigation, the East Pacific Rise seafloor-spreading center at 9°-10°N has been explored by marine geologists, geophysicists, chemists, and biologists, and has emerged as one of the best studied sections of the global midocean ridge. It is an example of a region for which there is now a great wealth of observational data, results, and data-driven theoretical studies. However, these have yet to be fully utilized, either by research scientists or educators. While the situation is improving, a large amount of data, results, and related theoretical models still exist either in an inert, noninteractive form (e.g., journal publications) or as unlinked and currently incompatible computer data or algorithms. Presented here is the prototype of a computational environment and toolset, called the Virtual Research Vessel, to improve the situation by providing marine scientists and educators with simultaneous access to data, maps, and numerical models. While infrastructure is desired and needed for ready access to data and the resulting maps via web GIS in order to link disparate data sets (data to data), it is argued that data must also be linked to models for better exploration of new relations between observables, refinement of numerical simulations, and the quantitative evaluation of scientific hypotheses. For widespread data access, web GIS is therefore only a preliminary step rather than a final solution, and the ongoing implementation of the Virtual Research Vessel (scheduled for final completion in 2004-2005) is a case study for the midocean ridge community to test the effectiveness of moving beyond the data-to-data mode towards data-to-models and data-to-interpretation.