Identifying Source Soils in Contemporary Estuarine Sediments: A New Compound-Specific Isotope Method

A new method is proposed for the identification and apportionment of contemporary source soils contributing to estuarine sediments. The method uses compound-specific isotopic analysis of naturally occurring biomarkers (fatty acids) derived from plants to link source soils to land use within a single catchment. For identification and apportionment of source soils in the estuarine samples, the method uses the isotopic mixing model, IsoSource. The feasible proportions obtained from IsoSource are then scaled to allow for the percent organic carbon in the source soils. With this approach, the estimation of each source soil contribution to a location in the estuary is independent of any degradation of the biomarkers through microbial or biogeochemical processes. Identification relies on the evaluation of the sediment sample relative to a “library” of reference source soils from different land use within the catchment. Selection of potential sources is geographically constrained by the requirement for a natural linkage between each source soil and the sediment site sampled. A case study, using this method, mapped the distribution of three main land use source soils (pasture, native forest, and pine forest) across the river delta in a small estuary fringed with mangroves. Rather than being uniformly distributed, the results indicated that the source soil contributions varied markedly across the delta, raising concerns about the validity of taking single cores to characterize the sediments of an estuary. Coupling the source apportionment results with land use data indicated that the mean percent contribution of pine forest soil in the river delta sediments was almost three times greater than the percent land use area of pine forest in the catchment. Furthermore, isotopic signatures indicated that most of the pine forest soil came from the much smaller areas exposed to erosion by clear cut harvesting and that the soil contribution from recently harvested areas of pine forest could be as much as 20 times greater than that land use area in this catchment. This is the first method that can identify and apportion, by land use on a catchment scale, the sources of soil contributing to the sediment at a location of an estuary. The results are given as a “best estimate”, within definable limits, of the proportional contribution of each potential source soil. Information obtained using this method will allow development of management strategies to alter land use practices to reduce the sediment load to rivers, and thus, the impact on the aquatic ecosystem downstream in estuaries.     

Light Requirements for Growth and Survival of Eelgrass (<Emphasis Type="Italic">Zostera marina</Emphasis> L.) in Pacific Northwest (USA) Estuaries

We developed light requirements for eelgrass in the Pacific Northwest, USA, to evaluate the effects of short- and long-term reductions in irradiance reaching eelgrass, especially related to turbidity and overwater structures. Photosynthesis-irradiance experiments and depth distribution field studies indicated that eelgrass productivity was maximum at a photosynthetic photon flux density (PPFD) of about 350–550 μmol quanta m⁻² s⁻¹. Winter plants had approximately threefold greater net apparent primary productivity rate at the same irradiance as summer plants. Growth studies using artificial shading as well as field monitoring of light and eelgrass growth indicated that long-term survival required at least 3 mol quanta m⁻² day⁻¹ on average during spring and summer (i.e., May-September), and that growth was saturated above about 7 mol quanta m⁻² day⁻¹. We conclude that non-light-limited growth of eelgrass in the Pacific Northwest requires an average of at least 7 mol quanta m⁻² day⁻¹ during spring and summer and that long-term survival requires a minimum average of 3 mol quanta m⁻² day⁻¹.     

Mineral assemblages of titanium-zirconium sands at the Central deposit,the east European platform

The study of mineral assemblages at the Central deposit allowed us to substantially refine the evolutionary model of the deposit and reveal the two main factors that control the variability of its mineralic space: (1) heterogeneity of the ore layer, consisting of a sublittoral bottom placer (the lower part) and a subaerial dune complex partly reworked in the course of a new transgression (the upper part), and (2) postore epigenetic alteration of the ore constituent of sands, which affected the quality of ore concentrates. The results obtained will be used in geological and technological mapping and development of the production program.     

Solar flares associated with singularities in the differential characteristics of magnetic fields

We analyze the flare events of October 28 and May 27, 2003 to examine the possible connection of solar flares with singularities in the differential characteristics of solar magnetic fields. The relation between flares and the behavior of F indicators (which describe structural singularities of the magnetic fields) introduced earlier is analyzed. It is shown that the generation of flares can be associated with self-intersections (or reconnection) of the F = 0 curves, where magnetic reconnection can occur. Consequently, the current sheets generated in such regions can be sources of energy release.     

Temporal variations in the position of the heliospheric equator

It is shown that the centroid of the heliospheric equator undergoes quasi-periodic oscillations. During the minimum of the 11-year cycle, the centroid shifts southwards (the so-called bashful-ballerina effect). The direction of the shift reverses during the solar maximum. The solar quadrupole is responsible for this effect. The shift is compared with the tilt of the heliospheric current sheet.     

Rotational and magnetic corrections to the g-mode frequencies in a one-dimensional MHD model

The central magnetic field and rotation of the solar radiative zone are responsible for corrections to the g-mode frequencies. Magnetogravitational spectra are calculated analytically in a simple one-dimensional MHD model that goes beyond the WKB approximation and avoid any cusp resonances that trap the wave within the radiative zone in the presence of a weak magnetic background. The calculations are compared with spacecraft observations of the 1% frequency shifts for candidate g-modes found in the SOHO GOLF experiment. The magnetic correction is the main contribution for a strong magnetic field satisfying the approximation used. It is shown that a constant magnetic field of 700 kG in the radiative zone provides the required frequency shift for the n = ?10 g-mode. The rotational correction, which is due to the Coriolis force in the one-dimensional model used, is much less than a percent (α_Ω ≤ 0.003).     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Geochemical Eccentricity of Ground Water Allied to Weathering of Basalts from the Deccan Volcanic Province,India: Insinuation on CO<Subscript>2</Subscript> Consumption

Analyses of 72 samples from Upper Panjhara basin in the northern part of Deccan Plateau, India, indicate that geochemical incongruity of groundwater is largely a function of mineral composition of the basaltic lithology. Higher proportion of alkaline earth elements to total cations and HCO₃>Cl + SO₄ reflect weathering of primary silicates as chief source of ions. Inputs of Cl, SO₄, and NO₃ are related to rainfall and localized anthropogenic factors. Groundwater from recharge area representing Ca + Mg–HCO₃ type progressively evolves to Ca + Na–HCO₃ and Na–Ca–HCO₃ class along flow direction replicates the role of cation exchange and precipitation processes. While the post-monsoon chemistry is controlled by silicate mineral dissolution + cation exchange reactions, pre-monsoon variability is attributable chiefly to precipitation reactions + anthropogenic factors. Positive correlations between Mg vs HCO₃ and Ca + Mg vs HCO₃ supports selective dissolution of olivine and pyroxene as dominant process in post-monsoon followed by dissolution of plagioclase feldspar and secondary carbonates. The pre-monsoon data however, points toward the dissolution of plagioclase and precipitation of CaCO₃ supported by improved correlation coefficients between Na + Ca vs HCO₃ and negative correlation of Ca vs HCO₃, respectively. It is proposed that the eccentricity in the composition of groundwater from the Panjhara basin is a function of selective dissolution of olivine > pyroxene followed by plagioclase feldspar. The data suggest siallitization (L < R and R _k) as dominant mechanism of chemical weathering of basalts, stimulating monosiallitic (kaolinite) and bisiallitic (montmorillonite) products. The chemical denudation rates for Panjhara basin worked out separately for the ground and surface water component range from 6.98 to 36.65 tons/km²/yr, respectively. The values of the CO₂ consumption rates range between 0.18 × 10⁶ mol//km²/yr (groundwater) and 0.9 × 10⁶ mol/km²/yr (surface water), which indicates that the groundwater forms a considerable fraction of CO₂ consumption, an inference, that is, not taken into contemplation in most of the studies.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.