Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 John Wiley & Sons, Ltd.     

Chemical abundances determined from meteor spectra: I. Ratios of the main chemical elements

Abstract— Relative chemical abundances of 13 meteoroids were determined by averaging the composition of the radiating gas along the fireball path that originated during their penetration into the Earth's atmosphere. Mg, Fe, Ni, Cr, Mn, and Co abundances, relative to Si, are similar to those reported for CI and CM carbonaceous chondrites and interplanetary dust particles. In contrast, relative abundances of Ca and Ti in meteor spectra indicate that these elements suffer incomplete evaporation processes. The chemical composition of all meteoroids studied in this work differs from that of 1P/Halley dust.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Peak metamorphic temperatures in type 6 ordinary chondrites: An evaluation of pyroxene and plagioclase geothermometry

Abstract— Quantifying the peak temperatures achieved during metamorphism is critical for understanding the thermal histories of ordinary chondrite parent bodies. Various geothermometers have been used to estimate equilibration temperatures for chondrites of the highest metamorphic grade (type 6), but results are inconsistent and span hundreds of degrees. Because different geothermometers and calibration models were used with different meteorites, it is unclear whether variations in peak temperatures represent actual ranges of metamorphic conditions within type 6 chondrites or differences in model calibrations. We addressed this problem by performing twopyroxene geothermometry, using QUILF95, on the same type 6 chondrites for which peak temperatures were estimated using the plagioclase geothermometer (Nakamuta and Motomura 1999). We also calculated temperatures for published pyroxene analyses from other type 6 H, L, and LL chondrites to determine the most representative peak metamorphic temperatures for ordinary chondrites. Pyroxenes record a narrow, overlapping range of temperatures in H6 (865–926 °C), L6 (812–934 °C), and LL6 (874–945 °C) chondrites. Plagioclase temperature estimates are 96–179 °C lower than pyroxenes in the same type 6 meteorites. Plagioclase estimates may not reflect peak metamorphic temperatures because chondrule glass probably recrystallized to plagioclase prior to reaching the metamorphic peak. The average temperature for H, L, and LL chondrites (～900 °C), which agrees with previously published oxygen isotope geothermometry, is at least 50 °C lower than the peak temperatures used in current asteroid thermal evolution models. This difference may require minor adjustments to thermal model calculations.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

Palaeoenvironmental significance of Late Permian palaeosols in the South‐Eastern Iberian Ranges,Spain

The Late Permian (Wuchiapingian) Alcotas Formation in the SE Iberian Ranges consists of one red alluvial succession where abundant soil profiles developed. Detailed petrographical and sedimentological studies in seven sections of the Alcotas Formation allow six different types of palaeosols, with distinctive characteristics and different palaeogeographical distribution, to be distinguished throughout the South‐eastern Iberian Basin. These characteristics are, in turn, related to topographic, climatic and tectonic controls. The vertical distribution of the palaeosols is used to differentiate the formation in three parts from bottom to top showing both drastic and gradual vertical upwards palaeoenvironmental changes in the sections. Reconstruction of palaeoenvironmental conditions based on palaeosols provides evidence for understanding the events that occurred during the Late Permian, some few millions of years before the well‐known Permian‐Triassic global crisis.     

Methane‐related microbial gypsum calcitization in stromatolites of a marine evaporative setting (Münder Formation,Upper Jurassic,Hils Syncline,north Germany)

Fossil stromatolites may reveal information about their hydrochemical palaeoenvironment, provided that assignment to a specific microbial community and a corresponding biogeochemical mechanism of formation can be made. Tithonian stromatolites of the Münder Formation at Thüste, north Germany, have traditionally been considered as formed by intertidal cyanobacterial communities. However, thin sections of the stromatolites show elongated angular traces of former gypsum crystals in a dense arrangement, but no algal or cyanobacterial filament traces. Moreover, high Fe²⁺ and Mn²⁺ contents, oxygen‐isotope and sulphur‐isotope ratios of carbonate‐bound sulphates, and sulphurized hydrocarbon biomarkers of the stromatolitic carbonate indicate that CaCO₃ precipitation occurred near the oxic–anoxic interface as a result of intensive bacterial sulphur cycling rather than photosynthetic activity. Furthermore, anaerobic oxidation of methane by Archaea may have driven CaCO₃ precipitation in deeper parts of the biofilm community, as reflected by high concentrations of squalane with a strongly negative δ¹³C in conjunction with evaporite pseudomorphs showing extremely low δ¹³C_Carb ratios. Consequently, the Thüste stromatolites are now interpreted as having initially formed by gypsum impregnation of biofilms. Subsequently, early Mg‐calcitic calcitization within the biofilms occurred because of combined bacterial iron, manganese and sulphate reduction, with an increasing contribution of anaerobic oxidation of methane with depth. This model plausibly explains the prominent preservation of signals derived from oxygen‐independent metabolic pathways, whereas virtually no geochemical record exists for an aerobic community that may, nevertheless, have prevailed at the stromatolite surface. Photic‐zone stromatolites with a prominent signal of anaerobic oxidation of methane may be common in, and indicative of, oxygen‐depleted sulphate‐bearing environments with high rates of methane production, conditions that possibly were fulfilled at the Archaean to Proterozoic transition.     

A GPR study of subsidence-creep-fault processes in Morelia, Michoacán, Mexico

In 1983, inhabitants of the City of Morelia, Michoacán, Mexico, began to observe a series of differential settlements causing damages to constructions along linear trends parallel to a system of regional faults. The same phenomenon occurs in others cities of the Mexican Volcanic Belt (MVB), such as Celaya, Aguascalientes, and Querétaro, and is linked to a structurally controlled subsidence, caused by groundwater withdrawal, and the presence of geological faults. We define this subsidence type as Subsidence-Creep-Fault Processes (SCFP), based on the necessary elements for their generation, and we studied them through geophysical and geotechnical techniques. In Morelia, the geophysical investigations have been carried out using ground-penetrating radar (GPR). GPR profiles, perpendicular to the axis of the surface fault generated by the SCFP were carried out. The common-offset single-fold profiling was used, with a central frequency of 50 MHz. In all cases it has been possible to visualize a fault plane dividing two blocks, the presence of synthetic and antithetic faults, influence zones from 20 m to 40 m, and a maximum “net throw” of 4 m. Exploration trenches followed the same direction of the profiles obtained with GPR (perpendicular to the axis of the surface fault). These trenches exposed a fault plane dividing two blocks with different lithology, generating a maximum “net throw” of 4.40 m; as well they help in the determination of influence zones that varied from 14 m to 40 m.     

Sedimentary petrology and geochemistry of siliciclastic rocks from the upper Jurassic Tordillo Formation (Neuquén Basin, western Argentina): Implications for provenance and tectonic setting

The Upper Jurassic Tordillo Formation is exposed along the western edge of the Neuquén Basin (west central Argentina) and consists of fluvial strata deposited under arid/semiarid conditions. The pebble composition of conglomerates, mineralogical composition of sandstones and pelitic rocks, and major- and trace-element geochemistry of sandstones, mudstones, and primary pyroclastic deposits are evaluated to determine the provenance and tectonic setting of the sedimentary basin. Conglomerates and sandstones derived almost exclusively from volcanic sources. The stratigraphic sections to the south show a clast population of conglomerates dominated by silicic volcanic fragments and a predominance of feldspathic litharenites. This framework composition records erosion of Triassic–Jurassic synrift volcaniclastic rocks and basement rocks from the Huincul arch, which was exhumed as a result of Late Jurassic inversion. In the northwestern part of the study area, conglomerates show a large proportion of mafic and acidic volcanic rock fragments, and sandstones are characterised by a high content of mafic volcanic rock fragments and plagioclase. These data suggest that the source of the sandstones and conglomerates was primarily the Andean magmatic arc, located west of the Neuquén Basin. The clay mineral assemblage is interpreted as the result of a complex set of factors, including source rock, climate, transport, and diagenesis. Postdepositional processes produced significant variations in the original compositions, especially the fine-grained deposits. The Tordillo sediments are characterised by moderate SiO₂ contents, variable abundances of K₂O and Na₂O, and a relatively high proportion of ferromagnesian elements. The degree of chemical weathering in the source area, expressed as the chemical index of alteration, is low to moderate. The major element geochemistry and Th/Sc, K/Rb, Co/Th, La/Sc, and Cr/Th values point to a significant input of detrital volcanic material of calcalkaline felsic and intermediate composition. However, major element geochemistry is not useful for interpreting the tectonic setting. Discrimination plots based on immobile trace elements, such as Ti, Zr, La, Sc, and Th, show that most data lie in the active continental margin field. Geochemical information is not sufficiently sensitive to differentiate the two different source areas recognized by petrographic and modal analyses of conglomerates and sandstones.