Stochastic analysis of unsaturated seepage through randomly heterogeneous earth embankments

Spatial variability of material properties is inherent in both natural soil deposits and earth structures, yet it is often ignored during geotechnical design. With the objective of developing novel methods for assessing the effects of soil variability on groundwater flow, this study presents a stochastic finite element model of seepage through a flood defense embankment with randomly heterogeneous material properties. Stochastic modeling is undertaken by means of a Monte Carlo simulation which involves a large number of finite element analyses, each with randomly varied porosity at element level, which leads to a corresponding random variation of both permeability and water retention properties across the embankment domain. This provides a statistical distribution of responses, such as total flow rate and time to reach steady state, instead of a single deterministic result as in conventional studies of seepage through unsaturated heterogeneous soils. As the degree of heterogeneity increases, water tends to flow along the most permeable paths inside the soil mass, resulting in an irregular shape of the predicted wetting fronts and pore pressure contours. The mean and standard deviation of the computed quantities strongly depend on the statistics of the input porosity field. Simulations are also conducted to compare the statistical variation of flow rate with and without dependency of the water retention curve on porosity. With recent growth in computer speed, stochastic finite element models based on the Monte Carlo approach can become a powerful design tool, especially if a quantitative assessment of geotechnical risks is required. Copyright © 2011 John Wiley & Sons, Ltd.     

Constraining the cosmological time variation of the fine - structure constant

The variation of the fine-structure constant α = e ² / ħc can be probed by comparing the wavelength of atomic transitions from the redshift of quasars in the Universe and laboratory over cosmological time scales t ~ 10¹⁰ yr. After a careful selection of pairs of lines, the Thong method with a derived analytical expression for the error analysis was applied to compute the α variation. We report a new constraint on the variation of the fine-structure constant based on the analysis of the C_IV, N_V, Mg_II, Al_III, and Si_IV doublet absorption lines. The weighted mean value of the variation in α derived from our analysis over the redshift range 0.4939 ≤ z ≤ 3.7 is = ( 0.09 ± 0.07)×10⁻⁵. This result is three orders of magnitude better than the results obtained by earlier analysis of the same data on the constraint on Δα/α .     

Moisture Movement Through Cracked Clay Soil Profiles

A continuum mechanics approach is used for the formulation of unsaturated hydraulic conductivity functions and the water storage functions for fractured or cracked clay soils in this parametric study. Suggested procedures are based on available research literature related to the behavior of cracked unsaturated porous media. The soil–water characteristic curve, hydraulic conductivity and water storage functions take on the character of bi-modal unsaturated soil property functions. The bimodal character arises out of the independent behavior of the cracks and the intact clay soil. Matric suction changes beneath a slab-on-grade foundation placed on a cracked clay soil profile are modeled for varied surface flux conditions using the proposed unsaturated hydraulic conductivity and water storage functions. The impact of various levels of surface cracking on soil suction distributions is discussed. Suction distribution patterns are dependent on the initial soil surface suction. In particular, the results are dependent upon whether the initial matric suction is less than or greater than the air entry of the cracked clay. There is an extremely wide range of possible conditions that could be modeled but the parametric study results presented in this paper are limited to a series of selected crack widths and densities for an exfiltration case and an infiltration case.     

The conditions of formation of sapphire and zircon in the areas of alkali-basaltoid volcanism in Central Vietnam

Study of the chemical composition of clinopyroxene and garnet megacrysts from the Dak Nong sapphire deposit and model calculations have shown that megacrysts originated from the crystallization of alkali basaltoid magma in a deep-seated intermediate chamber at 14–15 kbar, which is close to the Moho depth (50 km) in this part of southeastern Asia. The chamber was a source of heat and CO₂ fluids for the generation of crustal syenitic melts producing sapphires and zircons. The formation conditions of sapphires and zircons are significantly different. The presence of jadeite inclusions in placer zircons points to high pressures during their crystallization, which is confirmed by the ubiquitous decrepitation of CO₂-rich melt inclusions. Sapphires crystallized from iron-rich syenitic melt in the shallower Earth’s crust horizons with the participation of CO₂ and carbonate–H₂O–CO₂ fluids. The subsequent eruptions of alkali basalts favored the transportation of garnet and pyroxene megacrysts as well as sapphire and zircon xenocrysts to the surface. It is shown that sapphire deposits can be produced only during multistage basaltic volcanism with deep-seated intermediate chambers in the regions with thick continental crust. The widespread megacryst mineral assemblage (clinopyroxene, garnet, sanidine, ilmenite) and the presence of placer zircon megacrysts can be used as indicators for sapphire prospecting.     

The sedimentary, magmatic and tectonic evolution of the southwestern South China Sea revealed by seismic stratigraphic analysis

The southwestern South China Sea represents an area of continental crust frozen immediately before the onset of seafloor spreading. Here we compile a grid of multichannel seismic reflection data to characterize the continent-ocean transition just prior to full break-up. We identify a major continental block separated from the shelf margin by a basin of hyperextended crust. Oligocene-Early Miocene rifting was followed by mild compression and inversion prior to 16 Ma, probably linked to collision between the Dangerous Grounds, a continental block to the east of the study area, and Borneo. The timing of inversion supports models of seafloor spreading continuing until around 16 Ma, rather than becoming inactive at 20 Ma. The off-shelf banks experienced uplift prior to 16 Ma in an area, which had previously been a depocenter. The off-shelf banks continued to extend after this time when the rest of the region is in a phase of thermal subsidence. Post-rift magmatism is seen in the form of scattered seamounts (~5–10 km across) within or on the edge of the deeper basins, and are dated as Late Miocene and Pliocene. They are not clearly linked to any phase of tectonic activity. Further inversion of the off-shelf banks occurred in the Pliocene resulting in a major unconformity despite the lack of brittle faulting of that age. We speculate that this is part of a wider pattern of scattered magmatism throughout the South China Sea at this time. Prograding clinoforms are seen to build out from the shelf edge in the south of the study area during the Pliocene, after 5.3 Ma, and then more towards the north and east during the Pleistocene. At the same time a trough south of the off-shelf banks is filled with >1.35 km of mostly Pleistocene sediment. While we expect the bulk of the sediment to come from the Mekong River, we also suggest additional sediment supply from Borneo and the Malay Peninsula via the Molengraaff River and its predecessors.     

Theoretical and numerical studies of the Vishniac instability in supernova remnants

Vishniac instability has been theoretically studied in supernova remnants where it is supposed to explain the fragmentation of the interstellar medium. However its role is not fully demonstrated in these astrophysical objects. Conditions and assumptions required for the instability growth are explained in detail in the present paper. In addition the HYDRO-MUSCL hydrodynamic code has been used to simulate this instability in order to compare the numerical growth rate with the Vishniac analytical solution.     

A synthetic bandwidth measurement system for acoustic remotesensing of geoacoustic properties of bottom sediments

A series of narrowband acoustic transmissions may be combined into an equivalent broadband signal - with a corresponding increase in resolution in the time domain. This procedure is called synthetic bandwidth processing. If ray methods are used in combination with inverse theory to determine the otacoustic properties of bottom sediments, the increased time resolution aids in the identification of arrivals which are closely spaced in time. This, in turn, permits improved estimates of the spatial distribution of properties and more accurate otacoustic models. Signal processing techniques and guidelines for synthesizing multioctave broadband pulse-like signals are presented. Also described is the instrumentation used to implement the technique in real ocean environments     

The volcanic-hosted massive sulphide deposits of the Iberian Pyrite Belt Review and preface to the Thematic Issue

The Iberian Pyrite Belt (IPB) has, over the past decade, been an area of renewed mining activity and scientific research that has resulted in a wealth of new data and new geological and metallogenic concepts that are succinctly presented in this Thematic Issue. The reason for this interest in the IPB, which forms part of the Hercynian orogenic belt, is that its Late Devonian to Middle Carboniferous rocks host a huge quantity of volcanic-hosted massive sulphide (VMS) mineralization (1700 Mt of sulphides, totalling 14.6 Mt Cu, 13.0 Mt Pb, 34.9 Mt Zn, 46100 t Ag and 880 t Au). The mineralization and its environment display a number of typical signatures that can be related to the mineralogy and zoning of the sulphide orebodies, to the lead isotopes of the mineralization, to the geochemical and mineralogical variations in the hydrothermal alteration halos surrounding the orebodies, to the geochemical characteristics of the bimodal volcanics hosting the VMS, to the complex structural evolution during the Hercynian orogeny, to the presence of palaeofaults and synsedimentary structures that acted as channels and discharge traps for the metalliferous fluids, and to the gossans developed over VMS. Discriminant geological criteria have been deduced for each domain which can be helpful in mineral exploration, complementing the more traditional prospecting techniques. Although the question of the IPB's geodynamic setting is still under debate, any interpretation must now take into account some incontrovertible constraints: for example, the geochemical characteristics of a large part of the basic lavas are comparable to those of mantle-derived basalts emplaced in extensional tectonic settings, and the associated acidic rocks were produced by melting of a basic crustal protolith at low- to medium-pressures and a steep geothermal gradient, thus, the sulphide-bearing volcano-sedimentary sequence differs strongly from recent arc-related series. It is considered here that the tectonic setting was extensional and epicontinental and that it developed during the Hercynian plate convergence, that culminated in thin-skinned deformation and accretion of the South Portuguese terrane to the Iberian Paleozoic continental block.

---

Resumen (translated by E. Pascual) Durante la década pasada, la Faja Pirítica Ibérica (FPI) ha sido un área de actividad minera e investigación cientifica renovadas, lo que ha conducido a la obtención de nuevos datos y conceptos geológicos y metalogénicos, que se exponen sucintamente en este Número Especial. La razón de este interés en la FPI, que forma parte del cinturón orogénico hercínico, es que sus rocas, cuyas edades abarcan desde el Devónico tardío al Carbonífero Medio, albergan una enorme cantidad de mineralizaciones de sulfuros masivos ligados a vulcanismo (1700 millones de toneladas de sulfuros, que totalizan 14,6 Mt de Cu, 13,0 Mt de Pb, 34,9 Mt de Zn, 46100 toneladas de Ag y 880 toneladas de Au). Las mineralizaciones y su entorno muestran signaturas que se pueden relacionar con la mineralogía y la zonación de las masas de sulfuros, con los isótopos de plomo de la mineralización, con las variaciones en los halos de alteración hidrotermal alrededor de las mineralizaciones, con los caracteres geoquímicos de las rocas volcánicas bimodales que albergan los sulfuros masivos, con la compleja evolución tectónica del conjunto durante la orogenia hercínica, con la existencia de paleofallas y estructuras sinsedimentarias que actuaron como canales y trampas de descarga para los fluidos metalíferos y los gossans que se desarrollaron sobre los sulfuros. Se han deducido criterios geológicos discriminantes para cada área de conocimiento, que pueden ser útiles para la exploración minera, complementando las técnicas más tradicionales de prospección. Aunque la cuestión del entorno geodinámico de la FPI todavía es materia de debate, cualquier interpretación tiene que tener ahora en cuenta algunas restricciones incontrovertibles: por ejemplo, los caracteres geoquímicos de una gran parte de las rocas básicas son comparables a los de basaltos derivados del manto y emplazados en entornos tectónicos extensionales, y las rocas ácidas asociadas se produjeron a partir de un protolito cortical básico, a presiones bajas o intermedias y asociadas a un abrupto gradiente térmico. Por consiguiente, la secuencia vulcanosedimentaria que contiene los sulfuros masivos difiere claramente de las series recientes relacionadas con entornos de arco. Consideramos aquí que el entorno tectónico fue extensional y epicontinental y que tuvo lugar durante la convergencia de placas hercínica, que culminó en deformación “thin-skinned” y acreción del terreno constituído por la Zona Sudportuguesa al bloque continental paleozoico ibérico.

---

---

Received: 4 April 1996 / Accepted: 10 April 1997     

Explicit expression and regularization of the harmonic reproducing kernels for the earth's ellipsoid

The reproducing kernel functions and the integral formulae involving the solutions of the fixed gravimetric boundary value problems for the earth's ellipsoid are investigated in a suitable polar coordinate system defined on the boundary ellipsoid. The infinite series expressions of the reproducing kernel functions are represented explicitly up to the square of the first excentricity of the boundary surface. The obtained results show that the reproducing kernel functions for the earth's ellipsoid are, in contrast to the case of the spherical boundary, inhomogeneous and anisotropic. Moreover, the anisotropy is stronger near the pole. Subsequently the kernel functions are regularized in order to overcome the weak singularity of the integrals.