Multi-scale aquifer characterization and groundwater flow model parameterization using direct push technologies

Direct push (DP) technologies are typically used for cost-effective geotechnical characterization of unconsolidated soils and sediments. In more recent developments, DP technologies have been used for efficient hydraulic conductivity (K) characterization along vertical profiles with sampling resolutions of up to a few centimetres. Until date, however, only a limited number of studies document high-resolution in situ DP data for three-dimensional conceptual hydrogeological model development and groundwater flow model parameterization. This study demonstrates how DP technologies improve building of a conceptual hydrogeological model. We further evaluate the degree to which the DP-derived hydrogeological parameter K, measured across different spatial scales, improves performance of a regional groundwater flow model. The study area covers an area of ~60 km² with two overlying, mainly unconsolidated sand aquifers separated by a 5–7 m thick highly heterogeneous clay layer (in north-eastern Belgium). The hydrostratigraphy was obtained from an analysis of cored boreholes and about 265 cone penetration tests (CPTs). The hydrogeological parameter K was derived from a combined analysis of core and CPT data and also from hydraulic direct push tests. A total of 50 three-dimensional realizations of K were generated using a non-stationary multivariate geostatistical approach. To preserve the measured K values in the stochastic realizations, the groundwater model K realizations were conditioned on the borehole and direct push data. Optimization was performed to select the best performing model parameterization out of the 50 realizations. This model outperformed a previously developed reference model with homogeneous K fields for all hydrogeological layers. Comparison of particle tracking simulations, based either on the optimal heterogeneous or reference homogeneous groundwater model flow fields, demonstrate the impact DP-derived subsurface heterogeneity in K can have on groundwater flow and solute transport. We demonstrated that DP technologies, especially when calibrated with site-specific data, provide high-resolution 3D subsurface data for building more reliable conceptual models and increasing groundwater flow model performance.     

The Pleistocene evolution of the East Antarctic Ice Sheet in the Prydz bay region: stable isotopic evidence from ODP Site 1167

Ocean Drilling Program Leg 188, Prydz Bay, East Antarctica is part of a larger initiative to explore the Cenozoic history of the Antarctic Ice Sheet through direct drilling and sampling of the continental margins. In this paper, we present stable isotopic results from Ocean Drilling Program (ODP) Site 1167 located on the Prydz Channel Trough Mouth Fan (TMF), the first Antarctic TMF to be drilled. The foraminifer-based δ¹⁸O record is interpreted along with sedimentary and downhole logging evidence to reconstruct the Quaternary glacial history of Prydz Bay and the adjacent Lambert Glacier Amery Ice Shelf System (LGAISS). We report an electron spin resonance age date of 36.9±3.3 ka at 0.45 m below sea floor and correlate suspected glacial–interglacial cycles with the global isotopic stratigraphy to improve the chronology for Site 1167. The δ¹⁸O record based on planktonic (Neogloboquadrina pachyderma (s.)) and limited benthic results (Globocassidulina crassa), indicates a trend of ice sheet expansion that was interrupted by a period of reduced ice volume and possibly warmer conditions during the early–mid-Pleistocene (0.9–1.38 Ma). An increase in δ¹⁸O values after 900 ka appears to coincide with the mid-Pleistocene climate transition and the expansion of the northern hemisphere ice sheet. The δ¹⁸O record in the upper 50 m of the stratigraphic section indicates as few as three glacial–interglacial cycles, tentatively assigned as marine isotopic stages (MIS) 16–21, are preserved since the Brunhes/Matuyama paleomagnetic reversal (780 ka). This suggests that there is a large unconformity near the top of the section and/or that there may have been few extreme advances of the ice sheet since the mid-Pleistocene climate transition resulting in lowered sedimentation rates on the Prydz Channel TMF. The stable isotopic record from Site 1167 is one of the few available from the area south of the Antarctic Polar Front that has been linked with the global isotopic stratigraphy. Our results suggest the potential for the recovery of useful stable isotopic records in other TMFs.     

StalAge – An algorithm designed for construction of speleothem age models

Here we present a new algorithm (StalAge), which is designed to construct speleothem age models. The algorithm uses U-series ages and their corresponding age uncertainty for modelling and also includes stratigraphic information in order to further constrain and improve the age model. StalAge is applicable to problematic datasets that include outliers, age inversions, hiatuses and large changes in growth rate. Manual selection of potentially inaccurate ages prior to application is not required. StalAge can be applied by the general, non-expert user and has no adjustable free parameters. This offers the highest degree of reproducibility and comparability of speleothem records from different studies. StalAge consists of three major steps. Firstly, major outliers are identified. Secondly, age data are screened for minor outliers and age inversions, and the uncertainty of potential outliers is increased using an iterative procedure. Finally, the age model and corresponding 95%-confidence limits are calculated by a Monte-Carlo simulation fitting ensembles of straight lines to sub-sets of the age data.We apply StalAge to a synthetic stalagmite ’sample’ including several problematic features in order to test its performance and robustness. The true age is mostly within the 95%-confidence age limits of StalAge showing that the calculated age models are accurate even for very difficult samples. We also apply StalAge to three published speleothem datasets. One of those is annually laminated, and the lamina counting chronology agrees with the age model calculated by StalAge. For the other two speleothems the resulting age models are similar to the published age models, which are both based on smoothing splines. Calculated uncertainties are in the range of those calculated by combined application of Bayesian chronological ordering and a spline, showing that StalAge is efficient in using stratigraphic information in order to reduce age model uncertainty.The algorithm is written in the open source statistical software R and available from the authors or as an electronic supplement of this paper.     

Correlation between hydrogen isotope ratios of lipid biomarkers and sediment maturity

We investigated the influence of thermal maturity on the hydrogen isotope ratios of sedimentary hydrocarbons to prove that the isotope ratio of hydrocarbons mirrors paleoclimate signatures. δD values from n-alkanes and acyclic isoprenoids of two sediment sections (Kupferschiefer [KS], 258 Ma, and Posidonienschiefer [PS], 184 Ma) with different maturation history were investigated. Both covered thermal maturity from 0.48 to 1.3 R_c (vitrinite reflectance and reflectance calculated from MPI1). Sediment burial up to 4500 m caused thermal maturation of organic matter in the KS horizon from the Early Zechstein basin of Poland, whereas contact metamorphic thermal maturation originated in the Early Toarcian PS (Posidonienschiefer) of the North German Vlotho Massif. The δD values of the extracted n-alkanes positively correlate with thermal maturity in the KS (y = 56‰ × MPI1[x] − 160‰ [VSMOW]) and in the PS (y = 104‰ × MPI1[x] − 200‰ [VSMOW]). The δD values of isoprenoids (i.e., pristane, phytane) were even more enriched with increasing maturity (y = 179‰ × MPI1[x] − 341‰ [VSMOW] in the KS; y = 300‰ × MPI1[x] − 415‰ [VSMOW] in PS).These results explain why isotope ratios of n-alkanes and isoprenoids in mature sediments are generally enriched in D and do not have the expected isotopic difference between n-alkanes and isoprenoids of ∼190‰. Moreover, the correlation between sediment maturity parameters and δD values suggests that after correction the δD values of n-alkanes can be used to reconstruct climate and environment in the geological past.     

The aeolian dust accumulation curve

This article presents a simple physical concept of aeolian dust accumulation, based on the behaviour of the subprocesses of dust deposition and dust erosion. The concept is tested in an aeolian dust wind tunnel. The agreement between the accumulation curve predicted by the model and the accumulation curve obtained in the experiments is close to perfect and shows that it is necessary to discriminate between the processes of aeolian dust deposition and aeolian dust accumulation. Two important thresholds determine the accumulation process. For wind speeds below the deflation threshold, the aeolian accumulation of dust increases linearly with the wind speed. For wind velocities between the deflation threshold and the accumulation limit, the sedimentation balance is above unity and there is still accumulation, though it rapidly drops once the deflation threshold has been exceeded. At wind speeds beyond the accumulation limit, the sedimentation balance is below unity and there will no longer be an accumulation of dust. The thresholds have been determined in a wind tunnel test at friction velocity u_* = 0·34 m s^?1 (deflation threshold) and u_* = 0·43 m s^?1 (accumulation limit), but these values are only indicative since they depend heavily on the characteristics of the accumulation surface and of the airborne grains. Copyright © 2001 John Wiley & Sons, Ltd.     

Herbig-Haro objects—tracers of the formation of low-mass stars and sub-stellar objects

Herbig-Haro objects (HHOs) are caused by outflows from young objects. Since the outflow relies on mass accretion from a circumstellar disk, it indicates ongoing growth. Recent results of infrared observations yielded evidence for disks around brown dwarfs. This suggests that at least a certain fraction of brown dwarfs forms like stars. Thus, young sub-stellar objects might cause HHOs as well. We present selected results of a general survey for HHOs based on DSS-II plates and CCD images taken with the Tautenburg Schmidt telescope. Numerous young objects could be identified due to their association with newly detected HHOs. In some cases the luminosity is consistent with very low-mass stars or close to sub-stellar values. This holds for L1415-IRS and a few infrared sources embedded in other dark clouds (e.g., GF9, BHR111). The question on the minimum mass for outflow activity is addressed.     

Pulsation,rotation and sunspot cycle

A suggestion is made that the periods for solar pulsation, solar rotation, and sunspot cycle may be closely related one to another.     

The Effect of Surface Heterogeneity on the Temperature–Humidity Correlation and the Relative Transport Efficiency

We derive a conceptual model of the flow over heterogeneous terrain consisting of patches with contrasting Bowen ratios. Upward moving eddies are assumed to carry heterogeneous properties, whereas downward moving eddies carry homogeneous properties. This results in a decorrelation of temperature and humidity as the contrast between the patches increases. We show that this model is able to reproduce the relationship developed by Lamaud and Irvine (Boundary-Layer Meteorol. 120:87–109, 2006). Some details differ from their expression but are in accordance with data obtained over African savannah. We extend the conceptual model to a combination of any scalars, not necessarily linked through the surface energy balance (as is the case for temperature and humidity). To this end we introduce a new parameter that describes the surface heterogeneity in surface fluxes. The results of the current model can be used to predict the discrepancy between similarity relationships for different scalars over heterogeneous terrain.