The celestial mechanics approach: application to data of the GRACE mission

The celestial mechanics approach (CMA) has its roots in the Bernese GPS software and was extensively used for determining the orbits of high-orbiting satellites. The CMA was extended to determine the orbits of Low Earth Orbiting satellites (LEOs) equipped with GPS receivers and of constellations of LEOs equipped in addition with inter-satellite links. In recent years the CMA was further developed and used for gravity field determination. The CMA was developed by the Astronomical Institute of the University of Bern (AIUB). The CMA is presented from the theoretical perspective in (Beutler et al. 2010). The key elements of the CMA are illustrated here using data from 50 days of GPS, K-Band, and accelerometer observations gathered by the Gravity Recovery And Climate Experiment (GRACE) mission in 2007. We study in particular the impact of (1) analyzing different observables [Global Positioning System (GPS) observations only, inter-satellite measurements only], (2) analyzing a combination of observations of different types on the level of the normal equation systems (NEQs), (3) using accelerometer data, (4) different orbit parametrizations (short-arc, reduced-dynamic) by imposing different constraints on the stochastic orbit parameters, and (5) using either the inter-satellite ranges or their time derivatives. The so-called GRACE baseline, i.e., the achievable accuracy of the GRACE gravity field for a particular solution strategy, is established for the CMA.     

Modeling time series of ground water head fluctuations subjected to multiple stresses

The methods behind the predefined impulse response function in continuous time (PIRFICT) time series model are extended to cover more complex situations where multiple stresses influence ground water head fluctuations simultaneously. In comparison to autoregressive moving average (ARMA) time series models, the PIRFICT model is optimized for use on hydrologic problems. The objective of the paper is twofold. First, an approach is presented for handling multiple stresses in the model. Each stress has a specific parametric impulse response function. Appropriate impulse response functions for other stresses than precipitation are derived from analytical solutions of elementary hydrogeological problems. Furthermore, different stresses do not need to be connected in parallel in the model, as is the standard procedure in ARMA models. Second, general procedures are presented for modeling and interpretation of the results. The multiple-input PIRFICT model is applied to two real cases. In the first one, it is shown that this model can effectively decompose series of ground water head fluctuations into partial series, each representing the influence of an individual stress. The second application handles multiple observation wells. It is shown that elementary physical knowledge and the spatial coherence in the results of multiple wells in an area may be used to interpret and check the plausibility of the results. The methods presented can be used regardless of the hydrogeological setting. They are implemented in a computer package named Menyanthes (www.menyanthes.nl).     

Effect of channel deepening on tidal flow and sediment transport—part II: muddy channels

Natural tidal channels often need deepening for navigation purposes (larger vessels). The depth increase may lead to tidal amplification, salt intrusion over longer distances, and increasing sand and mud import. Increasing fine sediment import, in turn, may start a process in which the sediment concentration progressively increases until the river becomes hyper-turbid, which may lead to increased dredging volumes and to decreased ecological values. These effects can be modeled and studied using detailed 3D models. Reliable simplified models for a first quick engineering evaluation are however lacking. In this paper, we apply both simplified and detailed 3D models to analyze the effects of channel deepening in prismatic and weakly converging tidal channels with saturated mud flow. The objective is to gain quantitative understanding of the effects of channel deepening on mud transport. We developed a simplified tidal mud model describing most relevant processes and effects in saturated mud flows with only minor horizontal transport gradients (quasi uniform conditions). The simplified model is not valid for non-saturated mud flow conditions. This model can either be used in standalone mode or in post-processing mode with computed near-bed velocities from a 3D hydrodynamic model as an input. The standalone model has been compared to various field data sets. Mud transport processes in the mouth region of muddy tidal channels can be realistically represented by the simplified model, if sufficient salinity and sediment data are available for calibration. The simulation of tidal mud transport and the behavior of an estuarine turbidity maximum (ETM) in saturated and non-saturated mud flow conditions cannot be represented by the simplified model and requires the application of a detailed 3D model.     

Experimental verification of effective anisotropic crack theories in variable crack aspect ratio medium

Physical modelling of cracked/fractured media using downscaled laboratory experiments has been used with great success as a useful alternative for understanding the effect of anisotropy in the hydrocarbon reservoir characterization and in the crustal and mantle seismology. The main goal of this work was to experimentally verify the predictions of effective elastic parameters in anisotropic cracked media by Hudson and Eshelby–Cheng's effective medium models. For this purpose, we carried out ultrasonic measurements on synthetic anisotropic samples with low crack densities and different aspect ratios. Twelve samples were prepared with two different crack densities, 5% and 8%. Three samples for each crack density presented cracks with only one crack aspect ratio, whereas other three samples for each crack density presented cracks with three different aspect ratios in their composition. It results in samples with aspect ratio values varying from 0.13 to 0.26. All the cracked samples were simulated by penny‐shaped rubber inclusions in a homogeneous isotropic matrix made with epoxy resin. Moreover, an isotropic sample for reference was constructed with epoxy resin only. Regarding velocity predictions performed by the theoretical models, Eshelby–Cheng shows a better fit when compared with the experimental results for samples with single and mix crack aspect ratio (for both crack densities). From velocity values, our comparisons were also performed in terms of the ε, γ, and δ parameters (Thomsen parameters). The results show that Eshelby–Cheng effective medium model fits better with the measurements of ε and γ parameters for crack samples with only one type of crack aspect ratio.     

Image sequence based automatic multi-camera system calibration techniques

An image sequence-based, fully automatic and rather flexible procedure for the calibration of stationary multi-camera systems for 3-D observation of dynamic events is presented and analysed. While conventional close-range camera calibration techniques are either based on a stable point-field with known reference coordinates or on a temporarily stationary point-field with only approximately known 3-D coordinates, which is imaged from different locations and under different orientations with one single camera, the presented technique is based on stationary cameras and moving targets, making use of the image sequence acquisition nature of most solid state cameras. In the simplest version, only one easily detectable marker has to be tracked through image sequences of multiple pre-calibrated cameras, thus avoiding the necessity of homologous feature identification for the establishment of multi-view correspondences; 3-D coordinates of the marker position are not required. This single-marker method does not allow for the determination of the interior orientation. In an extended version allowing for full camera orientation and calibration, a reference bar of known length is moved through object space, with the problem of feature identification and establishment of multi-view correspondences being reduced to the tracking of two targets. The method can only be used with multi-camera systems and is most useful for 3-D motion analysis applications, but may be adapted to a wide range of other applications. The advantages of the method over conventional self-calibration techniques are the trivial establishment of multi-view correspondences, the fact that no temporarily stable target field has to be constructed, and the fact that each camera has to be set up only once. After an explanation of the technique, its performance is examined in detail based on extensive computer simulations, and the practical effectiveness is shown in a pilot study on industrial robot calibration. Based on these studies, recommendations are given concerning the number of reference bar locations, preferable reference bar orientation schemes and the achievable accuracy potential.     

The shoreface-connected ridges along the central Dutch coast — part 2: morphological modelling

The hydrodynamic and sediment transport processes on the Dutch inner-shelf were studied on the basis of short- and medium-term (days, month) field experiments during fair-weather and storm conditions. The primary goal of this study was to identify the mechanisms that may be responsible for the formation and maintenance of the shoreface-connected ridges in the study area. The hydrodynamic field observations do not consistently support existing theoretical models explaining the formation of such ridges. Instead, the water motion appears to be dominated by wind- and density-driven effects. A review of literature indicates that other field studies also fail to observe the flow response expected from theory. Sediment transport during fair-weather is very episodic, with bed load transport slightly dominant over suspended load transport. The sediment transport processes during storms are dominated by the mean fluxes, with waves acting as a stirring mechanism. The contribution of wave-oscillatory fluxes cannot be neglected and may be directed with and against the waves. Long-wave fluxes are present, but very small.     

AMG preconditioning for sedimentary basin simulations in Temis calculator

Basin modelling aims at reconstructing the time evolution of a sedimentary basin in order to make quantitative predictions of geological phenomena leading to hydrocarbon accumulations. It accounts for porous medium compaction, heat transfer, hydrocarbon generation and migration. These physical phenomena are modelled by partial differential equations [Schneider, F., Wolf, S., Faille, I., Pot, D. 2000. A 3D basin model for hydrocarbon potential evaluation: application to Congo offshore. Oil and Gas Science Technologie, Rev-IFP 55, 3–13] representing the mass balances of solid, fluids (water and oil) coupled with Darcy's law and a compaction law. These equations are discretized using a cell-centered Finite Volume method in space and an implicit Euler integration in time. At each time step, the resulting nonlinear system is solved using Newton's method ending up at each Newton iteration with the solution of a linear system which represents the most cpu-time-consuming part of the basin simulator.     

The role of syntectonic sedimentation in the evolution of doubly vergent thrust wedges and foreland folds

This paper shows the evolution of time-constrained two-dimensional scaled analogue models of doubly vergent thrust systems in the presence of syntectonic sedimentation. Two sets of experiments were considered: (1) the addition of a syntectonic layer composed of a polymer and overlying sand in the prowedge; and (2) the addition to the previous condition of a progradational sedimentary load. Results from the first set of experiments indicate that the foreland fold-and-thrust belt has a strong relationship with the competence of the syntectonic layers. When the competence is low, the deformation produces tight asymmetric detachment-folds. As the competence increases, the fold-and-thrust belt shows break-trough folds with longer and better defined foreland-vergence. Results from the second set of experiments indicate that structural vergence is determined by the sense of progradation of the syntectonic layers, and in the case of strong aggradation at the prowedge, extension and reactive diapirism form contemporaneously with the contraction. Three end-members are proposed for mountain front thrust systems formed in the presence of syntectonic polymer and sand sedimentation: (1) outcropping fold-and-thrust belt sequence, in the case of no syntectonic sedimentation; (2) long displacement blind-thrust sheets, in the case of under-filled basins and (3) short displacement blind-thrust sheets, in the case of over-filled basins. All results indicate that ductile units at the base of syntectonic layers increase the displacement of the underlying frontal thrusts at the prowedge, and reduce the critical taper. Results also indicate that at very high sedimentary rates and hyper-critical taper conditions the prowedge collapses. Conclusions drawn from this research may be applied as an analogue to foreland evolution and to evaluate hydrocarbon generation, migration, and entrapment in thrust belts in areas where seismic imaging is generally poor.