A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mathematical modeling of thermomechanical erosion beneath Proterozoic komatiitic basaltic sinuous rilles in the Cape Smith Belt,New Québec,Canada

Virtually all of the economic Ni–Cu–(platinum group element (PGE)) mineralization in the central part of the Cape Smith Belt of New Québec is hosted by thick olivine cumulate units in the Katinniq Member of the Raglan Formation at the base of the 1.9 Ga Chukotat Group. These units transgress underlying gabbros and pelitic metasediments, forming 50–200-m deep and 300–1,000-m wide V-shaped embayments and have been interpreted on the basis of surface geology, deep diamond core drilling, and magnetic inversion models to represent the remnants of one or more large, long (at least 20 km, possibly ≥50 km), sinuous, komatiitic basalt lava channels that formed by thermomechanical erosion of their substrates. We have used a mathematical model to test these hypotheses regarding komatiitic lava emplacement and erosion by lava. Our modeling predicts that an initially 10-m thick komatiitic basalt flow should have flowed turbulently near the vent and should have thermomechanically eroded unconsolidated pelitic sediment during emplacement to reach the observed degree of contamination of ≤10% at distances of ~30–60 km downstream from the source. Furthermore, our models predict that, at these distances downstream, a fully inflated 100-m thick komatiitic basalt flow would have had thermal erosion rates over consolidated gabbroic substrate of ~0.7–1.5 m/day, requiring ~70–140 days to incise a 100-m deep channel, depending on the initial temperature of the lava, the paleoslope, and the initial temperature and solidus temperature of the gabbro. These erosion rates would have been associated with volumetric flow rates of >10⁵–10⁶ m³/s and eruption volumes of >10³–10⁴ km³. Although these flow rates are orders of magnitude larger than those of most modern terrestrial basaltic flows, they are of the same order as those estimated for the largest terrestrial flood basalt flows and with those inferred for some of the largest extraterrestrial flows. Our predicted flow volumes are also of the same order as those of the largest terrestrial flood basalt units, consistent with the great thickness and widespread distribution of the Chukotat Group. Our modeling of thermomechanical erosion of gabbro by komatiitic basalt results in negligible contamination (<1%), and geochemical studies show that the spatially and petrogenetically related Chukotat basalts are uncontaminated, suggesting that the observed enrichments in U–Th–Light Rare Earth Elements (REE) > Middle REE–Heavy REE > Nb–Ta–Ti represent contamination by underlying Povungnituk semipelites. This result is consistent with present models for the genesis of the Ni–Cu–(PGE) mineralization in the Raglan Formation that involve thermomechanical erosion of unconsolidated, sulfidic semipelitic sediments, and decoupling of the miscible silicate and immiscible sulfide components.     

Optical characterization of CDOM in a marsh-influenced environment in the Changjiang (Yangtze River) Estuary

To elucidate the influence of tidal marshes on the amount and chemical structure of chromophoric dissolved organic matter (CDOM) in adjacent coastal waters, surface water samples were collected along three cross-marsh transects in tidal marshes within the Changjiang Estuary. In addition, three sediment cores were collected from the marshes and incubated for 10 h in the laboratory after adding overlying tidal waters. Nutrients, DOC, absorption, fluorescence excitation–emission matrix (EEM) spectroscopy, and other parameters were measured in these samples, which clearly illustrated releases by the tidal marshes of DOC, and CDOM absorbent and fluorescent molecules. To characterize the chemical nature of the marsh-derived CDOM, variations in a suite of optical indices (including molar absorptivity, absorption spectral slope, fluorescence quantum yield, and fluorescence ratios between certain excitation/emission locations) in these samples were plotted against the offshore distance of the stations as well as the incubation time. From the literature, it has been suggested that these variations are closely associated with molecular weight, aromaticity, and humic/fulvic or allochthonous/autochthonous ratios in the CDOM mixtures. This study showed that the CDOM fluorophores changed to more humic-like and were less protein-like from sea to land across the marsh gradients. Furthermore, subtractions between EEM spectra before and after laboratory incubations were applied to illustrate in greater detail the properties of the pure fluorescent compounds exchanged at the sediment–overlying water interface.     

New particle formation by ion-induced nucleation during dissipation stage of thunderstorm

A case of new particle formation observed during dissipation stage of a thunderstorm at a tropical station, Pune, India on 3 June 2008 is reported. The flash rate and rainfall intensity increased as high as 110 flashes per 5 minutes and 150 mm hour^− 1 respectively during the active stage of thunderstorm, and then gradually decreased during the dissipation stage. The number concentration of particles in the size range of 10–100 nm sharply increased from 350 particles cm^− 3 to ∼8000 particles cm^− 3 during the dissipation stage of a thunderstorm and grew to larger diameter subsequently. Observations suggest that the atmospheric conditions such as (i) reduced background aerosol concentration after heavy rain, (ii) high humidity condition, and (iii) increased ion concentration during the dissipation stage by corona discharges, favoured generation of new particles by ion-induced nucleation (IIN). Observations also suggest that generation of unipolar ions by corona discharges may be more favourable for IIN and subsequent growth of the particles.     

Atmospheric ultrafine aerosol number concentration and its correlation with vehicular flow at two sites in the western Himalayan region: Kullu-Manali,India

The concentration of ultrafine aerosol particles of aitken and nucleation mode having size in the range of 1–20 nm was monitored with water-based Condensation Particle Counter. The monitoring was carried out from midnight-to-midnight in every alternate day on a fortnightly basis to represent summer, monsoon and winter (autumn) seasons of 2008 at Mohal (1154 m amsl) and Kothi (2530 m amsl) in Kullu-Manali area of the northwestern Himalayan region of India. The results indicate that diurnal pattern has faint bimodal structure with two peaks, one in morning and the other in evening at both the sites but it is not as distinct as found in plains. There is rather a constant particle density pattern of large magnitude consistent with vehicular movement from morning till evening. The monthly 24 h average particle density gradually picks up from January, increases rapidly in summer months and then decreases in monsoon season at Mohal but at Kothi it keeps on rising from April to October with a slight more increase in September. The particle density is more in summer than in monsoon season at Mohal, a trend opposite to plains. It may be due to the development of warm thermal layer on valley floor while a cold layer develops along snowy hilltops in winter leading to convection of fine particle up the slopes of valley during daytime. At Kothi, the trend is same as it is in continental plains but opposite to Mohal. The relatively more value of particle density in September and October at both the sites may be due to month long International Kullu Dussehra fair in the valley. The vehicular survey conducted agrees well with entire study period averaged diurnal variations and monthly 24 h averaged value of fine particle density. The average value of ultrafine particle density at each hour of a day for entire study period is 20369 ± 1230 Ncm^− 3 and 14389 ± 1464 Ncm^− 3 at Mohal and Kothi sites, respectively. The comparison with earlier results shows a significant increase indicating impact of vehicular onslaught on pure air of this hilly region.     

Study on capillary rise from shallow groundwater and critical water table depth of a saline-sodic soil in western Songnen plain of China

The objective of this study was to develop an empirical equation for estimating the capillary rise in the saline-sodic soil area of Songnen Plain in China based on the Averianov formula. The capillary rise was observed under five controlled groundwater levels by lysimeters. Field experiment results indicated that capillary rise had close relationship with the groundwater table depth, soil moisture of 10–40 cm soil layer and leaf area index. These factors have been taken into account to develop the empirical equation for capillary rise simulation. The model parameters for Songnen Plain were derived by Levenberg–Marquardt and global optimization calculating method. The modeled capillary rise has a good agreement with the observed data (r ² = 0.875). With the simulation model, the critical water table depth was identified as 2.5 m, indicating that soil secondary salinization will not occur when the water table depth is deeper than 2.5 m. Therefore, in the irrigation areas, groundwater table depth should be controlled to be higher than 2.5 m to prevent the occurrence of soil secondary salinization. The results from this research will provide useful information for the water sources management and soil secondary salinization control in Songnen Plain of China, one of the most serious saline-affected areas in the world.     

Discovery of Ordovician–Silurian metamorphic monazite in garnet metapelites of the Alpine External Aiguilles Rouges Massif

The pre-Mesozoic, mainly Variscan metamorphic basement of the Col de Bérard area (Aiguilles Rouges Massif, External domain) consists of paragneisses and micaschists together with various orthogneisses and metabasites. Monazite in metapelites was analysed by the electron microprobe (EMPA-CHIME) age dating method. The monazites in garnet micaschists are dominantly of Variscan age (330–300 Ma). Garnet in these rocks displays well developed growth zonations in Fe–Mg–Ca–Mn and crystallized at maximal temperatures of 670°C/7 kbar to the west and 600°C/7–8 kbar to the east. In consequence the monazite is interpreted to date a slightly pressure-dominated Variscan amphibolite-facies evolution. In mylonitic garnet gneisses, large metamorphic monazite grains of Ordovician–Silurian (~440 Ma) age but also small monazite grains of Variscan (~300 Ma) age were discovered. Garnets in the mylonitic garnet gneisses display high-temperature homogenized Mg-rich profiles in their cores and crystallized near to ~800°C/6 kbar. The Ordovician–Silurian-age monazites can be assigned to a pre-Variscan high-temperature event recorded by the homogenised garnets. These monazite age data confirm Ordovician–Silurian and Devonian–Carboniferous metamorphic cycles which were already reported from other Alpine domains and further regions in the internal Variscides.     

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⁻¹.