Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessment of dissolved nutrients dispersal derived from offshore fish-farm using nitrogen stable isotope ratios (δ15N) in macroalgal bioassays

In this study, the dispersal of wastes from offshore fish farms was evaluated by analyzing nitrogen stable isotope ratios (δ¹⁵N) in macroalgae incubated in the water column at sites located at an increasing distance from the fish cages. Bioassays were performed at three fish farms situated in separate localities with different nutritional conditions (Canary Islands, Murcia and Catalonia) and varying in size, species of fish reared and annual production. Macroalgal bioassays were carried out in two different directions (DI and DII) and they were replicated at each distance in order to evaluate the effect of small-scale variability on the spatial extent of fish farm wastes. The results obtained with δ¹⁵N contribute to a better understanding of the application of nitrogen stable isotopes ratios in macroalgae as an effective bioindicator for tracing the dispersion of offshore fish farm wastes, and demonstrate that fish farm wastes can be traced even over distances of some km from the pollution source. In the Canary Islands, the maximum distance obtained for detection of fish farm wastes was between 450 and 700 m. Of the three installations studied, Murcia presented the greatest distance for detection of fish farm waste influence, ranging from between 1550 and 2450 m, whilst in Catalonia this distance was less than 120 m. In Catalonia, the results were masked by the influence of other sources of nitrogen, and thus fish farm wastes were detected at more reduced distances than expected. These results confirm that fish farm wastes can be traced using the nitrogen stable isotope ratios of macroalgae and that this method can also be useful for identifying areas of potential risk to some sensitive ecosystems, and as an early signal that changes in the community structure might occur.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.     

Provenance of Holocene beach sand in the Western Iberian margin: the use of the Kolmogorov–Smirnov test for the deciphering of sediment recycling in a modern coastal system

Detrital zircons from Holocene beach sand and igneous zircons from the Cretaceous syenite forming Cape Sines (Western Iberian margin) were dated using laser ablation – inductively coupled plasma – mass spectrometry. The U–Pb ages obtained were used for comparison with previous radiometric data from Carboniferous greywacke, Pliocene–Pleistocene sand and Cretaceous syenite forming the sea cliff at Cape Sines and the contiguous coast. New U–Pb dating of igneous morphologically simple and complex zircons from the syenite of the Sines pluton suggests that the history of zircon crystallization was more extensive (ca 87 to 74 Ma), in contrast to the findings of previous geochronology studies (ca 76 to 74 Ma). The U–Pb ages obtained in Holocene sand revealed a wide interval, ranging from the Cretaceous to the Archean, with predominance of Cretaceous (37%), Palaeozoic (35%) and Neoproterozoic (19%) detrital‐zircon ages. The paucity of round to sub‐rounded grains seems to indicate a short transportation history for most of the Cretaceous zircons (ca 95 to 73 Ma) which are more abundant in the beach sand that was sampled south of Cape Sines. Comparative analysis using the Kolmogorov–Smirnov statistical method, analysing sub‐populations separately, suggests that the zircon populations of the Carboniferous and Cretaceous rocks forming the sea cliff were reproduced faithfully in Quaternary sand, indicating sediment recycling. The similarity of the pre‐Cretaceous ages (>ca 280 Ma) of detrital zircons found in Holocene sand, as compared with Carboniferous greywacke and Pliocene–Pleistocene sand, provides support for the hypothesis that detritus was reworked into the beach from older sedimentary rocks exposed along the sea cliff. The largest percentage of Cretaceous zircons (<ca 95 Ma) found in Holocene sand, as compared with Pliocene–Pleistocene sand (secondary recycled source), suggests that the Sines pluton was the one of the primary sources that became progressively more exposed to erosion during Quaternary uplift. This work highlights the application of the Kolmogorov–Smirnov method in comparison of zircon age populations used to identify provenance and sediment recycling in modern and ancient detrital sedimentary sequences.     

Biogenic volatile compound emissions from a temperate forest,China: model simulation

Emissions of biogenic volatile organic compounds (BVOC) were measured using a relaxed eddy accumulation (REA) technique on an above-canopy tower in a temperate forest (Changbai Mountain, Jilin province, China) during the 2010 and 2011 summer seasons. Solar global radiation and photosynthetically active radiation (PAR) were also measured. Based on PAR energy dynamic balance, an empirical BVOC emission and PAR transfer model was developed that includes the processes of BVOC emissions and PAR transfer above the canopy level, including PAR absorption and consumption, and scattering by gases, liquids, and particles (GLPs). Simulated emissions of isoprene and monoterpenes were in agreement with observations. The averages of the relative estimator biases for the flux were 39.3 % for isoprene, and 27.1 % for monoterpenes in the 2010 and 2011 growing seasons, with NMSE (normalized mean square error) values of 0.133 and 0.101, respectively. The observed and simulated mean diurnal variations of isoprene and monoterpenes in the 2010 and 2011 growing seasons were evaluated for the validation of the empirical model. Under observed atmospheric conditions, the sensitivity analysis showed that emissions of isoprene and monoterpenes were more sensitive to changes in PAR than to water vapor content or to the magnitude of the scattering factor. The emissions of isoprene and monoterpenes in the 2010 and 2011 growing seasons (from June to September) were estimated using this empirical model along with hourly observational data, with mean hourly emissions of 1.71 and 1.55 mg m^?2 h^?1 for isoprene, and 0.48 and 0.47 mg m^?2 h^?1 for monoterpenes in 2010 and 2011, respectively. As formaldehyde (HCHO) is considered as the main oxidation product of isoprene and monoterpenes, it is necessary to investigate the link between HCHO and BVOC emissions. GOME-2 HCHO vertical column densities (VCDs) can be used to estimate BVOC emission fluxes in the Changbai Mountain temperate forest.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

An open source GIS‐based Planning Support System: Application to the land use plan of La Troncal,Ecuador

Planning Support Systems (PSS) comprise a wide variety of geo‐technological tools related to GIS and spatial modeling aimed at addressing land planning processes. This article describes the OpenRules system, a PSS based on a previous system called RULES. Among OpenRules new features are its architecture, based exclusively on free and open source software, and its applicability to all land use types, including rural and urban uses. In addition, OpenRules incorporates an unlimited number of land evaluation factors and a new objective in land use spatial allocation. OpenRules has been programmed in Java and implemented as a module of the free GIS software gvSIG, with full integration between the GIS and the decision support tools. Decision support tools include multicriteria evaluation, multiobjective linear programming and heuristic techniques, which support three basic stages of land use planning processes, namely land suitability evaluation, land use area optimization and land use spatial allocation. The application of OpenRules to the region of La Troncal, Ecuador, demonstrates its capability to generate alternative and coherent solutions through a scientific and justified procedure at low cost in terms of time and resources.     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.