Use of a parsimonious rainfall–run‐off model for predicting hydrological response in ungauged basins

A new parameter parsimonious rainfall–run‐off model, the Distance Distribution Dynamics (DDD) model, is used to simulate hydrological time series at ungauged sites in the Lygne basin in Norway. The model parameters were estimated as functions of catchment characteristics determined by geographical information system. The multiple regression equations relating catchment characteristics and model parameters were trained from 84 calibrated catchments located all over Norway, and all model parameters showed significant correlations with catchment characteristics. The significant correlation coefficients (with p‐value < 0.05) ranged from 0.22 to 0.55. The suitability of DDD for predictions in ungauged basins was tested for 17 catchments not used to estimate the multiple regression equations. For ten of the 17 catchments, deviations in Nash–Sutcliffe efficiency (NSE) criteria between the calibrated and regionalised model were less than 0.1, and for two calibrated catchments within the Lygne basin, the deviations were less than 0.08. The median NSE for the regionalized DDD for the 17 catchments for two time series was 0.66 and 0.72. Deviations in NSE between calibrated and regionalised models are well explained by the deviations between calibrated and regressed parameters describing spatial snow distribution and snowmelt respectively. The quality of the simulated run‐off series for the ungauged sites in the Lygne basin was assessed by comparing flow indices describing high, medium and low flow estimated from observed run‐off at the 17 catchments and for the simulated run‐off series. The indices estimated for the simulated series were generally well within the ranges defined by the 17 observed series. Copyright © 2014 John Wiley & Sons, Ltd.     

Multipole borehole acoustic theory: Source imbalances and the effects of an elastic logging tool

In recent years the emphasis in acoustic logging has been shifting from the wireline to the Logging While Drilling (LWD) environment, the latter being far different from the former in that both tool rigidity and tool radius are considerably greater. In this paper we present a generic mathematical formulation for the multipole borehole acoustic measurement (alternate and equal polarity case), including a detailed analysis on the effects of multipole source amplitude imbalances. It is shown that source imbalance induced mode contaminants have excitation amplitudes that are scaled by the sum of the relative source imbalances between diametrically opposed sources. Furthermore it is shown that for mode contaminants with odd modal number there can be a significant offset in the associated directivity pattern, even at low levels of source imbalance. However, it is also shown that source imbalance induced mode contaminants can be completely eliminated if the multipole source is accompanied with a ‘vertically’ (but not azimuthally) offset multipole receiver. Mathematically, it is demonstrated (for a centered tool) that a polarity weighted stack of these multipole receivers completely eliminates the source imbalance induced mode contaminants. Excitation amplitudes and phase slowness of borehole guided modes are presented for the most common excitation regimes, i.e., monopole, dipole, quadrupole and the hexapole excitation, the latter one showing to have the advantage of a higher formation shear cut-off frequency than the quadrupole mode. Special emphasis will be on the analysis of the dipole excitation and the differences that occur due to variations in tool rigidity, tool diameter and (isotropic) formation properties with the resulting conclusion that a formation flexural mode is not observable as a result of a LWD dipole excitation (this opposed to its dipole wireline counterpart). The guided mode excitation amplitudes are calculated as residues using a Laurent series expansion. This (unconventional) way of calculating the residue has the advantage that it is independent of the pole order, does not require the numerical evaluation of derivatives with respect to the vertical wavenumber and allows for an accurate and efficient FFT implementation.     

Magnetostratigraphy of deep drilling core SG-1 in the western Qaidam Basin (NE Tibetan Plateau) and its tectonic implications

The Qaidam Basin is the largest intermontane basin of the northeastern Tibetan Plateau and contains a continuous Cenozoic sequence of lacustrine sediments. A ~ 1000-m-deep drilling (SG-1) with an average core recovery of ~ 95% was carried out in the depocenter of the Chahansilatu playa (sub-depression) in the western Qaidam Basin, aimed to obtain a high-resolution record of the paleoenvironmental evolution and the erosion history. Stepwise alternating field and thermal demagnetization, together with rock magnetic results, revealed a stable remanent magnetization for most samples, carried by magnetite. The polarity sequence consisted of 16 normal and 15 reverse zones which can be correlated with chrons 1n to 2An of the global geomagnetic polarity time scale. Magnetostratigraphic results date the entire core SG-1 at ~ 2.77 Ma to ~ 0.1 Ma and yielded sediment accumulation rate (SAR) ranging from 26.1 cm/ka to 51.5 cm/ka. Maximum SARs occurred within the intervals of ~ 2.6–2.2 Ma and after ~ 0.8 Ma, indicating two episodes of erosion, which we relate to pulse tectonic uplift of the NE Tibetan Plateau with subsequent global cooling.     

Evaluation of thermal equations of state for CO2 in numerical simulations

Three commonly used thermal equations of state for carbon dioxide, as well as the ideal gas law, have been compared against a large number of measurement data taken from the literature. Complex equations of state reach a higher accuracy than simple ones. The inaccuracy of the density function can cause large errors in fluid property correlations, such as heat capacity or viscosity. The influence of this inaccuracy on the results of numerical simulations have been evaluated by two examples: The first one assumes isothermal gas expansion from a volume, while the second one considers heat transport along a fracture. For both examples, different equations of state have been utilized. The simulations have been performed on the scientific software platform OpenGeoSys. The difference among the particular simulation results is significant. Apparently small errors in the density function can cause considerably different results of otherwise identical simulation setups.     

Toward physical aspects affecting a possible leakage of geologically stored CO2 into the shallow subsurface

In geological formations, migration of CO₂ plume is very complex and irregular. To make CO₂ capture and storage technology feasible, it is important to quantify CO₂ amount associated with possible leakage through natural occurring faults and fractures in geologic medium. Present work examines the fracture aperture effect on CO₂ migration due to free convection. Numerical results reveal that fracture with larger-aperture intensify CO₂ leakage. Mathematical formulation and equations of state for the mixture are implemented within the object-oriented finite element code OpenGeoSys developed by the authors. The volume translated Peng–Robinson equation of state is used for material properties of CO₂ and water.     

Concept and workflow for 3D visualization of atmospheric data in a virtual reality environment for analytical approaches

In the future, climate change will strongly influence our environment and living conditions. Weather and Climate simulations that predict possible changes produce big data sets. The combination of various variables of climate models with spatial data from different sources helps to identify correlations and to study key processes. In this paper, the results of the Weather Research and Forecasting model are visualized for two regions. For this purpose, a continuous workflow that leads from the integration of heterogeneous raw data to 3D visualizations that can be displayed on a desktop computer or in an interactive virtual reality environment is developed. These easy-to-understand visualizations of complex data are the basis for scientific communication and for the evaluation and verification of models as well as for interdisciplinary discussions of the research results.     

A nested-cell approach for in situ remediation

We characterize the hydraulics of an extraction-injection well pair in arbitrarily oriented regional flow by the recirculation ratio, area, and average residence time in the recirculation zone. Erratic regional flow conditions may compromise the performance of the reactor between a single well pair. We propose an alternative four-well system: two downgradient extraction and two upgradient injection wells creating an inner cell nested within an outer cell. The outer cell protects the inner cell from the influence of regional flow. Compared to a two-well system, the proposed four-well system has several advantages: (1) the recirculation ratio within the nested inner cell is less sensitive to the regional flow direction; (2) a transitional recirculation zone between the inner and outer cells can capture flow leakage from the inner cell, minimizing the release of untreated contaminants; and (3) the size of the recirculation zone and residence times can be better controlled within the inner cell by changing the pumping rates. The system is applied at the Field Research Center in Oak Ridge, Tennessee, where experiments on microbial in situ reduction of uranium (VI) are under way.     

Analyzing bank filtration by deconvoluting time series of electric conductivity

Knowing the travel-time distributions from infiltrating rivers to pumping wells is important in the management of alluvial aquifers. Commonly, travel-time distributions are determined by releasing a tracer pulse into the river and measuring the breakthrough curve in the wells. As an alternative, one may measure signals of a time-varying natural tracer in the river and in adjacent wells and infer the travel-time distributions by deconvolution. Traditionally this is done by fitting a parametric function such as the solution of the one-dimensional advection-dispersion equation to the data. By choosing a certain parameterization, it is impossible to determine features of the travel-time distribution that do not follow the general shape of the parameterization, i.e., multiple peaks. We present a method to determine travel-time distributions by nonparametric deconvolution of electric-conductivity time series. Smoothness of the inferred transfer function is achieved by a geostatistical approach, in which the transfer function is assumed as a second-order intrinsic random time variable. Nonnegativity is enforced by the method of Lagrange multipliers. We present an approach to directly compute the best nonnegative estimate and to generate sets of plausible solutions. We show how the smoothness of the transfer function can be estimated from the data. The approach is applied to electric-conductivity measurements taken at River Thur, Switzerland, and five wells in the adjacent aquifer, but the method can also be applied to other time-varying natural tracers such as temperature. At our field site, electric-conductivity fluctuations appear to be an excellent natural tracer.     

Determination of first-order degradation rate constants from monitoring networks

In this article, different strategies for estimating first-order degradation rate constants from measured field data are compared by application to multiple, synthetic, contaminant plumes. The plumes were generated by numerical simulation of contaminant transport and degradation in virtual heterogeneous aquifers. These sites were then individually and independently investigated on the computer by installation of extensive networks of observation wells. From the data measured at the wells, that is, contaminant concentrations, hydraulic conductivities, and heads, first-order degradation rates were estimated by three 1D centerline methods, which use only measurements located on the plume axis, and a two-dimensional method, which uses all concentration measurements available downgradient from the contaminant source. Results for both strategies show that the true rate constant used for the numerical simulation of the plumes in general tends to be overestimated. Overestimation is stronger for narrow plumes from small source zones, with an average overestimation factor of about 5 and single values ranging from 0.5 to 20, decreasing for wider plumes, with an average overestimation factor of about 2 and similar spread. Reasons for this overestimation are identified in the velocity calculation, the dispersivity parameterization, and off-centerline measurements. For narrow plumes, the one- and the two-dimensional strategies show approximately the same amount of overestimation. For wider plumes, however, incorporation of all measurements in the two-dimensional approach reduces the estimation error. No significant relation between the number of observation wells in the monitoring network and the quality of the estimated rate constant is found for the two-dimensional approach.