A Model for the polar motion of the deformable Earth adequate for astrometric data

Refined analytical expressions for the frequencies corresponding to the Chandler motion of the pole and the diurnal rotation of the deformable Earth are derived. Numerical estimates of the period and amplitude of the polar oscillations are presented. The trajectory of the Chandler polar motion derived via numerical modeling is in qualitative and quantitative agreement with experimental data from the International Earth Rotation Service (IERS). An evolutionary model describing slow variations in the Earth’s rotation parameters under the action of the dissipative moments of the tidal gravitational forces on time scales considerably longer than the precession period of the Earth’s axis is constructed. The axis of the Earth’s figure tends to approach the angular momentum vector of the proper rotation.     

Groundwater abstraction and contamination in Lithuania as geoindicators of environmental change

Groundwater abstraction not only affects subsurface water, but also causes changes in the hydrosphere, lithosphere, and biosphere. There are certain symptoms of these environmental changes that might be designated as geoindicators. Geoindicators related to groundwater abstraction and contamination in Lithuania have been classified into three groups related to: (1) the decline/rise of groundwater levels, (2) the deterioration of groundwater quality, and (3) the intensification of water-rock interaction. On the basis of hydrogeological and hydrological investigations carried out over many years, different environmental impacts of intensive groundwater use are demonstrated. The most important among them are the formation of depression cones, intensification of groundwater pollution, saltwater encroachment, and the intensification of karst.     

Rare gas isotopes and parent trace elements in ultrabasic-alkaline-carbonatite complexes, Kola Peninsula: identification of lower mantle plume component

During the Devonian magmatism (370 Ma ago) ∼20 ultrabasic-alkaline-carbonatite complexes (UACC) were formed in the Kola Peninsula (north-east of the Baltic Shield). In order to understand mantle and crust sources and processes having set these complexes, rare gases were studied in ∼300 rocks and mineral separates from 9 UACC, and concentrations of parent Li, K, U, and Th were measured in ∼70 samples. ⁴He/³He ratios in He released by fusion vary from pure radiogenic values ∼10⁸ down to 6 × 10⁴. The cosmogenic and extraterrestrial sources as well as the radiogenic production are unable to account for the extremely high abundances of ³He, up to 4 × 10⁻⁹ cc/g, indicating a mantle-derived fluid in the Kola rocks. In some samples helium extracted by crushing shows quite low ⁴He/³He = 3 × 10⁴, well below the mean ratio in mid ocean ridge basalts (MORB), (8.9 ± 1.0) × 10⁴, indicating the contribution of ³He-rich plume component. Magnetites are principal carriers of this component. Trapped ³He is extracted from these minerals at high temperatures 1100°C to 1600°C which may correspond to decrepitation or annealing primary fluid inclusions, whereas radiogenic ⁴He is manly released at a temperature range of 500°C to 1200°C, probably corresponding to activation of ⁴He sites degraded by U, Th decay.Similar ⁴He/³He ratios were observed in Oligocene flood basalts from the Ethiopian plume. According to a paleo-plate-tectonic reconstruction, 450 Ma ago the Baltica (including the Kola Peninsula) continent drifted not far from the present-day site of that plume. It appears that both magmatic provinces could relate to one and the same deep-seated mantle source.The neon isotopic compositions confirm the occurrence of a plume component since, within a conventional ²⁰Ne/²²Ne versus ²¹Ne/²²Ne diagram, the regression line for Kola samples is indistinguishable from those typical of plumes, such as Loihi (Hawaii). ²⁰Ne/²²Ne ratios (up to 12.1) correlate well with ⁴⁰Ar/³⁶Ar ones, allowing to infer a source ⁴⁰Ar/³⁶Ar ratio of about 4000 for the mantle end-member, which is 10 times lower than that of the MORB source end-member. In (³He/²²Ne)_PRIM versus (⁴He/²¹Ne)_RAD plot the Kola samples are within array established for plume and MORB samples; almost constant production ratio of (⁴He/²¹Ne)_RAD ≅ 2 × 10⁷ is translated via this array into (³He/²²Ne)_PRIM ∼ 10. The latter value approaches the solar ratio implying the non-fractionated solar-like rare gas pattern in a plume source.The Kola UACC show systematic variations in the respective contributions of in situ-produced radiogenic isotopes and mantle-derived isotopes. Since these complexes were essentially plutonic, we propose that the depth of emplacement exerted a primary control on the retention of both trapped and radiogenic species, which is consistent with geological observations. The available data allow to infer the following sequence of processes for the emplacement and evolution of Kola Devonian UACC: 1) Ascent of the plume from the lower mantle to the subcontinental lithosphere; the plume triggered mantle metasomatism not later than ∼700 to 400 Ma ago. 2) Metasomatism of the lithosphere (beneath the central part of the Kola Peninsula), including enrichment in volatile (e.g., He, Ne) and in incompatible (e.g., U, Th) elements. 3) Multistage intrusions of parental melts, their degassing, and crystallisation differentiation ∼370 Ma ago. 4) Postcrystallisation migration of fluids, including loss of radiogenic and of trapped helium. Based on model compositions of the principle terrestrial reservoirs we estimate the contributions (by mass) of the plume material, the upper mantle material, and the atmosphere (air-saturated groundwater), into the source of parent melt at ∼2%, 97.95%, and ∼0.05%, respectively.     

Comparative Re-Os systematics of chondrites: Implications regarding early solar system processes

A suite of 47 carbonaceous, enstatite, and ordinary chondrites are examined for Re-Os isotopic systematics. There are significant differences in the ¹⁸⁷Re/¹⁸⁸Os and ¹⁸⁷Os/¹⁸⁸Os ratios of carbonaceous chondrites compared with ordinary and enstatite chondrites. The average ¹⁸⁷Re/¹⁸⁸Os for carbonaceous chondrites is 0.392 ± 0.015 (excluding the CK chondrite, Karoonda), compared with 0.422 ± 0.025 and 0.421 ± 0.013 for ordinary and enstatite chondrites (1σ standard deviations). These ratios, recast into elemental Re/Os ratios, are as follows: 0.0814 ± 0.0031, 0.0876 ± 0.0052 and 0.0874 ± 0.0027, respectively. Correspondingly, the ¹⁸⁷Os/¹⁸⁸Os ratios of carbonaceous chondrites average 0.1262 ± 0.0006 (excluding Karoonda), and ordinary and enstatite chondrites average 0.1283 ± 0.0017 and 0.1281 ± 0.0004, respectively (1σ standard deviations). The new results indicate that the Re/Os ratios of meteorites within each group are, in general, quite uniform. The minimal overlap between the isotopic compositions of ordinary and enstatite chondrites vs. carbonaceous chondrites indicates long-term differences in Re/Os for these materials, most likely reflecting chemical fractionation early in solar system history.A majority of the chondrites do not plot within analytical uncertainties of a 4.56-Ga reference isochron. Most of the deviations from the isochron are consistent with minor, relatively recent redistribution of Re and/or Os on a scale of millimeters to centimeters. Some instances of the redistribution may be attributed to terrestrial weathering; others are most likely the result of aqueous alteration or shock events on the parent body within the past 2 Ga.The ¹⁸⁷Os/¹⁸⁸Os ratio of Earth’s primitive upper mantle has been estimated to be 0.1296 ± 8. If this composition was set via addition of a late veneer of planetesimals after core formation, the composition suggests the veneer was dominated by materials that had Re/Os ratios most similar to ordinary and enstatite chondrites.     

A Comparison of Two-Phase Dynamic Upscaling Methods Based on Fluid Potentials

Although a number of methods for calculating dynamic pseudo-functions have been developed over the years, there is still a lack of understanding as to why a certain method will succeed in some cases but fail in others. In this paper, we describe the results of an assessment of several upscaling methods, namely the Kyte and Berry (KB) method, the Stone method, the Hewett and Archer (HA) method and the Transmissibility-Weighted (TW) method. We have analyzed the equations for deriving the methods and investigated the results of numerical simulations of gas displacing oil, in a variety of models to enable us to gain new insights into these, and related, upscaling methods. In particular, some novel observations on methods based on fluid potential are presented and the issue of using predicted fluid mobilities as a criterion of accuracy of an upscaling method is clarified.     

A decade of change in the Skidaway River estuary I. Hydrography and nutrients

The Skidaway River estuary is a tidally-dominated subtropical estuary in the southeastern USA surrounded by extensiveSpartina salt marshes. Weekly smapling at high and low tide began in 1986 for hydrography, nutrients, chlorophylla, particulate matter, and microbial and plankton biomass and composition; hydrographic and nutrient data during 1986–1996 are reported here. Salinity varied inversely with river discharge and exhibited variability at all time scales but with no long-term trend. Water temperature typically ranged over 25°C and was without apparent long-term frend. Seasonal cycles in concentrations of NO₃, NH₄, PO₄, Si(OH)₄, and DON were observed, with annual maxima generally occurring in late summer. Superimposed on seasonal cycles, all five nutrients exhibited steady increases in minimum, mean, and maximum concentrations; mean concentrations increased c. 50–150% during the decade. Nutrient concentrations were highly correlated with water temperature over the ten-year period, but weakly related to salinity and discharge. Nutrients were strongly correlated with one another, and the relative ratios among inorganic nutrients showed little long-term trend. Correlations among temperature and nutrient concentrations exhibited considerable inter-annual variability. Major spikes in organic and inorganic nutrient concentrations coincided with significant rainfall events; concentrations increased hyperbolically with rainfall. Although pristine compared to more heavily impacted waterways primarily outside the region, residential development and population density have been increasing rapidly during the past 15–20 years. Land use is apparently altering nutrient loading over the long-term (months-years), and superimposed on this are stochastic meteorological events that accelerate these changes over the short term (days-weeks).     

Numerical Homogenization of the Rigidity Tensor in Hooke's Law Using the Node-Based Finite Element Method

Combining a geological model with a geomechanical model, it generally turns out that the geomechanical model is built from units that are at least a 100 times larger in volume than the units of the geological model. To counter this mismatch in scales, the geological data model's heterogeneous fine-scale Young's moduli and Poisson's ratios have to be “upscaled” to one “equivalent homogeneous” coarse-scale rigidity. This coarse-scale rigidity relates the volume-averaged displacement, strain, stress, and energy to each other, in such a way that the equilibrium equation, Hooke's law, and the energy equation preserve their fine-scale form on the coarse scale. Under the simplifying assumption of spatial periodicity of the heterogeneous fine-scale rigidity, homogenization theory can be applied. However, even then the spatial variability is generally so complex that exact solutions cannot be found. Therefore, numerical approximation methods have to be applied. Here the node-based finite element method for the displacement as primary variable has been used. Three numerical examples showing the upper bound character of this finite element method are presented.     

Mathematical modeling of PCB bioaccumulation in Perna viridis

In the present work, we built a mathematical model of polychlorinated biphenyl (PCB) bioaccumulation in Perna viridis, namely, a one-compartment model with a time dependent incorporation rate R (μg g⁻¹ lipid per ppb water per day), with positive substrate cooperativity as the underlying physical mechanism. The temporal change of the PCB concentration Q (μg g⁻¹ lipid) in the soft tissues of the mussel depends on the competition of the input rate RW and the output rate kQ, where W is the concentration of PCB in water (ppb water) and k is the elimination rate (per day). From our experimental data, k=0.181±0.017 d⁻¹. The critical concentration in water W_c for positive substrate cooperativity was found to be 2.4 ppb. Below W_c, R is a constant. For a water concentration of 0.5 ppb Aroclor 1254, R=24.0±2.4 μg g⁻¹ lipid ppb⁻¹ d⁻¹. Above W_c, positive substrate cooperativity comes into effect and R becomes a function of time and dependent on the concentration Q in a form R=γQ/(Q+δ). This is the case for a water concentration of 5 ppb Aroclor 1254, where γ=15.1 μg g⁻¹ lipid ppb⁻¹ d¹ and δ≈200 μg g⁻¹ lipid. From this model, the uptake is exponentially increasing when the PCB concentration in the mussel is small compared to 200 μg g⁻¹ lipid, and hyperbolically increasing when the concentration is large compared to 200 μg g⁻¹ lipid, which are consistent with the experimental data. The model is useful for understanding the true processes taking place during the bioaccumulation and for risk assessment with higher confidence. Future experimental data which challenge the present model are anticipated and in fact desirable for improvement and perfection of the model.     

Determination of S-wave Site Response in Anchorage, Alaska in the 1–9 Hz Frequency Band

—?Site response was estimated at 19 sites in the Anchorage basin in south-central Alaska, using 15 local earthquakes recorded with good signal-to-noise ratio by a temporary weak motion network. The receiver-function-type horizontal-to-vertical spectral ratios (HVSR) were computed at 1–9 Hz frequency band and the resulting HVSR contour maps at 1, 5 and 9 Hz are presented here. The spatial site response distribution shows considerable variation from the foothills of the Chugach Mountains in the east to the western part of Anchorage. The site response increases by a factor of 3 and 3.5 at 1 and 5 Hz, respectively, from the area of older glacial deposits in the eastern part of the city to the area occupied by the Bootlegger Cove formation, particularly in the section adjoining Knik Arm. At 9 Hz, the variation of HVSR from the east to the west is smaller, approximately by a factor of 2. Moreover, the trend of the HVSR variation at 1 and 5 Hz shows good correlation with that of the soil class obtained from surface measurements of S-wave velocity in the 0–30?m depth range and available results on ground failure susceptibility of Anchorage.