Application of well log tomography to the Dundee and Rogers City Limestones, Michigan Basin, USA

The Michigan Basin is located in the upper Midwest region of the United States and is centered geographically over the Lower Peninsula of Michigan. It is filled primarily with Paleozoic carbonates and clastics, overlying Precambrian basement rocks and covered by Pleistocene glacial drift. In Michigan, more than 46,000 wells have been drilled in the basin, many producing significant quantities of oil and gas since the 1920s in addition to providing a wealth of data for subsurface visualization. Well log tomography, formerly log-curve amplitude slicing, is a visualization method recently developed at Michigan Technological University to correlate subsurface data by utilizing the high vertical resolution of well log curves. The well log tomography method was first successfully applied to the Middle Devonian Traverse Group within the Michigan Basin using gamma ray log curves. The purpose of this study is to prepare a digital data set for the Middle Devonian Dundee and Rogers City Limestones, apply the well log tomography method to this data and from this application, interpret paleogeographic trends in the natural radioactivity. Both the Dundee and Rogers City intervals directly underlie the Traverse Group and combined are the most prolific reservoir within the Michigan Basin. Differences between this study and the Traverse Group include increased well control and “slicing” of a more uniform lithology. Gamma ray log curves for the Dundee and Rogers City Limestones were obtained from 295 vertical wells distributed over the Lower Peninsula of Michigan, converted to Log ASCII Standard files, and input into the well log tomography program. The “slicing” contour results indicate that during the formation of the Dundee and Rogers City intervals, carbonates and evaporites with low natural radioactive signatures on gamma ray logs were deposited. This contrasts the higher gamma ray amplitudes from siliciclastic deltas that cyclically entered the basin during Traverse Group deposition. Additionally, a subtle north-south, low natural radioactive trend in the center of the basin may correlate with previously published Dundee facies tracts. Prominent trends associated with the distribution of limestone and dolomite are not observed because the regional range of gamma ray values for both carbonates are equivalent in the Michigan Basin and additional log curves are needed to separate these lithologies.

Evaluation of the performance and cost-effectiveness of pavement sections containing open-graded base courses

Moisture induced distresses have been the prevalent distress type affecting the deterioration of both asphalt and concrete pavement sections. While various surface techniques have been employed over the years to minimize the ingress of moisture into the pavement structural sections, subsurface drainage components like open-graded base courses remain the best alternative in minimizing the time the pavement structural sections are exposed to saturated conditions. This research therefore focuses on assessing the performance and cost-effectiveness of pavement sections containing both treated and untreated open-graded aggregate base materials. Three common roadway aggregates comprising of two virgin aggregates and one recycled aggregate were investigated using four open-ended gradations and two binder types. Laboratory tests were conducted to determine the hydraulic, mechanical and durability characteristics of treated and untreated open-graded mixes made from these three aggregate types. Results of the experimental program show that for the same gradation and mix design types, limestone samples have the greatest drainage capacity, stability to traffic loads and resistance to degradation from environmental conditions like freeze-thaw. However, depending on the gradation and mix design used, all three aggregate types namely limestone, natural gravel and recycled concrete can meet the minimum coefficient of hydraulic conductivity required for good drainage in most pavements. Tests results for both asphalt and cement treated open-graded samples indicate that a percent air void content within the range of 15-25 will produce a treated open-graded base course with sufficient drainage capacity and also long term stability under both traffic and environmental loads. Using the new Mechanistic and Empirical Design Guide software, computer simulations of pavement performance were conducted on pavement sections containing these open-graded base aggregate base materials to determine how the MEPDG predicted pavement performance is sensitive to drainage. Using three truck traffic levels and four climatic regions, results of the computer simulations indicate that the predicted performance was not sensitive to the drainage characteristics of the open-graded base course. Based on the result of the MEPDG predicted pavement performance, the cost-effectiveness of the pavement sections with open-graded base was computed on the assumption that the increase service life experienced by these sections was attributed to the positive effects of subsurface drainage. The two cost analyses used gave two contrasting results with the one indicating that the inclusion of open-graded base courses can lead to substantial savings.

Estimation of vertical groundwater fluxes into a streambed through continuous temperature profile monitoring and the relationship of groundwater fluxes to coaster brook trout spawning habitat

We hypothesized that the spatial distribution of groundwater inflows through river bottom sediments is a critical factor associated with the selection of coaster brook trout (a life history variant of Salvelinus fontinalis,) spawning sites. An 80-m reach of the Salmon Trout River, in the Huron Mountains of the upper peninsula of Michigan, was selected to test the hypothesis based on long-term documentation of coaster brook trout spawning at this site. Throughout this site, the river is relatively similar along its length with regard to stream channel and substrate features. A monitoring well system consisting of an array of 27 wells was installed to measure subsurface temperatures underneath the riverbed over a 13-month period. The monitoring well locations were separated into areas where spawning has and has not been observed. Over 200,000 total temperature measurements were collected from 5 depths within each of the 27 monitoring wells. Temperatures within the substrate at the spawning area were generally cooler and less variable than river temperatures. Substrate temperatures in the non-spawning area were generally warmer, more variable, and closely tracked temporal variations in river temperatures. Temperature data were inverted to obtain subsurface groundwater velocities using a numerical approximation of the heat transfer equation. Approximately 45,000 estimates of groundwater velocities were obtained. Estimated velocities in the spawning and non-spawning areas confirmed that groundwater velocities in the spawning area were primarily in the upward direction, and were generally greater in magnitude than velocities in the non-spawning area. In the non-spawning area there was a greater occurrence of velocities in the downward direction, and velocity estimates were generally lesser in magnitude than in the spawning area. Both the temperature and velocity results confirm the hypothesis that spawning sites correspond to areas of significant groundwater influx to the river bed.

Underground Stratigraphy of Montana

The purpose of this thesis is to study the subsurface stratigraphy of the state of Montana from information acquired from well logs, composite logs, and measured sections. Underground stratigraphy is important for the intelligent prospecting for oil, gas, or water. A knowledge of the strata beneath us can be gathered only by measurements of numerous outcrops or by deep drilling with careful sampling or logging of the formations passed through.

USING AUTO- AND CROSS-CORRELATIONS FROM SEISMIC NOISE TO MONITOR VELOCITY CHANGES AT VILLARRICA VOLCANO, CHILE

We used the Green's functions from auto-correlations and cross-correlations of seismic ambient noise to monitor temporal velocity changes in the subsurface at Villarrica volcano in the Southern Andes of Chile. Campaigns were conducted from March to October 2010 and February to April 2011 with 8 broadband and 6 short-period stations, respectively. We prepared the data by removing the instrument response, normalizing with a root-mean-square method, whitening the spectra, and filtering from 1 to 10 Hz. This frequency band was chosen based on the relatively high background noise level in that range. Hour-long auto- and cross-correlations were computed and the Green's functions stacked by day and total time. To track the temporal velocity changes we stretched a 24 hour moving window of correlation functions from 90% to 110% of the original and cross correlated them with the total stack. All of the stations' auto-correlations detected what is interpreted as an increase in velocity in 2010, with an average increase of 0.13%. Cross-correlations from station V01, near the summit, to the other stations show comparable changes that are also interpreted as increases in velocity. We attribute this change to the closing of cracks in the subsurface due either to seasonal snow loading or regional tectonics. In addition to the common increase in velocity across the stations, there are excursions in velocity on the same order lasting several days. Amplitude decreases as the station's distance from the vent increases suggesting these excursions may be attributed to changes within the volcanic edifice. In at least two occurrences the amplitudes at stations V06 and V07, the stations farthest from the vent, are smaller. Similar short temporal excursions were seen in the auto-correlations from 2011, however, there was little to no increase in the overall velocity.

A Survey of Image-based Relighting Techniques

Image-based Relighting (IBRL) has recently attracted a lot of research interest for its ability to relight real objects or scenes, from novel illuminations captured in natural/synthetic environments. Complex lighting effects such as subsurface scattering, interreflection, shadowing, mesostructural self-occlusion, refraction and other relevant phenomena can be generated using IBRL. The main advantage of image-based graphics is that the rendering time is independent of scene complexity as the rendering is actually a process of manipulating image pixels, instead of simulating light transport. The goal of this paper is to provide a complete and systematic overview of the research in Imagebased Relighting. We observe that essentially all IBRL techniques can be broadly classified into three categories (Fig. 9), based on how the scene/illumination information is captured: Reflectance function-based, Basis function-based and Plenoptic function-based. We discuss the characteristics of each of these categories and their representative methods. We also discuss about the sampling density and types of light source(s), relevant issues of IBRL.

LakeMIP Kivu: evaluating the representation of a large, deep tropical lake by a set of one-dimensional lake models

The African great lakes are of utmost importance for the local economy (fishing), as well as being essential to the survival of the local people. During the past decades, these lakes experienced fast changes in ecosystem structure and functioning, and their future evolution is a major concern. In this study, for the first time a set of one-dimensional lake models are evaluated for Lake Kivu (2.28°S; 28.98°E), East Africa. The unique limnology of this meromictic lake, with the importance of salinity and subsurface springs in a tropical high-altitude climate, presents a worthy challenge to the seven models involved in the Lake Model Intercomparison Project (LakeMIP). Meteorological observations from two automatic weather stations are used to drive the models, whereas a unique dataset, containing over 150 temperature profiles recorded since 2002, is used to assess the model’s performance. Simulations are performed over the freshwater layer only (60 m) and over the average lake depth (240 m), since salinity increases with depth below 60 m in Lake Kivu and some lake models do not account for the influence of salinity upon lake stratification. All models are able to reproduce the mixing seasonality in Lake Kivu, as well as the magnitude and seasonal cycle of the lake enthalpy change. Differences between the models can be ascribed to variations in the treatment of the radiative forcing and the computation of the turbulent heat fluxes. Fluctuations in wind velocity and solar radiation explain inter-annual variability of observed water column temperatures. The good agreement between the deep simulations and the observed meromictic stratification also shows that a subset of models is able to account for the salinity- and geothermal-induced effects upon deep-water stratification. Finally, based on the strengths and weaknesses discerned in this study, an informed choice of a one-dimensional lake model for a given research purpose becomes possible.

Lake Floor Morphology and Sediment Architecture of Lake Torneträsk, Northern Sweden

ABSTRACT. Here we present datasets from a hydroacoustic survey in July 2011 at Lake Torneträsk, northern Sweden. Our hydroacoustic data exhibit lake floor morphologies formed by glacial erosion and accumulation processes, insights into lacustrine sediment accumulation since the beginning of deglaciation, and information on seismic activity along the Pärvie Fault. Features of glacial scouring with a high-energy relief, steep slopes, and relative reliefs of more than 50 m are observed in the large W-basin. The remainder of the lacustrine subsurface appears to host a broad variety of well preserved formations from glacial accumulation related to the last retreat of the Fennoscandian ice sheet. Deposition of glaciolacustrine and lacustrine sediments is focused in areas situated in proximity to major inlets. Sediment accumulation in distal areas of the lake seldom exceeds 2 m or is not observable. We assume that lack of sediment deposition in the lake is a result of different factors, including low rates of erosion in the catchment, a previously high lake level leading to deposition of sediments in higher elevated paleodeltas, tributaries carrying low suspension loads as a result of sedimentation in upstream lakes, and an overall low productivity in the lake. A clear off-shore trace of the Pärvie Fault could not be detected from our hydroacoustic data. However, an absence of sediment disturbance in close proximity to the presumed fault trace implies minimal seismic activity since deposition of the glaciolacustrine and lacustrine sediments.

Comparison of iron isotope variations in modern and Ordovician siliceous Fe oxyhydroxide deposits

Formation pathways of ancient siliceous iron formations and related Fe isotopic fractionation are still not completely understood. Investigating these processes, however, is difficult as good modern analogues to ancient iron formations are scarce. Modern siliceous Fe oxyhydroxide deposits are found at marine hydrothermal vent sites, where they precipitate from diffuse, low temperature fluids along faults and fissures on the seafloor. These deposits exhibit textural and chemical features that are similar to some Phanerozoic iron formations, raising the question as to whether the latter could have precipitated from diffuse hydrothermal fluids rather than from hydrothermal plumes. In this study, we present the first data on modern Fe oxyhydroxide deposits from the Jan Mayen hydrothermal vent fields, Norwegian-Greenland Sea. The samples we investigated exhibited very low δ56Fe values between -2.09‰ and -0.66‰. Due to various degrees of partial oxidation, the Fe oxyhydroxides are with one exception either indistinguishable from low-temperature hydrothermal fluids from which they precipitated (-1.84‰ and -1.53‰ in δ56Fe) or are enriched in the heavy Fe isotopes. In addition, we investigated Fe isotope variations in Ordovician jasper beds from the Løkken ophiolite complex, Norway, which have been interpreted to represent diagenetic products of siliceous ferrihydrite precursors that precipitated in a hydrothermal plume, in order to compare different formation pathways of Fe oxyhydroxide deposits. Iron isotopes in the jasper samples have higher δ56Fe values (-0.38‰ to +0.89‰) relative to modern, high-temperature hydrothermal vent fluids (ca. -0.40‰ on average), supporting the fallout model. However, formation of the Ordovician jaspers by diffuse venting cannot be excluded, due to lithological differences of the subsurface of the two investigated vent systems. Our study shows that reliable interpretation of Fe isotope variations in modern and ancient marine Fe oxyhydroxide deposits depends on comprehensive knowledge of the geological context. Furthermore, we demonstrate that very negative δ56Fe values in such samples might not be the result of microbial dissimilatory iron reduction, but could be caused instead by inorganic reactions.

Combined effect of heterogeneity, anisotropy and saturation on steady state flow and transport: Structure recognition and numerical simulation

1 Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.

Upper ocean vertical supply: A neglected primary factor controlling the distribution of neodymium concentrations of open ocean surface waters?

Neodymium (Nd) isotopes are an important geochemical tool to trace the present and past water mass mixing as well as continental inputs. The distribution of Nd concentrations in open ocean surface waters (0�100 m) is generally assumed to be controlled by lateral mixing of Nd from coastal surface currents and by removal through reversible particle scavenging. However, using 228Ra activity as an indicator of coastal water mass influence, surface water Nd concentration data available on key oceanic transects as a whole do not support the above scenario. From a global compilation of available data, we find that more stratified regions are generally associated with low surface Nd concentrations. This implies that upper ocean vertical supply may be an as yet neglected primary factor in determining the basin-scale variations of surface water Nd concentrations. Similar to the mechanism of nutrients supply, it is likely that stratification inhibits vertical supply of Nd from the subsurface thermocline waters and thus the magnitude of Nd flux to the surface layer. Consistently, the estimated required input flux of Nd to the surface layer to maintain the observed concentrations could be nearly two orders of magnitudes larger than riverine/dust flux, and also larger than the model-based estimation on shelf-derived coastal flux. In addition, preliminary results from modeling experiments reveal that the input from shallow boundary sources, riverine input, and release from dust are actually not the primary factors controlling Nd concentrations most notably in the Pacific and Southern Ocean surface waters.

The Value of Adding Optics to Ecosystem Models: A Case Study

Many ecosystem models have been developed to study the ocean's biogeochemical properties, but most of these models use simple formulations to describe light penetration and spectral quality. Here, an optical model is coupled with a previously published ecosystem model that explicitly represents two phytoplankton (picoplankton and diatoms) and two zooplankton functional groups, as well as multiple nutrients and detritus. Surface ocean color fields and subsurface light fields are calculated by coupling the ecosystem model with an optical model that relates biogeochemical standing stocks with inherent optical properties (absorption, scattering); this provides input to a commercially available radiative transfer model (Ecolight). We apply this bio-optical model to the equatorial Pacific upwelling region, and find the model to be capable of reproducing many measured optical properties and key biogeochemical processes in this region. Our model results suggest that non-algal particles largely contribute to the total scattering or attenuation (> 50% at 660 nm) but have a much smaller contribution to particulate absorption (< 20% at 440 nm), while picoplankton dominate the total phytoplankton absorption (> 95% at 440 nm). These results are consistent with the field observations. In order to achieve such good agreement between data and model results, however, key model parameters, for which no field data are available, have to be constrained. Sensitivity analysis of the model results to optical parameters reveals a significant role played by colored dissolved organic matter through its influence on the quantity and quality of the ambient light. Coupling explicit optics to an ecosystem model provides advantages in generating: (1) a more accurate subsurface light-field, which is important for light sensitive biogeochemical processes such as photosynthesis and photo-oxidation, (2) additional constraints on model parameters that help to reduce uncertainties in ecosystem model simulations, and (3) model output which is comparable to basic remotely-sensed properties. In addition, the coupling of biogeochemical models and optics paves the road for future assimilation of ocean color and in-situ measured optical properties into the models.

Seasonal Climatology of Hydrographic Conditions in the Upwelling Region Off Northern Chile

Over 30 years of hydrographic data from the northern Chile (18 degreesS-24 degreesS) upwelling region are used to calculate the surface and subsurface seasonal climatology extending 400 km offshore. The data are interpolated to a grid with sufficient spatial resolution to preserve cross-shelf gradients and then presented as means within four seasons: austral winter (July-September), spring (October-December), summer (January-March), and fall (April-June). Climatological monthly wind forcing, surface temperature, and sea level from three coastal stations indicate equatorward (upwelling favorable) winds throughout the year, weakest in the north. Seasonal maximum alongshore wind stress is in late spring and summer (December-March). Major water masses of the region are identified in climatological T-S plots and their sources and implied circulation discussed. Surface fields and vertical transects of temperature and salinity confirm that upwelling occurs year-round, strongest in summer and weakest in winter, bringing relatively fresh water to the surface nearshore. Surface geostrophic flow nearshore is equatorward throughout the year. During summer, an anticyclonic circulation feature in the north which extends to at least 200 m depth is evident in geopotential anomaly and in both temperature and geopotential variance fields. Subsurface fields indicate generally poleward flow throughout the year, strongest in an undercurrent near the coast. This undercurrent is strongest in summer and most persistent and organized in the south (south of 21 degreesS), A subsurface oxygen minimum, centered at similar to 250 m, is strongest at lower latitudes. Low-salinity subsurface water intrudes into the study area near 100 m, predominantly in offshore regions, strongest during summer and fall and in the southernmost portion of the region. The climatological fields are compared to features off Baja within the somewhat analogous California Current and to measurements from higher latitudes within the Chile-Peru Current system.

Ocean (De)oxygenation Across the Last Deglaciation: Insights for the Future

Anthropogenic warming is expected to drive oxygen out of the ocean as the water temperature rises and the rate of exchange between subsurface waters and the atmosphere slows due to enhanced upper ocean density stratification. Observations from recent decades are tantalizingly consistent with this prediction, though these changes remain subtle in the face of natural variability. Earth system model projections unanimously predict a long-term decrease in the global ocean oxygen inventory, but show regional discrepancies, particularly in the most oxygen-depleted waters, owing to the complex interplay between oxygen supply pathways and oxygen consumption. The geological record provides an orthogonal perspective, showing how the oceanic oxygen content varied in response to prior episodes of climate change. These past changes were much slower than the current, anthropogenic change, but can help to appraise sensitivities, and point toward potentially dominant mechanisms of change. Consistent with the model projections, marine sediments recorded an overall expansion of low-oxygen waters in the upper ocean as it warmed at the end of the last ice age. This expansion was not linearly related with temperature, though, but reached a deoxygenation extreme midway through the warming. Meanwhile, the deep ocean became better oxygenated, opposite the general expectation. These observations require that significant changes in apparent oxygen utilization occurred, suggesting that they will also be important in the future.

Erosive tooth wear: a multifactorial condition of growing concern and increasing knowledge.

Dental erosion is often described solely as a surface phenomenon, unlike caries where it has been established that the destructive effects involve both the surface and the subsurface region. However, besides removal of the surface, erosion shows dissolution of mineral within the softened layer - beneath the surface. In order to distinguish this process from the carious process it is now called 'near surface demineralization'. Erosion occurs in low pH, but there is no fixed critical pH value concerning dental erosion. The critical pH value for enamel concerning caries (pH 5.5-5.7) has to be calculated from calcium and phosphate concentrations of plaque fluid. In the context of dental erosion, the critical pH value is calculated from the calcium and phosphate concentrations in the erosive solution itself. Thus, critical pH for enamel with regard to erosion will vary according to the erosive solution. Erosive tooth wear is becoming increasingly significant in the management of the long-term health of the dentition. What is considered as an acceptable amount of wear is dependent on the anticipated lifespan of the dentition and is, therefore, different for deciduous compared to permanent teeth. However, erosive damage to the teeth may compromise the patient's dentition for their entire lifetime and may require repeated and increasingly complex and expensive restorations. Therefore, it is important that diagnosis of the tooth wear process in children and adults is made early and that adequate preventive measures are undertaken. These measures can only be initiated when the risk factors are known and interactions between them are present.