Dual models of pore spaces

We present a model for pore spaces that consists of two parts related by duality: (1) a decomposition of an open polyhedral pore space into open contractible pore bodies separated by relatively open interfaces and (2) a pore network that is homotopy equivalent to the pore space. The dual model is unique and free of parameters, but it relies on regularity conditions for the pore space. We show how to approximate any pore space by the interior of a polyhedral complex such that the regularity conditions are fulfilled. Thus, we are able to calculate the dual model from synthetic porous media and images of real porous media. The pore bodies are unions of relatively open Delaunay cells with respect to the corners of the pore boundary, and the pore network consists of certain at most two-dimensional (2D) Voronoi cells with respect to the corners of the pore boundary. The pore network describes the neighborhood relations between the pore bodies. In particular, any relatively open 2D Delaunay face f separating two pore bodies has a unique (relatively open) dual network edge. In our model, f is a pore throat only if it is hit by its dual network edge. Thus, as opposed to widespread intuition, any pore throat is convex, and adjacent pore bodies are not necessarily separated by pore throats. Due to the duality between the pore network and the decomposition of the pore space into pore bodies it is straightforward to store the geometrical properties of the pore bodies [pore throats] as attributes of the dual network vertices [edges]. Such an attributed network is used to perform 2D drainage simulations. The results agree very well with those from a pore-morphology based modeling approach performed directly on the digital image of a porous medium. Contractibility of the pore bodies and homotopy equivalence of the pore space and the pore network is proven using discrete Morse theory and the nerve theorem from combinatorial topology.     

Evaluating the performance of ground-based and remotely sensed near real-time rainfall fields from a hydrological perspective

Abstract

The South African Weather Service (SAWS) issues routine experimental, near real-time rainfall maps from daily raingauge networks, radar networks and satellite images, as well as merged rainfall fields. These products are potentially useful for near real-time forecasting, especially in areas of fast hydrological response, and also to simulate the “now state” of various hydrological state variables such as soil moisture content, streamflow, and reservoir inflows. The purpose of this paper is to evaluate their skill as inputs to hydrological simulations and, in particular, the skill of the merged field in terms of better hydrological results relative to the individual products. Rainfall fields derived from raingauge, radar, satellite, conditioned satellite and the merged (gauge/radar/satellite) were evaluated for two selected days with relatively high amounts of rainfall, as well as for a continuous period of 90 days in the Mgeni catchment, South Africa. Streamflows simulated with the ACRU model indicate that the use of raingauge as well as merged fields of satellite/raingauge and satellite/radars/raingauge provides relatively realistic rainfall results, without much difference in their hydrological outputs, whereas the radar and raw satellite information by themselves cannot be used in operational hydrological application in their current status.

Citation Ghile, Y., Schulze, R. & Brown, C. (2010) Evaluating the performance of ground-based and remotely sensed near real-time rainfall fields from a hydrological perspective. Hydrol. Sci. J. 55(4), 497–511.     

Comparative hydrogeology – reference analysis of groundwater dynamics from neighbouring observation wells

ABSTRACT

Groundwater level fluctuations are caused by spatial and temporal superposition of processes within and outside the aquifer system. Most of the subsurface processes are usually observed on a small scale. Upscaling to the regional scale, as required for future climate change scenarios, is difficult due to data scarcity and increasing complexity. In contrast to the limited availability of system characteristics, high-resolution data records of groundwater hydrographs are more generally available. Exploiting the information contained in these records should thus be a priority for analysis of the chronical lack of data describing groundwater system characteristics. This study analyses the applicability of 63 indices derived from daily hydrographs to quantify different dynamics of groundwater levels in unconfined gravel aquifers from three groundwater regions (Bavaria, Germany). Based on the results of two different skill tests, the study aids index selection for different dynamic components of groundwater hydrographs.     

Analysing Debris-Flow Impact Models, Based on a Small Scale Modelling Approach

The objective of this study is to analyse adaptable debris-flow impact models, which are very important for mitigation measurements and buildings using their sphere of influence. For this reason, 16 debris-flow experiments, on a small-scale modelling approach, were performed. Impact forces were measured with a force plate panel, consisting of 24 aluminium devices, coaxially mounted with resistance strain gauges. Flow velocities, flow heights as well as horizontal impact forces were sampled with a frequency of 2.4 kHz. Sub datasets of sampled raw force data were defined by applying an average median filter, a low-pass filter routine. Further, estimated peak pressure values as well as empirical coefficients of hydraulic impact models were compared, and the influence of signal processing is discussed.     

Gravimetric monitoring of the first field‐wide steam injection in a fractured carbonate field in Oman – a feasibility study

Gas‐Oil Gravity Drainage is to be enhanced by steam injection in a highly fractured, low permeability carbonate field in Oman. Following a successful pilot, field‐wide steam injection is being implemented, first of this type in the world. A dedicated monitoring program has been designed to track changes in the reservoir. Various observations are to be acquired, including, surface deformation, temperature measurements, microseismic, well logs, pressure and saturation measurements to monitor the reservoir. Because significant changes in the reservoir density are expected, time‐lapse gravimetry is also being considered. In this paper we investigate the feasibility of gravimetric monitoring of the thermally enhanced gravity drainage process at the carbonate field in Oman. For this purpose, forward gravity modelling is performed. Based on field groundwater measurements, the estimates of the hydrological signal are considered and it is investigated under what conditions the groundwater influences can be minimized. Using regularized inversion of synthetic gravity data, we analyse the achievable accuracy of heat‐front position estimates. In case of large groundwater variations at the field, the gravity observations can be significantly affected and, consequently, the accuracy of heat‐front monitoring can be deteriorated. We show that, by applying gravity corrections based on local observations of groundwater, the hydrological influences can to a large extent be reduced and the accuracy of estimates can be improved. We conclude that gravimetric monitoring of the heat‐front evolution has a great potential.     

The burial efficiency of organic carbon in the sediments of Lake Kinneret

Even though lake sediments constitute a significant long-term carbon sink, studies on the regulation of carbon burial in lakes sediments have, to date, been surprisingly few. We investigated to what degree the organic carbon (OC) being deposited onto the bottom of Lake Kinneret (Israel) is buried in the sediment at four different sites with varying degrees of oxygenation and varying supply of allochthonous particles from the River Jordan. For estimation of the OC burial efficiency (OC BE), i.e., the ratio between buried and deposited OC, we calculated OC burial from dated sediment cores, and calculated OC deposition using three different approaches. Calculation of OC deposition from sediment trap-derived mass deposition rates multiplied with the OC content of surface sediment yielded OC BE values that were at odds with published values for sediments dominated by autochthonous OC sources. Calculation via sediment trap data on organic matter flux collected within the Lake Kinneret monitoring program, as well as calculation of OC deposition as the sum of OC burial plus OC mineralization, returned fairly congruent estimates of OC BE (range 10–41%), but only if the sediment trap data were corrected for the proportion of resuspended particles in the traps. Differences in OC BE between sites were small, indicating that OC source (common to all sites) was a more important regulator of OC BE in Lake Kinneret than oxygen exposure or mineral particles characteristics.     

Revisiting Force Model Error Modeling in GRACE Gravity Field Recovery

Surveys in Geophysics - The gravity field recovery from GRACE (Gravity Recovery and Climate Experiment) mission data is contaminated by both observation noise and dynamic force errors, especially...     

Interactive effects of environmental variability and human impacts on the long-term dynamics of an Amazonian floodplain lake and a South Atlantic coastal lagoon

Human activities are exposing freshwater ecosystems to a wide range of stressors, whose direct and indirect effects can be alleviated or exacerbated through interactive effects with dynamic environmental drivers. This study used long-term data from two Neotropical lacustrine freshwater systems (Batata Lake, an Amazonian floodplain lake and Imboassica lagoon, an Atlantic coastal lagoon) subjected to different kinds of environmental fluctuations (i.e., flood pulse and sandbar opening) and anthropogenic impacts (i.e., siltation and eutrophication). Our objective was to determine whether the effects of human perturbations are contingent on modifications of important biotic and abiotic characteristics through environmental variability. For both ecosystems, environmental variability consistently interacted with anthropogenic perturbations to alter most of the variables analyzed, such as nutrient dynamics, chlorophyll-a concentration, zooplankton and benthic invertebrate species richness, and temporal community stability, which indicates that interactive effects between environmental variability and anthropogenic perturbations may impact a myriad of ecosystem properties. Furthermore, the nature of these interactive effects was highly dependent on the variable considered and on the ecosystem analyzed. For example, at Imboassica lagoon, sandbar openings interacted synergistically with trophic state to increase the phosphorus concentration in the water column. At Batata Lake, flooding generally alleviated the negative effects of siltation on species richness by both diluting inorganic suspended material concentration and by promoting local recruitment from the regional species pool. Such results indicate that our ability to understand and predict the outcome of anthropogenic impacts on inland aquatic systems can be hampered if we consider human stressors as “static” phenomena disconnected from dynamic interactions with major local environmental drivers.     

Gravity field determination using boundary element methods

The Boundary Element Method (BEM), a numerical technique for solving boundary integral equations, is introduced to determine the earth's gravity field. After a short survey on its main principles, we apply this method to the fixed gravimetric boundary value problem (BVP), i.e. the determination of the earth's gravitational potential from measurements of the intensity of the gravity field in points on the earth's surface. We show how to linearize this nonlinear BVP using an implicit function theorem and how to transform the linearized BVP into a boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation using the single layer representation. A Galerkin method is used to transform the boundary integral equation into a linear system of equations. We discuss the major problems of this approach for setting up and solving the linear system. The BVP is numerically solved for a bounded part of the earth's surface using a high resolution reference gravity model, measured gravity values of high density, and a 50 50 m² digital terrain model to describe the earth's surface. We obtain a gravity field resolution of 1 1 km² with an accuracy of the order 10^–3 to 10^–4 in about 1 CPU-hour on a Siemens/Fujitsu SIMD vector pipeline machine using highly sophisticated numerical integration techniques and fast equation solvers. We conclude that BEM is a powerful numerical tool for solving boundary value problems and may be an alternative to classical geodetic techniques.