An alternative method for non-negative estimation of variance components

A typical problem of estimation principles of variance and covariance components is that they do not produce positive variances in general. This caveat is due, in particular, to a variety of reasons: (1) a badly chosen set of initial variance components, namely initial value problem (IVP), (2) low redundancy in functional model, (3) an improper stochastic model, and (4) data’s possibility of containing outliers. Accordingly, a lot of effort has been made in order to design non-negative estimates of variance components. However, the desires on non-negative and unbiased estimation can seldom be met simultaneously. Likewise, in order to search for a practical non-negative estimator, one has to give up the condition on unbiasedness, which implies that the estimator will be biased. On the other hand, unlike the variance components, the covariance components can be negative, so the methods for obtaining non-negative estimates of variance components are not applicable. This study presents an alternative method to non-negative estimation of variance components such that non-negativity of the variance components is automatically supported. The idea is based upon the use of the functions whose range is the set of all positive real numbers, namely positive-valued functions (PVFs), for unknown variance components in stochastic model instead of using variance components themselves. Using the PVF could eliminate the effect of IVP on the estimation process. This concept is reparameterized on the restricted maximum likelihood with no effect on the unbiasedness of the scheme. The numerical results show the successful estimation of non-negativity estimation of variance components (as positive values) as well as covariance components (as negative or positive values).     

New data on the Late Cenozoic basaltic volcanism in Syria,applied to its origin

New data on geology and 21 K–Ar dates of the Late Oligocene–Quaternary basalts in Syria, combined with analysis of the new and previous data are used to reconstruct the volcanic history and relations between it and tectonic events. Volcanism began at the end of Oligocene (26–24 Ma) and was concentrated in the Late Oligocene–Early Miocene along a N-trending band, which stretches from the Jebel Arab (Harrat Ash Shaam) up to Kurd Dagh and southern Turkey. Activity waned in the Middle Miocene (17–12 Ma), but was resumed in the same band in the Tortonian and increased in the Messinian and Early Pliocene (6.3–4 Ma), when volcanism spread to the Shin Plateau and its coastal extension. After a brief hiatus ～ 4–3.5 Ma, volcanism became still more intensive and spread from the N-trending band to the east into the northern margin of the Mesopotamian Foredeep and to the west into the Dead Sea Transform zone. Additional eruptions continued into the Holocene.Volcanism lasted > 25 million years in the Jebel Arab Highland and > 15 million years in the Aleppo Plateau. The long duration of volcanism in the same parts of the moving Arabian plate and absence of records of one-way migration of the activity mean that the magmatic sources moved together with the plate, i.e., they were situated within the lithosphere mantle. Coincidence of the tectonic and volcanic stages of the Arabian plate development proves that volcanic activity depended on the geodynamic situation, caused by the plate motion. Situated within the lithosphere, magmatic sources within this transverse band were possibly caused by thermal and deforming influences of the asthenospheric lateral flow, moved laterally from the Ethiopia–Afar deep superplume.     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

Pore scale study of permeability and tortuosity for flow through particulate media using Lattice Boltzmann method

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic permeability of granular media. The obtained results show that particles' size diameter has a two‐fold effect on the coefficient of permeability: one is through specific surface and the other is by tortuosity factor. Numerical study also reveals that tortuosity of granular soils decreases almost linearly with increasing the porosity. Copyright © 2010 John Wiley & Sons, Ltd.     

Oceanographic variability in shallow-water acoustics and the dualrole of the sea bottom

Acoustic propagation in shallow water is examined. Multipath propagation and extensive boundary interactions, which along with a host of other phenomena produce a highly variable and often unpredictable acoustic field, are discussed. The responsible mechanisms, and hence the acoustic effects, cover a wide range of temporal and spatial scales and are classified as either deterministic or random, although the two types often act in concert. Because of extensive interactions with the sound field, the bottom can severely degrade waterborne propagation, although the sea bottom (and subbottom) can provide a seismic path that not only is relatively stable, but exists even under environmental conditions that preclude an effective waterborne path. Propagation in the bottom is particularly significant at very low frequencies. These various aspects of shallow-water acoustics are illustrated using the results of experiments conducted in diverse geographic areas     

<Emphasis Type="Italic">3D</Emphasis> Modelling of Urban Area Using Synthetic Aperture Radar (SAR)

Synthetic aperture radar (SAR) is a newly-developed remote sensing technology that works in all weather and independent of daylight. Recent satellite designs such as TerraSAR-x, which have resolutions of a couple of meters and sub-meters, have provided appropriate data for modelling and monitoring of urban areas. Image classification and height information extraction is possible considering the nature of SAR data. In this paper, a proper classification method for high-resolution SAR images has been used in urban areas. This classifier is based on statistical models. First, statistical models that are well adapted to urban SAR images are selected. Initial labelling is performed using the maximum likelihood method. A method based on Markov random fields is applied to improve the results by considering neighbourhood information. Meanwhile, topographic information is extracted using the phase difference obtained from SAR interferometry. After classification and height extraction, the homogeneous regions consisting of locations with similar objects are determined. The homogeneous region adjacency graph are generated using vectors containing classification information, extracted objects, height of pixels forming each region, and information on the neighbouring areas. Height and classification information are then merged by assigning height conditions based on the nature of objects and optimizing an energy function. The results obtained, including buildings, streets, and corner reflectors, are easily recognizable. The overall accuracy is improved from 57% in the initial classification to 95% in the employed procedure. Moreover, the accuracy of height estimation is about 2.74 m, which is acceptable for height estimations of buildings with more than one floor.     

Impacts of Depositional Facies and Diagenesis on Reservoir Quality: A Case Study from the Rudist-bearing Sarvak Formation, Abadan Plain, SW Iran

An integrated geological-petrophysical analysis of the rudist-bearing sequence of the Cretaceous Sarvak Formation is given one giant oilfield, and provides an improved understanding of this main reservoir in the Abadan Plain, in the Zagros Basin, SW Iran. The main objective of this study is to evaluate reservoir potential of the Sarvak Formation, and then to utilize the calibrated well log signature to correlate reservoir potential in un-cored wells. Eight main facies are recognized and categori...     

The added value of joint PP–PS inversion for reservoir characterization: A case study using Jubarte PRM seismic data (offshore Brazil)

The added value of the joint pre-stack inversion of PP (incident P-wave and reflected P-wave) and PS (incident P-wave and reflected S-wave) seismic data for the time-lapse application is shown. We focus on the application of this technique to the time-lapse (four-dimensional) multicomponent Jubarte field permanent reservoir monitoring seismic data. The joint inversion results are less sensitive to noise in the input data and show a better match with the rock physics models calibrated for the field. Further, joint inversion improves S-impedance estimates and provides a more robust quantitative interpretation, allowing enhanced differentiation between pore pressure and fluid saturation changes, which will be extremely useful for reservoir management. Small changes in reservoir properties are expected in the short time between the time-lapse seismic acquisitions used in the Jubarte project (only 1 year apart). The attempt to recover subtle fourth-dimensional effects via elastic inversion is recurrent in reservoir characterization projects, either due to the small sensitivity of the reservoirs to fluid and pressure changes or the short interval between the acquisitions. Therefore, looking for methodologies that minimize the uncertainty of fourth-dimensional inversion outputs is of fundamental importance. Here, we also show the differences between PP only and joint PP–PS inversion workflows and parameterizations that can be applied in other projects. We show the impact of using multicomponent data as input for elastic seismic inversions in the analysis of the time-lapse differences of the elastic properties. The larger investment in the acquisition and processing of multicomponent seismic data is shown to be justified by the improved results from the fourth-dimensional joint inversion.     

Modelling of Soil Degradation in Semi-arid Area Using Remote Sensing and GIS Techniques,Southern Jordan As Case Study

Doklady Earth Sciences - This research dealt with the topic of modeling the soil lost in a semi-arid desert area in the Ma’an watershed and its surroundings in Ma’an Governorate using...