Stabilized single-point 4-node quadrilateral element for dynamic analysis of fluid saturated porous media

An accurate and efficient low-order quadrilateral mixed u?Cp element suitable for dynamic analysis of fluid saturated porous media is presented. The element uses physical hourglass stabilization to facilitate single-point integration for the solid phase, and non-residual stabilization of the fluid phase to circumvent instability in the incompressible-impermeable limit due to the use of equal-order interpolation for the displacement and pressure fields. Element behavior is verified and demonstrated through several numerical examples.     

Hydrochemical considerations and heavy metal variability in the middle Paraná River

The Paraná River is one of the largest drainage systems in the Americas. Its hydrology is characterized by an active teleconnection with the ENSO, and by a significant discharge increase trend, evident since the mid-1970s. An Eh–pH data set collected in the Paraná’s middle stretch suggests that large flood events, such as the one triggered by the 1982–1983 ENSO, are discernible in the plot, probably due to the influx of water draining flood plain water bodies. The total (particulate + dissolved) concentration of a set of heavy metals (Cr, Mn, Ni, Cu, Zn, As, Cd, and Pb) was determined in a downriver survey of the middle stretch. With the exception of Cu, Cd, and Pb, the metals exhibit a significantly increasing concentration trend towards the river mouth. The slopes of the regression lines imply that Zn and Ni, on one hand, and Mn and Cr, on the other would have common controlling sources. Another set of analyses were performed during the 1982–1983 flooding event; besides an increased variability observable during the flood arrival, most elements, with the only exception of Pb, did not show a variability coherent with the discharge series.     

Mercury in the Waters of the Jundiaí River, SP, Brazil: The Role of Dissolved Organic Matter

Many developing countries have regions of high demographic density, where untreated residuary waters from different sources are often discharged into rivers, streams and other water bodies. This paper discusses the reducing action of organic matter of anthropic origin on the mercury redox cycle in the Jundiaí River impacted by discharged wastes, and on the Piraí River, a non-impacted water body. The total mercury concentrations in these locations vary from 1.7 to 32?ng?L^?1 in the former and from 0.6 to 10.6?ng?L^?1 in the latter. Dissolved organic carbon concentrations of up to 68.3 and 6.5?mg?L^?1 were observed, confirming the higher impact on the Jundiaí River. It was found that an inverse correlation between the concentration of dissolved organic carbon and total mercury was stronger in the Jundiaí River, given that it receives higher organic loads, suggesting that organic matter exerts a reducing action on mercury, which is released as gas into the atmosphere. This correlation was not observed in the Piraí River, where the organic matter of natural origin is probably not sufficiently labile to act intensely upon the Hg redox cycle, favoring the metal transport.     

Obsidian Subsource Identification in the Sierra de Pachuca and Otumba Volcanic Regions,Central Mexico,by ICP‐MS and DBSCAN Statistical Analysis

Otumba and Sierra de Pachuca obsidian deposits in Central Mexico have been important sources of raw material since pre‐Hispanic times. Numerous archaeological investigations have suggested that the economical and political expansion of major Mesoamerican societies were linked to the control of obsidian sources and distribution of quarried material. Sierra de Pachuca contains several obsidian flows and numerous quarries throughout the region that were preferentially exploited by different cultures. The Otumba Volcanic Complex has four important obsidian domes, but three of them have not been studied in detail. A geochemical characterization of subsources from the Sierra de Pachuca and Otumba Volcanic Complex is an important step toward future sourcing of obsidian artifacts that would help provide insight into spheres of influence and trade by past cultures in Central Mexico. Having this purpose in mind, inductively coupled plasma mass spectrometry (ICP‐MS) was used to analyze obsidian samples collected from five separated locations at Sierra de Pachuca and four at Otumba, followed by statistical analysis (density‐based spatial clustering of applications with noise, DBSCAN). We were able to distinguish three chemically distinctive subsources in Sierra de Pachuca and three in Otumba. This study illustrates the importance of accurate characterization of obsidian raw material when attempting to define subsource usage.     

Vertical coupling and dynamical source for the intraseasonal variability in the deep Kuroshio Extension

Ocean Dynamics - In the power spectrum, the upper and deep parts of the Kuroshio Extension have distinctly different peaks. The former peaks around 200 days, while the latter is mainly at the...     

PETROLOGY AND CLASSIFICATION OF THE GARRAF,SPAIN CHONDRITE

Microscopic and electron microprobe studies indicate that the Garraf meteorite is a highly-recrystallized chondrite of petrologic type 6. Olivine (Fa24.7; PMD 1.1) and low-Ca pyroxene (Fs20.9; PMD 1.1) compositions indicate that it belongs to the L-group. Based on contents of noble gases, pervasive fracturing of silicates, common undulose extinction of olivine and plagioclase, and the lack of melt pockets and maskelynite, we place Garraf into shock facies b. We conclude that Garraf is a highly recrystallized L6b chondrite that, after recrystallization, was cataclased and comminuted by shock.     

On the non-gravitational forces in the 1986 return of P/Halley

In this paper we study the time evolution of the non-gravitational forces (NGF) along with the observational bias, during the 1986 return of comet Halley. First, evidence is presented which shows that a NGF model with constant coefficients along the radial and transverse directions (A1 andA2), and a constant coefficient radial-only observation bias model fails to represent the observations within the conservative optical measurement noise of one second of arc. Second, we present a stochastic approach to the problem, where coefficients for the radial, transverse and normal components of both the NGF and the observation bias, are allowed to vary in time as three statistically independent Gauss-Markov processes. Finally, an orbit is estimated with this model, to fit observations made during the last apparition, employing a stochastic optimal smoother.Results are given in terms of time histories of the coefficient estimates along with their smoother computed uncertainties. A plot of the observation residuals is also included, showing a uniform and unbiased behavior. The analysis of the results confirms some of the assumptions often made when modelling cometary motion, but questions others. In particular, the normal bias coefficient shows an unexpected pre-perihelion peak (–1200±250 km) which proves that a radial-only observation bias model may lead to biased orbital estimates and unrealistic computed uncertainties.     

Estimating the nucleus density of Comet 19P/Borrelly

The nucleus bulk density of Comet 19P/Borrelly has been estimated by modeling the sublimation-induced non-gravitational force acting upon the orbital motion, thereby reproducing the empirical perihelion advance (i.e., the shortening of the orbital period). The nucleus has been modeled as a prolate ellipsoid, covered by various surface activity maps which reproduce the observed water production rate. The theoretical water production rate of active areas has been obtained by applying a sophisticated thermophysical model. This model takes into account net sublimation of ice and thermal reradiation from the surface, solid state conductivity, sub-surface sublimation and recondensation, mass and heat transport by diffusing gas, layer absorption of solar energy, a full treatment of local time-dependent illumination conditions, and a detailed consideration of nucleus/coma interaction mechanisms. The outgassing properties of the modeled nucleus are physically consistent with the gas kinetic structure of the innermost coma since the molecular backflux and surface gas density required in the thermophysical model (as functions of the nucleus surface temperature and the sub-surface temperature profile) have been obtained from Direct Simulation Monte Carlo modeling of inelastic intermolecular collisions in the cometary Knudsen layer. The calculation of local normal forces acting on the nucleus due to outgassing has been made within the same framework—recoil and/or impact momentum transfer to the nucleus caused by sublimating molecules and by recondensing and/or scattered coma molecules is therefore evaluated in accordance with local nucleus/coma conditions. According to this model, the density is found to be 100-300 kg, depending on the applied spin axis orientation and surface activity map. This range can be narrowed down to 180-300 kg by also requiring that the empirical changes (per orbital revolution) of the argument of perihelion and the longitude of the ascending node are reproduced.     

Geometry of half-grabens containing a mid-level viscous décollement

In this work, we explore by means of analogue models how different basin-bounding fault geometries and thickness of a viscous layer within the otherwise brittle pre-rift sequence influence the deformation and sedimentary patterns of basins related to extension. The experimental device consists of a rigid wooden basement in the footwall to simulate a listric fault. The hangingwall consists of a sequence of pre-rift deposits, including the shallow interlayered viscous layer, and a syn-rift sequence deposited at constant intervals during extension. Two different geometries exist of listric normal faults, dip at 30 and 60° at surface. This imposes different geometries in the hangingwall anticlines and their associated sedimentary basins. A strong contrast exists between models with and without a viscous layer. With a viscous décollement, areas near the main basement fault show a wide normal drag and the hangingwall basin is gently synclinal, with dips in the fault side progressively shallowing upwards. A secondary roll-over structure appears in some of the models. Other structures are: (1) reverse faults dipping steeply towards the main fault, (2) antithetic faults in the footwall, appearing only in models with the 30° dipping fault and silicone-level thicknesses of 1 and 1.5 cm and (3) listric normal faults linked to the termination of the detachment level opposite to the main fault, with significant thickness changes in the syn-tectonic units. The experiments demonstrate the importance of detachment level in conditioning the geometry of extensional sedimentary basins and the possibility of syncline basin geometries associated with a main basement fault. Comparison with several basins with half-graben geometries containing a mid-level décollement supports the experimental results and constrains their interpretation.