Prediction of Remediation of a Heterogeneous Aquifer: A Case Study

Contaminant plumes whose characteristic length is smaller than the horizontal integral scale of the hydraulic conductivity, K, are abundant in shallow, phreatic aquifers. In such cases, the aquifer can be regarded as layered, with K being only a function of the vertical coordinate. The heterogeneity of K has a critical role upon the efficiency of remediation of such sites, for example, by Pump and Treat schemes. The expected efficiency is a random variable, with uncertainty. Quantifying this uncertainty can be of great importance to decision making. In this study, we focus on a case study in the coastal aquifer of Israel and compare two different approaches for constructing realizations of K: continuous and indicator. We observe a significant difference between the constructed realizations, which results in a considerable difference in the predicted remediation efficiency and its uncertainty. Furthermore, we study the effect of conditioning the realizations by a rather limited number of K data points. We find that the conditioning results in a major reduction of the uncertainty. In addition, we compare the results of the transport model to a simplified semi‐analytical solution that is based on assuming radial flow. We find a good agreement with the three‐dimensional numerical model. This result illustrates that the simplified solution can be used for prediction of the remediation efficiency when the flow at the plume vicinity can be regarded as radial.     

S-system theory applied to array-based GNSS ionospheric sensing

The GPS carrier-phase and code data have proven to be valuable sources of measuring the Earth’s ionospheric total electron content (TEC). With the development of new GNSSs with multi frequency data, many more ionosphere-sensing combinations of different precision can be formed as input of ionospheric modelling. We present the general way of interpreting such combinations through an application of S-system theory and address how their precision propagates into that of the unbiased TEC solution. Presenting the data relevant to TEC determination, we propose the usage of an array of GNSS antennas to improve the TEC precision and to expedite the rather long observational time-span required for high-precision TEC determination.     

Sparse Bayesian linearized amplitude-versus-angle inversion

The technique of seismic amplitude-versus-angle inversion has been widely used to estimate lithology and fluid properties in seismic exploration. The amplitude-versus-angle inversion problem is intrinsically ill-posed and generally stabilized by the use of L2-norm regularization methods but with drawback of smoothing important boundaries between adjacent layers. In this study, we propose a sparse Bayesian linearized solution for amplitude-versus-angle inversion problem to preserve the sharp geological interfaces. In this regard, a priori constraint term with two regularization functions is presented: the sparse constraint regularization and the low-frequency model information. In addition, to obtain high-resolution reflectivity estimation, the model parameters decorrelation technique combined with dipole decomposition method is employed. We validate the applicability of the presented method by both synthetic and real seismic data from the Gulf of Mexico. The accuracy improvement of the presented method is also confirmed by comparing the results with the commonly used Bayesian linearized amplitude-versus-angle inversion.     

A hybrid optimization scheme for self-potential measurements due to multiple sheet-like bodies in arbitrary 2D resistivity distributions

Self-potential is a passive geophysical method that can be applied in a straightforward manner with minimum requirements in the field. Nonetheless, interpretation of self-potential data is particularly challenging due to the inherited non-uniqueness present in all potential methods. Incorporating information regarding the target of interest can facilitate interpretation and increase the reliability of the final output. In the current paper, a novel method for detecting multiple sheet-like targets is presented. A numerical framework is initially described that simulates sheet-like bodies in an arbitrary 2D resistivity distribution. A scattered field formulation based on finite differences is employed that allows the edges of the sheet to be independent of the grid geometry. A novel analytical solution for two-layered models is derived and subsequently used to validate the accuracy of the proposed numerical scheme. Lastly, a hybrid optimization is proposed that couples linear least-squares with particle-swarm optimization in order to effectively locate the edges of multiple sheet-like bodies. Through numerical and real data, it is proven that the hybrid optimization overcomes local minimal that occurs in complex resistivity distributions and converges substantially faster compared to traditional particle-swarm optimization.     

Drag Optimization of A Planing Vessel Based on the Stability Criteria Limits

China Ocean Engineering - One of the important phases of designing a craft is stability analysis and optimization of the drag force in cruise speed. In this research the longitudinal and lateral...     

The Ohalo II prehistoric camp (19.5 Ky): New evidence for environmental and tectonic changes at the Sea of Galilee

Combined archaeological data, shore surveys, and aerial photos of submerged sediments in the Sea of Galilee provide new insights into environmental and tectonic events, their dating, and their impact on the Ohalo II prehistoric camp (ca. 19,500 yr B.P.) and its surroundings. The Ohalo II waterlogged campsite contains excellently preserved brush hut remains and other in situ features, all embedded in late Pleistocene lacustrine strata. The findings indicate relatively short occupation of the site, not more than months or several years at a time. The high quality in situ preservation of delicate organic materials, as well as the short occupation period, suggests a quick and gentle burial by fine sediments. The evident fast submergence (water level rise of the Sea of Galilee) could have been the result of climatic fluctuations towards the end of the last glaciation and/or small‐scale tectonic subsidence. The site is located on a tectonic block formed in the western fault belt of the Dead Sea Rift. We present new evidence of post‐occupational folding of the late Pleistocene strata and recent tilting and faulting. A westward tectonic tilt may have caused the blockage of the old Jordan River outlet after A. D. 1106. Excellent preservation of the fault traces to the east of the site is attributed to the young age of the displacement on the fault. The last displacement apparently post‐dates the blockage of the old Jordan River. © 2002 Wiley Periodicals, Inc.